Files

282 lines
21 KiB
XML
Raw Permalink Normal View History

<?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="schur_8cpp" kind="file" language="C++">
<compoundname>schur.cpp</compoundname>
<includes refid="schur_8hpp" local="yes">serif/solver/preconditioner/backend/mfem/gravity/schur.hpp</includes>
<includes local="no">cmath</includes>
<includes refid="solver_2preconditioner_2backend_2mfem_2inverse_2validation_8hpp" local="yes">serif/solver/preconditioner/backend/mfem/inverse/validation.hpp</includes>
<incdepgraph>
<node id="14">
<label>serif/meta/enums.hpp</label>
<link refid="enums_8hpp"/>
<childnode refid="15" relation="include">
</childnode>
<childnode refid="16" relation="include">
</childnode>
<childnode refid="17" relation="include">
</childnode>
</node>
<node id="2">
<label>serif/solver/preconditioner/backend/mfem/gravity/schur.hpp</label>
<link refid="schur_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>
</node>
<node id="6">
<label>serif/solver/preconditioner/backend/mfem/inverse.hpp</label>
<link refid="backend_2mfem_2inverse_8hpp"/>
<childnode refid="7" relation="include">
</childnode>
<childnode refid="8" relation="include">
</childnode>
<childnode refid="9" relation="include">
</childnode>
<childnode refid="10" relation="include">
</childnode>
</node>
<node id="19">
<label>serif/solver/preconditioner/backend/mfem/inverse/validation.hpp</label>
<link refid="solver_2preconditioner_2backend_2mfem_2inverse_2validation_8hpp"/>
<childnode refid="3" relation="include">
</childnode>
<childnode refid="20" relation="include">
</childnode>
<childnode refid="6" relation="include">
</childnode>
<childnode refid="21" relation="include">
</childnode>
</node>
<node id="9">
<label>serif/solver/preconditioner/contract.hpp</label>
<link refid="contract_8hpp"/>
</node>
<node id="10">
<label>serif/solver/preconditioner/failure.hpp</label>
<link refid="solver_2preconditioner_2failure_8hpp"/>
<childnode refid="11" relation="include">
</childnode>
<childnode refid="12" relation="include">
</childnode>
</node>
<node id="12">
<label>serif/utils/error/codes.hpp</label>
<link refid="codes_8hpp"/>
<childnode refid="13" relation="include">
</childnode>
<childnode refid="14" relation="include">
</childnode>
</node>
<node id="21">
<label>serif/utils/error/errors.hpp</label>
<link refid="errors_8hpp"/>
<childnode refid="20" relation="include">
</childnode>
<childnode refid="22" relation="include">
</childnode>
<childnode refid="23" relation="include">
</childnode>
<childnode refid="11" relation="include">
</childnode>
<childnode refid="13" relation="include">
</childnode>
<childnode refid="24" relation="include">
</childnode>
<childnode refid="25" relation="include">
</childnode>
<childnode refid="26" relation="include">
</childnode>
<childnode refid="30" relation="include">
</childnode>
<childnode refid="12" relation="include">
</childnode>
<childnode refid="34" relation="include">
</childnode>
</node>
<node id="30">
<label>serif/utils/misc/concepts/exceptions.hpp</label>
<link refid="utils_2misc_2concepts_2exceptions_8hpp"/>
<childnode refid="31" relation="include">
</childnode>
<childnode refid="16" relation="include">
</childnode>
<childnode refid="11" relation="include">
</childnode>
<childnode refid="32" relation="include">
</childnode>
</node>
<node id="32">
<label>serif/utils/misc/concepts/numeric.hpp</label>
<link refid="numeric_8hpp"/>
<childnode refid="33" relation="include">
</childnode>
<childnode refid="16" relation="include">
</childnode>
</node>
<node id="26">
<label>serif/utils/misc/terminal/colors.hpp</label>
<link refid="colors_8hpp"/>
<childnode refid="20" relation="include">
</childnode>
<childnode refid="23" relation="include">
</childnode>
<childnode refid="27" relation="include">
</childnode>
<childnode refid="28" relation="include">
</childnode>
<childnode refid="29" relation="include">
</childnode>
</node>
<node id="34">
<label>serif/utils/types/strings.hpp</label>
<link refid="types_2strings_8hpp"/>
<childnode refid="24" relation="include">
</childnode>
<childnode refid="35" relation="include">
</childnode>
<childnode refid="20" relation="include">
</childnode>
</node>
<node id="1">
<label>src/lib/serif/solver/preconditioner/backend/mfem/gravity/schur.cpp</label>
<link refid="schur_8cpp"/>
<childnode refid="2" relation="include">
</childnode>
<childnode refid="18" relation="include">
</childnode>
<childnode refid="19" relation="include">
</childnode>
</node>
<node id="24">
<label>algorithm</label>
</node>
<node id="18">
<label>cmath</label>
</node>
<node id="33">
<label>concepts</label>
</node>
<node id="13">
<label>cstdint</label>
</node>
<node id="23">
<label>cstdio</label>
</node>
<node id="27">
<label>cstdlib</label>
</node>
<node id="31">
<label>exception</label>
</node>
<node id="3">
<label>expected</label>
</node>
<node id="22">
<label>format</label>
</node>
<node id="35">
<label>iostream</label>
</node>
<node id="17">
<label>magic_enum/magic_enum_all.hpp</label>
</node>
<node id="7">
<label>memory</label>
</node>
<node id="8">
<label>mfem/linalg/operator.hpp</label>
</node>
<node id="5">
<label>mfem/linalg/sparsemat.hpp</label>
</node>
<node id="4">
<label>span</label>
</node>
<node id="11">
<label>string</label>
</node>
<node id="20">
<label>string_view</label>
</node>
<node id="16">
<label>type_traits</label>
</node>
<node id="29">
<label>unistd.h</label>
</node>
<node id="15">
<label>unordered_map</label>
</node>
<node id="28">
<label>utility</label>
</node>
<node id="25">
<label>version</label>
</node>
</incdepgraph>
<innernamespace refid="namespaceserif">serif</innernamespace>
<innernamespace refid="namespaceserif_1_1solver">serif::solver</innernamespace>
<innernamespace refid="namespaceserif_1_1solver_1_1preconditioner">serif::solver::preconditioner</innernamespace>
<innernamespace refid="namespaceserif_1_1solver_1_1preconditioner_1_1backend">serif::solver::preconditioner::backend</innernamespace>
<innernamespace refid="namespaceserif_1_1solver_1_1preconditioner_1_1backend_1_1mfem">serif::solver::preconditioner::backend::mfem</innernamespace>
<briefdescription>
</briefdescription>
<detaileddescription>
</detaileddescription>
<programlisting>
<codeline lineno="1"><highlight class="preprocessor">#include<sp/>&quot;<ref refid="schur_8hpp" kindref="compound">serif/solver/preconditioner/backend/mfem/gravity/schur.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;cmath&gt;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="4"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&quot;<ref refid="solver_2preconditioner_2backend_2mfem_2inverse_2validation_8hpp" kindref="compound">serif/solver/preconditioner/backend/mfem/inverse/validation.hpp</ref>&quot;</highlight><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_1solver_1_1preconditioner_1_1backend_1_1mfem" kindref="compound">serif::solver::preconditioner::backend::mfem</ref></highlight></codeline>
<codeline lineno="7"><highlight class="normal">{</highlight></codeline>
<codeline lineno="8"><highlight class="normal"><sp/><sp/><sp/><sp/>std::expected&lt;std::shared_ptr&lt;const<sp/>::mfem::SparseMatrix&gt;,<sp/>Failure&gt;<sp/><ref refid="namespaceserif_1_1solver_1_1preconditioner_1_1backend_1_1mfem_1aa1bd34ababff2c8a3fab4566f1e73b5e" kindref="member">try_gravity_schur</ref>(const<sp/>::mfem::SparseMatrix&amp;<sp/>divergence,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>std::span&lt;const<sp/>double&gt;<sp/>mapped_mass_diagonal)</highlight></codeline>
<codeline lineno="9"><highlight class="normal"><sp/><sp/><sp/><sp/>{</highlight></codeline>
<codeline lineno="10"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/></highlight><highlight class="normal">enum<sp/><ref refid="namespaceserif_1_1utils_1_1errors_1ab04484c13887f78337c2331fb2911346" kindref="member">utils::errors::SERiFErrorCode</ref>;</highlight></codeline>
<codeline lineno="11"><highlight class="normal"></highlight></codeline>
<codeline lineno="12"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1discretization_1_1backend_1_1mfem_1af75f98db90618153ea4ca169945733f7a1f064262ee7a1d63d490d9880aeb4aa0" kindref="member">divergence</ref>.Finalized()<sp/>||<sp/><ref refid="namespaceserif_1_1discretization_1_1backend_1_1mfem_1af75f98db90618153ea4ca169945733f7a1f064262ee7a1d63d490d9880aeb4aa0" kindref="member">divergence</ref>.Height()<sp/>&lt;=<sp/>0<sp/>||<sp/><ref refid="namespaceserif_1_1discretization_1_1backend_1_1mfem_1af75f98db90618153ea4ca169945733f7a1f064262ee7a1d63d490d9880aeb4aa0" kindref="member">divergence</ref>.Width()<sp/>&lt;=<sp/>0<sp/>||<sp/>mapped_mass_diagonal.size()<sp/>!=<sp/></highlight><highlight class="keyword">static_cast&lt;</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">&gt;</highlight><highlight class="normal">(<ref refid="namespaceserif_1_1discretization_1_1backend_1_1mfem_1af75f98db90618153ea4ca169945733f7a1f064262ee7a1d63d490d9880aeb4aa0" kindref="member">divergence</ref>.Width()))</highlight></codeline>
<codeline lineno="13"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>{</highlight></codeline>
<codeline lineno="14"><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_1solver_1_1preconditioner_1_1backend_1_1mfem_1_1detail_1a707023cbec2faf9488d727eb98373bbe" kindref="member">detail::failure&lt;PRECONDITIONER_SIZE_MISMATCH&gt;</ref>(</highlight><highlight class="stringliteral">&quot;The<sp/>gravity<sp/>Schur<sp/>surrogate<sp/>requires<sp/>finalized<sp/>computational<sp/>divergence<sp/>and<sp/>one<sp/>mapped<sp/>mass<sp/>entry<sp/>per<sp/>gradient<sp/>coordinate.&quot;</highlight><highlight class="normal">);</highlight></codeline>
<codeline lineno="15"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="16"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::Vector<sp/>reciprocal(<ref refid="namespaceserif_1_1discretization_1_1backend_1_1mfem_1af75f98db90618153ea4ca169945733f7a1f064262ee7a1d63d490d9880aeb4aa0" kindref="member">divergence</ref>.Width());</highlight></codeline>
<codeline lineno="17"><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/>index<sp/>=<sp/>0;<sp/><ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1ab297c6eedd8df878b6af2135b05e4616" kindref="member">index</ref><sp/>&lt;<sp/>reciprocal.Size();<sp/>++<ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1ab297c6eedd8df878b6af2135b05e4616" kindref="member">index</ref>)</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/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>entry<sp/>=<sp/>mapped_mass_diagonal[<ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1ab297c6eedd8df878b6af2135b05e4616" kindref="member">index</ref>];</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/>(!std::isfinite(entry)<sp/>||<sp/>entry<sp/>&lt;=<sp/>0.0<sp/>||<sp/>!std::isfinite(1.0<sp/>/<sp/>entry))</highlight></codeline>
<codeline lineno="21"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>{</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">return</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1solver_1_1preconditioner_1_1backend_1_1mfem_1_1detail_1a707023cbec2faf9488d727eb98373bbe" kindref="member">detail::failure&lt;PRECONDITIONER_INVALID_OPERATOR&gt;</ref>(</highlight><highlight class="stringliteral">&quot;Mapped<sp/>mass<sp/>entry<sp/>{}<sp/>must<sp/>be<sp/>positive,<sp/>finite<sp/>and<sp/>have<sp/>a<sp/>representable<sp/>reciprocal<sp/>before<sp/>forming<sp/>the<sp/>gravity<sp/>Schur<sp/>surrogate.&quot;</highlight><highlight class="normal">,<sp/>index);</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/>reciprocal[<ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1ab297c6eedd8df878b6af2135b05e4616" kindref="member">index</ref>]<sp/>=<sp/>1.0<sp/>/<sp/>entry;</highlight></codeline>
<codeline lineno="25"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="26"><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/>index<sp/>=<sp/>0;<sp/><ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1ab297c6eedd8df878b6af2135b05e4616" kindref="member">index</ref><sp/>&lt;<sp/><ref refid="namespaceserif_1_1discretization_1_1backend_1_1mfem_1af75f98db90618153ea4ca169945733f7a1f064262ee7a1d63d490d9880aeb4aa0" kindref="member">divergence</ref>.NumNonZeroElems();<sp/>++<ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1ab297c6eedd8df878b6af2135b05e4616" kindref="member">index</ref>)</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/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!std::isfinite(<ref refid="namespaceserif_1_1discretization_1_1backend_1_1mfem_1af75f98db90618153ea4ca169945733f7a1f064262ee7a1d63d490d9880aeb4aa0" kindref="member">divergence</ref>.GetData()[index]))</highlight></codeline>
<codeline lineno="29"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>{</highlight></codeline>
<codeline lineno="30"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1solver_1_1preconditioner_1_1backend_1_1mfem_1_1detail_1a707023cbec2faf9488d727eb98373bbe" kindref="member">detail::failure&lt;PRECONDITIONER_NONFINITE_VALUE&gt;</ref>(</highlight><highlight class="stringliteral">&quot;Computational<sp/>divergence<sp/>contains<sp/>a<sp/>nonfinite<sp/>coefficient.&quot;</highlight><highlight class="normal">);</highlight></codeline>
<codeline lineno="31"><highlight class="normal"><sp/><sp/><sp/><sp/><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></codeline>
<codeline lineno="33"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::SparseMatrix<sp/>weighted(divergence);</highlight></codeline>
<codeline lineno="34"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>weighted.ScaleColumns(reciprocal);</highlight></codeline>
<codeline lineno="35"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>std::unique_ptr&lt;::mfem::SparseMatrix&gt;<sp/>transpose(::mfem::Transpose(divergence));</highlight></codeline>
<codeline lineno="36"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>std::shared_ptr&lt;const<sp/>::mfem::SparseMatrix&gt;<sp/>schur(::mfem::Mult(weighted,<sp/>*transpose));</highlight></codeline>
<codeline lineno="37"><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/>index<sp/>=<sp/>0;<sp/><ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1ab297c6eedd8df878b6af2135b05e4616" kindref="member">index</ref><sp/>&lt;<sp/>schur-&gt;NumNonZeroElems();<sp/>++<ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1ab297c6eedd8df878b6af2135b05e4616" kindref="member">index</ref>)</highlight></codeline>
<codeline lineno="38"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>{</highlight></codeline>
<codeline lineno="39"><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/>(!std::isfinite(schur-&gt;GetData()[index]))</highlight></codeline>
<codeline lineno="40"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>{</highlight></codeline>
<codeline lineno="41"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1solver_1_1preconditioner_1_1backend_1_1mfem_1_1detail_1a707023cbec2faf9488d727eb98373bbe" kindref="member">detail::failure&lt;PRECONDITIONER_NONFINITE_VALUE&gt;</ref>(</highlight><highlight class="stringliteral">&quot;The<sp/>gravity<sp/>Schur<sp/>product<sp/>overflowed.<sp/>Rescale<sp/>the<sp/>field<sp/>coordinates<sp/>before<sp/>preparing<sp/>its<sp/>inverse.&quot;</highlight><highlight class="normal">);</highlight></codeline>
<codeline lineno="42"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="43"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="44"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>schur;</highlight></codeline>
<codeline lineno="45"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="46"><highlight class="normal">}<sp/></highlight><highlight class="comment">//<sp/>namespace<sp/>serif::solver::preconditioner::backend::mfem</highlight><highlight class="normal"></highlight></codeline>
</programlisting>
<location file="src/lib/serif/solver/preconditioner/backend/mfem/gravity/schur.cpp"/>
</compounddef>
</doxygen>