Commit 930a6951 authored by thomaskummer's avatar thomaskummer
Browse files

euler lifev update

parent 32fa5749
......@@ -1083,25 +1083,25 @@ 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) = dUdx[j];
// F(0,2) = dUdx[k];
// F(1,0) = dUdy[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) = dUdy[k];
// F(2,0) = dUdz[i];
// F(2,1) = dUdz[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) = dUdx[j];
F(0,2) = dUdx[k];
F(1,0) = dUdy[i];
F(1,1) = 1.0 + dUdy[j];
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];
f0(1) = f0_[j];
......@@ -1110,7 +1110,7 @@ EMSolver<Mesh, ElectroSolver>::computeDeformedFiberDirection (VectorEpetra& f_,
f0.normalize();
auto f = F * f0;
f_[i] = f(0);
f_[i] = 5.; //f(0);
f_[j] = f(1);
f_[k] = f(2);
}
......
......@@ -247,7 +247,7 @@ public:
UInt (0) );
m_exporter->addVariable ( ExporterData<RegionMesh<LinearTetra> >::VectorField,
"Fibers0",
"Fibers",
M_emSolver.structuralOperatorPtr()->dispFESpacePtr(),
M_emSolver.structuralOperatorPtr()->fPtr(),
//M_emSolver.structuralOperatorPtr()->EMMaterial()->fiberVectorPtr(),
......
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