MeanField_SERiF/docs/xml/evaluator_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="evaluator_8cpp" kind="file" language="C++">
    <compoundname>evaluator.cpp</compoundname>
    <includes refid="evaluator_8hpp" local="yes">serif/backend/mfem/mapping/evaluator.hpp</includes>
    <includes refid="errors_8hpp" local="yes">serif/utils/error/errors.hpp</includes>
    <incdepgraph>
      <node id="22">
        <label>serif/backend/mfem/mapping/compactification/compactification.hpp</label>
        <link refid="compactification_8hpp"/>
        <childnode refid="13" relation="include">
        </childnode>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
        <childnode refid="23" relation="include">
        </childnode>
        <childnode refid="24" relation="include">
        </childnode>
      </node>
      <node id="23">
        <label>serif/backend/mfem/mapping/compactification/input_result.hpp</label>
        <link refid="input__result_8hpp"/>
        <childnode refid="7" relation="include">
        </childnode>
        <childnode refid="11" relation="include">
        </childnode>
      </node>
      <node id="12">
        <label>serif/backend/mfem/mapping/compactification/kelvin.hpp</label>
        <link refid="backend_2mfem_2mapping_2compactification_2kelvin_8hpp"/>
        <childnode refid="13" relation="include">
        </childnode>
        <childnode refid="7" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="25" relation="include">
        </childnode>
      </node>
      <node id="24">
        <label>serif/backend/mfem/mapping/compactification/variation.hpp</label>
        <link refid="compactification_2variation_8hpp"/>
        <childnode refid="7" relation="include">
        </childnode>
        <childnode refid="11" relation="include">
        </childnode>
      </node>
      <node id="26">
        <label>serif/backend/mfem/mapping/context.hpp</label>
        <link refid="context_8hpp"/>
        <childnode refid="11" relation="include">
        </childnode>
        <childnode refid="7" relation="include">
        </childnode>
      </node>
      <node id="27">
        <label>serif/backend/mfem/mapping/data.hpp</label>
        <link refid="data_8hpp"/>
        <childnode refid="28" relation="include">
        </childnode>
        <childnode refid="29" relation="include">
        </childnode>
        <childnode refid="7" relation="include">
        </childnode>
        <childnode refid="11" relation="include">
        </childnode>
        <childnode refid="30" relation="include">
        </childnode>
        <childnode refid="16" relation="include">
        </childnode>
        <childnode refid="12" relation="include">
        </childnode>
        <childnode refid="26" relation="include">
        </childnode>
        <childnode refid="27" relation="include">
        </childnode>
      </node>
      <node id="9">
        <label>serif/backend/mfem/mapping/domain_mapper.hpp</label>
        <link refid="backend_2mfem_2mapping_2domain__mapper_8hpp"/>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="7" relation="include">
        </childnode>
        <childnode refid="11" relation="include">
        </childnode>
        <childnode refid="5" relation="include">
        </childnode>
        <childnode refid="6" relation="include">
        </childnode>
        <childnode refid="12" relation="include">
        </childnode>
        <childnode refid="26" relation="include">
        </childnode>
        <childnode refid="27" relation="include">
        </childnode>
        <childnode refid="31" relation="include">
        </childnode>
        <childnode refid="32" relation="include">
        </childnode>
      </node>
      <node id="2">
        <label>serif/backend/mfem/mapping/evaluator.hpp</label>
        <link refid="evaluator_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="7" relation="include">
        </childnode>
        <childnode refid="8" relation="include">
        </childnode>
        <childnode refid="9" relation="include">
        </childnode>
        <childnode refid="26" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
      </node>
      <node id="14">
        <label>serif/backend/mfem/mapping/status.hpp</label>
        <link refid="status_8hpp"/>
        <childnode refid="15" relation="include">
        </childnode>
        <childnode refid="16" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="18" relation="include">
        </childnode>
        <childnode refid="19" relation="include">
        </childnode>
      </node>
      <node id="31">
        <label>serif/backend/mfem/mapping/variation.hpp</label>
        <link refid="variation_8hpp"/>
        <childnode refid="11" relation="include">
        </childnode>
        <childnode refid="7" relation="include">
        </childnode>
      </node>
      <node id="19">
        <label>serif/meta/enums.hpp</label>
        <link refid="enums_8hpp"/>
        <childnode refid="15" relation="include">
        </childnode>
        <childnode refid="20" relation="include">
        </childnode>
        <childnode refid="21" relation="include">
        </childnode>
      </node>
      <node id="18">
        <label>serif/utils/error/codes.hpp</label>
        <link refid="codes_8hpp"/>
        <childnode refid="16" relation="include">
        </childnode>
        <childnode refid="19" relation="include">
        </childnode>
      </node>
      <node id="33">
        <label>serif/utils/error/errors.hpp</label>
        <link refid="errors_8hpp"/>
        <childnode refid="13" relation="include">
        </childnode>
        <childnode refid="34" relation="include">
        </childnode>
        <childnode refid="35" relation="include">
        </childnode>
        <childnode refid="36" relation="include">
        </childnode>
        <childnode refid="16" relation="include">
        </childnode>
        <childnode refid="37" relation="include">
        </childnode>
        <childnode refid="38" relation="include">
        </childnode>
        <childnode refid="39" relation="include">
        </childnode>
        <childnode refid="42" relation="include">
        </childnode>
        <childnode refid="18" relation="include">
        </childnode>
        <childnode refid="46" relation="include">
        </childnode>
      </node>
      <node id="42">
        <label>serif/utils/misc/concepts/exceptions.hpp</label>
        <link refid="utils_2misc_2concepts_2exceptions_8hpp"/>
        <childnode refid="43" relation="include">
        </childnode>
        <childnode refid="20" relation="include">
        </childnode>
        <childnode refid="44" relation="include">
        </childnode>
      </node>
      <node id="44">
        <label>serif/utils/misc/concepts/numeric.hpp</label>
        <link refid="numeric_8hpp"/>
        <childnode refid="45" relation="include">
        </childnode>
        <childnode refid="20" relation="include">
        </childnode>
      </node>
      <node id="39">
        <label>serif/utils/misc/terminal/colors.hpp</label>
        <link refid="colors_8hpp"/>
        <childnode refid="13" relation="include">
        </childnode>
        <childnode refid="35" relation="include">
        </childnode>
        <childnode refid="40" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="41" relation="include">
        </childnode>
      </node>
      <node id="46">
        <label>serif/utils/types/strings.hpp</label>
        <link refid="types_2strings_8hpp"/>
        <childnode refid="37" relation="include">
        </childnode>
        <childnode refid="47" relation="include">
        </childnode>
        <childnode refid="13" relation="include">
        </childnode>
      </node>
      <node id="25">
        <label>serif/utils/user_options/mapping/compactification/kelvin.hpp</label>
        <link refid="utils_2user__options_2mapping_2compactification_2kelvin_8hpp"/>
      </node>
      <node id="32">
        <label>serif/utils/user_options/mapping/domain_mapper.hpp</label>
        <link refid="utils_2user__options_2mapping_2domain__mapper_8hpp"/>
      </node>
      <node id="1">
        <label>src/lib/serif/backend/mfem/mapping/evaluator.cpp</label>
        <link refid="evaluator_8cpp"/>
        <childnode refid="2" relation="include">
        </childnode>
        <childnode refid="33" relation="include">
        </childnode>
      </node>
      <node id="37">
        <label>algorithm</label>
      </node>
      <node id="45">
        <label>concepts</label>
      </node>
      <node id="16">
        <label>cstdint</label>
      </node>
      <node id="35">
        <label>cstdio</label>
      </node>
      <node id="40">
        <label>cstdlib</label>
      </node>
      <node id="43">
        <label>exception</label>
      </node>
      <node id="10">
        <label>expected</label>
      </node>
      <node id="34">
        <label>format</label>
      </node>
      <node id="47">
        <label>iostream</label>
      </node>
      <node id="21">
        <label>magic_enum/magic_enum_all.hpp</label>
      </node>
      <node id="3">
        <label>memory</label>
      </node>
      <node id="5">
        <label>mfem/fem/eltrans.hpp</label>
      </node>
      <node id="28">
        <label>mfem/fem/fe/fe_base.hpp</label>
      </node>
      <node id="29">
        <label>mfem/fem/fespace.hpp</label>
      </node>
      <node id="4">
        <label>mfem/fem/gridfunc.hpp</label>
      </node>
      <node id="6">
        <label>mfem/fem/intrules.hpp</label>
      </node>
      <node id="8">
        <label>mfem/general/array.hpp</label>
      </node>
      <node id="11">
        <label>mfem/linalg/densemat.hpp</label>
      </node>
      <node id="7">
        <label>mfem/linalg/vector.hpp</label>
      </node>
      <node id="30">
        <label>mfem/mesh/mesh.hpp</label>
      </node>
      <node id="36">
        <label>string</label>
      </node>
      <node id="13">
        <label>string_view</label>
      </node>
      <node id="20">
        <label>type_traits</label>
      </node>
      <node id="41">
        <label>unistd.h</label>
      </node>
      <node id="15">
        <label>unordered_map</label>
      </node>
      <node id="17">
        <label>utility</label>
      </node>
      <node id="38">
        <label>version</label>
      </node>
    </incdepgraph>
    <innernamespace refid="namespaceserif">serif</innernamespace>
    <innernamespace refid="namespaceserif_1_1mapping">serif::mapping</innernamespace>
    <briefdescription>
    </briefdescription>
    <detaileddescription>
    </detaileddescription>
    <programlisting>
