binding.hpp tuple utility serif/model/computations/point_constraints.hpp serif/model/equilibrium_plan.hpp serif/model/backend/mfem/inputs.hpp serif/geometry/surface/backend/mfem/pressure_boundary.hpp serif/geometry/deformation/backend/mfem/inputs.hpp src/include/serif/model/backend/mfem/equilibrium.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/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/detail/reference_tracking.hpp serif/geometry/deformation/backend/mfem/inputs.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/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 src/include/serif/model/backend/mfem/binding.hpp serif/model/backend/mfem/deformation.hpp serif/model/backend/mfem/inputs.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/shape_derivatives.hpp serif/model/detail/surface_plan.hpp serif/model/equilibrium_equations.hpp serif/model/equilibrium_plan.hpp serif/model/plan.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/quadrature.hpp serif/operators/backend/mfem/rejection.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/plan.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/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 src/include/serif/model/backend/mfem/binding.hpp src/include/serif/model/backend/mfem/equilibrium.hpp serif::model::backend::mfem::BoundModel serif serif::model serif::model::backend serif::model::backend::mfem serif::model::backend::mfem::detail #pragmaonce #include<tuple> #include<utility> #include"serif/model/computations/point_constraints.hpp" #include"serif/model/equilibrium_plan.hpp" #include"serif/model/backend/mfem/inputs.hpp" #include"serif/geometry/surface/backend/mfem/pressure_boundary.hpp" #include"serif/geometry/deformation/backend/mfem/inputs.hpp" namespaceserif::model::backend::mfem { template<typenameModel,typenameDescriptionT,typenameDiscretization,typenameBoundary,typenameConstraintsT,typenamePlanT> structBoundModel { usingDescription=DescriptionT; usingConstraintBindings=ConstraintsT; usingPlan=PlanT; std::shared_ptr<constModel>model; Descriptiondefinition; std::shared_ptr<constDiscretization>discretization; Boundaryboundary; ConstraintBindingsconstraints; template<typenameRow> [[nodiscard]]autoequation_discretization()const { returndiscretization->templateequation_discretization<Row>(); } }; template<utils::dependencies::IsKeyListPrescribed=utils::dependencies::KeyList<>,typenameModel,typenameDiscretization,typenameCoefficients> [[nodiscard]]autobind_model_definition(std::shared_ptr<constModel>model,std::shared_ptr<constDiscretization>discretization, constCoefficients&coefficients) { usingnamespaceutils::errors; serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(model!=nullptr&&discretization!=nullptr,"Bindingastellarmodelrequiresbothanowningphysicaldescriptionandanowningdiscretization."); autodefinition=resolve_equilibrium_equations(*model); usingDescription=decltype(definition); usingAvailableFields=typenameCoefficients::Keys; usingPlan=EquilibriumPlan<typename Model::EquationOfState, Description, Prescribed, AvailableFields>; static_assert(Plan::Surface::supported, "Theisobaricconditionhasnouniquesupportedequation.ChecktheEOSrelations,availableboundaryfields,derivativesandsolvedequationrows.ForgenuinelyambiguousrowsanequationmaynameitsIsobaricSurfaceField."); detail::validate_coefficient_owners(*discretization,coefficients,typenamePlan::Coefficients{}); constautoboundary=[&]() { ifconstexpr(Plan::has_surface&&Plan::Surface::supported) { usingFormulation=typenamePlan::Surface::Formulation; returngeometry::surface::backend::mfem::make_pressure_boundary<Formulation>(coefficients); } else { returnNoSurfaceCondition{}; } }(); autoconstraints=std::apply([&](constauto&...equation) { returnstd::tuple{computations::point_constraint_discretization(equation,*discretization)...}; },definition.constraints); usingBinding=BoundModel<Model,Description,Discretization,decltype(boundary),decltype(constraints),Plan>; returnstd::make_shared<constBinding>(Binding{.model=std::move(model),.definition=std::move(definition), .discretization=std::move(discretization),.boundary=boundary,.constraints=std::move(constraints)}); } namespacedetail { template<typenameGeometry> [[nodiscard]]std::size_tgeometry_size(constGeometry&coefficients) { ifconstexpr(requires{coefficients.Size();}) { returnstatic_cast<std::size_t>(coefficients.Size()); } else { returncoefficients.size(); } } template<typenameBound,typenameCoefficients,std::size_t...Indices> [[nodiscard]]autoconstraint_inputs(constBound&bound,constCoefficients&coefficients,std::index_sequence<Indices...>) { returnutils::dependencies::make_dependency_context(operators::backend::mfem::make_point_border_binding_input< typenamestd::tuple_element_t<Indices,typenameBound::ConstraintBindings>::Equation::Border>(std::get<Indices>(bound.constraints),coefficients)...); } } template<typenameBound,typenameCoefficients,typenameGeometry> [[nodiscard]]autoequilibrium_inputs(constBound&bound,constCoefficients&coefficients,constGeometry&geometry_coefficients) { namespaceexecution=operators::backend::mfem; usingDescription=typenameBound::Description; usingPlan=typenameBound::Plan; constauto&discretization=*bound.discretization; detail::validate_coefficient_owners(discretization,coefficients,typenamePlan::Coefficients{}); constautocommon=detail::model_inputs(discretization,coefficients,typenameDescription::InteriorEquations::Equations{},typenamePlan::Coefficients{}); constautodeformation=[&]() { ifconstexpr(Plan::has_shape) { returngeometry::deformation::backend::mfem::make_deformation_inputs(discretization.deformation.model,discretization.deformation.reference,geometry_coefficients); } else { returnutils::dependencies::make_dependency_context(utils::dependencies::own<execution::DisplacementCoefficients>( execution::bind(discretization.templatebinding<discretization::blocks::displacement::vector>(),geometry_coefficients),discretization::backend::mfem::CoefficientTracking{})); } }(); constautoshape=[&]() { ifconstexpr(Plan::has_shape) { returnutils::dependencies::make_dependency_context(execution::make_volume_binding_input<typenameDescription::ShapeForce,discretization::blocks::surface_deformation::shape_equilibrium>( discretization.templateequation_discretization<typenameDescription::ShapeForce::Test>(),coefficients)); } else { returnutils::dependencies::make_dependency_context(); } }(); constautoboundary=[&]() { usingnamespaceutils::dependencies; ifconstexpr(Plan::has_surface) { constautoeos=borrow<geometry::surface::keys::EquationOfState>(bound.model->equation_of_state(),track_by([](constauto&value) { returnstd::addressof(value); })); returngeometry::surface::backend::mfem::make_pressure_boundary_metadata(bound.boundary,eos,bound.definition.surface_condition); } else { returnmake_dependency_context(); } }(); constautoconstraints=detail::constraint_inputs(bound,coefficients,std::make_index_sequence<std::tuple_size_v<typenameBound::ConstraintBindings>>{}); returnutils::dependencies::combine_dependency_contexts(common,deformation,shape,boundary,constraints); } template<typenamePlan,typenameDiscretization,typenameGeometry> [[nodiscard]]typenamePlan::Layoutequilibrium_layout(constDiscretization&discretization,constGeometry&geometry_coefficients) { constautosize=[&]<typenameTerm>()->std::size_t { ifconstexpr(std::same_as<Term,discretization::blocks::surface_deformation::coefficients>||std::same_as<Term,discretization::blocks::surface_deformation::shape_equilibrium>) { returndetail::geometry_size(geometry_coefficients); } else { returnstatic_cast<std::size_t>(discretization.templatebinding<Term>().size()); } }; constautosizes=[&]<typename...Terms>(utils::dependencies::KeyList<Terms...>) { returnstd::array<std::size_t,sizeof...(Terms)>{size.templateoperator()<Terms>()...}; }; return{.state_sizes=sizes(typenamePlan::Values{}),.residual_sizes=sizes(typenamePlan::Residuals{})}; } }