polytropic.hpp serif/solver/continuation/family.hpp serif/solver/continuation/stellar_parameters.hpp src/compile_time_checks/serif/solver/continuation.cpp src/include/serif/solver/continuation/continuation.hpp serif/dimensions/dimensions.hpp serif/dimensions/partial.hpp serif/dimensions/quantities.hpp serif/dimensions/type_alias.hpp serif/eos/models/polytropic.hpp serif/eos/relations/catalog.hpp serif/eos/relations/relations.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/continuation/family.hpp serif/solver/continuation/parameters.hpp src/include/serif/solver/continuation/polytropic.hpp serif/solver/continuation/stellar_parameters.hpp serif/utils/dependencies/binding.hpp serif/utils/dependencies/detail/type_lists.hpp serif/utils/dependencies/key_list.hpp serif/utils/dependencies/key_list_operations.hpp serif/utils/error/codes.hpp serif/utils/error/errors.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 iostream limits magic_enum/magic_enum_all.hpp optional ranges span string string_view tuple type_traits unistd.h unordered_map utility vector version src/compile_time_checks/serif/backend_independent.cpp src/compile_time_checks/serif/solver/continuation.cpp src/include/serif/solver/continuation/continuation.hpp src/include/serif/solver/continuation/polytropic.hpp src/include/serif/solver/solver.hpp serif::solver::continuation::PreservePressureAtReferenceDensity serif serif::solver serif::solver::continuation #pragmaonce #include"serif/solver/continuation/family.hpp" #include"serif/solver/continuation/stellar_parameters.hpp" namespaceserif::solver::continuation{ classPreservePressureAtReferenceDensityfinal{ public: usingOutputParameter=PolytropicConstant; usingDependencyParameters=ParameterList<PolytropicIndex>; explicitPreservePressureAtReferenceDensity(dimensions::DensityValuedensity):m_density(density){ utils::errors::serif_ensure<utils::errors::SERiFErrorCode::INVALID_ARGUMENT, std::invalid_argument>( std::isfinite(density.value())&&density.value()>0.0,"Pressure-referencedensitymustbefiniteandpositive."); } template<model::StellarModelTypeModel> requiresModel::templatecontains_specification<eos::models::Polytrope> [[nodiscard]]autobind(constModel&target)const{ constauto&eos=target.templatespecification<eos::models::Polytrope>(); constdoubleindex=eos.polytropic_index(); constdoubleconstant=eos.polytropic_constant(); constdoubledensity=m_density.value(); constdoublepressure=eos.evaluate(eos::relations::PressureFromDensity{},m_density).value(); utils::errors::serif_ensure<utils::errors::SERiFErrorCode::INVALID_ARGUMENT, std::invalid_argument>( std::isfinite(pressure)&&pressure>0.0,"Thetargetpressureatthereferencedensitymustbefiniteandpositive."); returnlaw(PolytropicConstant{},depends_on(PolytropicIndex{}),[index,constant,density](doublecandidate){ if(candidate==index){ returnconstant; } returnstd::exp(std::log(constant)+(1.0/index-1.0/candidate)*std::log(density)); }); } private: dimensions::DensityValuem_density; }; [[nodiscard]]inlineautopreserve_pressure_at_reference_density(dimensions::DensityValuedensity){ returnPreservePressureAtReferenceDensity{density}; } }