MeanField_SERiF/docs/xml/polytropic_8cpp.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="polytropic_8cpp" kind="file" language="C++">
    <compoundname>polytropic.cpp</compoundname>
    <includes refid="eos_2models_2polytropic_8hpp" local="yes">serif/eos/models/polytropic.hpp</includes>
    <includes refid="utils_2misc_2finite_8hpp" local="yes">serif/utils/misc/finite.hpp</includes>
    <includes refid="eos_2exceptions_8hpp" local="yes">serif/eos/exceptions.hpp</includes>
    <includes refid="codes_8hpp" local="yes">serif/utils/error/codes.hpp</includes>
    <incdepgraph>
      <node id="17">
        <label>serif/dimensions/dimensions.hpp</label>
        <link refid="dimensions_8hpp"/>
        <childnode refid="18" relation="include">
        </childnode>
        <childnode refid="12" relation="include">
        </childnode>
      </node>
      <node id="20">
        <label>serif/dimensions/partial.hpp</label>
        <link refid="partial_8hpp"/>
        <childnode refid="12" relation="include">
        </childnode>
        <childnode refid="16" relation="include">
        </childnode>
      </node>
      <node id="12">
        <label>serif/dimensions/quantities.hpp</label>
        <link refid="quantities_8hpp"/>
        <childnode refid="9" relation="include">
        </childnode>
        <childnode refid="13" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
        <childnode refid="15" relation="include">
        </childnode>
        <childnode refid="8" relation="include">
        </childnode>
        <childnode refid="16" relation="include">
        </childnode>
      </node>
      <node id="19">
        <label>serif/dimensions/type_alias.hpp</label>
        <link refid="type__alias_8hpp"/>
        <childnode refid="12" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
      </node>
      <node id="43">
        <label>serif/eos/exceptions.hpp</label>
        <link refid="eos_2exceptions_8hpp"/>
        <childnode refid="30" relation="include">
        </childnode>
        <childnode refid="23" relation="include">
        </childnode>
        <childnode refid="29" relation="include">
        </childnode>
      </node>
      <node id="2">
        <label>serif/eos/models/polytropic.hpp</label>
        <link refid="eos_2models_2polytropic_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="19" relation="include">
        </childnode>
        <childnode refid="20" relation="include">
        </childnode>
      </node>
      <node id="4">
        <label>serif/eos/relations/catalog.hpp</label>
        <link refid="catalog_8hpp"/>
        <childnode refid="5" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="12" relation="include">
        </childnode>
      </node>
      <node id="3">
        <label>serif/eos/relations/relations.hpp</label>
        <link refid="eos_2relations_2relations_8hpp"/>
        <childnode refid="4" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
      </node>
      <node id="38">
        <label>serif/meta/enums.hpp</label>
        <link refid="enums_8hpp"/>
        <childnode refid="39" relation="include">
        </childnode>
        <childnode refid="9" relation="include">
        </childnode>
        <childnode refid="40" relation="include">
        </childnode>
      </node>
      <node id="37">
        <label>serif/utils/error/codes.hpp</label>
        <link refid="codes_8hpp"/>
        <childnode refid="30" relation="include">
        </childnode>
        <childnode refid="38" relation="include">
        </childnode>
      </node>
      <node id="27">
        <label>serif/utils/error/errors.hpp</label>
        <link refid="errors_8hpp"/>
        <childnode refid="18" relation="include">
        </childnode>
        <childnode refid="24" relation="include">
        </childnode>
        <childnode refid="28" relation="include">
        </childnode>
        <childnode refid="29" relation="include">
        </childnode>
        <childnode refid="30" relation="include">
        </childnode>
        <childnode refid="31" relation="include">
        </childnode>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="33" relation="include">
        </childnode>
        <childnode refid="25" relation="include">
        </childnode>
        <childnode refid="37" relation="include">
        </childnode>
        <childnode refid="41" relation="include">
        </childnode>
      </node>
      <node id="25">
        <label>serif/utils/misc/concepts/exceptions.hpp</label>
        <link refid="utils_2misc_2concepts_2exceptions_8hpp"/>
        <childnode refid="26" relation="include">
        </childnode>
        <childnode refid="9" relation="include">
        </childnode>
        <childnode refid="16" relation="include">
        </childnode>
      </node>
      <node id="16">
        <label>serif/utils/misc/concepts/numeric.hpp</label>
        <link refid="numeric_8hpp"/>
        <childnode refid="14" relation="include">
        </childnode>
        <childnode refid="9" relation="include">
        </childnode>
      </node>
      <node id="5">
        <label>serif/utils/misc/concepts/type_uniqueness.hpp</label>
        <link refid="concepts_2type__uniqueness_8hpp"/>
        <childnode refid="6" relation="include">
        </childnode>
      </node>
      <node id="21">
        <label>serif/utils/misc/finite.hpp</label>
        <link refid="utils_2misc_2finite_8hpp"/>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="23" relation="include">
        </childnode>
        <childnode refid="24" relation="include">
        </childnode>
        <childnode refid="16" relation="include">
        </childnode>
        <childnode refid="25" relation="include">
        </childnode>
        <childnode refid="27" relation="include">
        </childnode>
      </node>
      <node id="8">
        <label>serif/utils/misc/std_helper/cleaning.hpp</label>
        <link refid="cleaning_8hpp"/>
        <childnode refid="9" relation="include">
        </childnode>
      </node>
      <node id="10">
        <label>serif/utils/misc/std_helper/variadic_cleaning.hpp</label>
        <link refid="variadic__cleaning_8hpp"/>
        <childnode refid="11" relation="include">
        </childnode>
      </node>
      <node id="33">
        <label>serif/utils/misc/terminal/colors.hpp</label>
        <link refid="colors_8hpp"/>
        <childnode refid="18" relation="include">
        </childnode>
        <childnode refid="28" relation="include">
        </childnode>
        <childnode refid="34" relation="include">
        </childnode>
        <childnode refid="35" relation="include">
        </childnode>
        <childnode refid="36" relation="include">
        </childnode>
      </node>
      <node id="6">
        <label>serif/utils/misc/traits/type_uniqueness.hpp</label>
        <link refid="traits_2type__uniqueness_8hpp"/>
        <childnode refid="7" relation="include">
        </childnode>
        <childnode refid="8" relation="include">
        </childnode>
      </node>
      <node id="41">
        <label>serif/utils/types/strings.hpp</label>
        <link refid="types_2strings_8hpp"/>
        <childnode refid="31" relation="include">
        </childnode>
        <childnode refid="42" relation="include">
        </childnode>
        <childnode refid="18" relation="include">
        </childnode>
      </node>
      <node id="1">
        <label>src/lib/serif/eos/models/polytropic.cpp</label>
        <link refid="polytropic_8cpp"/>
        <childnode refid="2" relation="include">
        </childnode>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="43" relation="include">
        </childnode>
        <childnode refid="37" relation="include">
        </childnode>
      </node>
      <node id="31">
        <label>algorithm</label>
      </node>
      <node id="7">
        <label>array</label>
      </node>
      <node id="22">
        <label>cmath</label>
      </node>
      <node id="13">
        <label>compare</label>
      </node>
      <node id="14">
        <label>concepts</label>
      </node>
      <node id="11">
        <label>cstddef</label>
      </node>
      <node id="30">
        <label>cstdint</label>
      </node>
      <node id="28">
        <label>cstdio</label>
      </node>
      <node id="34">
        <label>cstdlib</label>
      </node>
      <node id="26">
        <label>exception</label>
      </node>
      <node id="24">
        <label>format</label>
      </node>
      <node id="42">
        <label>iostream</label>
      </node>
      <node id="40">
        <label>magic_enum/magic_enum_all.hpp</label>
      </node>
      <node id="23">
        <label>stdexcept</label>
      </node>
      <node id="29">
        <label>string</label>
      </node>
      <node id="18">
        <label>string_view</label>
      </node>
      <node id="15">
        <label>tuple</label>
      </node>
      <node id="9">
        <label>type_traits</label>
      </node>
      <node id="36">
        <label>unistd.h</label>
      </node>
      <node id="39">
        <label>unordered_map</label>
      </node>
      <node id="35">
        <label>utility</label>
      </node>
      <node id="32">
        <label>version</label>
      </node>
    </incdepgraph>
    <innernamespace refid="namespaceserif">serif</innernamespace>
    <innernamespace refid="namespaceserif_1_1eos">serif::eos</innernamespace>
    <innernamespace refid="namespaceserif_1_1eos_1_1models">serif::eos::models</innernamespace>
    <briefdescription>
    </briefdescription>
    <detaileddescription>
    </detaileddescription>
    <programlisting>
