graph.hpp limits memory tuple utility serif/backend/mfem/model/computations/interior.hpp serif/backend/mfem/model/computations/surface_potential.hpp serif/backend/mfem/model/computations/potential_infinity.hpp serif/backend/mfem/model/computations/pressure_boundary.hpp serif/model/equilibrium_plan.hpp serif/backend/mfem/geometry/deformation/computations.hpp serif/backend/mfem/model/volume_displacement.hpp serif/backend/mfem/geometry/surface/boundary_computations.hpp serif/model/keys.hpp serif/model/plan.hpp src/include/serif/backend/mfem/model/equilibrium.hpp src/include/serif/backend/mfem/solver/admissibility/geometry.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/element_coordinates.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/displacement_mapping.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_field_block.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/transformations.hpp serif/backend/mfem/mapping/variation.hpp serif/backend/mfem/mapping/volume_cache.hpp serif/backend/mfem/model/computations/interior.hpp serif/backend/mfem/model/computations/potential_infinity.hpp serif/backend/mfem/model/computations/pressure_boundary.hpp serif/backend/mfem/model/computations/scalar_row.hpp serif/backend/mfem/model/computations/shape_derivatives.hpp serif/backend/mfem/model/computations/surface_potential.hpp serif/backend/mfem/model/deformation.hpp src/include/serif/backend/mfem/model/graph.hpp serif/backend/mfem/model/scalar_row.hpp serif/backend/mfem/model/surface_potential.hpp serif/backend/mfem/model/system.hpp serif/backend/mfem/model/volume.hpp serif/backend/mfem/model/volume_actions.hpp serif/backend/mfem/model/volume_displacement.hpp serif/backend/mfem/operators/actions.hpp serif/backend/mfem/operators/element_data.hpp serif/backend/mfem/operators/element_evaluation.hpp serif/backend/mfem/operators/element_preparation.hpp serif/backend/mfem/operators/element_transpose.hpp serif/backend/mfem/operators/geometry_derivatives.hpp serif/backend/mfem/operators/linearized_block.hpp serif/backend/mfem/operators/operand_evaluation.hpp serif/backend/mfem/operators/preparation.hpp serif/backend/mfem/operators/quadrature.hpp serif/backend/mfem/operators/rejection.hpp serif/backend/mfem/operators/replaced_rows.hpp serif/backend/mfem/operators/volume_application.hpp serif/backend/mfem/operators/volume_computations.hpp serif/backend/mfem/operators/volume_data.hpp serif/backend/mfem/operators/volume_field_block.hpp serif/backend/mfem/operators/volume_inputs.hpp serif/backend/mfem/operators/volume_stages.hpp serif/backend/mfem/operators/workspace.hpp serif/backend/mfem/shared/finite.hpp serif/backend/mfem/shared/parallel.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/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/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/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/key_list.hpp serif/utils/dependencies/key_list_operations.hpp serif/utils/dependencies/keys.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/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 atomic cmath compare concepts cstddef cstdint cstdio cstdlib exception expected format functional initializer_list iostream limits magic_enum/magic_enum_all.hpp 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/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/graph.hpp src/include/serif/backend/mfem/solver/admissibility/geometry.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 serif::model serif::model::detail #pragmaonce #include<limits> #include<memory> #include<tuple> #include<utility> #include"serif/backend/mfem/model/computations/interior.hpp" #include"serif/backend/mfem/model/computations/surface_potential.hpp" #include"serif/backend/mfem/model/computations/potential_infinity.hpp" #include"serif/backend/mfem/model/computations/pressure_boundary.hpp" #include"serif/model/equilibrium_plan.hpp" #include"serif/backend/mfem/geometry/deformation/computations.hpp" #include"serif/backend/mfem/model/volume_displacement.hpp" #include"serif/backend/mfem/geometry/surface/boundary_computations.hpp" #include"serif/model/keys.hpp" #include"serif/model/plan.hpp" namespaceserif::model{ namespacedetail{ template<typenameLayout,typenameOwner,typename...ContributionKeys> [[nodiscard]]autophysical_system_computation(Layoutlayout,std::shared_ptr<constOwner>owner, utils::dependencies::KeyList<ContributionKeys...>){ usingnamespaceutils::dependencies; usingnamespacediscretization::blocks; usingnamespaceutils::errors; conststd::size_trows=layout.templatetotal_size<RepresentedValueKind::residual>(); conststd::size_tcolumns=layout.templatetotal_size<RepresentedValueKind::value>(); serif_ensure<SERiFErrorCode::INVALID_DOF_COUNT>(rows<=std::numeric_limits<int>::max()&&columns<=std::numeric_limits<int>::max(),"Stellarsystemdimensionsexceedthebackend'sintegerrange."); returnmake_computation<keys::PhysicalSystem,InputList<ContributionKeys...>>([layout,owner=std::move(owner)](constauto&...contributions){ computations::Contributionscombined; combined.actions.emplace(); constautoappend=[&](constcomputations::Contributions&part){ combined.residual.insert(combined.residual.end(),part.residual.begin(),part.residual.end()); combined.jacobian.insert(combined.jacobian.end(),part.jacobian.begin(),part.jacobian.end()); combined.replacements.insert(combined.replacements.end(),part.replacements.begin(),part.replacements.end()); constauto&actions=part.actions?*part.actions:part.jacobian; combined.actions->insert(combined.actions->end(),actions.begin(),actions.end()); }; (append(contributions),...); returnbackend::mfem::make_physical_system( static_cast<int>(layout.templatetotal_size<RepresentedValueKind::residual>()), static_cast<int>(layout.templatetotal_size<RepresentedValueKind::value>()), std::move(combined.residual),std::move(combined.jacobian),std::move(combined.replacements),owner,std::move(combined.actions)); }); } template<typename...Equations,typename...Additional> [[nodiscard]]constevalautocontribution_keys(utils::types::TypeList<Equations...>,utils::dependencies::KeyList<Additional...>){ returnutils::dependencies::KeyList<computations::InteriorContribution<typename Equations::Test>...,Additional...>{}; } template<typenameRegistry,typenameDiscretization,typenameLayout,typename...Equations> [[nodiscard]]autointerior_computations(conststd::shared_ptr<constRegistry>&registry,constDiscretization&discretization, constLayout&layout,operators::backend::mfem::VolumeOptionsoptions,utils::types::TypeList<Equations...>){ returnstd::tuple_cat(operators::backend::mfem::make_volume_equation_computations<typename Equations::Test>( std::shared_ptr<constEquations>(registry,&registry->templateget<Equations>()), discretization.templateequation_discretization<typenameEquations::Test>(),options)..., computations::interior_computations<Equations>(layout)...); } } template<ExecutableInteriorEquationsRegistry,typenameDiscretization,discretization::blocks::SystemLayoutTypeLayout, typename...AdditionalComputations,typename...AdditionalContributions> [[nodiscard]]automake_model_graph(std::shared_ptr<constRegistry>registry,std::shared_ptr<constDiscretization>discretization, constLayout&layout,std::tuple<AdditionalComputations...>additional={}, utils::dependencies::KeyList<AdditionalContributions...>contributions={},operators::backend::mfem::VolumeOptionsoptions={}){ usingnamespaceutils::errors; usingEquations=typenameRegistry::Equations; serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(registry!=nullptr&&discretization!=nullptr,"Modelgraphconstructionrequiresanowningequationregistryanddiscretizationbinding."); constautointerior=detail::interior_computations(registry,*discretization,layout,options,Equations{}); constautophysical=detail::physical_system_computation(layout,discretization,detail::contribution_keys(Equations{},contributions)); returnstd::apply([](auto...computations){ returnutils::dependencies::make_computation_graph(std::move(computations)...); },std::tuple_cat(interior,std::move(additional),std::tuple{physical})); } template<ExecutableInteriorEquationsRegistry,typenameDiscretization,discretization::blocks::SystemLayoutTypeLayout, typename...AdditionalComputations,typename...AdditionalContributions> [[nodiscard]]automake_model_graph(Registryregistry,std::shared_ptr<constDiscretization>discretization, constLayout&layout,std::tuple<AdditionalComputations...>additional={}, utils::dependencies::KeyList<AdditionalContributions...>contributions={},operators::backend::mfem::VolumeOptionsoptions={}){ returnmake_model_graph(std::make_shared<constRegistry>(std::move(registry)),std::move(discretization),layout,std::move(additional),contributions,std::move(options)); } namespacedetail{ template<typenameBound,typenameLayout> [[nodiscard]]autoshape_computations(conststd::shared_ptr<constBound>&bound,constLayout&layout){ namespacedeformation=geometry::deformation::backend::mfem; namespaceexecution=operators::backend::mfem; ifconstexpr(Bound::Plan::has_shape){ usingSurface=std::remove_cvref_t<decltype(bound->discretization->deformation.model.surface)>; usingMetadata=typenameexecution::volume_keys::Scoped<discretization::blocks::gravity::potential>::Bindings; constauto&displacement=bound->discretization->templatebinding<discretization::blocks::displacement::vector>(); returnstd::tuple_cat(deformation::make_deformation_computations<Surface>(), deformation::make_volume_displacement_computations<Metadata>(displacement), std::tuple{computations::shape_derivatives(displacement), computations::surface_potential_contribution(bound->definition.shape_force,layout)}); } else{ returnstd::tuple{}; } } template<typenameBound,typenameLayout> [[nodiscard]]autoboundary_computations(constBound&,constLayout&layout){ usingPlan=typenameBound::Plan; ifconstexpr(Plan::has_surface){ usingFormulation=typenamePlan::Surface::Formulation; returnstd::tuple_cat(geometry::surface::backend::mfem::make_pressure_boundary_computations<Formulation>(), std::tuple{computations::pressure_boundary_contribution<Formulation>(layout)}); } else{ returnstd::tuple{}; } } template<typenamePlan> [[nodiscard]]constevalautoequilibrium_contribution_keys(){ usingnamespaceutils::dependencies; usingShape=std::conditional_t<Plan::has_shape, KeyList<computations::ShapeContribution>,KeyList<>>; usingSurface=std::conditional_t<Plan::has_surface,KeyList<computations::PressureBoundaryContribution>,KeyList<>>; usingInfinity=std::conditional_t<Plan::has_primal_gravity,KeyList<computations::PotentialInfinityContribution>,KeyList<>>; returnCombinedKeyList<Shape,Surface,Infinity>{}; } } template<typenameBound,typenameLayout,typename...NumericalComputations> [[nodiscard]]automake_equilibrium_graph(std::shared_ptr<constBound>bound,constLayout&layout, std::tuple<NumericalComputations...>numerical={},operators::backend::mfem::VolumeOptionsoptions={}){ usingnamespaceutils::errors; usingDescription=typenameBound::Description; usingRegistry=typenameDescription::InteriorEquations; serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(bound!=nullptr,"Anequilibriumgraphrequiresaboundmodelowner.Bindthephysicaldescriptionanddiscretizationbeforeconstructingthegraph."); constautoshape=detail::shape_computations(bound,layout); constautosurface=detail::boundary_computations(*bound,layout); constautoinfinity=[&](){ ifconstexpr(Bound::Plan::has_primal_gravity){ returnstd::tuple{computations::potential_infinity_contribution(layout, bound->discretization->templatebinding<discretization::blocks::gravity::potential>())}; }else{ returnstd::tuple{}; } }(); constautocontributions=detail::equilibrium_contribution_keys<typename Bound::Plan>(); returnmake_model_graph(std::shared_ptr<constRegistry>(bound,&bound->definition.interior_equations),bound,layout, std::tuple_cat(shape,surface,infinity,std::move(numerical)),contributions,options); } }