Commit 4a82bfc9 authored by GILLES Sebastien's avatar GILLES Sebastien
Browse files

#820 Add new output files for comparison with Freefem (solid residual is...

#820 Add new output files for comparison with Freefem (solid residual is currently not quite right).
parent d14c0dd9
......@@ -229,18 +229,26 @@ namespace HappyHeart
decltype(auto) variable_holder = this->GetNonCstVariableHolder();
variable_holder.solid_newton_index = this->GetSnes().GetSnesIteration(__FILE__, __LINE__);
auto& vectors_and_matrices = this->GetNonCstVectorsAndMatrices();
const auto& displacement_previous_time_iteration =
vectors_and_matrices.GetDisplacementPreviousTimeIteration();
variable_holder.template DevPrint<DevPhase::solid_newton>(vectors_and_matrices.GetEvaluationState(),
"dispSr_before_solid_solve");
variable_holder.template DevPrint<DevPhase::solid_newton>(displacement_previous_time_iteration,
"dispSr_old_before_solid_solve");
const auto& numbering_subset = GetNumberingSubset();
auto& vectors_and_matrices = this->GetNonCstVectorsAndMatrices();
auto& rhs = this->GetNonCstSystemRhs(numbering_subset);
rhs.ZeroEntries(__FILE__, __LINE__);
HyperelasticityNS::Private::ComputeDynamicTimeSchemeResidualContribution(vectors_and_matrices, rhs);
const auto& displacement_previous_time_iteration =
vectors_and_matrices.GetDisplacementPreviousTimeIteration();
const auto& velocity_previous_time_iteration =
vectors_and_matrices.GetVelocityPreviousTimeIteration();
......
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