Commit 2a124391 authored by thomaskummer's avatar thomaskummer
Browse files

intrepid lifev update

parent c0f6b90c
......@@ -376,6 +376,8 @@ int main (int argc, char** argv)
VectorSmall<2> VCirc, VCircNew, VCircPert, VFe, VFeNew, VFePert, R, dp;
MatrixSmall<2,2> JFe, JCirc, JR;
VectorSmall<2> AvgWorkVent;
UInt iter (0);
Real t (0);
......@@ -893,17 +895,17 @@ int main (int argc, char** argv)
//============================================
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);
Real dispNormInf;
disp.normInf(dispNormInf);
//Real patchPower1 = heartSolver.externalPower(disp, dispPre, dETFESpace, p("lv"), dt_mechanics, 900);
//Real patchPower2 = heartSolver.externalPower(disp, dispPre, dETFESpace, p("rv"), dt_mechanics, 901);
AvgWorkVent(0) += leftVentPower * dt_mechanics;
AvgWorkVent(1) += rightVentPower * dt_mechanics;
if ( 0 == comm->MyPID() )
{
std::cout << "\n******************************************";
std::cout << "\nLeft ventricular power is " << leftVentPower;// << " " << p("lv");
std::cout << "\nRight ventricular power is " << dispNormInf;
std::cout << "\nInstantaneous left and right ventricular power is " << leftVentPower << " / " << rightVentPower;
std::cout << "\nAveraged left and right ventricular power is " << AvgWorkVent(0) / t << " / " << AvgWorkVent(1) / t;
std::cout << "\n******************************************\n\n";
}
......
......@@ -399,7 +399,7 @@ public:
integrate ( boundary ( dispETFESpace->mesh(), bdFlag),
myBDQR,
dispETFESpace,
p * J * dot( FmT * Nface, phi_i)
value(-1.0) * p * J * dot( FmT * Nface, phi_i)
//p * J * dot( FmT * Nface, phi_i)
//value(-1.0) * J * dot (vE1, FmT * Nface) * phi_i) >> intergral
) >> traction;
......
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