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

#1292 Print to understand.

parent 84383186
......@@ -164,6 +164,14 @@ namespace MoReFEM
std::cerr << "Bug in SNES implementation: third argument of snes function implemention is expected "
"to be the residual, which is born in VariationalFormulationT by the attribute RHS(). It is not the case "
"here.\n" << desc << std::endl;
variational_formulation.GetSystemRhs(displacement_numbering_subset).View(variational_formulation.GetMpi(), __FILE__, __LINE__);
std::cout << "Other" << std::endl;
VecView(a_residual, PETSC_VIEWER_STDOUT_(variational_formulation.GetMpi().GetCommunicator()));
assert(false);
}
}
......@@ -465,10 +473,10 @@ namespace MoReFEM
auto& system_matrix = GetNonCstSystemMatrix(displacement_numbering_subset, displacement_numbering_subset);
system_matrix.ZeroEntries(__FILE__, __LINE__);
#ifndef NDEBUG
AssertSameNumberingSubset(GetMatrixTangentStiffness(), system_matrix);
AssertSameNumberingSubset(GetMatrixMassPerSquareTimeStep(), system_matrix);
#endif // NDEBUG
#ifndef NDEBUG
AssertSameNumberingSubset(GetMatrixTangentStiffness(), system_matrix);
AssertSameNumberingSubset(GetMatrixMassPerSquareTimeStep(), system_matrix);
#endif // NDEBUG
Wrappers::Petsc::AXPY<NonZeroPattern::same>(1.,
GetMatrixMassPerSquareTimeStep(),
......
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