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

euler lifev update

parent d34f9c68
......@@ -1071,14 +1071,24 @@ EMSolver<Mesh, ElectroSolver>::computeDeformedFiberDirection (VectorEpetra& f_,
j = dUdx.blockMap().GID (p + n);
k = dUdx.blockMap().GID (p + 2 * n);
// F(0,0) = 1.0 + dUdx[i];
// F(0,1) = dUdy[i];
// F(0,2) = dUdz[i];
// F(1,0) = dUdx[j];
// F(1,1) = 1.0 + dUdy[j];
// F(1,2) = dUdz[j];
// F(2,0) = dUdx[k];
// F(2,1) = dUdy[k];
// F(2,2) = 1.0 + dUdz[k];
F(0,0) = 1.0 + dUdx[i];
F(0,1) = dUdy[i];
F(0,2) = dUdz[i];
F(1,0) = dUdx[j];
F(0,1) = dUdx[j];
F(0,2) = dUdx[k];
F(1,0) = dUdy[i];
F(1,1) = 1.0 + dUdy[j];
F(1,2) = dUdz[j];
F(2,0) = dUdx[k];
F(2,1) = dUdy[k];
F(1,2) = dUdy[k];
F(2,0) = dUdz[i];
F(2,1) = dUdz[j];
F(2,2) = 1.0 + dUdz[k];
f0(0) = f0_[i];
......
......@@ -249,8 +249,8 @@ public:
m_exporter->addVariable ( ExporterData<RegionMesh<LinearTetra> >::VectorField,
"Fibers",
M_emSolver.structuralOperatorPtr()->dispFESpacePtr(),
//M_emSolver.structuralOperatorPtr()->fPtr(),
M_emSolver.structuralOperatorPtr()->EMMaterial()->fiberVectorPtr(),
M_emSolver.structuralOperatorPtr()->fPtr(),
//M_emSolver.structuralOperatorPtr()->EMMaterial()->fiberVectorPtr(),
UInt (0) );
m_exporter->addVariable ( ExporterData<RegionMesh<LinearTetra> >::VectorField,
......
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