Commit 061b5c69 authored by thomaskummer's avatar thomaskummer
Browse files

intrepid lifev update

parent fda806a7
......@@ -386,10 +386,10 @@ public:
return M_activationModelPtr;
}
// vectorPtr_Type activationTimePtr()
// {
// return M_activationTimePtr;
// }
vectorPtr_Type activationTimePtr()
{
return M_activationTimePtr;
}
void saveSolution (Real time, const bool& restart = 0);
......@@ -502,7 +502,7 @@ public:
meshPtr_Type M_localMeshPtr;
meshPtr_Type M_fullMeshPtr;
//vectorPtr_Type M_activationTimePtr;
vectorPtr_Type M_activationTimePtr;
bool M_oneWayCoupling;
......@@ -531,7 +531,7 @@ EMSolver<Mesh, ElectroSolver>::EMSolver(commPtr_Type comm) :
M_mechanicsExporterPtr ( ),
M_localMeshPtr ( ),
M_fullMeshPtr ( ),
//M_activationTimePtr ( ),
M_activationTimePtr ( ),
M_oneWayCoupling (true),
M_wteTotal ( ),
// M_wtePassive ( ),
......@@ -555,7 +555,7 @@ EMSolver<Mesh, ElectroSolver>::EMSolver (const EMSolver& solver) :
M_mechanicsExporterPtr ( solver.M_mechanicsExporterPtr),
M_localMeshPtr ( solver.M_localMeshPtr),
M_fullMeshPtr ( solver.M_fullMeshPtr),
//M_activationTimePtr ( solver.M_activationTimePtr),
M_activationTimePtr ( solver.M_activationTimePtr),
M_oneWayCoupling ( solver.M_oneWayCoupling),
M_wteTotal (solver.M_wteTotal),
// M_wtePassive (solver.M_wtePassive),
......@@ -630,6 +630,9 @@ EMSolver<Mesh, ElectroSolver>::setupElectroSolver ( GetPot& dataFile )
M_electroSolverPtr -> setup (dataFile, ionicModelPtr->Size() );
// }
M_activationTimePtr.reset (new vector_Type ( M_electroSolverPtr->potentialPtr() -> map() ));
*M_activationTimePtr = -1.0;
if (M_commPtr -> MyPID() == 0)
{
std::cout << "\nEMSolver: setupElectroSolver - done";
......@@ -974,7 +977,7 @@ EMSolver<Mesh, ElectroSolver>::solveElectrophysiology (function_Type& stimulus,
M_electroSolverPtr -> solveOneICIStep();
//M_electroSolverPtr -> registerActivationTime (*M_activationTimePtr, time, 0.9);
M_electroSolverPtr -> registerActivationTime (*M_activationTimePtr, time, 0.9);
if (M_commPtr -> MyPID() == 0)
{
std::cout << "\nEMSolver: solveElectrophysiology - done";
......
......@@ -284,6 +284,12 @@ public:
M_emSolver.electroSolverPtr()->feSpacePtr(),
M_emSolver.activationModelPtr()->fiberActivationPtr(),
UInt (0) );
m_exporter -> addVariable ( ExporterData<RegionMesh<LinearTetra> >::ScalarField,
"Activation Time",
M_emSolver.electroSolverPtr() -> feSpacePtr(),
M_activationTimePtr,
UInt (0) );
for (int i = 0; i < M_emSolver.electroSolverPtr()->ionicModelPtr()->Size(); ++i)
{
......
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