393 lines
46 KiB
XML
393 lines
46 KiB
XML
|
<?xml version='1.0' encoding='UTF-8' standalone='no'?>
|
|||
|
|
<doxygen xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="compound.xsd" version="1.13.2" xml:lang="en-US">
|
|||
|
|
<compounddef id="scaling_8cpp" kind="file" language="C++">
|
|||
|
|
<compoundname>scaling.cpp</compoundname>
|
|||
|
|
<includes refid="scaling_8hpp" local="yes">serif/utils/laneemden/scaling.hpp</includes>
|
|||
|
|
<includes local="no">algorithm</includes>
|
|||
|
|
<includes local="no">cmath</includes>
|
|||
|
|
<includes local="no">numbers</includes>
|
|||
|
|
<incdepgraph>
|
|||
|
|
<node id="6">
|
|||
|
|
<label>serif/dimensions/dimensions.hpp</label>
|
|||
|
|
<link refid="dimensions_8hpp"/>
|
|||
|
|
<childnode refid="7" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="8" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
</node>
|
|||
|
|
<node id="8">
|
|||
|
|
<label>serif/dimensions/quantities.hpp</label>
|
|||
|
|
<link refid="quantities_8hpp"/>
|
|||
|
|
<childnode refid="9" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="10" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="11" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="12" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="13" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="14" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
</node>
|
|||
|
|
<node id="31">
|
|||
|
|
<label>serif/meta/enums.hpp</label>
|
|||
|
|
<link refid="enums_8hpp"/>
|
|||
|
|
<childnode refid="32" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="9" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="33" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
</node>
|
|||
|
|
<node id="30">
|
|||
|
|
<label>serif/utils/error/codes.hpp</label>
|
|||
|
|
<link refid="codes_8hpp"/>
|
|||
|
|
<childnode refid="21" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="31" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
</node>
|
|||
|
|
<node id="19">
|
|||
|
|
<label>serif/utils/error/errors.hpp</label>
|
|||
|
|
<link refid="errors_8hpp"/>
|
|||
|
|
<childnode refid="7" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="18" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="20" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="5" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="21" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="22" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="23" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="24" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="28" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="30" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="34" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
</node>
|
|||
|
|
<node id="2">
|
|||
|
|
<label>serif/utils/laneemden/scaling.hpp</label>
|
|||
|
|
<link refid="scaling_8hpp"/>
|
|||
|
|
<childnode refid="3" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="4" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="5" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="6" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="15" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
</node>
|
|||
|
|
<node id="15">
|
|||
|
|
<label>serif/utils/laneemden/solution.hpp</label>
|
|||
|
|
<link refid="solution_8hpp"/>
|
|||
|
|
<childnode refid="16" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="3" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="4" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="5" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="7" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="17" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="18" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="19" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
</node>
|
|||
|
|
<node id="28">
|
|||
|
|
<label>serif/utils/misc/concepts/exceptions.hpp</label>
|
|||
|
|
<link refid="utils_2misc_2concepts_2exceptions_8hpp"/>
|
|||
|
|
<childnode refid="29" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="9" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="14" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
</node>
|
|||
|
|
<node id="14">
|
|||
|
|
<label>serif/utils/misc/concepts/numeric.hpp</label>
|
|||
|
|
<link refid="numeric_8hpp"/>
|
|||
|
|
<childnode refid="11" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="9" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
</node>
|
|||
|
|
<node id="13">
|
|||
|
|
<label>serif/utils/misc/std_helper/cleaning.hpp</label>
|
|||
|
|
<link refid="cleaning_8hpp"/>
|
|||
|
|
<childnode refid="9" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
</node>
|
|||
|
|
<node id="24">
|
|||
|
|
<label>serif/utils/misc/terminal/colors.hpp</label>
|
|||
|
|
<link refid="colors_8hpp"/>
|
|||
|
|
<childnode refid="7" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="20" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="25" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="26" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="27" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
</node>
|
|||
|
|
<node id="34">
|
|||
|
|
<label>serif/utils/types/strings.hpp</label>
|
|||
|
|
<link refid="types_2strings_8hpp"/>
|
|||
|
|
<childnode refid="22" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="35" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="7" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
</node>
|
|||
|
|
<node id="1">
|
|||
|
|
<label>src/lib/serif/utils/laneemden/scaling.cpp</label>
|
|||
|
|
<link refid="scaling_8cpp"/>
|
|||
|
|
<childnode refid="2" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="22" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="36" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
<childnode refid="37" relation="include">
|
|||
|
|
</childnode>
|
|||
|
|
</node>
|
|||
|
|
<node id="22">
|
|||
|
|
<label>algorithm</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="36">
|
|||
|
|
<label>cmath</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="10">
|
|||
|
|
<label>compare</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="11">
|
|||
|
|
<label>concepts</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="16">
|
|||
|
|
<label>cstddef</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="21">
|
|||
|
|
<label>cstdint</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="20">
|
|||
|
|
<label>cstdio</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="25">
|
|||
|
|
<label>cstdlib</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="29">
|
|||
|
|
<label>exception</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="3">
|
|||
|
|
<label>expected</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="18">
|
|||
|
|
<label>format</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="35">
|
|||
|
|
<label>iostream</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="33">
|
|||
|
|
<label>magic_enum/magic_enum_all.hpp</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="37">
|
|||
|
|
<label>numbers</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="4">
|
|||
|
|
<label>optional</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="5">
|
|||
|
|
<label>string</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="7">
|
|||
|
|
<label>string_view</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="12">
|
|||
|
|
<label>tuple</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="9">
|
|||
|
|
<label>type_traits</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="27">
|
|||
|
|
<label>unistd.h</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="32">
|
|||
|
|
<label>unordered_map</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="26">
|
|||
|
|
<label>utility</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="17">
|
|||
|
|
<label>vector</label>
|
|||
|
|
</node>
|
|||
|
|
<node id="23">
|
|||
|
|
<label>version</label>
|
|||
|
|
</node>
|
|||
|
|
</incdepgraph>
|
|||
|
|
<innernamespace refid="namespaceserif">serif</innernamespace>
|
|||
|
|
<innernamespace refid="namespaceserif_1_1utils">serif::utils</innernamespace>
|
|||
|
|
<innernamespace refid="namespaceserif_1_1utils_1_1laneemden">serif::utils::laneemden</innernamespace>
|
|||
|
|
<briefdescription>
|
|||
|
|
</briefdescription>
|
|||
|
|
<detaileddescription>
|
|||
|
|
</detaileddescription>
|
|||
|
|
<programlisting>
|
|||
|
|
<codeline lineno="1"><highlight class="preprocessor">#include<sp/>"<ref refid="scaling_8hpp" kindref="compound">serif/utils/laneemden/scaling.hpp</ref>"</highlight><highlight class="normal"></highlight></codeline>
|
|||
|
|
<codeline lineno="2"><highlight class="normal"></highlight></codeline>
|
|||
|
|
<codeline lineno="3"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/><algorithm></highlight><highlight class="normal"></highlight></codeline>
|
|||
|
|
<codeline lineno="4"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/><cmath></highlight><highlight class="normal"></highlight></codeline>
|
|||
|
|
<codeline lineno="5"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/><numbers></highlight><highlight class="normal"></highlight></codeline>
|
|||
|
|
<codeline lineno="6"><highlight class="normal"></highlight></codeline>
|
|||
|
|
<codeline lineno="7"><highlight class="normal"></highlight><highlight class="keyword">namespace<sp/></highlight><highlight class="normal"><ref refid="namespaceserif_1_1utils_1_1laneemden" kindref="compound">serif::utils::laneemden</ref><sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="8"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">namespace<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
|||
|
|
<codeline lineno="9"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">constexpr</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">long</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>four_pi<sp/>=<sp/>4.0L<sp/>*<sp/>std::numbers::pi_v<long<sp/>double>;</highlight></codeline>
|
|||
|
|
<codeline lineno="10"><highlight class="normal"></highlight></codeline>
|
|||
|
|
<codeline lineno="12"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">bool</highlight><highlight class="normal"><sp/>positive_finite(</highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>value)<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="13"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::isfinite(value)<sp/>&&<sp/><ref refid="namespaceserif_1_1discretization_1_1backend_1_1mfem_1af75f98db90618153ea4ca169945733f7a2063c1608d6e0baf80249c42e2be5804" kindref="member">value</ref><sp/>><sp/>0.0;</highlight></codeline>
|
|||
|
|
<codeline lineno="14"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="15"><highlight class="normal"></highlight></codeline>
|
|||
|
|
<codeline lineno="17"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">bool</highlight><highlight class="normal"><sp/>finite_surface(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="structserif_1_1utils_1_1laneemden_1_1_solution" kindref="compound">Solution</ref><sp/>&solution)<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="18"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::isfinite(solution.n)<sp/>&&<sp/>solution.n<sp/>>=<sp/>0.0<sp/>&&<sp/>solution.surface<sp/>&&<sp/>positive_finite(solution.surface->xi)<sp/>&&</highlight></codeline>
|
|||
|
|
<codeline lineno="19"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::isfinite(solution.surface->derivative)<sp/>&&<sp/>solution.surface->derivative<sp/><<sp/>0.0;</highlight></codeline>
|
|||
|
|
<codeline lineno="20"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="21"><highlight class="normal"></highlight></codeline>
|
|||
|
|
<codeline lineno="23"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::expected<PhysicalScales,<sp/>std::string><sp/>make_scales(</highlight></codeline>
|
|||
|
|
<codeline lineno="24"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="structserif_1_1utils_1_1laneemden_1_1_solution" kindref="compound">Solution</ref><sp/>&solution,<sp/></highlight><highlight class="keywordtype">long</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>radial_scale,<sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>central_density,<sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>gravitational_constant)<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="25"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">long</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>surface_xi<sp/>=<sp/>solution.surface->xi;</highlight></codeline>
|
|||
|
|
<codeline lineno="26"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">long</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>dimensionless_mass<sp/>=<sp/>-surface_xi<sp/>*<sp/>surface_xi<sp/>*<sp/>solution.surface->derivative;</highlight></codeline>
|
|||
|
|
<codeline lineno="27"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>radius<sp/>=<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(radial_scale<sp/>*<sp/>surface_xi);</highlight></codeline>
|
|||
|
|
<codeline lineno="28"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>mass<sp/>=<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(four_pi<sp/>*<sp/>central_density<sp/>*<sp/>radial_scale<sp/>*<sp/>radial_scale<sp/>*<sp/>radial_scale<sp/>*<sp/>dimensionless_mass);</highlight></codeline>
|
|||
|
|
<codeline lineno="29"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>enthalpy<sp/>=<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(four_pi<sp/>*<sp/>gravitational_constant<sp/>*<sp/>central_density<sp/>*<sp/>radial_scale<sp/>*<sp/>radial_scale);</highlight></codeline>
|
|||
|
|
<codeline lineno="30"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>alpha<sp/>=<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(radial_scale);</highlight></codeline>
|
|||
|
|
<codeline lineno="31"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!positive_finite(alpha)<sp/>||<sp/>!positive_finite(radius)<sp/>||<sp/>!positive_finite(mass)<sp/>||<sp/>!positive_finite(enthalpy))<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="32"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(</highlight><highlight class="stringliteral">"Lane–Emden<sp/>physical<sp/>scales<sp/>overflow<sp/>or<sp/>underflow;<sp/>choose<sp/>representable<sp/>consistent<sp/>units."</highlight><highlight class="normal">);</highlight></codeline>
|
|||
|
|
<codeline lineno="33"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="34"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::optional<double><sp/>polytropic_constant;</highlight></codeline>
|
|||
|
|
<codeline lineno="35"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(solution.n<sp/>><sp/>0.0)<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="36"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>constant<sp/>=<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(enthalpy<sp/>/<sp/>((solution.n<sp/>+<sp/>1.0L)<sp/>*<sp/>std::pow(</highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">long</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(central_density),<sp/>1.0L<sp/>/<sp/>solution.n)));</highlight></codeline>
|
|||
|
|
<codeline lineno="37"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(positive_finite(constant))<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="38"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>polytropic_constant<sp/>=<sp/>constant;</highlight></codeline>
|
|||
|
|
<codeline lineno="39"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="40"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="41"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="structserif_1_1utils_1_1laneemden_1_1_physical_scales" kindref="compound">PhysicalScales</ref>{</highlight></codeline>
|
|||
|
|
<codeline lineno="42"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.n<sp/>=<sp/>solution.n,</highlight></codeline>
|
|||
|
|
<codeline lineno="43"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.radial_scale<sp/>=<sp/>alpha,</highlight></codeline>
|
|||
|
|
<codeline lineno="44"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.mass<sp/>=<sp/>dimensions::QuantityValue<dimensions::Mass>{mass},</highlight></codeline>
|
|||
|
|
<codeline lineno="45"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.radius<sp/>=<sp/>dimensions::QuantityValue<dimensions::Length>{radius},</highlight></codeline>
|
|||
|
|
<codeline lineno="46"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.central_density<sp/>=<sp/>dimensions::QuantityValue<dimensions::MassDensity>{central_density},</highlight></codeline>
|
|||
|
|
<codeline lineno="47"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.central_enthalpy<sp/>=<sp/>dimensions::QuantityValue<dimensions::SpecificEnthalpy>{enthalpy},</highlight></codeline>
|
|||
|
|
<codeline lineno="48"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.gravitational_constant<sp/>=<sp/>dimensions::QuantityValue<dimensions::GravitationalConstant>{gravitational_constant},</highlight></codeline>
|
|||
|
|
<codeline lineno="49"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.polytropic_constant<sp/>=<sp/>polytropic_constant</highlight></codeline>
|
|||
|
|
<codeline lineno="50"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
|
|||
|
|
<codeline lineno="51"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="52"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="53"><highlight class="normal"></highlight></codeline>
|
|||
|
|
<codeline lineno="54"><highlight class="normal"><sp/><sp/><sp/><sp/>std::expected<PhysicalScales,<sp/>std::string><sp/><ref refid="namespaceserif_1_1utils_1_1laneemden_1a4c3ea7fb533fd5a9306e993593eb20fc" kindref="member">scale_from_mass</ref>(</highlight></codeline>
|
|||
|
|
<codeline lineno="55"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="structserif_1_1utils_1_1laneemden_1_1_solution" kindref="compound">Solution</ref><sp/>&solution,</highlight></codeline>
|
|||
|
|
<codeline lineno="56"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>dimensions::QuantityValue<dimensions::Mass><sp/>mass,</highlight></codeline>
|
|||
|
|
<codeline lineno="57"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>dimensions::QuantityValue<dimensions::MassDensity><sp/>central_density,</highlight></codeline>
|
|||
|
|
<codeline lineno="58"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>dimensions::QuantityValue<dimensions::GravitationalConstant><sp/>gravitational_constant)<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="59"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!finite_surface(solution)<sp/>||<sp/>!positive_finite(mass.value())<sp/>||<sp/>!positive_finite(central_density.value())<sp/>||<sp/>!positive_finite(gravitational_constant.value()))<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="60"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(</highlight><highlight class="stringliteral">"Lane–Emden<sp/>mass<sp/>scaling<sp/>needs<sp/>a<sp/>finite<sp/>surface<sp/>and<sp/>positive<sp/>finite<sp/>mass,<sp/>central<sp/>density,<sp/>and<sp/>G."</highlight><highlight class="normal">);</highlight></codeline>
|
|||
|
|
<codeline lineno="61"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="62"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">long</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>xi<sp/>=<sp/>solution.surface->xi;</highlight></codeline>
|
|||
|
|
<codeline lineno="63"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">long</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>dimensionless_mass<sp/>=<sp/>-xi<sp/>*<sp/>xi<sp/>*<sp/>solution.surface->derivative;</highlight></codeline>
|
|||
|
|
<codeline lineno="64"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">long</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>alpha<sp/>=<sp/>std::cbrt(mass.value()<sp/>/<sp/>(four_pi<sp/>*<sp/>central_density.value()<sp/>*<sp/>dimensionless_mass));</highlight></codeline>
|
|||
|
|
<codeline lineno="65"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>make_scales(solution,<sp/>alpha,<sp/>central_density.value(),<sp/>gravitational_constant.value());</highlight></codeline>
|
|||
|
|
<codeline lineno="66"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="67"><highlight class="normal"></highlight></codeline>
|
|||
|
|
<codeline lineno="68"><highlight class="normal"><sp/><sp/><sp/><sp/>std::expected<PhysicalScales,<sp/>std::string><sp/><ref refid="namespaceserif_1_1utils_1_1laneemden_1ab57cad62436b0fcf5cb9df68347b196b" kindref="member">scale_from_enthalpy</ref>(</highlight></codeline>
|
|||
|
|
<codeline lineno="69"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="structserif_1_1utils_1_1laneemden_1_1_solution" kindref="compound">Solution</ref><sp/>&solution,</highlight></codeline>
|
|||
|
|
<codeline lineno="70"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>dimensions::QuantityValue<dimensions::SpecificEnthalpy><sp/>central_enthalpy,</highlight></codeline>
|
|||
|
|
<codeline lineno="71"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>dimensions::QuantityValue<dimensions::MassDensity><sp/>central_density,</highlight></codeline>
|
|||
|
|
<codeline lineno="72"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>dimensions::QuantityValue<dimensions::GravitationalConstant><sp/>gravitational_constant)<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="73"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!finite_surface(solution)<sp/>||<sp/>!positive_finite(central_enthalpy.value())<sp/>||<sp/>!positive_finite(central_density.value())<sp/>||<sp/>!positive_finite(gravitational_constant.value()))<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="74"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(</highlight><highlight class="stringliteral">"Lane–Emden<sp/>enthalpy<sp/>scaling<sp/>needs<sp/>a<sp/>finite<sp/>surface<sp/>and<sp/>positive<sp/>finite<sp/>enthalpy,<sp/>central<sp/>density,<sp/>and<sp/>G."</highlight><highlight class="normal">);</highlight></codeline>
|
|||
|
|
<codeline lineno="75"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="76"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">long</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>alpha<sp/>=<sp/>std::sqrt(central_enthalpy.value()<sp/>/<sp/>(four_pi<sp/>*<sp/>gravitational_constant.value()<sp/>*<sp/>central_density.value()));</highlight></codeline>
|
|||
|
|
<codeline lineno="77"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>make_scales(solution,<sp/>alpha,<sp/>central_density.value(),<sp/>gravitational_constant.value());</highlight></codeline>
|
|||
|
|
<codeline lineno="78"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="79"><highlight class="normal"></highlight></codeline>
|
|||
|
|
<codeline lineno="80"><highlight class="normal"><sp/><sp/><sp/><sp/>std::expected<PhysicalScales,<sp/>std::string><sp/><ref refid="namespaceserif_1_1utils_1_1laneemden_1adf057a3c9420e9b4e92c20c89126d3fd" kindref="member">scale_from_polytropic_constant</ref>(</highlight></codeline>
|
|||
|
|
<codeline lineno="81"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="structserif_1_1utils_1_1laneemden_1_1_solution" kindref="compound">Solution</ref><sp/>&solution,<sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>polytropic_constant,</highlight></codeline>
|
|||
|
|
<codeline lineno="82"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>dimensions::QuantityValue<dimensions::MassDensity><sp/>central_density,</highlight></codeline>
|
|||
|
|
<codeline lineno="83"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>dimensions::QuantityValue<dimensions::GravitationalConstant><sp/>gravitational_constant)<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="84"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!finite_surface(solution)<sp/>||<sp/>solution.n<sp/><=<sp/>0.0<sp/>||<sp/>!positive_finite(polytropic_constant)<sp/>||<sp/>!positive_finite(central_density.value())<sp/>||</highlight></codeline>
|
|||
|
|
<codeline lineno="85"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>!positive_finite(gravitational_constant.value()))<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="86"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(</highlight><highlight class="stringliteral">"Lane–Emden<sp/>EOS<sp/>scaling<sp/>needs<sp/>n<sp/>><sp/>0,<sp/>a<sp/>finite<sp/>surface,<sp/>and<sp/>positive<sp/>finite<sp/>K,<sp/>central<sp/>density,<sp/>and<sp/>G."</highlight><highlight class="normal">);</highlight></codeline>
|
|||
|
|
<codeline lineno="87"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="88"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>enthalpy<sp/>=<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">((solution.n<sp/>+<sp/>1.0L)<sp/>*<sp/>polytropic_constant<sp/>*<sp/>std::pow(</highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">long</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(central_density.value()),<sp/>1.0L<sp/>/<sp/>solution.n));</highlight></codeline>
|
|||
|
|
<codeline lineno="89"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>scales<sp/>=<sp/><ref refid="namespaceserif_1_1utils_1_1laneemden_1ab57cad62436b0fcf5cb9df68347b196b" kindref="member">scale_from_enthalpy</ref>(solution,<sp/>dimensions::QuantityValue<dimensions::SpecificEnthalpy>{enthalpy},<sp/>central_density,<sp/>gravitational_constant);</highlight></codeline>
|
|||
|
|
<codeline lineno="90"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(scales)<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="91"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>scales->polytropic_constant<sp/>=<sp/>polytropic_constant;</highlight></codeline>
|
|||
|
|
<codeline lineno="92"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="93"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>scales;</highlight></codeline>
|
|||
|
|
<codeline lineno="94"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="95"><highlight class="normal"></highlight></codeline>
|
|||
|
|
<codeline lineno="96"><highlight class="normal"><sp/><sp/><sp/><sp/>std::expected<PhysicalSample,<sp/>std::string><sp/><ref refid="namespaceserif_1_1utils_1_1laneemden_1aec2842a8929dc7a2b5dac2adf8a067b6" kindref="member">evaluate_physical</ref>(</highlight></codeline>
|
|||
|
|
<codeline lineno="97"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="structserif_1_1utils_1_1laneemden_1_1_solution" kindref="compound">Solution</ref><sp/>&solution,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="structserif_1_1utils_1_1laneemden_1_1_physical_scales" kindref="compound">PhysicalScales</ref><sp/>&scales,<sp/>dimensions::QuantityValue<dimensions::Length><sp/>radius)<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="98"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!finite_surface(solution)<sp/>||<sp/>solution.n<sp/>!=<sp/>scales.n<sp/>||<sp/>!std::isfinite(radius.value())<sp/>||<sp/>radius.value()<sp/><<sp/>0.0<sp/>||</highlight></codeline>
|
|||
|
|
<codeline lineno="99"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>!positive_finite(scales.radial_scale)<sp/>||<sp/>!positive_finite(scales.radius.value())<sp/>||<sp/>!positive_finite(scales.mass.value())<sp/>||</highlight></codeline>
|
|||
|
|
<codeline lineno="100"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>!positive_finite(scales.central_density.value())<sp/>||<sp/>!positive_finite(scales.central_enthalpy.value())<sp/>||<sp/>!positive_finite(scales.gravitational_constant.value()))<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="101"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(</highlight><highlight class="stringliteral">"Lane–Emden<sp/>physical<sp/>evaluation<sp/>requires<sp/>matching<sp/>finite-surface<sp/>scales<sp/>and<sp/>a<sp/>finite<sp/>nonnegative<sp/>radius."</highlight><highlight class="normal">);</highlight></codeline>
|
|||
|
|
<codeline lineno="102"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="103"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">long</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>gravitational_mass<sp/>=<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">long</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(scales.gravitational_constant.value())<sp/>*<sp/>scales.mass.value();</highlight></codeline>
|
|||
|
|
<codeline lineno="104"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>density<sp/>=<sp/>0.0;</highlight></codeline>
|
|||
|
|
<codeline lineno="105"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>pressure<sp/>=<sp/>0.0;</highlight></codeline>
|
|||
|
|
<codeline lineno="106"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>enthalpy<sp/>=<sp/>0.0;</highlight></codeline>
|
|||
|
|
<codeline lineno="107"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>mass<sp/>=<sp/>scales.mass.value();</highlight></codeline>
|
|||
|
|
<codeline lineno="108"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>potential;</highlight></codeline>
|
|||
|
|
<codeline lineno="109"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>gravity;</highlight></codeline>
|
|||
|
|
<codeline lineno="110"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(radius.value()<sp/>>=<sp/>scales.radius.value())<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="111"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>potential<sp/>=<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(-gravitational_mass<sp/>/<sp/>radius.value());</highlight></codeline>
|
|||
|
|
<codeline lineno="112"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>gravity<sp/>=<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(gravitational_mass<sp/>/<sp/>radius.value()<sp/>/<sp/>radius.value());</highlight></codeline>
|
|||
|
|
<codeline lineno="113"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}<sp/></highlight><highlight class="keywordflow">else</highlight><highlight class="normal"><sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="114"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>xi<sp/>=<sp/>std::min(radius.value()<sp/>/<sp/>scales.radial_scale,<sp/>solution.surface->xi);</highlight></codeline>
|
|||
|
|
<codeline lineno="115"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>sample<sp/>=<sp/><ref refid="namespaceserif_1_1utils_1_1laneemden_1ac58bd2b83384f8454be7fd0b11cd9332" kindref="member">evaluate</ref>(solution,<sp/>xi);</highlight></codeline>
|
|||
|
|
<codeline lineno="116"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!sample)<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="117"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(sample.error());</highlight></codeline>
|
|||
|
|
<codeline lineno="118"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="119"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>theta<sp/>=<sp/>std::max(0.0,<sp/>sample->theta);</highlight></codeline>
|
|||
|
|
<codeline lineno="120"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>density<sp/>=<sp/>scales.central_density.value()<sp/>*<sp/>std::pow(theta,<sp/>solution.n);</highlight></codeline>
|
|||
|
|
<codeline lineno="121"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>enthalpy<sp/>=<sp/>scales.central_enthalpy.value()<sp/>*<sp/>theta;</highlight></codeline>
|
|||
|
|
<codeline lineno="122"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>pressure<sp/>=<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(</highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">long</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(density)<sp/>*<sp/>enthalpy<sp/>/<sp/>(solution.n<sp/>+<sp/>1.0L));</highlight></codeline>
|
|||
|
|
<codeline lineno="123"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>potential<sp/>=<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(-gravitational_mass<sp/>/<sp/>scales.radius.value()<sp/>-<sp/>enthalpy);</highlight></codeline>
|
|||
|
|
<codeline lineno="124"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>gravity<sp/>=<sp/>-scales.central_enthalpy.value()<sp/>/<sp/>scales.radial_scale<sp/>*<sp/>sample->derivative;</highlight></codeline>
|
|||
|
|
<codeline lineno="125"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mass<sp/>=<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(-four_pi<sp/>*<sp/>scales.central_density.value()<sp/>*<sp/>std::pow(</highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">long</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(scales.radial_scale),<sp/>3)<sp/>*<sp/>xi<sp/>*<sp/>xi<sp/>*<sp/>sample->derivative);</highlight></codeline>
|
|||
|
|
<codeline lineno="126"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="127"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!std::isfinite(density)<sp/>||<sp/>!std::isfinite(pressure)<sp/>||<sp/>!std::isfinite(enthalpy)<sp/>||<sp/>!std::isfinite(potential)<sp/>||<sp/>!std::isfinite(gravity)<sp/>||<sp/>!std::isfinite(mass))<sp/>{</highlight></codeline>
|
|||
|
|
<codeline lineno="128"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(</highlight><highlight class="stringliteral">"Lane–Emden<sp/>physical<sp/>fields<sp/>overflow<sp/>in<sp/>the<sp/>supplied<sp/>units."</highlight><highlight class="normal">);</highlight></codeline>
|
|||
|
|
<codeline lineno="129"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="130"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="structserif_1_1utils_1_1laneemden_1_1_physical_sample" kindref="compound">PhysicalSample</ref>{</highlight></codeline>
|
|||
|
|
<codeline lineno="131"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.density<sp/>=<sp/>dimensions::QuantityValue<dimensions::MassDensity>{density},</highlight></codeline>
|
|||
|
|
<codeline lineno="132"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.pressure<sp/>=<sp/>dimensions::QuantityValue<dimensions::Pressure>{pressure},</highlight></codeline>
|
|||
|
|
<codeline lineno="133"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.enthalpy<sp/>=<sp/>dimensions::QuantityValue<dimensions::SpecificEnthalpy>{enthalpy},</highlight></codeline>
|
|||
|
|
<codeline lineno="134"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.potential<sp/>=<sp/>dimensions::QuantityValue<dimensions::GravitationalPotential>{potential},</highlight></codeline>
|
|||
|
|
<codeline lineno="135"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.enclosed_mass<sp/>=<sp/>dimensions::QuantityValue<dimensions::Mass>{mass},</highlight></codeline>
|
|||
|
|
<codeline lineno="136"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.gravity<sp/>=<sp/>dimensions::QuantityValue<dimensions::Acceleration>{gravity}</highlight></codeline>
|
|||
|
|
<codeline lineno="137"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
|
|||
|
|
<codeline lineno="138"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
|||
|
|
<codeline lineno="139"><highlight class="normal">}</highlight></codeline>
|
|||
|
|
</programlisting>
|
|||
|
|
<location file="src/lib/serif/utils/laneemden/scaling.cpp"/>
|
|||
|
|
</compounddef>
|
|||
|
|
</doxygen>
|