state.hpp array limits serif/model/plan.hpp serif/model/state.hpp serif/backend/mfem/discretization/term_coefficients.hpp src/include/serif/backend/mfem/model/equilibrium_evaluator.hpp serif/backend/mfem/discretization/dof_map.hpp serif/backend/mfem/discretization/reconstruction.hpp serif/backend/mfem/discretization/term_binding.hpp serif/backend/mfem/discretization/term_coefficients.hpp src/include/serif/backend/mfem/model/state.hpp serif/dimensions/dimensions.hpp serif/dimensions/quantities.hpp serif/discretization/blocks/base.hpp serif/discretization/blocks/concepts.hpp serif/discretization/blocks/fields.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/traits.hpp serif/discretization/domain/types.hpp serif/discretization/forms/forms.hpp serif/discretization/forms/operands.hpp serif/discretization/forms/value.hpp serif/discretization/functional_space/space.hpp serif/discretization/reconstruction/surface_power.hpp serif/equations/concepts.hpp serif/equations/contributions.hpp serif/equations/point.hpp serif/equations/registry.hpp serif/meta/enums.hpp serif/meta/names.hpp serif/model/plan.hpp serif/model/state.hpp serif/operators/concepts.hpp serif/operators/plan.hpp serif/utils/dependencies/binding.hpp serif/utils/dependencies/dependency_context.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/misc/concepts/concat.hpp serif/utils/misc/concepts/exceptions.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 algorithm array cmath compare concepts cstddef cstdint cstdio cstdlib exception format functional iostream limits magic_enum/magic_enum_all.hpp memory mfem/fem/fe/fe_base.hpp mfem/fem/fespace.hpp mfem/fem/pfespace.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/state.cpp src/include/serif/backend/mfem/model/equilibrium.hpp src/include/serif/backend/mfem/model/equilibrium_evaluator.hpp src/include/serif/backend/mfem/model/state.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::detail::StateCoefficients< Discretization, Layout, Prescribed, utils::dependencies::KeyList< Fields... > > serif serif::model serif::model::backend serif::model::backend::mfem serif::model::backend::mfem::detail #pragmaonce #include<array> #include<limits> #include"serif/model/plan.hpp" #include"serif/model/state.hpp" #include"serif/backend/mfem/discretization/term_coefficients.hpp" namespaceserif::model::backend::mfem::detail{ template<typenameInputs,typenameTerm> conceptCoefficientInput=requires(constInputs&inputs){ typenameInputs::Keys; {inputs.templateget<Term>().values}->std::convertible_to<const::mfem::Vector&>; requiresstd::same_as<typenamestd::remove_cvref_t<decltype(inputs.templateget<Term>().binding)>::Term,Term>; }; template<typenameInputs,typename...Terms> [[nodiscard]]constevalboolcoefficient_inputs(utils::dependencies::KeyList<Terms...>){ return(CoefficientInput<Inputs, Terms>&&...); } template<typenamePlan,typenameInputs> conceptPrescribedInputsFor=utils::dependencies::IsDependencyContext<Inputs>&& (utils::dependencies::key_count_v<utils::dependencies::CombinedKeyList<typename Plan::Values, typename Inputs::Keys>>== utils::dependencies::key_count_v<typename Plan::Values>+utils::dependencies::key_count_v<typename Inputs::Keys>)&& coefficient_inputs<Inputs>(typenamemodel::detail::RemoveKeys<typename Plan::Coefficients, typename Plan::Values>::Type{}); template<typenameDiscretization,typenameLayout,typenamePrescribed,typenameFields> classStateCoefficients; template<typenameDiscretization,typenameLayout,typenamePrescribed,typename...Fields> classStateCoefficients<Discretization,Layout,Prescribed,utils::dependencies::KeyList<Fields...>>final{ public: usingKeys=utils::dependencies::CombinedKeyList<utils::dependencies::KeyList<Fields...>,typenamePrescribed::Keys>; static_assert(utils::dependencies::key_count_v<Keys>==sizeof...(Fields)+utils::dependencies::key_count_v<typename Prescribed::Keys>, "Aphysicalfieldcannotbebothsolvedandprescribed.Removeitsprescribedinputoritssolvedcoordinate."); StateCoefficients(constDiscretization&discretization,constStellarState<Layout>&state,constPrescribed&prescribed): m_discretization(discretization),m_prescribed(prescribed),m_fields{field_view(state.templateget<Fields>())...}{ } [[nodiscard]]autocontext()const&{ usingnamespaceutils::dependencies; constautofields=make_dependency_context(own<Fields>(operators::backend::mfem::bind( m_discretization.templatebinding<Fields>(),m_fields[key_index<Fields,KeyList<Fields...>>()]))...); returncombine_dependency_contexts(fields,m_prescribed); } autocontext()const&&=delete; private: [[nodiscard]]static::mfem::Vectorfield_view(std::span<constdouble>values){ usingnamespaceutils::errors; serif_ensure<SERiFErrorCode::INVALID_DOF_COUNT>(values.size()<=static_cast<std::size_t>(std::numeric_limits<int>::max()), "AstellarfieldexceedsMFEM'scoefficientindexrange.Supplyarepresentablediscretization."); return::mfem::Vector(const_cast<double*>(values.data()),static_cast<int>(values.size())); } constDiscretization&m_discretization; constPrescribed&m_prescribed; std::array<::mfem::Vector,sizeof...(Fields)>m_fields; }; template<typenameDiscretization,typenameLayout,typenamePrescribed> [[nodiscard]]autostate_coefficients(constDiscretization&discretization,constStellarState<Layout>&state,constPrescribed&prescribed){ usingFields=typenamemodel::detail::RemoveKeys<typenameLayout::Values, utils::dependencies::KeyList<discretization::blocks::surface_deformation::coefficients>>::Type; returnStateCoefficients<Discretization, Layout, Prescribed, Fields>(discretization,state,prescribed); } }