Commit 32fa5749 authored by thomaskummer's avatar thomaskummer
Browse files

euler lifev update

parent 807cc3c7
......@@ -1063,7 +1063,7 @@ template<typename Mesh , typename ElectroSolver>
void
EMSolver<Mesh, ElectroSolver>::computeDeformedFiberDirection (VectorEpetra& f_, VectorEpetra& f0_, VectorEpetra& disp, solidFESpacePtr_Type feSpacePtr)
{
f_ = VectorEpetra(disp.map(), Unique);
// f_ = VectorEpetra(disp.map(), Unique);
VectorEpetra dUdx (disp);
VectorEpetra dUdy (disp);
......@@ -1110,9 +1110,9 @@ EMSolver<Mesh, ElectroSolver>::computeDeformedFiberDirection (VectorEpetra& f_,
f0.normalize();
auto f = F * f0;
f_[i] = f0_[i]; //f(0);
f_[j] = f0_[j]; //f(1);
f_[k] = f0_[k]; //f(2);
f_[i] = f(0);
f_[j] = f(1);
f_[k] = f(2);
}
}
......
Markdown is supported
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