Commit 43da6fbc authored by thomaskummer's avatar thomaskummer
Browse files

euler lifev update

parent f0ce555b
......@@ -1063,17 +1063,17 @@ 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, Unique);
VectorEpetra dUdx (disp.map(), Repeated);
VectorEpetra dUdy (disp.map(), Repeated);
VectorEpetra dUdz (disp.map(), Repeated);
VectorEpetra dUdx (Unique);
VectorEpetra dUdy (disp.map());
VectorEpetra dUdz (disp.map());
dUdx = GradientRecovery::ZZGradient (feSpacePtr, disp, 0);
dUdy = GradientRecovery::ZZGradient (feSpacePtr, disp, 1);
dUdz = GradientRecovery::ZZGradient (feSpacePtr, disp, 2);
f_ = dUdx;
f_.reset(new vector_Type (dUdx ));
return;
int n = f_.epetraVector().MyLength() / 3;
......
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