fitted.hpp algorithm cmath limits serif/backend/mfem/geometry/deformation/computations.hpp serif/solver/seed/fitted.hpp serif/utils/dependencies/computation_cache.hpp src/include/serif/backend/mfem/solver/seed/state.hpp src/include/serif/solver/seed/state.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/shared/reference_tracking.hpp src/include/serif/backend/mfem/solver/seed/fitted.hpp serif/dimensions/dimensions.hpp serif/dimensions/partial.hpp serif/dimensions/quantities.hpp serif/dimensions/type_alias.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/eos/models/polytropic.hpp serif/eos/relations/catalog.hpp serif/eos/relations/relations.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/meta/enums.hpp serif/model/detail/specification_list.hpp serif/model/equilibrium_equations.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/stellar_model.hpp serif/solver/seed/fitted.hpp serif/solver/seed/polytrope.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/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/types/strings.hpp serif/utils/types/type_list.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/fespace.hpp mfem/linalg/sparsemat.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/fitted.hpp src/include/serif/backend/mfem/solver/seed/state.hpp src/include/serif/solver/seed/state.hpp src/include/serif/solver/solver.hpp serif serif::solver serif::solver::seed serif::solver::seed::backend serif::solver::seed::backend::mfem serif::solver::seed::backend::mfem::detail #pragmaonce #include<algorithm> #include<cmath> #include<limits> #include"serif/backend/mfem/geometry/deformation/computations.hpp" #include"serif/solver/seed/fitted.hpp" #include"serif/utils/dependencies/computation_cache.hpp" namespaceserif::solver::seed::backend::mfem{ namespacedetail{ template<typenameDiscretization> [[nodiscard]]autoradial_seed_surface(constDiscretization&discretization){ namespacedeformation=geometry::deformation; namespacedependencies=utils::dependencies; usingSurface=std::remove_cvref_t<decltype(discretization.deformation.model.surface)>; conststd::vector<double>shape; constautoinputs=deformation::backend::mfem::make_deformation_inputs(discretization.deformation.model,discretization.deformation.reference,shape); autocache=dependencies::make_computation_cache<deformation::DeformationFailure>(deformation::backend::mfem::make_deformation_graph<Surface>(),inputs); usingOwner=std::remove_cvref_t<decltype(cache.templateget<deformation::keys::SurfaceMap>())>; usingResult=std::expected<Owner,std::string>; constautoprepared=cache.try_prepare(inputs,dependencies::KeyList<deformation::keys::SurfaceMap>{}); if(!prepared){ returnResult(std::unexpected(prepared.error().error.message)); } returnResult(cache.templateget<deformation::keys::SurfaceMap>()); } } template<typenameDiscretization> [[nodiscard]]std::expected<FittedRadialSeed,std::string>fit_radial_seed(constDiscretization&discretization, constRadialSeed&seed,doublespherical_relative_tolerance=5.0e-4){ if(!std::isfinite(spherical_relative_tolerance)||spherical_relative_tolerance<0.0){ returnstd::unexpected("Radial-guessfittingrequiresafinitenonnegativesphericaltolerance."); } constautosurface=detail::radial_seed_surface(discretization); if(!surface){ returnstd::unexpected(surface.error()); } doubleminimum=std::numeric_limits<double>::infinity(); doublemaximum=0.0; for(constauto&point:(*surface)->reference_positions()){ constdoubleradius=std::hypot(point[0]-seed.center[0],point[1]-seed.center[1],point[2]-seed.center[2]); if(!std::isfinite(radius)||radius<=0.0){ returnstd::unexpected("Radial-guessfittingrequiresfinitepositivereference-surfaceradiiabouttheseedcenter."); } minimum=std::min(minimum,radius); maximum=std::max(maximum,radius); } constdoubleradius=minimum+0.5*(maximum-minimum); if(!std::isfinite(radius)||radius<=0.0||(maximum-minimum)/radius>spherical_relative_tolerance){ returnstd::unexpected("Radial-guessfittingrequiresanonemptysphericalreferencesurfacewithintheselectedrelativetolerance."); } returnserif::solver::seed::fit_radial_seed(seed,dimensions::LengthValue{radius}); } }