Commit 9d9ce1e6 authored by thomaskummer's avatar thomaskummer
Browse files

euler lifev update

parent 2ff25d0a
......@@ -1027,14 +1027,24 @@ EMSolver<Mesh, ElectroSolver>::computeI4f (VectorEpetra& i4f, VectorEpetra& f0_,
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];
......@@ -1053,6 +1063,8 @@ template<typename Mesh , typename ElectroSolver>
void
EMSolver<Mesh, ElectroSolver>::computeDeformedFiberDirection (VectorEpetra& f_, VectorEpetra& f0_, VectorEpetra& disp, solidFESpacePtr_Type feSpacePtr)
{
f_ = VectorEpetra(disp.map(), Unique);
VectorEpetra dUdx (disp);
VectorEpetra dUdy (disp);
VectorEpetra dUdz (disp);
......
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