Commit c56fd2a4 authored by thomaskummer's avatar thomaskummer
Browse files

cervin lifev update

parent b0fa36da
......@@ -267,7 +267,7 @@ int main (int argc, char** argv)
//============================================
// Building Matrices
//============================================
solver.oneWayCoupling();
solver.twoWayCoupling();
solver.structuralOperatorPtr()->setNewtonParameters(dataFile);
solver.buildSystem();
......@@ -604,7 +604,7 @@ int main (int argc, char** argv)
chronoSave.start();
solver.saveSolution (-1.0);
heartSolver.exporter()->postProcess(-1.0);
heartSolver.postProcess(-1.0);
if ( 0 == comm->MyPID() )
{
......@@ -634,7 +634,7 @@ int main (int argc, char** argv)
solver.bcInterfacePtr() -> updatePhysicalSolverVariables();
solver.solveMechanics();
//solver.saveSolution (i-1);
//heartSolver.exporter()->postProcess(i-1);
//heartSolver.postProcess(i-1);
}
auto maxI4fValue ( solver.activationModelPtr()->I4f().maxValue() );
......@@ -680,7 +680,7 @@ int main (int argc, char** argv)
if ( ! restart )
{
solver.saveSolution(t);
heartSolver.exporter()->postProcess(t);
heartSolver.postProcess(t);
circulationSolver.exportSolution( circulationOutputFile );
}
......@@ -933,7 +933,7 @@ int main (int argc, char** argv)
if ( save )
{
solver.saveSolution(t);
heartSolver.exporter()->postProcess(t);
heartSolver.postProcess(t);
}
}
......
......@@ -254,17 +254,36 @@ public:
m_exporter.reset (new exporter_Type());
setupExporter<mesh_Type>(*m_exporter, M_emSolver.localMeshPtr(), M_emSolver.comm(), outputFileName, problemFolder);
m_exporter->addVariable ( ExporterData<RegionMesh<LinearTetra> >::VectorField,
"displacement",
M_emSolver.structuralOperatorPtr()->dispFESpacePtr(),
M_emSolver.structuralOperatorPtr()->displacementPtr(),
UInt (0) );
m_exporter->addVariable ( ExporterData<RegionMesh<LinearTetra> >::VectorField,
"displacement",
M_emSolver.structuralOperatorPtr()->dispFESpacePtr(),
M_emSolver.structuralOperatorPtr()->displacementPtr(),
UInt (0) );
m_exporter->addVariable ( ExporterData<RegionMesh<LinearTetra> >::ScalarField,
"Von Mises Stress",
M_emSolver.electroSolverPtr()->feSpacePtr(),
M_emSolver.tensionEstimator().vonMisesStressPtr(),
UInt (0) );
m_exporter->addVariable ( ExporterData<RegionMesh<LinearTetra> >::VectorField,
"fibers",
M_emSolver.structuralOperatorPtr()->dispFESpacePtr(),
M_emSolver.structuralOperatorPtr()->EMMaterial()->fiberVectorPtr(),
UInt (0) );
m_exporter->addVariable ( ExporterData<RegionMesh<LinearTetra> >::VectorField,
"fibers",
M_emSolver.structuralOperatorPtr()->dispFESpacePtr(),
M_emSolver.structuralOperatorPtr()->EMMaterial()->sheetVectorPtr(),
UInt (0) );
m_exporter->addVariable ( ExporterData<RegionMesh<LinearTetra> >::ScalarField,
"Activation",
M_emSolver.electroSolverPtr()->feSpacePtr(),
M_emSolver.activationModelPtr()->fiberActivationPtr(),
UInt (0) );
m_exporter->addVariable ( ExporterData<RegionMesh<LinearTetra> >::ScalarField,
"Von Mises Stress",
M_emSolver.electroSolverPtr()->feSpacePtr(),
M_emSolver.tensionEstimator().vonMisesStressPtr(),
UInt (0) );
}
......@@ -284,7 +303,7 @@ public:
void postProcess(const Real& time)
{
M_emSolver.tensionEstimator().setDisplacement ( M_emSolver.structuralOperatorPtr()-> displacement() );
M_emSolver.tensionEstimator().setDisplacement ( M_emSolver.structuralOperatorPtr()->displacement() );
M_emSolver.tensionEstimator().analyzeTensionsRecoveryVonMisesStress();
m_exporter->postProcess(time);
......
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