surface_power.hpp cmath span string_view serif/discretization/domain/physical_domains.hpp serif/utils/error/errors.hpp src/include/serif/backend/mfem/discretization/reconstruction.hpp serif/discretization/domain/concepts.hpp serif/discretization/domain/physical_domains.hpp serif/discretization/domain/types.hpp src/include/serif/discretization/reconstruction/surface_power.hpp serif/meta/enums.hpp serif/utils/error/codes.hpp serif/utils/error/errors.hpp serif/utils/misc/concepts/exceptions.hpp serif/utils/misc/concepts/numeric.hpp serif/utils/misc/std_helper/cleaning.hpp serif/utils/misc/std_helper/variadic_cleaning.hpp serif/utils/misc/terminal/colors.hpp serif/utils/types/strings.hpp algorithm array cmath concepts cstddef cstdint cstdio cstdlib exception format iostream magic_enum/magic_enum_all.hpp numeric span string string_view type_traits unistd.h unordered_map utility version src/compile_time_checks/serif/geometry/surface_boundary.cpp src/compile_time_checks/serif/model/finite_elements.cpp src/compile_time_checks/serif/model/scalar_row.cpp src/compile_time_checks/serif/model/state.cpp src/compile_time_checks/serif/operators/computation_owners.cpp src/compile_time_checks/serif/operators/plan.cpp src/include/serif/backend/mfem/discretization/basis_evaluation.hpp src/include/serif/backend/mfem/discretization/boundary_trace.hpp src/include/serif/backend/mfem/discretization/dof_map.hpp src/include/serif/backend/mfem/discretization/dof_selection.hpp src/include/serif/backend/mfem/discretization/element_coordinates.hpp src/include/serif/backend/mfem/discretization/field_values.hpp src/include/serif/backend/mfem/discretization/mass_diagonal.hpp src/include/serif/backend/mfem/discretization/reconstruction.hpp src/include/serif/backend/mfem/discretization/tensor_basis.hpp src/include/serif/backend/mfem/discretization/term_binding.hpp src/include/serif/backend/mfem/discretization/term_coefficients.hpp src/include/serif/backend/mfem/geometry/deformation/displacement_mapping.hpp src/include/serif/backend/mfem/geometry/surface/boundary_computations.hpp src/include/serif/backend/mfem/geometry/surface/boundary_data.hpp src/include/serif/backend/mfem/geometry/surface/boundary_inputs.hpp src/include/serif/backend/mfem/geometry/surface/boundary_rows.hpp src/include/serif/backend/mfem/geometry/surface/pressure_boundary.hpp src/include/serif/backend/mfem/geometry/surface/pressure_field_block.hpp src/include/serif/backend/mfem/model/binding.hpp src/include/serif/backend/mfem/model/computations/interior.hpp src/include/serif/backend/mfem/model/computations/potential_infinity.hpp src/include/serif/backend/mfem/model/computations/pressure_boundary.hpp src/include/serif/backend/mfem/model/computations/scalar_row.hpp src/include/serif/backend/mfem/model/computations/surface_potential.hpp src/include/serif/backend/mfem/model/discretization.hpp src/include/serif/backend/mfem/model/equilibrium.hpp src/include/serif/backend/mfem/model/equilibrium_evaluator.hpp src/include/serif/backend/mfem/model/finite_elements.hpp src/include/serif/backend/mfem/model/graph.hpp src/include/serif/backend/mfem/model/inputs.hpp src/include/serif/backend/mfem/model/normalization/computations.hpp src/include/serif/backend/mfem/model/state.hpp src/include/serif/backend/mfem/model/surface_potential.hpp src/include/serif/backend/mfem/model/volume.hpp src/include/serif/backend/mfem/model/volume_actions.hpp src/include/serif/backend/mfem/model/volume_displacement.hpp src/include/serif/backend/mfem/operators/element_data.hpp src/include/serif/backend/mfem/operators/element_evaluation.hpp src/include/serif/backend/mfem/operators/element_preparation.hpp src/include/serif/backend/mfem/operators/element_transpose.hpp src/include/serif/backend/mfem/operators/geometry_derivatives.hpp src/include/serif/backend/mfem/operators/jacobian_operator.hpp src/include/serif/backend/mfem/operators/operand_evaluation.hpp src/include/serif/backend/mfem/operators/volume_application.hpp src/include/serif/backend/mfem/operators/volume_computations.hpp src/include/serif/backend/mfem/operators/volume_data.hpp src/include/serif/backend/mfem/operators/volume_field_block.hpp src/include/serif/backend/mfem/operators/volume_inputs.hpp src/include/serif/backend/mfem/operators/volume_stages.hpp src/include/serif/backend/mfem/operators/workspace.hpp src/include/serif/backend/mfem/serialize/configuration.hpp src/include/serif/backend/mfem/serialize/discretization.hpp src/include/serif/backend/mfem/serialize/options.hpp src/include/serif/backend/mfem/serialize/resources.hpp src/include/serif/backend/mfem/serialize/top.hpp 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/backend/mfem/solver/normalization/metric_computations.hpp src/include/serif/backend/mfem/solver/seed/state.hpp src/include/serif/discretization/reconstruction/surface_power.hpp src/include/serif/model/finite_elements.hpp src/include/serif/solver/continuation/equilibrium.hpp src/include/serif/solver/equilibrium.hpp src/include/serif/solver/seed/state.hpp src/include/serif/solver/solver.hpp src/lib/serif/backend/mfem/discretization/basis_evaluation.cpp src/lib/serif/backend/mfem/discretization/boundary_trace.cpp src/lib/serif/backend/mfem/discretization/element_coordinates.cpp src/lib/serif/backend/mfem/discretization/mass_diagonal.cpp src/lib/serif/backend/mfem/discretization/reconstruction.cpp src/lib/serif/backend/mfem/model/deformation.cpp src/lib/serif/backend/mfem/model/finite_elements.cpp src/lib/serif/backend/mfem/model/surface_potential.cpp src/lib/serif/backend/mfem/model/system_selection.cpp src/lib/serif/backend/mfem/serialize/resources.cpp src/lib/serif/backend/mfem/serialize/top.cpp serif::discretization::reconstruction::SurfaceRegion serif::discretization::reconstruction::SurfaceAdjacent serif::discretization::reconstruction::SurfacePower serif::discretization::reconstruction::RadialProfile serif::discretization::reconstruction::SurfacePowerProfile serif serif::discretization serif::discretization::reconstruction #pragmaonce #include<cmath> #include<span> #include<string_view> #include"serif/discretization/domain/physical_domains.hpp" #include"serif/utils/error/errors.hpp" namespaceserif::discretization::reconstruction{ structSurfaceRegion{ std::string_viewboundary_name; }; template<typenameBoundaryT=domain::StellarSurfaceBoundary> structSurfaceAdjacent:SurfaceRegion{ constexprSurfaceAdjacent():SurfaceRegion{BoundaryT::name}{} }; template<typenameTermT> structSurfacePower{ usingTerm=TermT; doubledensity_exponent{1.0}; SurfaceRegionregion{SurfaceAdjacent<>{}}; }; classRadialProfile{ public: virtual~RadialProfile()=default; [[nodiscard]]virtualstd::string_viewname()constnoexcept=0; virtualvoidevaluate(std::span<constdouble>nodes,doubleinward_coordinate, std::span<double>values,std::span<double>derivatives)const=0; }; classSurfacePowerProfilefinal:publicRadialProfile{ public: explicitSurfacePowerProfile(doubledensity_exponent):m_exponent(density_exponent){ usingnamespaceutils::errors; serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(std::isfinite(m_exponent)&&m_exponent>=1.0, "SurfacePowerrequiresafinitedensityexponent>=1;received{}.Choosethepolytropicindexofavacuummodelinthesupportedrange.",m_exponent); } [[nodiscard]]std::string_viewname()constnoexceptoverride{return"surface_power";} [[nodiscard]]doubledensity_exponent()constnoexcept{returnm_exponent;} voidevaluate(std::span<constdouble>nodes,doubleinward_coordinate, std::span<double>values,std::span<double>derivatives)constoverride{ usingnamespaceutils::errors; serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(nodes.size()>=2&&values.size()==nodes.size()&&derivatives.size()==nodes.size()&& nodes.front()==0.0&&nodes.back()==1.0&&std::isfinite(inward_coordinate)&&inward_coordinate>=0.0&&inward_coordinate<=1.0, "SurfacePowerevaluationrequiresclosedinwardnodeson[0,1]andmatchingoutputsizes;receivedcoordinate{}and{}nodes.",inward_coordinate,nodes.size()); for(std::size_tnode=1;node<nodes.size();++node){ serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(std::isfinite(nodes[node])&&nodes[node]>nodes[node-1], "SurfacePowerrequiresstrictlyincreasingradialnodes."); } values[0]=1.0; derivatives[0]=0.0; for(std::size_tnode=1;node<nodes.size();++node){ doublecardinal=1.0; doublecardinal_derivative=0.0; for(std::size_tother=1;other<nodes.size();++other){ if(other==node){continue;} constdoubledenominator=nodes[node]-nodes[other]; cardinal_derivative=cardinal_derivative*(inward_coordinate-nodes[other])/denominator+cardinal/denominator; cardinal*=(inward_coordinate-nodes[other])/denominator; } constdoubleratio=inward_coordinate/nodes[node]; constdoublepower=std::pow(ratio,m_exponent); constdoublepower_derivative=m_exponent/nodes[node]*std::pow(ratio,m_exponent-1.0); values[node]=power*cardinal; derivatives[node]=power_derivative*cardinal+power*cardinal_derivative; serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(std::isfinite(values[node])&&std::isfinite(derivatives[node]), "SurfacePowerexponent{}cannotbeevaluatedfinitelyontheseradialnodes.",m_exponent); values[0]-=values[node]; derivatives[0]-=derivatives[node]; } } private: doublem_exponent; }; }