<codeline lineno="1"><highlight class="preprocessor">#include<sp/>&quot;<ref refid="eos_2models_2polytropic_8hpp" kindref="compound">serif/eos/models/polytropic.hpp</ref>&quot;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="2"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&quot;<ref refid="utils_2misc_2finite_8hpp" kindref="compound">serif/utils/misc/finite.hpp</ref>&quot;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="3"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&quot;<ref refid="eos_2exceptions_8hpp" kindref="compound">serif/eos/exceptions.hpp</ref>&quot;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="4"><highlight class="normal"></highlight></codeline>
<codeline lineno="5"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&quot;<ref refid="codes_8hpp" kindref="compound">serif/utils/error/codes.hpp</ref>&quot;</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_1eos_1_1models" kindref="compound">serif::eos::models</ref><sp/>{</highlight></codeline>
<codeline lineno="8"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/></highlight><highlight class="normal"><ref refid="namespaceserif_1_1utils_1_1misc_1a81150dfc59c33d2e767a9e572cb15fce" kindref="member">utils::misc::validate_finite</ref>;</highlight></codeline>
<codeline lineno="9"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/></highlight><highlight class="normal"><ref refid="namespaceserif_1_1utils_1_1misc_1a21655a5b8f48ae49f385e0a7ea836f3a" kindref="member">utils::misc::validate_greater_than</ref>;</highlight></codeline>
<codeline lineno="10"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/></highlight><highlight class="normal"><ref refid="namespaceserif_1_1utils_1_1misc_1ade362937d8732832b60a4d829b15d611" kindref="member">utils::misc::validate_less_than</ref>;</highlight></codeline>
<codeline lineno="11"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/></highlight><highlight class="normal"><ref refid="namespaceserif_1_1utils_1_1misc_1a7f0bf6f20695cc8bb3c052c8a8f16dff" kindref="member">utils::misc::validate_nonnegativity</ref>;</highlight></codeline>
<codeline lineno="12"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/></highlight><highlight class="normal"><ref refid="namespaceserif_1_1utils_1_1errors_1ab04484c13887f78337c2331fb2911346" kindref="member">utils::errors::SERiFErrorCode</ref>;</highlight></codeline>
<codeline lineno="13"><highlight class="normal"></highlight></codeline>
<codeline lineno="14"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/>namespace<sp/></highlight><highlight class="normal">dimensions;</highlight></codeline>
<codeline lineno="15"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/>namespace<sp/></highlight><highlight class="normal">relations;</highlight></codeline>
<codeline lineno="16"><highlight class="normal"></highlight></codeline>
<codeline lineno="17"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="classserif_1_1eos_1_1models_1_1_polytrope_1a01becccced132cf0db4681296ccf6104" kindref="member">Polytrope::Polytrope</ref>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="structserif_1_1eos_1_1models_1_1_polytrope_parameters" kindref="compound">PolytropeParameters</ref>&amp;<sp/>parameters)<sp/>:<sp/><ref refid="classserif_1_1eos_1_1models_1_1_polytrope" kindref="compound">Polytrope</ref>(parameters.n,<sp/>parameters.K)<sp/>{}</highlight></codeline>
<codeline lineno="18"><highlight class="normal"></highlight></codeline>
<codeline lineno="19"><highlight class="normal"><sp/><sp/><sp/><sp/>Polytrope::Polytrope(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>n,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>K):</highlight></codeline>
<codeline lineno="20"><highlight class="normal"><sp/><sp/><sp/><sp/>m_polytropic_index(n),</highlight></codeline>
<codeline lineno="21"><highlight class="normal"><sp/><sp/><sp/><sp/>m_polytropic_constant(K),</highlight></codeline>
<codeline lineno="22"><highlight class="normal"><sp/><sp/><sp/><sp/>m_enthalpy_scale((n<sp/>+<sp/>1.0)<sp/>*<sp/>K)<sp/>{</highlight></codeline>
<codeline lineno="23"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>validate_finite&lt;EOSEvaluationError,<sp/>SERiFErrorCode::EOS_OUTSIDE_DOMAIN&gt;(m_polytropic_index);</highlight></codeline>
<codeline lineno="24"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>validate_greater_than&lt;EOSEvaluationError,<sp/>SERiFErrorCode::EOS_OUTSIDE_DOMAIN&gt;(m_polytropic_index,<sp/>1.0);</highlight></codeline>
<codeline lineno="25"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>validate_finite&lt;EOSEvaluationError,<sp/>SERiFErrorCode::EOS_OUTSIDE_DOMAIN&gt;(m_polytropic_constant);</highlight></codeline>
<codeline lineno="26"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>utils::errors::serif_ensure&lt;SERiFErrorCode::EOS_OUTSIDE_DOMAIN,<sp/>EOSEvaluationError&gt;(m_polytropic_constant<sp/>&gt;<sp/>0.0,<sp/></highlight><highlight class="stringliteral">&quot;The<sp/>polytropic<sp/>constant<sp/>must<sp/>be<sp/>positive.<sp/>Received<sp/>{}.&quot;</highlight><highlight class="normal">,<sp/>m_polytropic_constant);</highlight></codeline>
<codeline lineno="27"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="28"><highlight class="normal"></highlight></codeline>
<codeline lineno="29"><highlight class="normal"><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>Polytrope::polytropic_index()<sp/>const<sp/>noexcept<sp/>{</highlight></codeline>
<codeline lineno="30"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>m_polytropic_index;</highlight></codeline>
<codeline lineno="31"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="32"><highlight class="normal"></highlight></codeline>
<codeline lineno="33"><highlight class="normal"><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>Polytrope::polytropic_constant()<sp/>const<sp/>noexcept<sp/>{</highlight></codeline>
<codeline lineno="34"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>m_polytropic_constant;</highlight></codeline>
<codeline lineno="35"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="36"><highlight class="normal"></highlight></codeline>
<codeline lineno="37"><highlight class="normal"><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>Polytrope::enthalpy_scale()<sp/>const<sp/>noexcept<sp/>{</highlight></codeline>
<codeline lineno="38"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>m_enthalpy_scale;</highlight></codeline>
<codeline lineno="39"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="40"><highlight class="normal"></highlight></codeline>
<codeline lineno="41"><highlight class="normal"><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/><ref refid="namespaceserif_1_1dimensions_1a60ad5bcd33f38e5e0e272d9df76008ec" kindref="member">PressureValue</ref><sp/>Polytrope::evaluate(PressureFromDensity,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>DensityValue<sp/>density)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="42"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>validate_nonnegativity&lt;EOSEvaluationError,<sp/>SERiFErrorCode::EOS_INVALID_THERMODYNAMIC_INPUT&gt;(density.value());</highlight></codeline>
<codeline lineno="43"><highlight class="normal"></highlight></codeline>
<codeline lineno="44"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(density.value()<sp/>==<sp/>0.0)<sp/>{</highlight></codeline>
<codeline lineno="45"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>P<sp/>=<sp/>0</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="46"><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="namespaceserif_1_1dimensions_1a60ad5bcd33f38e5e0e272d9df76008ec" kindref="member">PressureValue</ref>{0.0};</highlight></codeline>
<codeline lineno="47"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="48"><highlight class="normal"></highlight></codeline>
<codeline lineno="49"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>P<sp/>=<sp/>K<sp/>*<sp/>(ρ)^(1<sp/>+<sp/>1/n)</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="50"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1dimensions_1a60ad5bcd33f38e5e0e272d9df76008ec" kindref="member">PressureValue</ref>{</highlight></codeline>
<codeline lineno="51"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_polytropic_constant<sp/>*<sp/>std::pow(density.value(),<sp/>1.0<sp/>+<sp/>1.0<sp/>/<sp/>m_polytropic_index)</highlight></codeline>
<codeline lineno="52"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="53"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="54"><highlight class="normal"></highlight></codeline>
<codeline lineno="55"><highlight class="normal"><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/><ref refid="namespaceserif_1_1operators_1a84d7e8ff528dc3361769cc6f08e6cd49" kindref="member">DensityValue</ref><sp/>Polytrope::evaluate(DensityFromPressure,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>PressureValue<sp/>pressure)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="56"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>validate_nonnegativity&lt;EOSEvaluationError,<sp/>SERiFErrorCode::EOS_INVALID_THERMODYNAMIC_INPUT&gt;(pressure.value());</highlight></codeline>
<codeline lineno="57"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(pressure.value()<sp/>==<sp/>0.0)<sp/>{</highlight></codeline>
<codeline lineno="58"><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="namespaceserif_1_1operators_1a84d7e8ff528dc3361769cc6f08e6cd49" kindref="member">DensityValue</ref>{0.0};</highlight></codeline>
<codeline lineno="59"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="60"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1operators_1a84d7e8ff528dc3361769cc6f08e6cd49" kindref="member">DensityValue</ref>{std::pow(pressure.value()<sp/>/<sp/>m_polytropic_constant,<sp/>m_polytropic_index<sp/>/<sp/>(m_polytropic_index<sp/>+<sp/>1.0))};</highlight></codeline>
<codeline lineno="61"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="62"><highlight class="normal"></highlight></codeline>
<codeline lineno="63"><highlight class="normal"><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/>SecondPartialDerivative&lt;Pressure,<sp/>MassDensity,<sp/>MassDensity&gt;<sp/>Polytrope::second_partial_derivative(</highlight></codeline>
<codeline lineno="64"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>PressureFromDensity,<sp/>WithRespectTo&lt;MassDensity&gt;,<sp/>WithRespectTo&lt;MassDensity&gt;,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>DensityValue<sp/>density</highlight></codeline>
<codeline lineno="65"><highlight class="normal"><sp/><sp/><sp/><sp/>)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="66"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>validate_nonnegativity&lt;EOSEvaluationError,<sp/>SERiFErrorCode::EOS_INVALID_THERMODYNAMIC_INPUT&gt;(density.value());</highlight></codeline>
<codeline lineno="67"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>utils::errors::serif_ensure&lt;SERiFErrorCode::EOS_INVALID_THERMODYNAMIC_INPUT,<sp/>EOSEvaluationError&gt;(</highlight></codeline>
<codeline lineno="68"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>density.value()<sp/>&gt;<sp/>0.0<sp/>||<sp/>m_polytropic_index<sp/>==<sp/>1.0,</highlight></codeline>
<codeline lineno="69"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="stringliteral">&quot;The<sp/>second<sp/>pressure<sp/>derivative<sp/>is<sp/>singular<sp/>at<sp/>zero<sp/>density<sp/>for<sp/>polytropic<sp/>index<sp/>{}.<sp/>Apply<sp/>the<sp/>vacuum<sp/>density<sp/>trace<sp/>directly<sp/>and<sp/>evaluate<sp/>bulk<sp/>Newton<sp/>coefficients<sp/>only<sp/>at<sp/>positive<sp/>density.&quot;</highlight><highlight class="normal">,</highlight></codeline>
<codeline lineno="70"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_polytropic_index);</highlight></codeline>
<codeline lineno="71"><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/>exponent<sp/>=<sp/>1.0<sp/>/<sp/>m_polytropic_index;</highlight></codeline>
<codeline lineno="72"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>SecondPartialDerivative&lt;Pressure,<sp/>MassDensity,<sp/>MassDensity&gt;{</highlight></codeline>
<codeline lineno="73"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_polytropic_constant<sp/>*<sp/>(1.0<sp/>+<sp/>exponent)<sp/>*<sp/>exponent<sp/>*<sp/>std::pow(density.value(),<sp/>exponent<sp/>-<sp/>1.0)};</highlight></codeline>
<codeline lineno="74"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="75"><highlight class="normal"></highlight></codeline>
<codeline lineno="76"><highlight class="normal"><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/><ref refid="namespaceserif_1_1dimensions_1a9d34f8c867a871e5ebe50e99b54d0d16" kindref="member">SpecificEnthalpyValue</ref><sp/>Polytrope::evaluate(SpecificEnthalpyFromDensity,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>DensityValue<sp/>density)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="77"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>validate_nonnegativity&lt;EOSEvaluationError,<sp/>SERiFErrorCode::EOS_INVALID_THERMODYNAMIC_INPUT&gt;(density.value());</highlight></codeline>
<codeline lineno="78"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(density.value()<sp/>==<sp/>0.0)<sp/>{</highlight></codeline>
<codeline lineno="79"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>h<sp/>=<sp/>0</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="80"><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="namespaceserif_1_1dimensions_1a9d34f8c867a871e5ebe50e99b54d0d16" kindref="member">SpecificEnthalpyValue</ref>{0.0};</highlight></codeline>
<codeline lineno="81"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="82"><highlight class="normal"></highlight></codeline>
<codeline lineno="83"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>h<sp/>=<sp/>hs<sp/>*<sp/>(ρ)^(1/n)</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="84"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1dimensions_1a9d34f8c867a871e5ebe50e99b54d0d16" kindref="member">SpecificEnthalpyValue</ref>{</highlight></codeline>
<codeline lineno="85"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_enthalpy_scale<sp/>*<sp/>std::pow(density.value(),<sp/>1.0<sp/>/<sp/>m_polytropic_index)</highlight></codeline>
<codeline lineno="86"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="87"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="88"><highlight class="normal"></highlight></codeline>
<codeline lineno="90"><highlight class="normal"><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/><ref refid="namespaceserif_1_1dimensions_1a9c52fa3afb2594825dd52879fa155bd3" kindref="member">BarotropicPressurePotentialValue</ref><sp/>Polytrope::evaluate(PressurePotentialFromDensity,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>DensityValue<sp/>density)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="91"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1dimensions_1a9c52fa3afb2594825dd52879fa155bd3" kindref="member">BarotropicPressurePotentialValue</ref>{<ref refid="namespaceserif_1_1eos_1ade79763ef51f539b78059d8898a2d2a2" kindref="member">evaluate</ref>(SpecificEnthalpyFromDensity{},<sp/>density).<ref refid="namespaceserif_1_1discretization_1_1backend_1_1mfem_1af75f98db90618153ea4ca169945733f7a2063c1608d6e0baf80249c42e2be5804" kindref="member">value</ref>()};</highlight></codeline>
<codeline lineno="92"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="93"><highlight class="normal"></highlight></codeline>
<codeline lineno="94"><highlight class="normal"><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/><ref refid="namespaceserif_1_1operators_1a84d7e8ff528dc3361769cc6f08e6cd49" kindref="member">DensityValue</ref><sp/>Polytrope::evaluate(DensityFromSpecificEnthalpy,<sp/>SpecificEnthalpyValue<sp/>specific_enthalpy)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="95"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>validate_finite&lt;EOSEvaluationError,<sp/>SERiFErrorCode::EOS_INVALID_THERMODYNAMIC_INPUT&gt;(specific_enthalpy.value());</highlight></codeline>
<codeline lineno="96"><highlight class="normal"></highlight></codeline>
<codeline lineno="97"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(specific_enthalpy.value()<sp/>&lt;=<sp/>0.0)<sp/>{</highlight></codeline>
<codeline lineno="98"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>ρ<sp/>=<sp/>0</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="99"><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="namespaceserif_1_1operators_1a84d7e8ff528dc3361769cc6f08e6cd49" kindref="member">DensityValue</ref>{0.0};</highlight></codeline>
<codeline lineno="100"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="101"><highlight class="normal"></highlight></codeline>
<codeline lineno="102"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>ρ<sp/>=<sp/>(h/hs)^n</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="103"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1operators_1a84d7e8ff528dc3361769cc6f08e6cd49" kindref="member">DensityValue</ref>{std::pow(specific_enthalpy.value()<sp/>/<sp/>m_enthalpy_scale,<sp/>m_polytropic_index)};</highlight></codeline>
<codeline lineno="104"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="105"><highlight class="normal"></highlight></codeline>
<codeline lineno="106"><highlight class="normal"><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/><ref refid="namespaceserif_1_1dimensions_1a60ad5bcd33f38e5e0e272d9df76008ec" kindref="member">PressureValue</ref><sp/>Polytrope::evaluate(PressureFromSpecificEnthalpy,<sp/>SpecificEnthalpyValue<sp/>specific_enthalpy)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="107"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1operators_1a84d7e8ff528dc3361769cc6f08e6cd49" kindref="member">DensityValue</ref><sp/>density<sp/>=<sp/><ref refid="namespaceserif_1_1eos_1ade79763ef51f539b78059d8898a2d2a2" kindref="member">evaluate</ref>(DensityFromSpecificEnthalpy{},<sp/>specific_enthalpy);</highlight></codeline>
<codeline lineno="108"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(specific_enthalpy.value()<sp/>&lt;=<sp/>0.0)<sp/>{</highlight></codeline>
<codeline lineno="109"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>P<sp/>=<sp/>0</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="110"><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="namespaceserif_1_1dimensions_1a60ad5bcd33f38e5e0e272d9df76008ec" kindref="member">PressureValue</ref>{0.0};</highlight></codeline>
<codeline lineno="111"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="112"><highlight class="normal"></highlight></codeline>
<codeline lineno="113"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>P<sp/>=<sp/>ρ<sp/>*<sp/>h<sp/>/<sp/>(n<sp/>+<sp/>1)</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="114"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1dimensions_1a60ad5bcd33f38e5e0e272d9df76008ec" kindref="member">PressureValue</ref>{density.value()<sp/>*<sp/>specific_enthalpy.value()/<sp/>(m_polytropic_index<sp/>+<sp/>1.0)};</highlight></codeline>
<codeline lineno="115"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="116"><highlight class="normal"></highlight></codeline>
<codeline lineno="117"><highlight class="normal"><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/><ref refid="namespaceserif_1_1dimensions_1a9d34f8c867a871e5ebe50e99b54d0d16" kindref="member">SpecificEnthalpyValue</ref><sp/>Polytrope::evaluate(SpecificEnthalpyFromPressure,<sp/>PressureValue<sp/>pressure)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="118"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>validate_nonnegativity&lt;EOSEvaluationError,<sp/>SERiFErrorCode::EOS_INVALID_THERMODYNAMIC_INPUT&gt;(pressure.value());</highlight></codeline>
<codeline lineno="119"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(pressure.value()<sp/>==<sp/>0.0)<sp/>{</highlight></codeline>
<codeline lineno="120"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>h<sp/>=<sp/>0</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="121"><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="namespaceserif_1_1dimensions_1a9d34f8c867a871e5ebe50e99b54d0d16" kindref="member">SpecificEnthalpyValue</ref>{0.0};</highlight></codeline>
<codeline lineno="122"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="123"><highlight class="normal"></highlight></codeline>
<codeline lineno="124"><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/>indexPlusOne<sp/>=<sp/>m_polytropic_index<sp/>+<sp/>1.0;</highlight></codeline>
<codeline lineno="125"><highlight class="normal"></highlight></codeline>
<codeline lineno="126"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>h<sp/>=<sp/>(n+1)<sp/>*<sp/>(K^(n/(n+1)))<sp/>*<sp/>(P^(1/(n+1)))</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="127"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1dimensions_1a9d34f8c867a871e5ebe50e99b54d0d16" kindref="member">SpecificEnthalpyValue</ref>{</highlight></codeline>
<codeline lineno="128"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>indexPlusOne<sp/>*<sp/>std::pow(m_polytropic_constant,<sp/>m_polytropic_index<sp/>/<sp/>indexPlusOne)<sp/>*</highlight></codeline>
<codeline lineno="129"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::pow(pressure.value(),<sp/>1.0<sp/>/<sp/>indexPlusOne)</highlight></codeline>
<codeline lineno="130"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="131"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="132"><highlight class="normal"></highlight></codeline>
<codeline lineno="133"><highlight class="normal"><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/>PartialDerivative&lt;MassDensity,<sp/>SpecificEnthalpy&gt;<sp/>Polytrope::partial_derivative(DensityFromSpecificEnthalpy,<sp/>WithRespectTo&lt;SpecificEnthalpy&gt;,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>SpecificEnthalpyValue<sp/>specific_enthalpy)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="134"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>validate_finite&lt;EOSEvaluationError,<sp/>SERiFErrorCode::EOS_INVALID_THERMODYNAMIC_INPUT&gt;(specific_enthalpy.value());</highlight></codeline>
<codeline lineno="135"><highlight class="normal"></highlight></codeline>
<codeline lineno="136"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(specific_enthalpy.value()<sp/>&lt;<sp/>0.0)<sp/>{</highlight></codeline>
<codeline lineno="137"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>dρ/dh<sp/>=<sp/>0</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="138"><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/>PartialDerivative&lt;MassDensity,<sp/>SpecificEnthalpy&gt;{0.0};</highlight></codeline>
<codeline lineno="139"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="140"><highlight class="normal"></highlight></codeline>
<codeline lineno="141"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(specific_enthalpy.value()<sp/>==<sp/>0.0)<sp/>{</highlight></codeline>
<codeline lineno="142"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>dρ/dh<sp/>=<sp/>n<sp/>==<sp/>1<sp/>-&gt;<sp/>1/hs</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="143"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>dρ/dh<sp/>=<sp/>else<sp/><sp/><sp/>-&gt;<sp/>0</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="144"><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/>PartialDerivative&lt;MassDensity,<sp/>SpecificEnthalpy&gt;{</highlight></codeline>
<codeline lineno="145"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_polytropic_index<sp/>==<sp/>1.0<sp/>?<sp/>1.0<sp/>/<sp/>m_enthalpy_scale<sp/>:<sp/>0.0</highlight></codeline>
<codeline lineno="146"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="147"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="148"><highlight class="normal"></highlight></codeline>
<codeline lineno="149"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>dρ/dh<sp/>=<sp/>(n/hs)<sp/>*<sp/>(h/hs)^(n-1)</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="150"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>PartialDerivative&lt;MassDensity,<sp/>SpecificEnthalpy&gt;{</highlight></codeline>
<codeline lineno="151"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_polytropic_index<sp/>/<sp/>m_enthalpy_scale<sp/>*<sp/>std::pow(specific_enthalpy.value()<sp/>/<sp/>m_enthalpy_scale,<sp/>m_polytropic_index<sp/>-<sp/>1.0)</highlight></codeline>
<codeline lineno="152"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="153"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="154"><highlight class="normal"></highlight></codeline>
<codeline lineno="155"><highlight class="normal"><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/>PartialDerivative&lt;Pressure,<sp/>SpecificEnthalpy&gt;<sp/>Polytrope::partial_derivative(PressureFromSpecificEnthalpy,<sp/>WithRespectTo&lt;SpecificEnthalpy&gt;,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>SpecificEnthalpyValue<sp/>specific_enthalpy)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="156"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1operators_1a84d7e8ff528dc3361769cc6f08e6cd49" kindref="member">DensityValue</ref><sp/>density<sp/>=<sp/><ref refid="namespaceserif_1_1eos_1ade79763ef51f539b78059d8898a2d2a2" kindref="member">evaluate</ref>(DensityFromSpecificEnthalpy{},<sp/>specific_enthalpy);</highlight></codeline>
<codeline lineno="157"><highlight class="normal"></highlight></codeline>
<codeline lineno="158"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>dP/dh<sp/>=<sp/>ρ</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="159"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>PartialDerivative&lt;Pressure,<sp/>SpecificEnthalpy&gt;{density.value()};</highlight></codeline>
<codeline lineno="160"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="161"><highlight class="normal"></highlight></codeline>
<codeline lineno="162"><highlight class="normal"><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/>PartialDerivative&lt;Pressure,<sp/>MassDensity&gt;<sp/>Polytrope::partial_derivative(PressureFromDensity,<sp/>WithRespectTo&lt;MassDensity&gt;,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>DensityValue<sp/>density)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="163"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>validate_nonnegativity&lt;EOSEvaluationError,<sp/>SERiFErrorCode::EOS_INVALID_THERMODYNAMIC_INPUT&gt;(density.value());</highlight></codeline>
<codeline lineno="164"><highlight class="normal"></highlight></codeline>
<codeline lineno="165"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(density.value()<sp/>==<sp/>0.0)<sp/>{</highlight></codeline>
<codeline lineno="166"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>dP/dρ<sp/>=<sp/>0</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="167"><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/>PartialDerivative&lt;Pressure,<sp/>MassDensity&gt;{0.0};</highlight></codeline>
<codeline lineno="168"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="169"><highlight class="normal"></highlight></codeline>
<codeline lineno="170"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>dP/dρ<sp/>=<sp/>K<sp/>*<sp/>(1<sp/>+<sp/>1/n)<sp/>*<sp/>(ρ)^(1/n)</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="171"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>PartialDerivative&lt;Pressure,<sp/>MassDensity&gt;{</highlight></codeline>
<codeline lineno="172"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_polytropic_constant<sp/>*<sp/>(1.0<sp/>+<sp/>1.0<sp/>/<sp/>m_polytropic_index)<sp/>*<sp/>std::pow(density.value(),<sp/>1.0<sp/>/<sp/>m_polytropic_index)</highlight></codeline>
<codeline lineno="173"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="174"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="175"><highlight class="normal">}</highlight></codeline>
    </programlisting>
    <location file="src/lib/serif/eos/models/polytropic.cpp"/>
  </compounddef>
</doxygen>