Commit 6bb8543b authored by Kristof Sarosi's avatar Kristof Sarosi
Browse files

introduced forced restart

parent 9934331e
...@@ -399,8 +399,8 @@ int main (int argc, char** argv) ...@@ -399,8 +399,8 @@ int main (int argc, char** argv)
std::string restartInput = ""; //command_line.follow ("noRestart", 2, "-r", "--restart"); std::string restartInput = ""; //command_line.follow ("noRestart", 2, "-r", "--restart");
const bool restart ( restartInput != "noRestart" ); const bool restart ( restartInput != "noRestart" );
if ( restart ) //if ( restart )
{ //{
LifeChrono chronoRestart; LifeChrono chronoRestart;
chronoRestart.start(); chronoRestart.start();
...@@ -483,7 +483,7 @@ int main (int argc, char** argv) ...@@ -483,7 +483,7 @@ int main (int argc, char** argv)
std::cout << " Restart data at TIME = " << t_ << " imported after " << chronoRestart.diff() << " s"; std::cout << " Restart data at TIME = " << t_ << " imported after " << chronoRestart.diff() << " s";
std::cout << "\n\n*****************************************************************\n"; std::cout << "\n\n*****************************************************************\n";
} }
} }
// Circulation export and/or restart // Circulation export and/or restart
...@@ -525,14 +525,14 @@ int main (int argc, char** argv) ...@@ -525,14 +525,14 @@ int main (int argc, char** argv)
modifyPressureBC(bcValues); modifyPressureBC(bcValues);
solver.bcInterfacePtr() -> updatePhysicalSolverVariables(); solver.bcInterfacePtr() -> updatePhysicalSolverVariables();
solver.solveMechanics(); solver.solveMechanics();
} //}
//============================================ //============================================
// Preload // Preload
//============================================ //============================================
if ( ! restart ) /*if ( ! restart )
{ {
solver.structuralOperatorPtr() -> data() -> dataTime() -> setTime(0.0); solver.structuralOperatorPtr() -> data() -> dataTime() -> setTime(0.0);
...@@ -620,7 +620,7 @@ int main (int argc, char** argv) ...@@ -620,7 +620,7 @@ int main (int argc, char** argv)
std::cout << "\n*****************************************************************\n"; std::cout << "\n*****************************************************************\n";
} }
} }*/
//============================================ //============================================
...@@ -650,11 +650,11 @@ int main (int argc, char** argv) ...@@ -650,11 +650,11 @@ int main (int argc, char** argv)
return p; return p;
}; };
if ( ! restart ) /*if ( ! restart )
{ {
heartSolver.postProcess(t); heartSolver.postProcess(t);
circulationSolver.exportSolution( circulationOutputFile ); circulationSolver.exportSolution( circulationOutputFile );
} }*/
LifeChrono chronoExport; LifeChrono chronoExport;
chronoExport.start(); chronoExport.start();
...@@ -722,7 +722,7 @@ int main (int argc, char** argv) ...@@ -722,7 +722,7 @@ int main (int argc, char** argv)
// Iterate mechanics / circulation // Iterate mechanics / circulation
//============================================ //============================================
if ( makeMechanicsCirculationCoupling ) /*if ( makeMechanicsCirculationCoupling )
{ {
iter = 0; iter = 0;
const double dt_circulation ( dt_mechanics / 1000 ); const double dt_circulation ( dt_mechanics / 1000 );
...@@ -925,9 +925,11 @@ int main (int argc, char** argv) ...@@ -925,9 +925,11 @@ int main (int argc, char** argv)
dispPre = disp; dispPre = disp;
} }*/
VFeNew[0] = LV.volume(disp, dETFESpace, - 1);
VFeNew[1] = RV.volume(disp, dETFESpace, 1);
/*
//============================================ //============================================
// Export FE-solution // Export FE-solution
//============================================ //============================================
...@@ -947,7 +949,7 @@ int main (int argc, char** argv) ...@@ -947,7 +949,7 @@ int main (int argc, char** argv)
std::cout << "\nPrevious " << dt_save << " ms computed in " << chronoDiffToLastSave << " s"; std::cout << "\nPrevious " << dt_save << " ms computed in " << chronoDiffToLastSave << " s";
std::cout << "\n<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<\n\n"; std::cout << "\n<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<\n\n";
} }
} }*/
} }
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment