Adapt existing stellar preparation to Newton without giving acceptance to the model. More...

#include <stellar_problem.hpp>

Collaboration diagram for serif::solver::backend::mfem::StellarProblem< Evaluator, Inputs >:

Public Types
using	Vector = serif::Vector

using	System = std::remove_const_t<typename Evaluator::Result::element_type>

using	Model = std::remove_cvref_t<decltype(std::declval<const Evaluator&>().model())>

using	Layout = typename Evaluator::Layout

using	Failure = typename Evaluator::Cache::Failure

Public Member Functions
	StellarProblem (Evaluator evaluator, Inputs prescribed=Inputs{}, ::mfem::Vector displacement={})
	Retain one cache and all external inputs throughout a solve.

auto	model_owner () const
	Retain the immutable physical definition in the resulting checkpoint.

const auto &	binding () const
	Inspect typed FE resources for explicitly selected numerical strategies.

const Layout &	layout () const
	Inspect the coordinate identities used by every physical state in this solve.

Evaluator &	evaluator ()
	Access the shared physical preparation service for numerical extensions.

auto	evaluate (const model::StellarState< Layout > &state)
	Evaluate the complete equations, preserving all previously published systems.

auto	prepare (const model::StellarState< Layout > &state)
	Prepare the inverse belonging to exactly this physical candidate.

linear::Report	solve_linear (const auto &prepared, const Vector &rhs, Vector &correction, const auto &options)
	Solve with the operator and inverse retained by this candidate's preparation.

auto	geometry_limit (const model::StellarState< Layout > &state, std::span< const double > direction, const double maximum_step, const double safety, const bool enabled)
	Snapshot the accepted geometry direction for step restriction and backtracking.

double	next_step (const double step, const Backtracking &options, const auto &evaluation) const
	Skip unsafe backtracking rungs after a structured orientation failure.

Static Public Member Functions
static std::string	failure_message (const Failure &failure)
	Preserve the original failure message without classifying by its text.

static std::string	failure_message (const geometry::deformation::DeformationFailure &failure)
	Geometry preflight already reports an ordinary deformation failure.

static bool	recoverable (const Failure &failure)
	Backtracking may repair a numerical candidate, never an incompatible discretization.

Private Member Functions
auto	invoke (auto callback)
	Keep fixed displacement outside the unknown vector and omit it for moving stars.

Static Private Member Functions
static bool	inverted (const Failure &failure)
	Recognize orientation failure from the structured physical rejection.

Private Attributes
Evaluator	m_evaluator

Inputs	m_prescribed

::mfem::Vector	m_fixed_displacement

admissibility::backend::mfem::GeometryDirection	m_geometry

Detailed Description

template<typename Evaluator, typename Inputs = utils::dependencies::DependencyContext<>>
class serif::solver::backend::mfem::StellarProblem< Evaluator, Inputs >

Adapt existing stellar preparation to Newton without giving acceptance to the model.

Note: The prescribed input context and a copy of fixed geometry are retained for the entire solve. Explicitly borrowed resources keep their existing lifetime requirements. The public solver factory constructs this object; astronomers normally work with a discretization and a state.