Detailed Description

template<typename ShapeFunction>
class ApplicationUtils::ComputeNaturalCoordsSolverImplementation< ShapeFunction >

Definition at line 42 of file ComputeNaturalCoordsSolver.h.

#include <ComputeNaturalCoordsSolver.h>

Inheritance diagram for ApplicationUtils::ComputeNaturalCoordsSolverImplementation< ShapeFunction >:

Collaboration diagram for ApplicationUtils::ComputeNaturalCoordsSolverImplementation< ShapeFunction >:

Public Member Functions
Eigen::Vector3d	solve (MeshLib::Element const &e, Eigen::Vector3d const &real_coords, int const max_iter, double const real_coords_tolerance) const override

Public Member Functions inherited from ApplicationUtils::ComputeNaturalCoordsSolverInterface
virtual	~ComputeNaturalCoordsSolverInterface ()=default

Private Types
using	ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, ElementDim>

Static Private Attributes
static constexpr int	ElementDim = ShapeFunction::DIM

Member Typedef Documentation

◆ ShapeMatricesType

template<typename ShapeFunction >

using ApplicationUtils::ComputeNaturalCoordsSolverImplementation< ShapeFunction >::ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, ElementDim>

private

Definition at line 46 of file ComputeNaturalCoordsSolver.h.

Member Function Documentation

◆ solve()

template<typename ShapeFunction >

Eigen::Vector3d ApplicationUtils::ComputeNaturalCoordsSolverImplementation< ShapeFunction >::solve	(	MeshLib::Element const &	e,
		Eigen::Vector3d const &	real_coords,
		int const	max_iter,
		double const	real_coords_tolerance ) const

inlineoverridevirtual

Implements ApplicationUtils::ComputeNaturalCoordsSolverInterface.

Definition at line 49 of file ComputeNaturalCoordsSolver.h.

    {
        using Problem = ComputeNaturalCoordsRootFindingProblem<ShapeFunction>;
        using LJM = typename Problem::LocalJacobianMatrix;
        using LRV = typename Problem::LocalResidualVector;
 
        Problem problem(e, Eigen::Map<const LRV>(real_coords.data()));
 
        Eigen::PartialPivLU<LJM> linear_solver(ElementDim);
        LJM jacobian;
 
        const double increment_tolerance = 0;
        auto const newton_solver = NumLib::makeNewtonRaphson(
            linear_solver,
            [&problem](LJM& J) { problem.updateJacobian(J); },
            [&problem](LRV& res) { problem.updateResidual(res); },
            [&problem](auto& delta_r) { problem.updateSolution(delta_r); },
            {max_iter, real_coords_tolerance, increment_tolerance});
 
        auto const opt_iter = newton_solver.solve(jacobian);
 
        if (opt_iter)
        {
            DBUG("Newton solver succeeded after {} iterations", *opt_iter);
 
            Eigen::Vector3d natural_coords = Eigen::Vector3d::Zero();
            natural_coords.head<ElementDim>() = problem.getNaturalCoordinates();
            return natural_coords;
        }
 
        OGS_FATAL(
            "Newton solver failed. Please consider increasing the error "
            "tolerance or the max. number of Newton iterations.");
    }