<codeline lineno="1"><highlight class="preprocessor">#include<sp/>&quot;<ref refid="evaluator_8hpp" kindref="compound">serif/backend/mfem/mapping/evaluator.hpp</ref>&quot;</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/>&quot;<ref refid="errors_8hpp" kindref="compound">serif/utils/error/errors.hpp</ref>&quot;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="4"><highlight class="normal"></highlight></codeline>
<codeline lineno="5"><highlight class="normal"></highlight></codeline>
<codeline lineno="6"><highlight class="normal"></highlight><highlight class="keyword">namespace<sp/></highlight><highlight class="normal"><ref refid="namespaceserif_1_1mapping" kindref="compound">serif::mapping</ref><sp/>{</highlight></codeline>
<codeline lineno="7"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/>namespace<sp/></highlight><highlight class="normal">utils::errors;</highlight></codeline>
<codeline lineno="8"><highlight class="normal"></highlight></codeline>
<codeline lineno="9"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="classserif_1_1mapping_1_1_grid_function_mapping_evaluator_1ad8973a6103c042d0e31929cf67685fd6" kindref="member">GridFunctionMappingEvaluator::GridFunctionMappingEvaluator</ref>(</highlight></codeline>
<codeline lineno="10"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="classserif_1_1mapping_1_1_domain_mapper" kindref="compound">DomainMapper</ref><sp/>&amp;mapper,</highlight></codeline>
<codeline lineno="11"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::GridFunction<sp/>&amp;displacement,</highlight></codeline>
<codeline lineno="12"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::GridFunction<sp/>&amp;compactification_coordinate</highlight></codeline>
<codeline lineno="13"><highlight class="normal"><sp/><sp/><sp/><sp/>)<sp/>:</highlight></codeline>
<codeline lineno="14"><highlight class="normal"><sp/><sp/><sp/><sp/>m_mapper(mapper),</highlight></codeline>
<codeline lineno="15"><highlight class="normal"><sp/><sp/><sp/><sp/>m_displacement(displacement),</highlight></codeline>
<codeline lineno="16"><highlight class="normal"><sp/><sp/><sp/><sp/>m_compactification_coordinate(compactification_coordinate),</highlight></codeline>
<codeline lineno="17"><highlight class="normal"><sp/><sp/><sp/><sp/>m_displacement_space(displacement.FESpace()),</highlight></codeline>
<codeline lineno="18"><highlight class="normal"><sp/><sp/><sp/><sp/>m_compactification_space(compactification_coordinate.FESpace()),</highlight></codeline>
<codeline lineno="19"><highlight class="normal"><sp/><sp/><sp/><sp/>m_displacement_space_sequence(m_displacement_space<sp/>!=<sp/>nullptr<sp/>?<sp/>m_displacement_space-&gt;GetSequence()<sp/>:<sp/>-1),</highlight></codeline>
<codeline lineno="20"><highlight class="normal"><sp/><sp/><sp/><sp/>m_compactification_space_sequence(m_compactification_space<sp/>!=<sp/>nullptr<sp/>?<sp/>m_compactification_space-&gt;GetSequence()<sp/>:<sp/>-1),</highlight></codeline>
<codeline lineno="21"><highlight class="normal"><sp/><sp/><sp/><sp/>m_workspace(mapper.get_dimension())<sp/>{</highlight></codeline>
<codeline lineno="22"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>serif_ensure&lt;SERiFErrorCode::MISSING_DISPLACEMENT_FINITE_ELEMENT_SPACE&gt;(m_displacement_space<sp/>!=<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">,<sp/></highlight><highlight class="stringliteral">&quot;Displacement<sp/>grid<sp/>function<sp/>is<sp/>not<sp/>associated<sp/>with<sp/>a<sp/>finite-element<sp/>space.&quot;</highlight><highlight class="normal">);</highlight></codeline>
<codeline lineno="23"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>serif_ensure&lt;SERiFErrorCode::MISSING_DISPLACEMENT_FINITE_ELEMENT_SPACE&gt;(m_compactification_space<sp/>!=<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">,<sp/></highlight><highlight class="stringliteral">&quot;Compactification-coordinate<sp/>grid<sp/>function<sp/>is<sp/>not<sp/>associated<sp/>with<sp/>a<sp/>finite-element<sp/>space.&quot;</highlight><highlight class="normal">);</highlight></codeline>
<codeline lineno="24"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>serif_ensure&lt;SERiFErrorCode::MISMATCHED_MESHING_MAPS&gt;(m_displacement_space-&gt;GetMesh()<sp/>==<sp/>m_compactification_space-&gt;GetMesh(),<sp/></highlight><highlight class="stringliteral">&quot;Displacement<sp/>and<sp/>compactification-coordinate<sp/>grid<sp/>functions<sp/>are<sp/>not<sp/>associated<sp/>with<sp/>the<sp/>same<sp/>mesh.&quot;</highlight><highlight class="normal">);</highlight></codeline>
<codeline lineno="25"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>serif_ensure&lt;SERiFErrorCode::INVALID_DIMENSION&gt;(m_displacement.VectorDim()<sp/>==<sp/>mapper.get_dimension(),<sp/></highlight><highlight class="stringliteral">&quot;Displacement<sp/>grid<sp/>function<sp/>vector<sp/>dimension<sp/>does<sp/>not<sp/>match<sp/>the<sp/>mapper<sp/>dimension.<sp/>({}<sp/>vs.<sp/>{})&quot;</highlight><highlight class="normal">,<sp/>m_displacement.VectorDim(),<sp/>mapper.get_dimension());</highlight></codeline>
<codeline lineno="26"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>serif_ensure&lt;SERiFErrorCode::INVALID_COMPACTIFICATION_RANK&gt;(m_compactification_coordinate.VectorDim()<sp/>==<sp/>1,<sp/></highlight><highlight class="stringliteral">&quot;Compactification-coordinate<sp/>grid<sp/>function<sp/>must<sp/>be<sp/>a<sp/>scalar<sp/>field.<sp/>Instead<sp/>it<sp/>has<sp/>vector<sp/>dimension<sp/>{}.&quot;</highlight><highlight class="normal">,<sp/>m_compactification_coordinate.VectorDim());</highlight></codeline>
<codeline lineno="27"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>serif_ensure&lt;SERiFErrorCode::INVALID_DIMENSION&gt;(m_displacement_space-&gt;GetMesh()-&gt;SpaceDimension()<sp/>==<sp/>mapper.get_dimension(),<sp/></highlight><highlight class="stringliteral">&quot;Displacement<sp/>finite-element<sp/>space<sp/>mesh<sp/>dimension<sp/>does<sp/>not<sp/>match<sp/>the<sp/>mapper<sp/>dimension.<sp/>({}<sp/>vs.<sp/>{})&quot;</highlight><highlight class="normal">,<sp/>m_displacement_space-&gt;GetMesh()-&gt;SpaceDimension(),<sp/>mapper.get_dimension());</highlight></codeline>
<codeline lineno="28"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="29"><highlight class="normal"></highlight></codeline>
<codeline lineno="30"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">void</highlight><highlight class="normal"><sp/>GridFunctionMappingEvaluator::invalidate_cache()<sp/>noexcept<sp/>{</highlight></codeline>
<codeline lineno="31"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_displacement_data.reset();</highlight></codeline>
<codeline lineno="32"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_compactification_data.reset();</highlight></codeline>
<codeline lineno="33"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_cached_element_id<sp/>=<sp/>-1;</highlight></codeline>
<codeline lineno="34"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="35"><highlight class="normal"></highlight></codeline>
<codeline lineno="36"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">void</highlight><highlight class="normal"><sp/>GridFunctionMappingEvaluator::validate_field_bindings()</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="37"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>serif_ensure&lt;SERiFErrorCode::CACHE_FE_BINDING_ERROR&gt;(m_displacement.FESpace()<sp/>==<sp/>m_displacement_space,<sp/></highlight><highlight class="stringliteral">&quot;Displacement<sp/>grid<sp/>function<sp/>is<sp/>no<sp/>longer<sp/>associated<sp/>with<sp/>the<sp/>cached<sp/>finite-element<sp/>space.&quot;</highlight><highlight class="normal">);</highlight></codeline>
<codeline lineno="38"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>serif_ensure&lt;SERiFErrorCode::CACHE_FE_BINDING_ERROR&gt;(m_compactification_coordinate.FESpace()<sp/>==<sp/>m_compactification_space,<sp/></highlight><highlight class="stringliteral">&quot;Compactification-coordinate<sp/>grid<sp/>function<sp/>is<sp/>no<sp/>longer<sp/>associated<sp/>with<sp/>the<sp/>cached<sp/>finite-element<sp/>space.&quot;</highlight><highlight class="normal">);</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/></highlight><highlight class="keywordtype">bool</highlight><highlight class="normal"><sp/>GridFunctionMappingEvaluator::invalidate_for_changed_spaces()<sp/>{</highlight></codeline>
<codeline lineno="42"><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/>displacement_sequence<sp/>=<sp/>m_displacement_space-&gt;GetSequence();</highlight></codeline>
<codeline lineno="43"><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/>compactification_sequence<sp/>=<sp/>m_compactification_space-&gt;GetSequence();</highlight></codeline>
<codeline lineno="44"><highlight class="normal"></highlight></codeline>
<codeline lineno="45"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(displacement_sequence<sp/>==<sp/>m_displacement_space_sequence<sp/>&amp;&amp;<sp/>compactification_sequence<sp/>==<sp/>m_compactification_space_sequence)<sp/>{</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/></highlight><highlight class="keyword">false</highlight><highlight class="normal">;</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/>invalidate_cache();</highlight></codeline>
<codeline lineno="50"><highlight class="normal"></highlight></codeline>
<codeline lineno="51"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_displacement_space_sequence<sp/>=<sp/>displacement_sequence;</highlight></codeline>
<codeline lineno="52"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_compactification_space_sequence<sp/>=<sp/>compactification_sequence;</highlight></codeline>
<codeline lineno="53"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">true</highlight><highlight class="normal">;</highlight></codeline>
<codeline lineno="54"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="55"><highlight class="normal"></highlight></codeline>
<codeline lineno="56"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">void</highlight><highlight class="normal"><sp/>GridFunctionMappingEvaluator::refresh()<sp/>{</highlight></codeline>
<codeline lineno="57"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>validate_field_bindings();</highlight></codeline>
<codeline lineno="58"><highlight class="normal"></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/>(invalidate_for_changed_spaces())<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">;</highlight></codeline>
<codeline lineno="61"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="62"><highlight class="normal"></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">int</highlight><highlight class="normal"><sp/>element_id<sp/>=<sp/>m_cached_element_id;</highlight></codeline>
<codeline lineno="64"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>invalidate_cache();</highlight></codeline>
<codeline lineno="65"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(element_id<sp/>&gt;=<sp/>0)<sp/>{</highlight></codeline>
<codeline lineno="66"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>load_element(element_id);</highlight></codeline>
<codeline lineno="67"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="68"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="69"><highlight class="normal"></highlight></codeline>
<codeline lineno="70"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">void</highlight><highlight class="normal"><sp/>GridFunctionMappingEvaluator::load_element(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>element_id)<sp/>{</highlight></codeline>
<codeline lineno="71"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>validate_field_bindings();</highlight></codeline>
<codeline lineno="72"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>(void)invalidate_for_changed_spaces();</highlight></codeline>
<codeline lineno="73"><highlight class="normal"></highlight></codeline>
<codeline lineno="74"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(element_id<sp/>==<sp/>m_cached_element_id)<sp/>{</highlight></codeline>
<codeline lineno="75"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal">;</highlight></codeline>
<codeline lineno="76"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="77"><highlight class="normal"></highlight></codeline>
<codeline lineno="78"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::FiniteElementSpace&amp;<sp/>displacement_space<sp/>=<sp/>*m_displacement_space;</highlight></codeline>
<codeline lineno="79"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::FiniteElementSpace&amp;<sp/>compactification_space<sp/>=<sp/>*m_compactification_space;</highlight></codeline>
<codeline lineno="80"><highlight class="normal"></highlight></codeline>
<codeline lineno="81"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>serif_ensure&lt;SERiFErrorCode::INVALID_ELEMENT_ID_ERROR&gt;(element_id<sp/>&gt;=<sp/>0,<sp/></highlight><highlight class="stringliteral">&quot;Element<sp/>ID<sp/>must<sp/>be<sp/>non-negative.<sp/>Instead<sp/>{}<sp/>was<sp/>provided.&quot;</highlight><highlight class="normal">,<sp/>element_id);</highlight></codeline>
<codeline lineno="82"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>serif_ensure&lt;SERiFErrorCode::INVALID_ELEMENT_ID_ERROR&gt;(element_id<sp/>&lt;<sp/>displacement_space.GetMesh()-&gt;GetNE(),<sp/></highlight><highlight class="stringliteral">&quot;Element<sp/>ID<sp/>must<sp/>be<sp/>less<sp/>than<sp/>the<sp/>number<sp/>of<sp/>elements<sp/>in<sp/>the<sp/>mesh<sp/>({}).<sp/>Instead<sp/>{}<sp/>was<sp/>provided.&quot;</highlight><highlight class="normal">,<sp/>displacement_space.GetMesh()-&gt;GetNE(),<sp/>element_id);</highlight></codeline>
<codeline lineno="83"><highlight class="normal"></highlight></codeline>
<codeline lineno="84"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::DofTransformation<sp/>*displacement_transformation<sp/>=<sp/>displacement_space.GetElementVDofs(element_id,<sp/>m_displacement_dofs);</highlight></codeline>
<codeline lineno="85"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>compactification_space.GetElementDofs(element_id,<sp/>m_compactification_dofs);</highlight></codeline>
<codeline lineno="86"><highlight class="normal"></highlight></codeline>
<codeline lineno="87"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_displacement.GetSubVector(m_displacement_dofs,<sp/>m_element_displacement);</highlight></codeline>
<codeline lineno="88"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_compactification_coordinate.GetSubVector(m_compactification_dofs,<sp/>m_element_compactification);</highlight></codeline>
<codeline lineno="89"><highlight class="normal"></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/>(displacement_transformation<sp/>!=<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">)<sp/>{</highlight></codeline>
<codeline lineno="91"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>displacement_transformation-&gt;InvTransformPrimal(m_element_displacement);</highlight></codeline>
<codeline lineno="92"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="93"><highlight class="normal"></highlight></codeline>
<codeline lineno="94"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::FiniteElement<sp/>&amp;displacement_element<sp/>=<sp/>*displacement_space.GetFE(element_id);</highlight></codeline>
<codeline lineno="95"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::FiniteElement<sp/>&amp;compactification_element<sp/>=<sp/>*compactification_space.GetFE(element_id);</highlight></codeline>
<codeline lineno="96"><highlight class="normal"></highlight></codeline>
<codeline lineno="97"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_displacement_data<sp/>=<sp/>std::make_unique&lt;ElementDisplacementData&gt;(<ref refid="namespaceserif_1_1mapping_1a89488756ce19ff5351dfbeb4235849c3" kindref="member">element_displacement_data_from_element_vdofs</ref>(displacement_element,<sp/>m_element_displacement));</highlight></codeline>
<codeline lineno="98"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_compactification_data<sp/>=<sp/>std::make_unique&lt;ElementCompactificationData&gt;(compactification_element,<sp/>m_element_compactification);</highlight></codeline>
<codeline lineno="99"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_cached_element_id<sp/>=<sp/>element_id;</highlight></codeline>
<codeline lineno="100"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="101"><highlight class="normal"></highlight></codeline>
<codeline lineno="102"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>GridFunctionMappingEvaluator::evaluate_point(mfem::ElementTransformation<sp/>&amp;transformation,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint<sp/>&amp;integration_point,<sp/>MappingPointContext<sp/>&amp;context)<sp/>{</highlight></codeline>
<codeline lineno="103"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>load_element(transformation.ElementNo);</highlight></codeline>
<codeline lineno="104"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementMappingData<sp/>data{</highlight></codeline>
<codeline lineno="105"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.displacement<sp/>=<sp/>*m_displacement_data,</highlight></codeline>
<codeline lineno="106"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.compactification<sp/>=<sp/>*m_compactification_data</highlight></codeline>
<codeline lineno="107"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="108"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>m_mapper.evaluate_point(data,<sp/>transformation,<sp/>integration_point,<sp/>m_workspace,<sp/>context);</highlight></codeline>
<codeline lineno="109"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="110"><highlight class="normal"></highlight></codeline>
<codeline lineno="111"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>GridFunctionMappingEvaluator::evaluate_volume(mfem::ElementTransformation<sp/>&amp;transformation,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint<sp/>&amp;integration_point,<sp/>VolumeMappingContext<sp/>&amp;context)<sp/>{</highlight></codeline>
<codeline lineno="112"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>load_element(transformation.ElementNo);</highlight></codeline>
<codeline lineno="113"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementMappingData<sp/>data{</highlight></codeline>
<codeline lineno="114"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.displacement<sp/>=<sp/>*m_displacement_data,</highlight></codeline>
<codeline lineno="115"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.compactification<sp/>=<sp/>*m_compactification_data</highlight></codeline>
<codeline lineno="116"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="117"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>m_mapper.evaluate_volume(data,<sp/>transformation,<sp/>integration_point,<sp/>m_workspace,<sp/>context);</highlight></codeline>
<codeline lineno="118"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="119"><highlight class="normal"></highlight></codeline>
<codeline lineno="120"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>GridFunctionMappingEvaluator::evaluate_face(mfem::FaceElementTransformations<sp/>&amp;transformations,<sp/>FaceElementSide<sp/>side,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint<sp/>&amp;integration_point,<sp/>FaceMappingContext<sp/>&amp;context)<sp/>{</highlight></codeline>
<codeline lineno="121"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::ElementTransformation<sp/>*transformation<sp/>=<sp/>side<sp/>==<sp/>FaceElementSide::ELEMENT_1<sp/>?<sp/>transformations.Elem1<sp/>:<sp/>transformations.Elem2;</highlight></codeline>
<codeline lineno="122"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>serif_ensure&lt;SERiFErrorCode::ELEMENT_NULLPTR_ERROR&gt;(transformation<sp/>!=<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">,<sp/></highlight><highlight class="stringliteral">&quot;Face<sp/>evaluation<sp/>requested<sp/>for<sp/>a<sp/>side<sp/>that<sp/>does<sp/>not<sp/>exist.&quot;</highlight><highlight class="normal">);</highlight></codeline>
<codeline lineno="123"><highlight class="normal"></highlight></codeline>
<codeline lineno="124"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>load_element(transformation-&gt;ElementNo);</highlight></codeline>
<codeline lineno="125"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementMappingData<sp/>data{</highlight></codeline>
<codeline lineno="126"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.displacement<sp/>=<sp/>*m_displacement_data,</highlight></codeline>
<codeline lineno="127"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.compactification<sp/>=<sp/>*m_compactification_data</highlight></codeline>
<codeline lineno="128"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="129"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>m_mapper.evaluate_face(data,<sp/>transformations,<sp/>side,<sp/>integration_point,<sp/>m_workspace,<sp/>context);</highlight></codeline>
<codeline lineno="130"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="131"><highlight class="normal">}</highlight></codeline>
    </programlisting>
    <location file="src/lib/serif/backend/mfem/mapping/evaluator.cpp"/>
  </compounddef>
</doxygen>