state.hpp algorithm cmath expected format span string vector serif/backend/mfem/discretization/dof_selection.hpp serif/backend/mfem/model/binding.hpp serif/model/state.hpp serif/backend/mfem/solver/seed/fitted.hpp serif/backend/mfem/mapping/evaluator.hpp serif/backend/mfem/projection/radial.hpp src/include/serif/solver/seed/state.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/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/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_boundary.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/evaluator.hpp serif/backend/mfem/mapping/status.hpp serif/backend/mfem/mapping/variation.hpp serif/backend/mfem/mapping/volume_cache.hpp serif/backend/mfem/model/binding.hpp serif/backend/mfem/model/inputs.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/volume_inputs.hpp serif/backend/mfem/operators/workspace.hpp serif/backend/mfem/projection/radial.hpp serif/backend/mfem/shared/finite.hpp serif/backend/mfem/shared/reference_tracking.hpp serif/backend/mfem/solver/seed/fitted.hpp src/include/serif/backend/mfem/solver/seed/state.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/polytropic.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/specification_list.hpp serif/model/detail/surface_plan.hpp serif/model/equilibrium_equations.hpp serif/model/equilibrium_plan.hpp serif/model/plan.hpp serif/model/specification.hpp serif/model/specifications/fixed_angular_momentum.hpp serif/model/specifications/fixed_total_mass.hpp serif/model/specifications/physical_names.hpp serif/model/state.hpp serif/model/stellar_model.hpp serif/operators/concepts.hpp serif/operators/plan.hpp serif/solver/seed/fitted.hpp serif/solver/seed/polytrope.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/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/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/laneemden/laneemden.hpp serif/utils/laneemden/scaling.hpp serif/utils/laneemden/solution.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/projection/radial.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/general/array.hpp mfem/linalg/densemat.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/seed/state.hpp src/include/serif/solver/seed/state.hpp src/include/serif/solver/solver.hpp serif::solver::seed::backend::mfem::ProjectionOptions serif serif::solver serif::solver::seed serif::solver::seed::backend serif::solver::seed::backend::mfem #pragmaonce #include<algorithm> #include<cmath> #include<expected> #include<format> #include<span> #include<string> #include<vector> #include"serif/backend/mfem/discretization/dof_selection.hpp" #include"serif/backend/mfem/model/binding.hpp" #include"serif/model/state.hpp" #include"serif/backend/mfem/solver/seed/fitted.hpp" #include"serif/backend/mfem/mapping/evaluator.hpp" #include"serif/backend/mfem/projection/radial.hpp" namespaceserif::solver::seed::backend::mfem{ structProjectionOptionsfinal{ doublesurface_radius_relative_tolerance{5.0e-4}; }; template<model::StellarModelTypeModel,typenameDiscretization,typenameProfile> requires(std::same_as<Profile,RadialSeed>||std::same_as<Profile,FittedRadialSeed>)&&requires(constModel&model,constDiscretization&discretization){ requires(!std::same_as<typenamestd::remove_cvref_t<decltype(model::resolve_equilibrium_equations(model))>::ShapeForce,model::NoShapeForce>); discretization.deformation.model.surface; discretization.deformation.reference; } [[nodiscard]]automake_polytropic_state(constModel&model,constDiscretization&discretization,constProfile&seed,constProjectionOptionsoptions={}){ namespaceblocks=serif::discretization::blocks; namespacedependencies=utils::dependencies; namespaceprojection=utils::projection::backend::mfem; usingDescription=std::remove_cvref_t<decltype(model::resolve_equilibrium_equations(model))>; usingPlan=model::EquilibriumPlan<typename Model::EquationOfState, Description, dependencies::KeyList<blocks::displacement::vector>>; usingState=model::StellarState<typename Plan::Layout>; usingResult=std::expected<State,std::string>; constRadialSeed&source=serif::solver::seed::detail::radial_source(seed); constdimensions::LengthValueseed_radius=serif::solver::seed::detail::radial_radius(seed); if(!std::isfinite(options.surface_radius_relative_tolerance)||options.surface_radius_relative_tolerance<0.0){ returnResult(std::unexpected("Theradialseedsurface-radiustolerancemustbefiniteandnonnegative.")); } ifconstexpr(std::same_as<typenameDescription::SurfaceCondition,geometry::surface::Isobaric>){ if(model::resolve_equilibrium_equations(model).surface_condition.target_pressure().value()!=0.0){ returnResult(std::unexpected("AcompleteLane-Emdenseedrequiresazero-pressureisobaricsurface.")); } } std::vector<double>shape; ifconstexpr(Plan::has_shape){ constautoprepared=detail::radial_seed_surface(discretization); if(!prepared){ returnResult(std::unexpected(prepared.error())); } constauto&surface=**prepared; for(constauto&point:surface.reference_positions()){ constdoubleradius=std::hypot(point[0]-source.center[0],point[1]-source.center[1],point[2]-source.center[2]); constdoublescale=std::max(radius,seed_radius.value()); if(std::abs(radius-seed_radius.value())>options.surface_radius_relative_tolerance*scale){ returnResult(std::unexpected(std::format( "Thereferencesurfaceradius{}differsfromtheradialseedradius{}(relativeerror{},allowed{}).Chooseasphericalreferencemeshattheseedradius,resolvingitsgeometrysufficiently.", radius,seed_radius.value(),std::abs(radius-seed_radius.value())/scale,options.surface_radius_relative_tolerance))); } } shape.resize(surface.input_size(),0.0); } Statestate(model::backend::mfem::equilibrium_layout<Plan>(discretization,shape)); ::mfem::GridFunctiondisplacement(&discretization.templatebinding<blocks::displacement::vector>().space()); displacement=0.0; mapping::GridFunctionMappingEvaluatormapping(discretization.mapper,displacement,discretization.compactification); constdoublegravitational_mass=source.scales.gravitational_constant.value()*serif::solver::seed::detail::radial_mass(seed).value(); constprojection::RadialProjectionOptionsprojection_options{ .center=source.center, .mapping=&mapping, .scalar_at_infinity=[](const::mfem::Vector&)->utils::projection::RadialValue{return0.0;} }; constautophysical=[&](doubleradius){ returnevaluate_radial_seed(seed,dimensions::LengthValue{radius}); }; constautoproject=[&]<typenameTerm>(automember)->std::expected<void,std::string>{ constauto&binding=discretization.templatebinding<Term>(); constautoprofile=[&](doubleradius)->utils::projection::RadialValue{ returnphysical(radius).transform([member](constutils::laneemden::PhysicalSample&sample){return(sample.*member).value();}); }; autofield=projection::project_scalar(binding.space(),profile,projection_options); if(!field){ returnstd::unexpected(field.error()); } const::mfem::Vectorreduced=binding.gather(*field); std::copy_n(reduced.GetData(),reduced.Size(),state.templateget<Term>().begin()); return{}; }; std::expected<void,std::string>status; constautoassign=[&]<typenameTerm>(automember){ ifconstexpr(dependencies::ContainsKey<Term,typenamePlan::Values>){ if(status){ status=project.templateoperator()<Term>(member); } } }; assign.templateoperator()<blocks::density::mass>(&utils::laneemden::PhysicalSample::density); assign.templateoperator()<blocks::gravity::potential>(&utils::laneemden::PhysicalSample::potential); if(!status){ returnResult(std::unexpected(status.error())); } ifconstexpr(std::same_as<typenameModel::EquationOfState,eos::models::Polytrope>&& std::same_as<typenameDescription::SurfaceCondition,geometry::surface::Isobaric>){ usingFormulation=typenamePlan::Surface::Formulation; ifconstexpr(!std::same_as<Formulation,void>){ usingCarrier=typenameFormulation::Test; if(model::resolve_equilibrium_equations(model).surface_condition.target_pressure().value()==0.0){ //Projectionroundoffattheanalyticsurfacemustnotspoiltheexactprescribedtrace. constautoselected=serif::discretization::backend::mfem::select_boundary_dofs<typename Formulation::Support>( discretization.templatebinding<Carrier>()); for(constintdof:selected.reduced_dofs()){ state.templateget<Carrier>()[dof]=0.0; } } } } ifconstexpr(dependencies::ContainsKey<blocks::barotropic_constant::mass_normalization,typenamePlan::Values>){ //Usethenonrotatingsphericalvalue-GM/RastheinitialguessforeithersurfacelevelorBernoulliconstant. state.templateget<blocks::barotropic_constant::mass_normalization>()[0]=-gravitational_mass/seed_radius.value(); } ifconstexpr(dependencies::ContainsKey<blocks::rotation::angular_velocity,typenamePlan::Values>){ state.templateget<blocks::rotation::angular_velocity>()[0]=seed.angular_velocity; } returnResult(std::move(state)); } }