Commit eb9444c0 authored by GILLES Sebastien's avatar GILLES Sebastien
Browse files

#820 Delta solid solve yields the expected result.

parent f6841aa5
......@@ -417,30 +417,30 @@ namespace HappyHeart
delta_residual.Scale(-1., __FILE__, __LINE__);
}
// decltype(auto) tangent_varf = GetTangentSolidVarf();
//
// delta_solid_varf_solver_->SolveLinear(tangent_varf,
// tangent_varf,
// GetDeltaResidual(),
// GetNonCstDeltaSolution(),
// __FILE__, __LINE__);
//
//
//
// {
// decltype(auto) delta_sol = GetDeltaSolution();
// const auto filename = parent::GetOutputDirectory(numbering_subset)
// + "/"
// + DifferentialPreffix(differential::yes)
// + "solid_varf_sol_" + GetIterationTag() + ".m";
//
// delta_sol.View(parent::MpiHappyHeart(),
// filename,
// __FILE__, __LINE__,
// PETSC_VIEWER_ASCII_MATLAB);
//
// }
//
decltype(auto) tangent_varf = GetTangentSolidVarf();
delta_solid_varf_solver_->SolveLinear(tangent_varf,
tangent_varf,
GetDeltaResidual(),
GetNonCstDeltaSolution(),
__FILE__, __LINE__);
{
decltype(auto) delta_sol = GetDeltaSolution();
const auto filename = parent::GetOutputDirectory(numbering_subset)
+ "/"
+ DifferentialPreffix(differential::yes)
+ "solid_varf_sol_" + GetIterationTag() + ".m";
delta_sol.View(parent::MpiHappyHeart(),
filename,
__FILE__, __LINE__,
PETSC_VIEWER_ASCII_MATLAB);
}
// \todo #820 At the moment, I accept my norm is not the same as Bruno's; see what Petsc gives when
......
......@@ -226,7 +226,7 @@ namespace HappyHeart
#ifndef NDEBUG
invariant_holder.Reset();
#endif // NDEBUG
#endif // NDEBUG
}
......
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