geometry.hpp algorithm memory tuple vector serif/backend/mfem/geometry/deformation/geometry_step.hpp serif/backend/mfem/mapping/domain_mapper.hpp serif/backend/mfem/model/graph.hpp serif/model/state.hpp serif/discretization/blocks/vector_view.hpp serif/solver/nonlinear/options.hpp src/include/serif/backend/mfem/solver/nonlinear/stellar_problem.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/geometry_step.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 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 src/include/serif/backend/mfem/solver/admissibility/geometry.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/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/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/geometry_step.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/model/state.hpp serif/operators/concepts.hpp serif/operators/contributions.hpp serif/operators/plan.hpp serif/solver/linear/linear.hpp serif/solver/nonlinear/options.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/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::solver::admissibility::backend::mfem::GeometrySamplingRules serif::solver::admissibility::backend::mfem::GeometryDirection serif serif::solver serif::solver::admissibility serif::solver::admissibility::backend serif::solver::admissibility::backend::mfem serif::solver::admissibility::backend::mfem::detail #pragmaonce #include<algorithm> #include<memory> #include<tuple> #include<vector> #include"serif/backend/mfem/geometry/deformation/geometry_step.hpp" #include"serif/backend/mfem/mapping/domain_mapper.hpp" #include"serif/backend/mfem/model/graph.hpp" #include"serif/model/state.hpp" #include"serif/discretization/blocks/vector_view.hpp" #include"serif/solver/nonlinear/options.hpp" namespaceserif::solver::admissibility::backend::mfem{ structGeometrySamplingRules{}; namespacedetail{ template<typenameEquation,typenameScope=typenameEquation::Test> [[nodiscard]]autoequation_rule_keys(){ return[]<std::size_t...Index>(std::index_sequence<Index...>){ returnutils::dependencies::KeyList<typename operators::backend::mfem::volume_keys::Scoped<Scope>::templateBasis<Index>...>{}; }(std::make_index_sequence<std::tuple_size_v<equations::ContributionsOf<Equation>>>{}); } [[nodiscard]]inlinestd::shared_ptr<conststd::vector<mapping::DMFixedMeshGeometry>>sample_reference_geometry( ::mfem::Mesh&mesh,constintelement,const::mfem::IntegrationRule&rule){ ::mfem::ElementTransformation&transformation=*mesh.GetElementTransformation(element); autosamples=std::make_shared<std::vector<mapping::DMFixedMeshGeometry>>(); samples->reserve(static_cast<std::size_t>(rule.GetNPoints())); for(intpoint=0;point<rule.GetNPoints();++point){ autogeometry=mapping::sample_fixed_mesh_geometry(transformation,rule.IntPoint(point)); if(!geometry){ returnnullptr; } samples->push_back(std::move(*geometry)); } returnsamples; } } template<typenameBound> [[nodiscard]]automake_geometry_computations(conststd::shared_ptr<constBound>&bound){ usingnamespaceutils::dependencies; ifconstexpr(!Bound::Plan::has_shape){returnstd::tuple{};} else{ //Thesurfacepotentialhasnovolumebasisbatch.Thegeneralgeometryrulesandeveryactive //interiorequation'srulesremainintheinspectionunion. usingInputs=decltype([]<typename...Equation>(utils::types::TypeList<Equation...>){ returnCombinedKeyList<decltype(detail::equation_rule_keys<Equation>())...>{}; }(typenameBound::Description::InteriorEquations::Equations{})); returnstd::tuple{make_computation<GeometrySamplingRules,Inputs>([bound](constauto&...basis){ std::vector<geometry::deformation::backend::mfem::GeometrySampleRule>rules; constauto&space=bound->discretization->templatebinding<discretization::blocks::displacement::vector>().space(); //Thereferencemeshisfixedforthelifeoftheserules,soitsgeometryissampledhereonce //ratherthanineveryNewtonsteplimit.EquationrulessharethesamplestheirBasisalreadyowns. for(intelement=0;element<space.GetNE();++element){ constauto*finite_element=space.GetFE(element); constintorder=std::max(finite_element->GetOrder()+2,2*space.GetMesh()->Dimension()*finite_element->GetOrder()); constauto&rule=::mfem::IntRules.Get(finite_element->GetGeomType(),order); rules.push_back({.element=element,.integration_rule=&rule, .reference_geometry=detail::sample_reference_geometry(*space.GetMesh(),element,rule)}); } constautoappend=[&](constauto&batch){ for(constauto&element:batch){ rules.push_back({.element=element->element,.integration_rule=element->rule, .reference_geometry=std::shared_ptr<conststd::vector<mapping::DMFixedMeshGeometry>>(element,&element->reference_geometry)}); } }; (append(basis),...); std::ranges::sort(rules,[](constauto&left,constauto&right){ if(left.element!=right.element){returnleft.element<right.element;} returnstd::less<const::mfem::IntegrationRule*>{}(left.integration_rule,right.integration_rule); }); constautoduplicate=std::ranges::unique(rules,[](constauto&left,constauto&right){ returnleft.element==right.element&&left.integration_rule==right.integration_rule; }); rules.erase(duplicate.begin(),duplicate.end()); returnrules; })}; } } classGeometryDirection{ public: template<typenameEvaluator,typenameLayout> [[nodiscard]]autolimit(constEvaluator&evaluator,constmodel::StellarState<Layout>&state,std::span<constdouble>direction, constdoublemaximum_step,constdoublesafety,constboolenabled) ->std::expected<geometry::deformation::GeometryStepLimit,geometry::deformation::DeformationFailure>{ geometry::deformation::GeometryStepLimitunrestricted{.step=maximum_step,.boundary_step=maximum_step}; m_geometry_ready=false; ifconstexpr(Evaluator::Plan::has_shape){ constauto&derivatives=evaluator.cache().templateget<model::computations::ShapeDerivatives>(); usingMap=std::remove_cvref_t<decltype(*derivatives.deformation->map)>; ifconstexpr(geometry::deformation::is_linear_deformation<Map>){ usingShape=discretization::blocks::surface_deformation::coefficients; usingDisplacement=discretization::blocks::displacement::vector; constauto&displacement=evaluator.binding().discretization->templatebinding<Displacement>(); constautoshape_direction=discretization::blocks::StateView<Layout, const double>(state.layout(),direction).templateget<Shape>(); ::mfem::Vectorshape(const_cast<double*>(shape_direction.data()),static_cast<int>(shape_direction.size())); ::mfem::Vectorreduced(displacement.size()); derivatives.tangent.action()->Mult(shape,reduced); m_displacement_direction=displacement.scatter_true(reduced); m_accepted_displacement=displacement.scatter_true(evaluator.cache().templateget<geometry::deformation::backend::mfem::volume_displacement_keys::Coefficients>()); m_rules=evaluator.cache().templateshare<GeometrySamplingRules>(); m_geometry_ready=true; if(enabled){ returnestimate(evaluator,maximum_step,safety); } } } returnunrestricted; } [[nodiscard]]doublenext(constauto&evaluator,constdoublestep,constBacktracking&options,constboolinverted)const{ if(!inverted||!m_geometry_ready){ returnstep*options.contraction; } constautogeometry=estimate(evaluator,step,options.fraction_to_boundary); if(!geometry){ returnstep*options.contraction; } returnnext_backtracking_step(step,geometry->minimum_determinant_at_accepted, geometry->minimum_determinant_at_maximum,true,options); } private: [[nodiscard]]autoestimate(constauto&evaluator,constdoublestep,constdoublesafety)const{ constauto&discretization=*evaluator.binding().discretization; returngeometry::deformation::backend::mfem::estimate_geometry_step_limit(discretization.mapper, discretization.templatebinding<discretization::blocks::displacement::vector>().space(),discretization.compactification, m_accepted_displacement,m_displacement_direction,*m_rules,{.maximum_step=step,.fraction_to_boundary=safety}); } ::mfem::Vectorm_accepted_displacement; ::mfem::Vectorm_displacement_direction; std::shared_ptr<conststd::vector<geometry::deformation::backend::mfem::GeometrySampleRule>>m_rules; boolm_geometry_ready{false}; }; }