save flow field only when it is newly computed.

LLNL · Mar 18, 2024 · 73c32ea · 73c32ea
1 parent 272c85c
commit 73c32ea
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/src/advdiff_solver.cpp b/src/advdiff_solver.cpp
@@ -243,8 +243,12 @@ void AdvDiffSolver::GetFlowField(ParameterizedProblem *flow_problem)
    stokes_solver->InitVariables();
    stokes_solver->SetSolutionSaveMode(save_flow);
 
+   bool flow_loaded = false;
    if (load_flow && FileExists(flow_file))
+   {
       stokes_solver->LoadSolution(flow_file);
+      flow_loaded = true;
+   }
    else
    {
       stokes_solver->SetParameterizedProblem(flow_problem);
@@ -255,7 +259,7 @@ void AdvDiffSolver::GetFlowField(ParameterizedProblem *flow_problem)
       stokes_solver->Solve();
    }
 
-   if (save_flow)
+   if (save_flow && (!flow_loaded))
       stokes_solver->SaveSolution(flow_file);
 
    DeletePointers(flow_coeffs);