Commit c4e8250c authored by thomaskummer's avatar thomaskummer
Browse files

intrepid lifev update

parent 80cce8b1
......@@ -618,6 +618,8 @@ int main (int argc, char** argv)
VCirc = VFe;
VectorEpetra dispCurrent ( disp );
VectorEpetra dispPre ( disp );
ID bdPowerFlag = dataFile ( ("solid/boundary_conditions/LVEndo/flag") , 0 );
printCoupling("Initial values");
......@@ -887,6 +889,25 @@ int main (int argc, char** argv)
}
//============================================
// Power computations
//============================================
Real leftVentPower = heartSolver.externalPower(disp, dispPre, dETFESpace, p("lv"), dt_mechanics, 454);
Real rightVentPower = heartSolver.externalPower(disp, dispPre, dETFESpace, p("rv"), dt_mechanics, 455);
//Real patchPower1 = heartSolver.externalPower(disp, dispPre, dETFESpace, p("lv"), dt_mechanics, 454);
//Real patchPower2 = heartSolver.externalPower(disp, dispPre, dETFESpace, p("lv"), dt_mechanics, 454);
if ( 0 == comm->MyPID() )
{
std::cout << "\n******************************************";
std::cout << "\nLeft ventricular power is " << leftVentPower;
std::cout << "\nRight ventricular power is " << rightVentPower;
std::cout << "\n******************************************\n\n";
}
dispPre = disp;
//============================================
// Export FE-solution
//============================================
......
......@@ -663,9 +663,9 @@ EMSolver<Mesh, ElectroSolver>::setupMechanicalSolver ( GetPot& dataFile)
M_commPtr) );
std::string data_file_name = dataFile.get (0, "NO_DATA_FILENAME_FOUND");
dFESpace->setQuadRule(quadRuleTetra1pt);
if (M_commPtr -> MyPID() == 0) std::cout << "\nPolynomial degree: " << dFESpace->polynomialDegree();
if (M_commPtr -> MyPID() == 0) std::cout << "\nDof: " << dFESpace->dof() << std::endl;
dFESpace->setQuadRule(quadRuleTetra4pt);
//if (M_commPtr -> MyPID() == 0) std::cout << "\nPolynomial degree: " << dFESpace->polynomialDegree();
//if (M_commPtr -> MyPID() == 0) std::cout << "\nDof: " << dFESpace->dof() << std::endl;
setupMechanicalBC (data_file_name, "solid", dFESpace);
M_EMStructuralOperatorPtr.reset (new structuralOperator_Type() );
......
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