equilibrium_evaluator.hpp expected span serif/backend/mfem/model/binding.hpp serif/backend/mfem/model/failure.hpp serif/backend/mfem/model/state.hpp serif/backend/mfem/model/preconditioning/computations.hpp serif/backend/mfem/shared/finite.hpp serif/utils/dependencies/computation_cache.hpp src/include/serif/backend/mfem/model/equilibrium.hpp serif/backend/mfem/discretization/basis_evaluation.hpp serif/backend/mfem/discretization/boundary_trace.hpp serif/backend/mfem/discretization/coefficient_tracking.hpp serif/backend/mfem/discretization/dof_map.hpp serif/backend/mfem/discretization/dof_selection.hpp serif/backend/mfem/discretization/quadrature/concepts.hpp serif/backend/mfem/discretization/quadrature/configure.hpp serif/backend/mfem/discretization/quadrature/mfem_resolver.hpp serif/backend/mfem/discretization/quadrature/mfem_static_order.hpp serif/backend/mfem/discretization/quadrature/order.hpp serif/backend/mfem/discretization/reconstruction.hpp serif/backend/mfem/discretization/tensor_basis.hpp serif/backend/mfem/discretization/term_binding.hpp serif/backend/mfem/discretization/term_coefficients.hpp serif/backend/mfem/geometry/deformation/computations.hpp serif/backend/mfem/geometry/deformation/constrained_radial_extension.hpp serif/backend/mfem/geometry/deformation/inputs.hpp serif/backend/mfem/geometry/deformation/radial_interpolation.hpp serif/backend/mfem/geometry/deformation/reference_geometry.hpp serif/backend/mfem/geometry/deformation/surface_metric.hpp serif/backend/mfem/geometry/surface/boundary_computations.hpp serif/backend/mfem/geometry/surface/boundary_data.hpp serif/backend/mfem/geometry/surface/boundary_inputs.hpp serif/backend/mfem/geometry/surface/boundary_rows.hpp serif/backend/mfem/geometry/surface/pressure_boundary.hpp serif/backend/mfem/mapping/compactification/compactification.hpp serif/backend/mfem/mapping/compactification/input_result.hpp serif/backend/mfem/mapping/compactification/kelvin.hpp serif/backend/mfem/mapping/compactification/variation.hpp serif/backend/mfem/mapping/context.hpp serif/backend/mfem/mapping/data.hpp serif/backend/mfem/mapping/domain_mapper.hpp serif/backend/mfem/mapping/status.hpp serif/backend/mfem/mapping/variation.hpp serif/backend/mfem/mapping/volume_cache.hpp serif/backend/mfem/model/binding.hpp serif/backend/mfem/model/computations/physical_blocks.hpp serif/backend/mfem/model/computations/solver_layout.hpp src/include/serif/backend/mfem/model/equilibrium_evaluator.hpp serif/backend/mfem/model/failure.hpp serif/backend/mfem/model/inputs.hpp serif/backend/mfem/model/normalized_system.hpp serif/backend/mfem/model/preconditioned_system.hpp serif/backend/mfem/model/preconditioning/computations.hpp serif/backend/mfem/model/state.hpp serif/backend/mfem/model/system.hpp serif/backend/mfem/operators/actions.hpp serif/backend/mfem/operators/element_data.hpp serif/backend/mfem/operators/linearized_block.hpp serif/backend/mfem/operators/quadrature.hpp serif/backend/mfem/operators/rejection.hpp serif/backend/mfem/operators/volume_data.hpp serif/backend/mfem/operators/volume_inputs.hpp serif/backend/mfem/operators/workspace.hpp serif/backend/mfem/shared/finite.hpp serif/backend/mfem/shared/reference_tracking.hpp serif/backend/mfem/solver/normalization/operators.hpp serif/backend/mfem/solver/preconditioner/border.hpp serif/backend/mfem/solver/preconditioner/computations/coupling.hpp serif/backend/mfem/solver/preconditioner/computations/inverse.hpp serif/backend/mfem/solver/preconditioner/computations/ordering.hpp serif/backend/mfem/solver/preconditioner/computations/physical_blocks.hpp serif/backend/mfem/solver/preconditioner/coupling.hpp serif/backend/mfem/solver/preconditioner/graph.hpp serif/backend/mfem/solver/preconditioner/inverse.hpp serif/backend/mfem/solver/preconditioner/inverse/amg.hpp serif/backend/mfem/solver/preconditioner/inverse/dense.hpp serif/backend/mfem/solver/preconditioner/inverse/diagonal.hpp serif/backend/mfem/solver/preconditioner/inverse/sparse_direct.hpp serif/backend/mfem/solver/preconditioner/solver.hpp serif/dimensions/dimensions.hpp serif/dimensions/partial.hpp serif/dimensions/quantities.hpp serif/dimensions/type_alias.hpp serif/discretization/blocks/base.hpp serif/discretization/blocks/concepts.hpp serif/discretization/blocks/fields.hpp serif/discretization/blocks/keys.hpp serif/discretization/blocks/runtime.hpp serif/discretization/blocks/system_layout.hpp serif/discretization/blocks/traits.hpp serif/discretization/blocks/vector_view.hpp serif/discretization/domain/concepts.hpp serif/discretization/domain/ids/boundary.hpp serif/discretization/domain/ids/concepts.hpp serif/discretization/domain/ids/domain.hpp serif/discretization/domain/ids/lists/concepts.hpp serif/discretization/domain/ids/lists/lists.hpp serif/discretization/domain/ids/lists/traits.hpp serif/discretization/domain/ids/validation/concepts.hpp serif/discretization/domain/ids/validation/unique_ids.hpp serif/discretization/domain/physical_domains.hpp serif/discretization/domain/relation/concepts.hpp serif/discretization/domain/relation/lists/concepts.hpp serif/discretization/domain/relation/lists/relation_list.hpp serif/discretization/domain/relation/lists/traits.hpp serif/discretization/domain/relation/relations.hpp serif/discretization/domain/relation/validation/concepts.hpp serif/discretization/domain/relation/validation/validation.hpp serif/discretization/domain/resolver.hpp serif/discretization/domain/schema/concepts.hpp serif/discretization/domain/schema/domain_schema.hpp serif/discretization/domain/schema/schemas.hpp serif/discretization/domain/schema/traits.hpp serif/discretization/domain/types.hpp serif/discretization/forms/concept.hpp serif/discretization/forms/form_traits.hpp serif/discretization/forms/forms.hpp serif/discretization/forms/operands.hpp serif/discretization/forms/runtime.hpp serif/discretization/forms/value.hpp serif/discretization/functional_space/space.hpp serif/discretization/quadrature/policy.hpp serif/discretization/quadrature/query.hpp serif/discretization/quadrature/rules.hpp serif/discretization/reconstruction/surface_power.hpp serif/eos/evaluation.hpp serif/eos/models/concepts.hpp serif/eos/models/traits.hpp serif/eos/relations/catalog.hpp serif/eos/relations/concepts.hpp serif/eos/relations/relations.hpp serif/eos/relations/traits.hpp serif/equations/companion.hpp serif/equations/concepts.hpp serif/equations/contributions.hpp serif/equations/point.hpp serif/equations/primal_gravity.hpp serif/equations/registry.hpp serif/equations/rigid_rotation.hpp serif/equations/surface_potential.hpp serif/geometry/deformation/actions.hpp serif/geometry/deformation/composition.hpp serif/geometry/deformation/computations.hpp serif/geometry/deformation/detail/action_storage.hpp serif/geometry/deformation/failure.hpp serif/geometry/deformation/linearization.hpp serif/geometry/deformation/models.hpp serif/geometry/deformation/radial_extension.hpp serif/geometry/deformation/radial_interpolation.hpp serif/geometry/deformation/reference_surface.hpp serif/geometry/deformation/spaces.hpp serif/geometry/deformation/surface_embedding.hpp serif/geometry/deformation/workspace.hpp serif/geometry/surface/automatic.hpp serif/geometry/surface/compilation.hpp serif/geometry/surface/computations.hpp serif/geometry/surface/detail/automatic_relations.hpp serif/geometry/surface/detail/batches.hpp serif/geometry/surface/detail/evaluation.hpp serif/geometry/surface/detail/relations.hpp serif/geometry/surface/isobaric.hpp serif/geometry/surface/linearization.hpp serif/geometry/surface/samples.hpp serif/meta/enums.hpp serif/meta/names.hpp serif/model/detail/surface_plan.hpp serif/model/equilibrium_equations.hpp serif/model/equilibrium_plan.hpp serif/model/keys.hpp serif/model/plan.hpp serif/model/state.hpp serif/operators/concepts.hpp serif/operators/linearized_system_view.hpp serif/operators/plan.hpp serif/solver/normalization/diagonal.hpp serif/solver/normalization/scales.hpp serif/solver/preconditioner/contract.hpp serif/solver/preconditioner/failure.hpp serif/solver/preconditioner/keys.hpp serif/solver/preconditioner/recipes/block.hpp serif/solver/preconditioner/recipes/bordered.hpp serif/solver/preconditioner/recipes/coupled.hpp serif/solver/preconditioner/strategies/coupling/order.hpp serif/solver/preconditioner/strategies/coupling/policies.hpp serif/solver/preconditioner/strategies/inverse/dense.hpp serif/solver/preconditioner/strategies/inverse/diagonal.hpp serif/solver/preconditioner/strategies/inverse/fixed_amg.hpp serif/solver/preconditioner/strategies/inverse/sparse_direct.hpp serif/utils/checks/finite.hpp serif/utils/dependencies/binding.hpp serif/utils/dependencies/computation.hpp serif/utils/dependencies/computation_cache.hpp serif/utils/dependencies/computation_graph.hpp serif/utils/dependencies/dependency_context.hpp serif/utils/dependencies/detail/cache_storage.hpp serif/utils/dependencies/detail/cache_traits.hpp serif/utils/dependencies/detail/graph.hpp serif/utils/dependencies/detail/result.hpp serif/utils/dependencies/detail/schema.hpp serif/utils/dependencies/detail/type_lists.hpp serif/utils/dependencies/key_list.hpp serif/utils/dependencies/key_list_operations.hpp serif/utils/dependencies/observers.hpp serif/utils/dependencies/preparation.hpp serif/utils/dependencies/resource.hpp serif/utils/dependencies/snapshot.hpp serif/utils/error/codes.hpp serif/utils/error/errors.hpp serif/utils/error/nonfinite.hpp serif/utils/format/strings.hpp serif/utils/misc/concepts/concat.hpp serif/utils/misc/concepts/exceptions.hpp serif/utils/misc/concepts/filter.hpp serif/utils/misc/concepts/numeric.hpp serif/utils/misc/concepts/type_uniqueness.hpp serif/utils/misc/std_helper/cleaning.hpp serif/utils/misc/std_helper/variadic_cleaning.hpp serif/utils/misc/terminal/colors.hpp serif/utils/misc/traits/type_uniqueness.hpp serif/utils/types/strings.hpp serif/utils/types/type_list.hpp serif/utils/user_options/mapping/compactification/kelvin.hpp serif/utils/user_options/mapping/domain_mapper.hpp algorithm array cmath compare concepts cstddef cstdint cstdio cstdlib exception expected format functional initializer_list iostream limits magic_enum/magic_enum_all.hpp memory mfem/fem/eltrans.hpp mfem/fem/fe/fe_base.hpp mfem/fem/fe/fe_h1.hpp mfem/fem/fe/fe_l2.hpp mfem/fem/fespace.hpp mfem/fem/geom.hpp mfem/fem/gridfunc.hpp mfem/fem/intrules.hpp mfem/fem/pfespace.hpp mfem/linalg/densemat.hpp mfem/linalg/operator.hpp mfem/linalg/sparsemat.hpp mfem/linalg/vector.hpp mfem/mesh/mesh.hpp nameof.hpp numeric optional ranges span stdexcept string string_view tuple type_traits unistd.h unordered_map utility variant vector version src/compile_time_checks/serif/model/state.cpp src/include/serif/backend/mfem/model/equilibrium.hpp src/include/serif/backend/mfem/model/equilibrium_evaluator.hpp src/include/serif/backend/mfem/solver/equilibrium.hpp src/include/serif/backend/mfem/solver/nonlinear/stellar_problem.hpp src/include/serif/solver/continuation/equilibrium.hpp src/include/serif/solver/equilibrium.hpp src/include/serif/solver/solver.hpp serif::model::backend::mfem::EquilibriumEvaluator serif serif::model serif::model::backend serif::model::backend::mfem serif::model::backend::mfem::detail #pragmaonce #include<expected> #include<span> #include"serif/backend/mfem/model/binding.hpp" #include"serif/backend/mfem/model/failure.hpp" #include"serif/backend/mfem/model/state.hpp" #include"serif/backend/mfem/model/preconditioning/computations.hpp" #include"serif/backend/mfem/shared/finite.hpp" #include"serif/utils/dependencies/computation_cache.hpp" namespaceserif::model::backend::mfem{ namespacedetail{ template<typenameBound,typenameCoefficients,typenameGeometry> [[nodiscard]]autocandidate_inputs(constBound&bound,constCoefficients&coefficients,constGeometry&geometry){ ifconstexpr(Bound::Plan::has_shape){ returnequilibrium_inputs(bound,coefficients,std::span<constdouble>{geometry}); }else{ returnequilibrium_inputs(bound,coefficients,geometry); } } } template<typenameBound,typenameGraph,typenameInputs,typenameAdditionalInputs=utils::dependencies::DependencyContext<>, typenameResultKey=keys::PhysicalSystem, utils::dependencies::IsKeyListEvaluationServices=utils::dependencies::KeyList<>, utils::dependencies::IsKeyListLinearServices=utils::dependencies::KeyList<>, typenameObservers=utils::dependencies::Observers<>> classEquilibriumEvaluatorfinal{ public: usingPlan=typenameBound::Plan; usingLayout=typenamePlan::Layout; usingSurfaceFormulation=typenamePlan::Surface::Formulation; usingCache=utils::dependencies::ComputationCache<Graph, Inputs, EquilibriumFailure, Observers>; usingResult=std::remove_cvref_t<decltype(std::declval<constCache&>().templateget<ResultKey>())>; usingEvaluation=std::expected<Result,typenameCache::Failure>; template<typenameRequestedKey> usingEvaluationFor=std::expected<std::remove_cvref_t<decltype(std::declval<constCache&>().templateget<RequestedKey>())>,typenameCache::Failure>; EquilibriumEvaluator(std::shared_ptr<constBound>bound,Graphgraph,Layoutlayout,AdditionalInputsadditional,Observersobservers=Observers{}): m_bound(std::move(bound)),m_cache(std::move(graph),std::move(observers)),m_layout(std::move(layout)),m_additional(std::move(additional)){ } [[nodiscard]]constLayout&layout()constnoexcept{ returnm_layout; } [[nodiscard]]constCache&cache()constnoexcept{ returnm_cache; } [[nodiscard]]constauto&model()constnoexcept{ return*m_bound->model; } [[nodiscard]]constBound&binding()constnoexcept{ return*m_bound; } private: template<typenameRequestedKey,typenameCoefficients,typenameGeometry> requires(detail::coefficient_inputs<Coefficients>(typenamePlan::Coefficients{})) [[nodiscard]]autoevaluate_result(constCoefficients&coefficients,constGeometry&geometry, utils::dependencies::PreparationModemode=utils::dependencies::PreparationMode::reuse_unchanged)->EvaluationFor<RequestedKey>{ m_cache.clear_preparation(); constboolfinite=[&]<typename...Terms>(utils::dependencies::KeyList<Terms...>){ return(utils::checks::is_finite(coefficients.templateget<Terms>().values)&&...); }(typenamePlan::Coefficients{}); if(!finite){ returnstd::unexpected(typenameCache::Failure{.computation_index=Graph::size, .error=EquilibriumFailure{operators::backend::mfem::PreparationRejection{ .reason=operators::backend::mfem::PreparationRejectionReason::NON_FINITE, .message="Astellarcandidatecontainsanon-finitephysicalcoefficient.Supplyfinitefieldandscalarvaluesbeforeevaluatingit." }}}); } constautoinputs=utils::dependencies::combine_dependency_contexts(detail::candidate_inputs(*m_bound,coefficients,geometry),m_additional); usingRequested=utils::dependencies::CombinedKeyList<utils::dependencies::KeyList<RequestedKey>,EvaluationServices, std::conditional_t<std::same_as<RequestedKey,preconditioning::PreparedSystem>,LinearServices,utils::dependencies::KeyList<>>>; constautopreparation=m_cache.try_prepare(inputs,Requested{},mode); if(!preparation){ returnstd::unexpected(preparation.error()); } returnm_cache.templateget<RequestedKey>(); } template<typenameRequestedKey,utils::dependencies::IsDependencyContextPrescribed=utils::dependencies::DependencyContext<>> requiresPlan::has_shape&&detail::PrescribedInputsFor<Plan, Prescribed> [[nodiscard]]autoevaluate_result(constStellarState<Layout>&state,constPrescribed&prescribed=Prescribed{}, utils::dependencies::PreparationModemode=utils::dependencies::PreparationMode::reuse_unchanged)->EvaluationFor<RequestedKey>{ validate_state(state); constautofields=detail::state_coefficients(*m_bound->discretization,state,prescribed); constautocoefficients=fields.context(); returnevaluate_result<RequestedKey>(coefficients,state.templateget<discretization::blocks::surface_deformation::coefficients>(),mode); } template<typenameRequestedKey,utils::dependencies::IsDependencyContextPrescribed=utils::dependencies::DependencyContext<>> requires(!Plan::has_shape)&&detail::PrescribedInputsFor<Plan, Prescribed> [[nodiscard]]autoevaluate_result(constStellarState<Layout>&state,const::mfem::Vector&displacement,constPrescribed&prescribed=Prescribed{}, utils::dependencies::PreparationModemode=utils::dependencies::PreparationMode::reuse_unchanged)->EvaluationFor<RequestedKey>{ validate_state(state); constautofields=detail::state_coefficients(*m_bound->discretization,state,prescribed); constautocoefficients=fields.context(); returnevaluate_result<RequestedKey>(coefficients,displacement,mode); } voidvalidate_state(constStellarState<Layout>&state){ usingnamespaceutils::errors; m_cache.clear_preparation(); serif_ensure<SERiFErrorCode::INVALID_DOF_COUNT>(state.layout()==m_layout, "Thestellarstate'sfieldorequationsizesdifferfromthisevaluator'slayout.Constructthestatewiththematchingdiscretizationlayout."); } public: template<typename...InputsT> requiresrequires(EquilibriumEvaluator&evaluator,constInputsT&...inputs){ evaluator.templateevaluate_result<ResultKey>(inputs...); } [[nodiscard]]Evaluationevaluate(constInputsT&...inputs){ returnevaluate_result<ResultKey>(inputs...); } template<typename...InputsT> requiresutils::dependencies::ContainsKey<preconditioning::PreparedSystem, typename Graph::Results>&& requires(EquilibriumEvaluator&evaluator,constInputsT&...inputs){ evaluator.templateevaluate_result<preconditioning::PreparedSystem>(inputs...); } [[nodiscard]]autoprepare_linear_system(constInputsT&...inputs)->decltype(evaluate_result<preconditioning::PreparedSystem>(inputs...)){ returnevaluate_result<preconditioning::PreparedSystem>(inputs...); } private: std::shared_ptr<constBound>m_bound; Cachem_cache; Layoutm_layout; AdditionalInputsm_additional; }; }