equilibrium.hpp expected serif/model/backend/mfem/binding.hpp serif/model/backend/mfem/configured_discretization.hpp serif/model/normalization/computations.hpp serif/model/backend/mfem/discretization.hpp serif/model/backend/mfem/failure.hpp serif/model/backend/mfem/problem.hpp serif/model/graph.hpp serif/model/execution_validation.hpp serif/model/stellar_model.hpp serif/utils/checks/finite.hpp serif/dimensions/dimensions.hpp serif/dimensions/partial.hpp serif/dimensions/quantities.hpp serif/dimensions/type_alias.hpp serif/discretization/backend/mfem/basis_evaluation.hpp serif/discretization/backend/mfem/boundary_trace.hpp serif/discretization/backend/mfem/coefficient_tracking.hpp serif/discretization/backend/mfem/dof_map.hpp serif/discretization/backend/mfem/dof_selection.hpp serif/discretization/backend/mfem/mass_diagonal.hpp serif/discretization/backend/mfem/point_dofs.hpp serif/discretization/backend/mfem/tensor_basis.hpp serif/discretization/backend/mfem/term_binding.hpp serif/discretization/blocks/base.hpp serif/discretization/blocks/concepts.hpp serif/discretization/blocks/fields.hpp serif/discretization/blocks/runtime.hpp serif/discretization/blocks/system_layout.hpp serif/discretization/blocks/traits.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/backend/mfem/concepts.hpp serif/discretization/quadrature/backend/mfem/configure.hpp serif/discretization/quadrature/backend/mfem/mfem_resolver.hpp serif/discretization/quadrature/backend/mfem/mfem_static_order.hpp serif/discretization/quadrature/backend/order.hpp serif/discretization/quadrature/policy.hpp serif/discretization/quadrature/query.hpp serif/discretization/quadrature/rules.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/concepts.hpp serif/equations/contributions.hpp serif/equations/point.hpp serif/equations/registry.hpp serif/geometry/deformation/actions.hpp serif/geometry/deformation/backend/mfem/computations.hpp serif/geometry/deformation/backend/mfem/detail/reference_tracking.hpp serif/geometry/deformation/backend/mfem/inputs.hpp serif/geometry/deformation/backend/mfem/radial_interpolation.hpp serif/geometry/deformation/backend/mfem/reference_geometry.hpp serif/geometry/deformation/backend/mfem/surface_metric.hpp serif/geometry/deformation/backend/mfem/volume_displacement.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/backend/mfem/automatic.hpp serif/geometry/surface/backend/mfem/boundary_computations.hpp serif/geometry/surface/backend/mfem/boundary_data.hpp serif/geometry/surface/backend/mfem/boundary_inputs.hpp serif/geometry/surface/backend/mfem/boundary_rows.hpp serif/geometry/surface/backend/mfem/detail/automatic_relations.hpp serif/geometry/surface/backend/mfem/pressure_boundary.hpp serif/geometry/surface/compilation.hpp serif/geometry/surface/computations.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/mapping/compactification/compactification.hpp serif/mapping/compactification/input_result.hpp serif/mapping/compactification/kelvin.hpp serif/mapping/compactification/variation.hpp serif/mapping/context.hpp serif/mapping/data.hpp serif/mapping/domain_mapper.hpp serif/mapping/status.hpp serif/mapping/transformations.hpp serif/mapping/variation.hpp serif/mapping/volume_cache.hpp serif/meta/enums.hpp serif/meta/names.hpp serif/model/backend/mfem/actions.hpp serif/model/backend/mfem/binding.hpp serif/model/backend/mfem/configured_discretization.hpp serif/model/backend/mfem/deformation.hpp serif/model/backend/mfem/discretization.hpp src/include/serif/model/backend/mfem/equilibrium.hpp serif/model/backend/mfem/failure.hpp serif/model/backend/mfem/inputs.hpp serif/model/backend/mfem/normalized_system.hpp serif/model/backend/mfem/problem.hpp serif/model/backend/mfem/system.hpp serif/model/backend/mfem/volume.hpp serif/model/backend/mfem/volume_actions.hpp serif/model/computations/interior.hpp serif/model/computations/point_constraints.hpp serif/model/computations/pressure_boundary.hpp serif/model/computations/shape.hpp serif/model/computations/shape_derivatives.hpp serif/model/description.hpp serif/model/detail/specification_endpoints.hpp serif/model/detail/specification_list.hpp serif/model/detail/surface_plan.hpp serif/model/equilibrium_equations.hpp serif/model/equilibrium_plan.hpp serif/model/execution_validation.hpp serif/model/graph.hpp serif/model/keys.hpp serif/model/normalization/computations.hpp serif/model/normalization/plan.hpp serif/model/normalization/reference.hpp serif/model/plan.hpp serif/model/specification.hpp serif/model/specifications/fixed_angular_momentum.hpp serif/model/specifications/fixed_central_density.hpp serif/model/specifications/fixed_total_mass.hpp serif/model/specifications/physical_names.hpp serif/model/stellar_model.hpp serif/model/validation.hpp serif/operators/backend/mfem/element_coordinates.hpp serif/operators/backend/mfem/element_data.hpp serif/operators/backend/mfem/element_evaluation.hpp serif/operators/backend/mfem/element_preparation.hpp serif/operators/backend/mfem/element_transpose.hpp serif/operators/backend/mfem/linearized_block.hpp serif/operators/backend/mfem/nodal_computations.hpp serif/operators/backend/mfem/nodal_data.hpp serif/operators/backend/mfem/nodal_evaluation.hpp serif/operators/backend/mfem/nodal_inputs.hpp serif/operators/backend/mfem/nodal_rows.hpp serif/operators/backend/mfem/operand_evaluation.hpp serif/operators/backend/mfem/point_border.hpp serif/operators/backend/mfem/preparation.hpp serif/operators/backend/mfem/pressure_field_block.hpp serif/operators/backend/mfem/quadrature.hpp serif/operators/backend/mfem/rejection.hpp serif/operators/backend/mfem/replaced_rows.hpp serif/operators/backend/mfem/term_coefficients.hpp serif/operators/backend/mfem/volume_application.hpp serif/operators/backend/mfem/volume_computations.hpp serif/operators/backend/mfem/volume_data.hpp serif/operators/backend/mfem/volume_field_block.hpp serif/operators/backend/mfem/volume_inputs.hpp serif/operators/backend/mfem/volume_stages.hpp serif/operators/backend/mfem/workspace.hpp serif/operators/concepts.hpp serif/operators/contributions.hpp serif/operators/linearized_system_view.hpp serif/operators/plan.hpp serif/solver/normalization/backend/mfem/metric_computations.hpp serif/solver/normalization/backend/mfem/operators.hpp serif/solver/normalization/computations.hpp serif/solver/normalization/configuration.hpp serif/solver/normalization/detail/factor_computations.hpp serif/solver/normalization/detail/packing.hpp serif/solver/normalization/detail/validation.hpp serif/solver/normalization/diagonal.hpp serif/solver/normalization/dimensionless.hpp serif/solver/normalization/keys.hpp serif/solver/normalization/layout.hpp serif/solver/normalization/metrics.hpp serif/solver/normalization/physical_coordinates.hpp serif/solver/normalization/plan.hpp serif/solver/normalization/point_rows.hpp serif/solver/normalization/scales.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/computation_store.hpp serif/utils/dependencies/dependencies.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/incremental_computation.hpp serif/utils/dependencies/key_list.hpp serif/utils/dependencies/key_list_operations.hpp serif/utils/dependencies/keys.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/format/strings.hpp serif/utils/laneemden/laneemden.hpp serif/utils/laneemden/scaling.hpp serif/utils/laneemden/solution.hpp serif/utils/misc/compensated_sum.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 map memory mfem/fem/doftrans.hpp mfem/fem/eltrans.hpp mfem/fem/fe/fe_base.hpp mfem/fem/fe/fe_h1.hpp mfem/fem/fe/fe_l2.hpp mfem/fem/fe/fe_rt.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 span stdexcept string string_view tuple type_traits unistd.h unordered_map utility variant vector version serif::model::backend::mfem::StellarEquilibrium serif serif::model serif::model::backend serif::model::backend::mfem #pragmaonce #include<expected> #include"serif/model/backend/mfem/binding.hpp" #include"serif/model/backend/mfem/configured_discretization.hpp" #include"serif/model/normalization/computations.hpp" #include"serif/model/backend/mfem/discretization.hpp" #include"serif/model/backend/mfem/failure.hpp" #include"serif/model/backend/mfem/problem.hpp" #include"serif/model/graph.hpp" #include"serif/model/execution_validation.hpp" #include"serif/model/stellar_model.hpp" #include"serif/utils/checks/finite.hpp" namespaceserif::model::backend::mfem{ template<typenameBound,typenameProblem,typenameAdditionalInputs=utils::dependencies::DependencyContext<>,typenamePreparedKey=keys::PhysicalSystem> classStellarEquilibriumfinal{ public: usingPlan=typenameBound::Plan; usingSurfaceFormulation=typenamePlan::Surface::Formulation; usingPrepared=std::remove_cvref_t<decltype(std::declval<consttypenameProblem::Cache&>().templateget<PreparedKey>())>; StellarEquilibrium(std::shared_ptr<constBound>bound,Problemproblem,AdditionalInputsadditional=AdditionalInputs{}): m_bound(std::move(bound)),m_problem(std::move(problem)),m_additional(std::move(additional)){ } template<typenameCoefficients,typenameGeometry> [[nodiscard]]autotry_prepare(constCoefficients&coefficients,constGeometry&geometry, utils::dependencies::PreparationModemode=utils::dependencies::PreparationMode::reuse_unchanged){ constautoresult=m_problem.try_prepare_inputs([&]()->std::expected<typenameProblem::InputContext,EquilibriumFailure>{ constboolfinite=[&]<typename...Terms>(utils::dependencies::KeyList<Terms...>){ return(utils::checks::is_finite(coefficients.templateget<Terms>().values)&&...); }(typenamePlan::Coefficients{}); if(!finite){ returnstd::unexpected(EquilibriumFailure{operators::backend::mfem::PreparationRejection{ .reason=operators::backend::mfem::PreparationRejectionReason::NON_FINITE, .message="Astellarcandidatecontainsanon-finitephysicalcoefficient.Supplyfinitefieldandscalarvaluesbeforepreparingit." }}); } returnutils::dependencies::combine_dependency_contexts(equilibrium_inputs(*m_bound,coefficients,geometry),m_additional); },utils::dependencies::KeyList<PreparedKey>{},mode); if(result){ m_prepared=m_problem.cache().templateget<PreparedKey>(); } returnresult; } [[nodiscard]]conststd::shared_ptr<constAcceptedSystem>&accepted()const{ returnm_problem.accepted(); } [[nodiscard]]constPrepared&prepared()const{ usingnamespaceutils::errors; serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(m_prepared!=nullptr,"Thestellarequilibriumhasnocompletepreparation.Successfullyprepareacandidatebeforerequestingsolvercoordinates."); returnm_prepared; } [[nodiscard]]constauto&layout()constnoexcept{ returnm_problem.layout(); } [[nodiscard]]constauto&cache()constnoexcept{ returnm_problem.cache(); } [[nodiscard]]constauto&model()constnoexcept{ return*m_bound->model; } private: std::shared_ptr<constBound>m_bound; Problemm_problem; AdditionalInputsm_additional; Preparedm_prepared{}; }; template<utils::dependencies::IsKeyListPrescribed=utils::dependencies::KeyList<>, utils::dependencies::IsBindingsSpecificationMap=DefaultSpecificationBindings, utils::dependencies::IsKeyListExplicitZeroCouplings=utils::dependencies::KeyList<>,StellarModelTypeModel, typenameDiscretization,typenameCoefficients,typenameGeometry> [[nodiscard]]autobind_stellar_model(Modelmodel,Discretizationdiscretization,constCoefficients&coefficients, constGeometry&geometry,operators::backend::mfem::VolumeOptionsoptions={}){ constautomodel_owner=std::make_shared<constModel>(std::move(model)); constautodiscretization_owner=std::make_shared<constDiscretization>(std::move(discretization)); constautobound=bind_model_definition<Prescribed>(model_owner,discretization_owner,coefficients); usingBound=typenamedecltype(bound)::element_type; usingPlan=typenameBound::Plan; static_assert(executable_model_topology<Model, typename Bound::Description, Plan, SpecificationMap, ExplicitZeroCouplings>(), "Theselectedmodelbackendneedsthedeclaredprescriptioncoordinates,constraintequations,andexplicitgenerated-coordinateeffects.Supplytheirequationsorauthor-registerknownzeroeffects;thischeckdoesnotprovemathematicalclosure."); constautolayout=equilibrium_layout<Plan>(*discretization_owner,geometry); autograph=make_equilibrium_graph(bound,layout,std::tuple{},options); constautoinputs=equilibrium_inputs(*bound,coefficients,geometry); autoproblem=make_stellar_problem<EquilibriumFailure>(std::move(graph),inputs,layout); returnStellarEquilibrium<std::remove_const_t<Bound>,decltype(problem)>(bound,std::move(problem)); } template<utils::dependencies::IsKeyListPrescribed=utils::dependencies::KeyList<>, utils::dependencies::IsBindingsSpecificationMap=DefaultSpecificationBindings, utils::dependencies::IsKeyListExplicitZeroCouplings=utils::dependencies::KeyList<>,StellarModelTypeModel, typenameDiscretization,typenameOptions,typenameCoefficients,typenameGeometry> [[nodiscard]]autobind_stellar_model(ConfiguredDiscretization<Model,Discretization,Options>configuration, constCoefficients&coefficients,constGeometry&geometry,operators::backend::mfem::VolumeOptionsoptions={}){ constautomodel_owner=std::make_shared<constModel>(std::move(configuration.model)); constautodiscretization_owner=std::make_shared<constDiscretization>(std::move(configuration.finite_elements)); constautobound=bind_model_definition<Prescribed>(model_owner,discretization_owner,coefficients); usingBound=typenamedecltype(bound)::element_type; usingPlan=typenameBound::Plan; usingNormalizationPlan=normalization::PlanFor<decltype(configuration.normalization.plan),decltype(configuration.normalization.coordinates),Bound>; static_assert(executable_model_topology<Model, typename Bound::Description, Plan, SpecificationMap, ExplicitZeroCouplings>(), "Theconfiguredmodelneedscompleteequationandprescriptionbindings,includingexplicitlyknownzeroeffects."); constautolayout=equilibrium_layout<Plan>(*discretization_owner,geometry); static_assert(solver::normalization::NormalizationPreparable<NormalizationPlan,decltype(layout)>, "Thenormalizationplanmustprovideonecompatibleruleforeveryunknownandresidualrow."); constautonumerical_options=std::make_shared<constOptions>(std::move(configuration.normalization)); constautometrics=normalization::MetricsFor<NormalizationPlan>{}; constautoadditional=normalization::make_inputs(numerical_options,*discretization_owner,metrics); constautonumerical=normalization::make_computations<NormalizationPlan>(*bound,*numerical_options,layout,metrics); autograph=make_equilibrium_graph(bound,layout,numerical,options); constautoinputs=utils::dependencies::combine_dependency_contexts(equilibrium_inputs(*bound,coefficients,geometry),additional); autoproblem=make_stellar_problem<EquilibriumFailure>(std::move(graph),inputs,layout); returnStellarEquilibrium<std::remove_const_t<Bound>,decltype(problem),std::remove_cvref_t<decltype(additional)>,normalization::PreparedSystem>( bound,std::move(problem),additional); } }