Commit 499b154d authored by thomaskummer's avatar thomaskummer
Browse files

cervin lifev update

parent 0479a425
......@@ -231,8 +231,8 @@ monodomain_xml_file = ParamListE.xml
[../time_discretization]
initialtime = 0.
endtime = 100000.
timestep = 2.0
dt_loadstep = 1.0
timestep = 1.0
dt_loadstep = 0.5
activation_limit_loadstep = -0.1
theta = 0.35
zeta = 0.75
......@@ -269,7 +269,7 @@ monodomain_xml_file = ParamListE.xml
[../newton]
maxiter = 50
reltol = 1.e-7
abstol = 1e-5
abstol = 0.8e-6
NonLinearLineSearch = 0
[../solver]
......
......@@ -69,6 +69,7 @@ protected:
M_activationLimit_loadstep = M_datafile ( "solid/time_discretization/activation_limit_loadstep", 0.0 );
M_dt_mechanics = M_datafile ("solid/time_discretization/timestep", 1.0 );
M_dt_save = M_datafile ( "exporter/save", 10. );
M_safePreload = M_datafile ( "exporter/savePreload", false );
M_endtime = M_datafile ("solid/time_discretization/endtime", 100000);
M_mechanicsLoadstepIter = static_cast<UInt>( M_dt_loadstep / M_dt_activation );
M_mechanicsCouplingIter = static_cast<UInt>( M_dt_mechanics / M_dt_activation );
......@@ -105,6 +106,7 @@ protected:
Real M_activationLimit_loadstep;
Real M_dt_mechanics;
Real M_dt_save;
bool M_safePreload;
Real M_endtime;
UInt M_mechanicsLoadstepIter;
UInt M_mechanicsCouplingIter;
......@@ -187,7 +189,9 @@ public:
// Solve mechanics
M_emSolver.bcInterfacePtr() -> updatePhysicalSolverVariables();
M_emSolver.solveMechanics();
//solver.saveSolution (i-1);
// Safe preload steps
if ( M_safePreload ) solver.saveSolution (i-1);
}
if ( 0 == M_emSolver.comm()->MyPID() )
......
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