solution.hpp cstddef expected optional string string_view vector format serif/utils/error/errors.hpp src/include/serif/utils/laneemden/laneemden.hpp src/include/serif/utils/laneemden/scaling.hpp src/lib/serif/utils/laneemden/solution.cpp serif/meta/enums.hpp serif/utils/error/codes.hpp serif/utils/error/errors.hpp src/include/serif/utils/laneemden/solution.hpp serif/utils/misc/concepts/exceptions.hpp serif/utils/misc/concepts/numeric.hpp serif/utils/misc/terminal/colors.hpp serif/utils/types/strings.hpp algorithm concepts cstddef cstdint cstdio cstdlib exception expected format iostream magic_enum/magic_enum_all.hpp optional string string_view type_traits unistd.h unordered_map utility vector version src/compile_time_checks/serif/backend_independent.cpp src/compile_time_checks/serif/model/barotropic_density.cpp src/compile_time_checks/serif/model/density_primal.cpp src/compile_time_checks/serif/model/normalization.cpp src/compile_time_checks/serif/model/state.cpp src/compile_time_checks/serif/solver/public_contracts.cpp src/compile_time_checks/serif/utils/laneemden.cpp src/include/serif/backend/mfem/model/equilibrium.hpp src/include/serif/backend/mfem/model/normalization/computations.hpp src/include/serif/backend/mfem/serialize/configuration.hpp src/include/serif/backend/mfem/serialize/options.hpp src/include/serif/backend/mfem/solver/equilibrium.hpp src/include/serif/backend/mfem/solver/nonlinear/stellar_problem.hpp src/include/serif/backend/mfem/solver/seed/fitted.hpp src/include/serif/backend/mfem/solver/seed/state.hpp src/include/serif/model/configured_discretization.hpp src/include/serif/model/normalization/plan.hpp src/include/serif/model/normalization/reference.hpp src/include/serif/solver/continuation/equilibrium.hpp src/include/serif/solver/equilibrium.hpp src/include/serif/solver/normalization/configuration.hpp src/include/serif/solver/normalization/normalization.hpp src/include/serif/solver/seed/fitted.hpp src/include/serif/solver/seed/polytrope.hpp src/include/serif/solver/seed/state.hpp src/include/serif/solver/solver.hpp src/include/serif/utils/laneemden/laneemden.hpp src/include/serif/utils/laneemden/scaling.hpp src/include/serif/utils/laneemden/solution.hpp src/include/serif/utils/serialize/numerics.hpp src/include/serif/utils/serialize/serialize.hpp src/include/serif/utils/serialize/solver.hpp src/lib/serif/solver/seed/fitted.cpp src/lib/serif/solver/seed/polytrope.cpp src/lib/serif/utils/laneemden/scaling.cpp src/lib/serif/utils/laneemden/solution.cpp serif::utils::laneemden::IntegrationOptions serif::utils::laneemden::Sample serif::utils::laneemden::Solution std::formatter< serif::utils::laneemden::IntegrationOptions > std::formatter< serif::utils::laneemden::Sample > std::formatter< serif::utils::laneemden::Solution > serif serif::utils serif::utils::laneemden #pragmaonce #include<cstddef> #include<expected> #include<optional> #include<string> #include<string_view> #include<vector> #include<format> #include"serif/utils/error/errors.hpp" namespaceserif::utils::laneemden{ structIntegrationOptionsfinal{ doublestep=0.01; doublemaximum_radius=100.0; std::size_tmaximum_steps=100000; doublesurface_tolerance=1e-10; std::size_tcentral_series_terms=12; friendbooloperator==(constIntegrationOptions&,constIntegrationOptions&)=default; }; structSamplefinal{ doublexi; doubletheta; doublederivative; }; structSolutionfinal{ doublen; std::vector<Sample>samples; std::optional<Sample>surface; }; [[nodiscard]]std::expected<Solution,std::string>solve(doublen,IntegrationOptionsoptions={}); [[nodiscard]]std::expected<Sample,std::string>evaluate(constSolution&solution,doublexi); } template<> structstd::formatter<serif::utils::laneemden::IntegrationOptions>{ template<typenameFormatContext> autoformat(constserif::utils::laneemden::IntegrationOptions&options,FormatContext&ctx)const{ returnstd::format_to(ctx.out(),"step={},maximum_radius={},maximum_steps={},surface_tolerance={},central_series_terms={}", options.step,options.maximum_radius,options.maximum_steps,options.surface_tolerance,options.central_series_terms); } }; template<> structstd::formatter<serif::utils::laneemden::Sample>{ enumclassComponent{all,xi,theta,derivative}; Componentcomponent{Component::all}; constexprautoparse(std::format_parse_context&ctx){ usingnamespaceserif::utils::errors; component=Component::all; autoend=ctx.begin(); while(end!=ctx.end()&&*end!='}'){ ++end; } if(end==ctx.begin()){ returnend; } conststd::string_viewselection(ctx.begin(),end); if(selection=="all"){ component=Component::all; }elseif(selection=="xi"){ component=Component::xi; }elseif(selection=="theta"){ component=Component::theta; }elseif(selection=="derivative"){ component=Component::derivative; }else{ serif_error<std::format_error, SERiFErrorCode::INVALID_ARGUMENT>( "UnknownLane–Emdensampleselector'{}'.Useanemptyspecifier,'all','xi','theta',or'derivative'.",selection); } returnend; } template<typenameFormatContext> autoformat(constserif::utils::laneemden::Sample&sample,FormatContext&ctx)const{ switch(component){ caseComponent::xi: returnstd::format_to(ctx.out(),"{}",sample.xi); caseComponent::theta: returnstd::format_to(ctx.out(),"{}",sample.theta); caseComponent::derivative: returnstd::format_to(ctx.out(),"{}",sample.derivative); default: returnstd::format_to(ctx.out(),"(ξ={},θ={},dθ/dξ={})",sample.xi,sample.theta,sample.derivative); } } }; template<> structstd::formatter<serif::utils::laneemden::Solution>{ std::formatter<serif::utils::laneemden::Sample>sample_formatter; constexprautoparse(std::format_parse_context&ctx){ returnsample_formatter.parse(ctx); } template<typenameFormatContext> autoformat(constserif::utils::laneemden::Solution&solution,FormatContext&ctx)const{ ctx.advance_to(std::format_to(ctx.out(),"Solution(n={},samples=[",solution.n)); for(std::size_tindex=0;index<solution.samples.size();++index){ if(index!=0){ ctx.advance_to(std::format_to(ctx.out(),",")); } ctx.advance_to(sample_formatter.format(solution.samples[index],ctx)); } ctx.advance_to(std::format_to(ctx.out(),"],surface=")); if(solution.surface){ ctx.advance_to(sample_formatter.format(*solution.surface,ctx)); }else{ ctx.advance_to(std::format_to(ctx.out(),"none")); } returnstd::format_to(ctx.out(),")"); } };