finite_elements.hpp serif/backend/mfem/discretization/resources.hpp serif/backend/mfem/discretization/term_binding.hpp serif/backend/mfem/model/discretization.hpp serif/model/equilibrium_plan.hpp serif/model/mesh_options.hpp serif/model/stellar_model.hpp src/include/serif/model/finite_elements.hpp src/lib/serif/backend/mfem/model/finite_elements.cpp serif/backend/mfem/discretization/basis_evaluation.hpp serif/backend/mfem/discretization/dof_map.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/resources.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/inputs.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/discretization.hpp src/include/serif/backend/mfem/model/finite_elements.hpp serif/backend/mfem/operators/element_data.hpp serif/backend/mfem/operators/quadrature.hpp serif/backend/mfem/operators/volume_data.hpp serif/backend/mfem/operators/workspace.hpp serif/backend/mfem/shared/reference_tracking.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/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/policy.hpp serif/discretization/quadrature/query.hpp serif/discretization/quadrature/rules.hpp serif/discretization/reconstruction/surface_power.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/spaces.hpp serif/geometry/deformation/workspace.hpp serif/geometry/surface/automatic.hpp serif/geometry/surface/compilation.hpp serif/geometry/surface/detail/automatic_relations.hpp serif/geometry/surface/detail/relations.hpp serif/geometry/surface/isobaric.hpp serif/meta/enums.hpp serif/meta/names.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/mesh_options.hpp serif/model/plan.hpp serif/model/specification.hpp serif/model/stellar_model.hpp serif/operators/concepts.hpp serif/operators/plan.hpp serif/utils/dependencies/binding.hpp serif/utils/dependencies/computation.hpp serif/utils/dependencies/computation_graph.hpp serif/utils/dependencies/dependency_context.hpp serif/utils/dependencies/detail/graph.hpp serif/utils/dependencies/detail/type_lists.hpp serif/utils/dependencies/key_list.hpp serif/utils/dependencies/key_list_operations.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/vector.hpp mfem/mesh/mesh.hpp nameof.hpp numeric optional ranges span stdexcept string string_view tuple type_traits unistd.h unordered_map utility vector version src/compile_time_checks/serif/model/finite_elements.cpp src/include/serif/backend/mfem/model/finite_elements.hpp src/include/serif/model/finite_elements.hpp src/include/serif/solver/continuation/equilibrium.hpp src/include/serif/solver/equilibrium.hpp src/include/serif/solver/solver.hpp src/lib/serif/backend/mfem/model/finite_elements.cpp serif serif::model serif::model::backend serif::model::backend::mfem serif::model::backend::mfem::detail #pragmaonce #include"serif/backend/mfem/discretization/resources.hpp" #include"serif/backend/mfem/discretization/term_binding.hpp" #include"serif/backend/mfem/model/discretization.hpp" #include"serif/model/equilibrium_plan.hpp" #include"serif/model/mesh_options.hpp" #include"serif/model/stellar_model.hpp" namespaceserif::model::backend::mfem{ [[nodiscard]]discretization::backend::mfem::Resourcesmake_stellar_resources(constStellarMeshOptions&options); [[nodiscard]]::mfem::FiniteElementSpace&make_space(discretization::backend::mfem::Resources&resources, std::unique_ptr<::mfem::FiniteElementCollection>collection,intcomponents,intordering=::mfem::Ordering::byNODES); namespacedetail{ template<typenameSpecification> conceptReconstructionSpecification=requires{typenameSpecification::Term;}; template<typenameOwner,typename...Terms,ReconstructionSpecificationSpecification> voidapply_reconstruction(Owner&owner,utils::dependencies::KeyList<Terms...>,constSpecification&specification){ usingnamespaceutils::errors; usingTerm=typenameSpecification::Term; usingKeys=utils::dependencies::KeyList<Terms...>; ifconstexpr(!(std::same_as<Term,Terms>||...)){ serif_error<std::invalid_argument,SERiFErrorCode::INVALID_ARGUMENT>("Thereconstruction'sfieldisnotpresentinthisstellarmodel."); }else{ auto&binding=std::get<utils::dependencies::key_index<Term,Keys>()>(owner.bindings); usingdiscretization::backend::mfem::make_reconstruction; ifconstexpr(requires{make_reconstruction(specification,*binding);}){ autoplan=make_reconstruction(specification,*binding); ifconstexpr(discretization::blocks::SpatialTerm<Term>){ usingBinding=typenameOwner::templateBinding<Term>; binding=std::make_unique<Binding>(binding->space(),std::move(plan)); }else{ serif_error<std::invalid_argument,SERiFErrorCode::INVALID_ARGUMENT>("Spatialreconstructioncannotbeattachedtoaglobalscalarfield."); } }else{ serif_error<std::invalid_argument,SERiFErrorCode::INVALID_ARGUMENT>("Thisreconstructionformhasnopreparationfunctionfortheselectedbackendandfield."); } } } template<ReconstructionSpecification...Specifications> voidvalidate_reconstruction_requests(){ usingnamespaceutils::errors; constexprboolunique=(([]<typenameSpecification>(){ return(0+...+int(std::same_as<typenameSpecification::Term,typenameSpecifications::Term>))==1; }.templateoperator()<Specifications>())&&...); serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(unique,"Specifyatmostonereconstructionforeachfield."); } template<discretization::blocks::SpatialTermTerm> [[nodiscard]]::mfem::FiniteElementSpace&term_space(discretization::backend::mfem::Resources&resources){ usingTraits=discretization::blocks::SpatialTermTraits<Term>; usingenumdiscretization::functional_space::Family; static_assert(Traits::rank<=1&&(Traits::family==h1||Traits::family==l2), "Theselectedbackendsupportsscalar/vectorH1andL2terms."); constexprintcomponents=Traits::rank==0?1:3; ifconstexpr(Traits::family==h1){ returnmake_space(resources,std::make_unique<::mfem::H1_FECollection>(Traits::order,3),components,components==1?::mfem::Ordering::byNODES:::mfem::Ordering::byVDIM); }else{ returnmake_space(resources,std::make_unique<::mfem::L2_FECollection>(Traits::order,3),components,components==1?::mfem::Ordering::byNODES:::mfem::Ordering::byVDIM); } } template<typename...Terms,ReconstructionSpecification...Specifications> [[nodiscard]]automake_bindings(constStellarMeshOptions&options,utils::dependencies::KeyList<Terms...>keys,constSpecifications&...specifications){ usingnamespaceutils::dependencies; usingSchema=discretization::domain::schema::CoreEnvelopeVacuumDomainSchema; usingContext=decltype(make_dependency_context(borrow<Terms>(std::declval<discretization::backend::mfem::TermBinding<Terms,Schema>&>())...)); usingOwner=discretization::backend::mfem::OwnedBindings<Context>; validate_reconstruction_requests<Specifications...>(); autoowner=std::make_shared<Owner>(); owner->resources=make_stellar_resources(options); ([&]{ usingBinding=typenameOwner::templateBinding<Terms>; auto&binding=std::get<key_index<Terms,KeyList<Terms...>>()>(owner->bindings); ifconstexpr(discretization::blocks::SpatialTerm<Terms>){ binding=std::make_unique<Binding>(term_space<Terms>(owner->resources)); }else{ binding=std::make_unique<Binding>(*owner->resources.meshes.front()); } }(),...); (apply_reconstruction(*owner,keys,specifications),...); returnowner; } } template<StellarModelTypeModel,detail::ReconstructionSpecification...Specifications> [[nodiscard]]automake_stellar_discretization(constModel&,constStellarMeshOptions&options={},constSpecifications&...specifications){ usingDescription=decltype(resolve_equilibrium_equations(std::declval<constModel&>())); usingPlan=EquilibriumPlan<typenameModel::EquationOfState,Description>; usingTerms=utils::dependencies::CombinedKeyList<typenamePlan::Coefficients, utils::dependencies::KeyList<discretization::blocks::displacement::vector>>; autoowner=detail::make_bindings(options,Terms{},specifications...); returnmake_stellar_discretization(owner,owner->context(),*owner->resources.compactification,*owner->resources.mapper); } template<StellarModelTypeModel,typenameSurface,typenameInterior,typenameVacuum,detail::ReconstructionSpecification...Specifications> [[nodiscard]]automake_stellar_discretization(constModel&model,constStellarMeshOptions&options, geometry::deformation::DomainDeformation<Surface,Interior,Vacuum>deformation,constSpecifications&...specifications){ usingDescription=decltype(resolve_equilibrium_equations(std::declval<constModel&>())); usingPlan=EquilibriumPlan<typenameModel::EquationOfState,Description>; usingTerms=utils::dependencies::CombinedKeyList<typenamePlan::Coefficients, utils::dependencies::KeyList<discretization::blocks::displacement::vector>>; autoowner=detail::make_bindings(options,Terms{},specifications...); usingDisplacement=discretization::blocks::displacement::vector; constautocontext=owner->context(); constauto*collection=context.templateget<Displacement>().space().FEColl(); auto&scalar=make_space(owner->resources,std::unique_ptr<::mfem::FiniteElementCollection>{::mfem::FiniteElementCollection::New(collection->Name())},1); auto&reference_displacement=make_space(owner->resources, std::unique_ptr<::mfem::FiniteElementCollection>{::mfem::FiniteElementCollection::New(collection->Name())},3,::mfem::Ordering::byNODES); constautoreference=geometry::deformation::backend::mfem::deformation_discretization( scalar,reference_displacement,*owner->resources.meshes[1]); returnwith_deformation(make_stellar_discretization(owner,context,*owner->resources.compactification,*owner->resources.mapper), std::move(deformation),reference); } }