MeanField_SERiF/docs/xml/backend_2mfem_2discretization_2topology_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="backend_2mfem_2discretization_2topology_8cpp" kind="file" language="C++">
    <compoundname>topology.cpp</compoundname>
    <includes refid="backend_2mfem_2discretization_2topology_8hpp" local="yes">serif/backend/mfem/discretization/topology.hpp</includes>
    <includes local="no">utility</includes>
    <includes local="no">vector</includes>
    <includes local="yes">mfem/mesh/mesh.hpp</includes>
    <incdepgraph>
      <node id="2">
        <label>serif/backend/mfem/discretization/topology.hpp</label>
        <link refid="backend_2mfem_2discretization_2topology_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
      </node>
      <node id="3">
        <label>serif/discretization/domain/mesh/topology.hpp</label>
        <link refid="discretization_2domain_2mesh_2topology_8hpp"/>
        <childnode refid="4" relation="include">
        </childnode>
        <childnode refid="5" relation="include">
        </childnode>
      </node>
      <node id="1">
        <label>src/lib/serif/backend/mfem/discretization/topology.cpp</label>
        <link refid="backend_2mfem_2discretization_2topology_8cpp"/>
        <childnode refid="2" relation="include">
        </childnode>
        <childnode refid="6" relation="include">
        </childnode>
        <childnode refid="5" relation="include">
        </childnode>
        <childnode refid="7" relation="include">
        </childnode>
      </node>
      <node id="4">
        <label>cstddef</label>
      </node>
      <node id="7">
        <label>mfem/mesh/mesh.hpp</label>
      </node>
      <node id="6">
        <label>utility</label>
      </node>
      <node id="5">
        <label>vector</label>
      </node>
    </incdepgraph>
    <innernamespace refid="namespaceserif">serif</innernamespace>
    <innernamespace refid="namespaceserif_1_1backend">serif::backend</innernamespace>
    <innernamespace refid="namespaceserif_1_1backend_1_1mfem">serif::backend::mfem</innernamespace>
    <innernamespace refid="namespaceserif_1_1backend_1_1mfem_1_1discretization">serif::backend::mfem::discretization</innernamespace>
    <briefdescription>
    </briefdescription>
    <detaileddescription>
    </detaileddescription>
    <programlisting>
<codeline lineno="1"><highlight class="preprocessor">#include<sp/>&quot;<ref refid="backend_2mfem_2discretization_2topology_8hpp" kindref="compound">serif/backend/mfem/discretization/topology.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/>&lt;utility&gt;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="4"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&lt;vector&gt;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="5"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&quot;mfem/mesh/mesh.hpp&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_1backend_1_1mfem_1_1discretization" kindref="compound">serif::backend::mfem::discretization</ref><sp/>{</highlight></codeline>
<codeline lineno="8"><highlight class="normal"><sp/><sp/><sp/><sp/>::serif::discretization::domain::mesh::MeshTopology<sp/><ref refid="namespaceserif_1_1backend_1_1mfem_1_1discretization_1ad41b20dccb7312d80aa8e6f5282b8d8d" kindref="member">make_mesh_topology</ref>(const<sp/>::mfem::Mesh&amp;<sp/>mesh)<sp/>{</highlight></codeline>
<codeline lineno="9"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/>namespace<sp/></highlight><highlight class="normal">::serif::discretization::domain::mesh;</highlight></codeline>
<codeline lineno="10"><highlight class="normal"></highlight></codeline>
<codeline lineno="11"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>MeshTopologyData<sp/>data;</highlight></codeline>
<codeline lineno="12"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>data.elementDomainIDs.resize(mesh.GetNE());</highlight></codeline>
<codeline lineno="13"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::vector&lt;int&gt;<sp/>active_faces(mesh.GetNumFaces(),<sp/>-1);</highlight></codeline>
<codeline lineno="14"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>data.boundaryElementBoundaryIDs.resize(mesh.GetNBE());</highlight></codeline>
<codeline lineno="15"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>data.boundaryElementFaceIDs.resize(mesh.GetNBE());</highlight></codeline>
<codeline lineno="16"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>elementID<sp/>=<sp/>0;<sp/>elementID<sp/>&lt;<sp/>mesh.GetNE();<sp/>++elementID)<sp/>{</highlight></codeline>
<codeline lineno="17"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>data.elementDomainIDs[elementID]<sp/>=<sp/>mesh.GetAttribute(elementID);</highlight></codeline>
<codeline lineno="18"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="19"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>faceID<sp/>=<sp/>0;<sp/>faceID<sp/>&lt;<sp/>mesh.GetNumFaces();<sp/>++faceID)<sp/>{</highlight></codeline>
<codeline lineno="20"><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/>(mesh.Nonconforming()<sp/>&amp;&amp;<sp/>mesh.GetFaceInformation(faceID).tag<sp/>==<sp/>::mfem::Mesh::FaceInfoTag::MasterNonconforming)<sp/>{</highlight></codeline>
<codeline lineno="21"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Fine<sp/>faces<sp/>describe<sp/>the<sp/>actual<sp/>adjacency.<sp/>The<sp/>coarse<sp/>bookkeeping<sp/>face<sp/>is<sp/>not<sp/>exterior.</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="22"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">continue</highlight><highlight class="normal">;</highlight></codeline>
<codeline lineno="23"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="24"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>active_faces[faceID]<sp/>=<sp/></highlight><highlight class="keyword">static_cast&lt;</highlight><highlight class="keywordtype">int</highlight><highlight class="keyword">&gt;</highlight><highlight class="normal">(data.faceElements.size());</highlight></codeline>
<codeline lineno="25"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>data.faceElements.emplace_back();</highlight></codeline>
<codeline lineno="26"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mesh.GetFaceElements(faceID,<sp/>&amp;data.faceElements.back().firstElementID,<sp/>&amp;data.faceElements.back().secondElementID);</highlight></codeline>
<codeline lineno="27"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="28"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>boundaryElementID<sp/>=<sp/>0;<sp/>boundaryElementID<sp/>&lt;<sp/>mesh.GetNBE();<sp/>++boundaryElementID)<sp/>{</highlight></codeline>
<codeline lineno="29"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>data.boundaryElementBoundaryIDs[boundaryElementID]<sp/>=<sp/>mesh.GetBdrAttribute(boundaryElementID);</highlight></codeline>
<codeline lineno="30"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>data.boundaryElementFaceIDs[boundaryElementID]<sp/>=<sp/>active_faces[mesh.GetBdrElementFaceIndex(boundaryElementID)];</highlight></codeline>
<codeline lineno="31"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="32"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>MeshTopology(std::move(data));</highlight></codeline>
<codeline lineno="33"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="34"><highlight class="normal">}</highlight></codeline>
    </programlisting>
    <location file="src/lib/serif/backend/mfem/discretization/topology.cpp"/>
  </compounddef>
</doxygen>