Commit 03096b07 authored by GILLES Sebastien's avatar GILLES Sebastien
Browse files

#1292 Remove last mentions of Seldon.

parent 423a3c0f
......@@ -419,9 +419,9 @@ namespace MoReFEM::Advanced::LocalVariationalOperatorNS
xt::noalias(hydrostatic_tangent) = xt::linalg::outer(dI3dC, dI3dC);
xt::noalias(FOO) = (pressure_at_quad_point * std::pow(inv3, - 1.5), hydrostatic_tangent);
xt::noalias(hydrostatic_tangent) *= (pressure_at_quad_point * std::pow(inv3, - 1.5));
Seldon::Add(- 2 * sqrt_inv3 * pressure_at_quad_point, d2I3dCdC, hydrostatic_tangent);
xt::noalias(hydrostatic_tangent) -= 2 * sqrt_inv3 * pressure_at_quad_point * d2I3dCdC;
Internal::LocalVariationalOperatorNS::SecondPiolaKirchhoffStressTensorNS
::ComputeLinearPart(derivative_green_lagrange_at_quad_point,
......@@ -436,7 +436,7 @@ namespace MoReFEM::Advanced::LocalVariationalOperatorNS
Internal::LocalVariationalOperatorNS::SecondPiolaKirchhoffStressTensorNS::ComputeNonLinearPart<FeltSpaceDimensionT>(hydrostatic_stress, tangent_matrix_disp_disp);
Seldon::Add(1., linear_part, tangent_matrix_disp_disp);
xt::noalias(tangent_matrix_disp_disp) += linear_part;
//std::cout << "tangent_matrix_disp_disp"<< std::endl;
//tangent_matrix_disp_disp.Print();
......
......@@ -51,10 +51,7 @@ namespace MoReFEM
* \brief Transpose of Gradient2Strain where Gradient2Strain is defined so that
* strain = Gradient2Strain x gradient(displacement).
*
* \internal <b><tt>[internal]</tt></b> This is stored here as it is cheap to store and might be computed
* many times... and considering Seldon Transpose is not efficient it is best to avoid it as much as possible.
* \endinternal
*
* \todo #1491 Can probably be avoided entirely with Xtensor.
* \tparam DimensionT Dimension of the mesh in which Parameter is defined.
*
* \return Transpose of Gradient2Strain where Gradient2Strain is defined so that
......
......@@ -35,12 +35,6 @@
std::cerr << "Exception caught in the Petsc Snes definition: " << e.what() << std::endl; \
return PETSC_ERR_MIN_VALUE; \
} \
catch(Seldon::Error& e) \
{ \
std::cerr << "Exception caught in the Petsc Snes definition: " << e.What() << std::endl; \
return PETSC_ERR_MIN_VALUE; \
}
# endif // MOREFEM_NO_TRAP_SNES_EXCEPTION
......
......@@ -91,7 +91,7 @@ namespace MoReFEM
/*!
* \brief Check whether there is a nan or inf value in a Seldon matrix.
* \brief Check whether there is a nan or inf value in a \a LocalMatrix.
*
* \param[in] matrix Matrix being scrutinized.
*
......@@ -101,7 +101,7 @@ namespace MoReFEM
/*!
* \brief Check whether there is a nan or inf value in a Seldon vector.
* \brief Check whether there is a nan or inf value in a \a LocalVector.
*
* \param[in] vector Vector being scrutinized.
*
......
......@@ -57,7 +57,7 @@ namespace MoReFEM
* exception was thrown.
* \param[in] exception_message Explanation message encapsulated within the exception. Usually the output
* of \a what() method for exceptions that derive from std::exception, but some exceptions from third-party
* libraries do not and define their own method for that (\a What() for Seldon for instance).
* libraries do not and define their own method for that (currently no such library is in use within MoReFEM).
*/
void PrintAndAbort(const Wrappers::Mpi& mpi,
const std::string& exception_message);
......
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