MeanField_SERiF/docs/xml/constrained__radial__extension_8hpp.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="constrained__radial__extension_8hpp" kind="file" language="C++">
    <compoundname>constrained_radial_extension.hpp</compoundname>
    <includes local="no">memory</includes>
    <includes local="no">span</includes>
    <includes local="no">utility</includes>
    <includes local="yes">mfem/fem/fespace.hpp</includes>
    <includes local="yes">mfem/linalg/sparsemat.hpp</includes>
    <includes refid="radial__extension_8hpp" local="yes">serif/geometry/deformation/radial_extension.hpp</includes>
    <includes refid="action__storage_8hpp" local="yes">serif/geometry/deformation/detail/action_storage.hpp</includes>
    <includedby refid="backend_2mfem_2geometry_2deformation_2computations_8hpp" local="yes">src/include/serif/backend/mfem/geometry/deformation/computations.hpp</includedby>
    <incdepgraph>
      <node id="1">
        <label>src/include/serif/backend/mfem/geometry/deformation/constrained_radial_extension.hpp</label>
        <link refid="constrained__radial__extension_8hpp"/>
        <childnode refid="2" relation="include">
        </childnode>
        <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="18" relation="include">
        </childnode>
      </node>
      <node id="18">
        <label>serif/geometry/deformation/detail/action_storage.hpp</label>
        <link refid="action__storage_8hpp"/>
        <childnode refid="19" relation="include">
        </childnode>
        <childnode refid="20" relation="include">
        </childnode>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
      </node>
      <node id="10">
        <label>serif/geometry/deformation/failure.hpp</label>
        <link refid="geometry_2deformation_2failure_8hpp"/>
        <childnode refid="11" relation="include">
        </childnode>
        <childnode refid="12" relation="include">
        </childnode>
      </node>
      <node id="13">
        <label>serif/geometry/deformation/models.hpp</label>
        <link refid="models_8hpp"/>
        <childnode refid="14" relation="include">
        </childnode>
        <childnode refid="15" relation="include">
        </childnode>
      </node>
      <node id="7">
        <label>serif/geometry/deformation/radial_extension.hpp</label>
        <link refid="radial__extension_8hpp"/>
        <childnode refid="2" relation="include">
        </childnode>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="8" relation="include">
        </childnode>
        <childnode refid="9" relation="include">
        </childnode>
        <childnode refid="13" relation="include">
        </childnode>
        <childnode refid="16" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
      </node>
      <node id="9">
        <label>serif/geometry/deformation/radial_interpolation.hpp</label>
        <link refid="geometry_2deformation_2radial__interpolation_8hpp"/>
        <childnode refid="4" relation="include">
        </childnode>
        <childnode refid="8" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
      </node>
      <node id="16">
        <label>serif/geometry/deformation/spaces.hpp</label>
        <link refid="spaces_8hpp"/>
      </node>
      <node id="17">
        <label>serif/geometry/deformation/workspace.hpp</label>
        <link refid="geometry_2deformation_2workspace_8hpp"/>
      </node>
      <node id="19">
        <label>algorithm</label>
      </node>
      <node id="14">
        <label>array</label>
      </node>
      <node id="20">
        <label>cmath</label>
      </node>
      <node id="11">
        <label>expected</label>
      </node>
      <node id="21">
        <label>functional</label>
      </node>
      <node id="22">
        <label>initializer_list</label>
      </node>
      <node id="2">
        <label>memory</label>
      </node>
      <node id="5">
        <label>mfem/fem/fespace.hpp</label>
      </node>
      <node id="6">
        <label>mfem/linalg/sparsemat.hpp</label>
      </node>
      <node id="15">
        <label>optional</label>
      </node>
      <node id="3">
        <label>span</label>
      </node>
      <node id="12">
        <label>string</label>
      </node>
      <node id="4">
        <label>utility</label>
      </node>
      <node id="8">
        <label>vector</label>
      </node>
    </incdepgraph>
    <invincdepgraph>
      <node id="5">
        <label>src/compile_time_checks/serif/model/state.cpp</label>
        <link refid="state_8cpp"/>
      </node>
      <node id="2">
        <label>src/include/serif/backend/mfem/geometry/deformation/computations.hpp</label>
        <link refid="backend_2mfem_2geometry_2deformation_2computations_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="12" relation="include">
        </childnode>
        <childnode refid="13" relation="include">
        </childnode>
        <childnode refid="15" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="20" relation="include">
        </childnode>
      </node>
      <node id="1">
        <label>src/include/serif/backend/mfem/geometry/deformation/constrained_radial_extension.hpp</label>
        <link refid="constrained__radial__extension_8hpp"/>
        <childnode refid="2" relation="include">
        </childnode>
      </node>
      <node id="4">
        <label>src/include/serif/backend/mfem/model/equilibrium.hpp</label>
        <link refid="backend_2mfem_2model_2equilibrium_8hpp"/>
        <childnode refid="5" relation="include">
        </childnode>
        <childnode refid="6" relation="include">
        </childnode>
      </node>
      <node id="14">
        <label>src/include/serif/backend/mfem/model/equilibrium_evaluator.hpp</label>
        <link refid="equilibrium__evaluator_8hpp"/>
        <childnode refid="4" relation="include">
        </childnode>
      </node>
      <node id="3">
        <label>src/include/serif/backend/mfem/model/graph.hpp</label>
        <link refid="backend_2mfem_2model_2graph_8hpp"/>
        <childnode refid="4" relation="include">
        </childnode>
        <childnode refid="11" relation="include">
        </childnode>
      </node>
      <node id="12">
        <label>src/include/serif/backend/mfem/model/normalization/computations.hpp</label>
        <link refid="backend_2mfem_2model_2normalization_2computations_8hpp"/>
        <childnode refid="4" relation="include">
        </childnode>
      </node>
      <node id="13">
        <label>src/include/serif/backend/mfem/model/preconditioning/computations.hpp</label>
        <link refid="backend_2mfem_2model_2preconditioning_2computations_8hpp"/>
        <childnode refid="14" relation="include">
        </childnode>
      </node>
      <node id="15">
        <label>src/include/serif/backend/mfem/serialize/top.hpp</label>
        <link refid="backend_2mfem_2serialize_2top_8hpp"/>
        <childnode refid="16" relation="include">
        </childnode>
      </node>
      <node id="11">
        <label>src/include/serif/backend/mfem/solver/admissibility/geometry.hpp</label>
        <link refid="geometry_8hpp"/>
        <childnode refid="6" relation="include">
        </childnode>
      </node>
      <node id="7">
        <label>src/include/serif/backend/mfem/solver/equilibrium.hpp</label>
        <link refid="backend_2mfem_2solver_2equilibrium_8hpp"/>
        <childnode refid="8" relation="include">
        </childnode>
      </node>
      <node id="6">
        <label>src/include/serif/backend/mfem/solver/nonlinear/stellar_problem.hpp</label>
        <link refid="stellar__problem_8hpp"/>
        <childnode refid="7" relation="include">
        </childnode>
      </node>
      <node id="17">
        <label>src/include/serif/backend/mfem/solver/seed/fitted.hpp</label>
        <link refid="backend_2mfem_2solver_2seed_2fitted_8hpp"/>
        <childnode refid="18" relation="include">
        </childnode>
        <childnode refid="19" relation="include">
        </childnode>
      </node>
      <node id="18">
        <label>src/include/serif/backend/mfem/solver/seed/state.hpp</label>
        <link refid="backend_2mfem_2solver_2seed_2state_8hpp"/>
        <childnode refid="19" relation="include">
        </childnode>
      </node>
      <node id="20">
        <label>src/include/serif/geometry/deformation/finite_elements.hpp</label>
        <link refid="geometry_2deformation_2finite__elements_8hpp"/>
      </node>
      <node id="9">
        <label>src/include/serif/solver/continuation/equilibrium.hpp</label>
        <link refid="solver_2continuation_2equilibrium_8hpp"/>
        <childnode refid="10" relation="include">
        </childnode>
      </node>
      <node id="8">
        <label>src/include/serif/solver/equilibrium.hpp</label>
        <link refid="solver_2equilibrium_8hpp"/>
        <childnode refid="9" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
      </node>
      <node id="19">
        <label>src/include/serif/solver/seed/state.hpp</label>
        <link refid="solver_2seed_2state_8hpp"/>
        <childnode refid="10" relation="include">
        </childnode>
      </node>
      <node id="10">
        <label>src/include/serif/solver/solver.hpp</label>
        <link refid="solver_2solver_8hpp"/>
      </node>
      <node id="16">
        <label>src/lib/serif/backend/mfem/serialize/top.cpp</label>
        <link refid="backend_2mfem_2serialize_2top_8cpp"/>
      </node>
    </invincdepgraph>
    <innerclass refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension" prot="public">serif::geometry::deformation::backend::mfem::ConstrainedRadialExtension</innerclass>
    <innerclass refid="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1_1_workspace" prot="public">serif::geometry::deformation::backend::mfem::ConstrainedRadialExtension::Workspace</innerclass>
    <innernamespace refid="namespaceserif">serif</innernamespace>
    <innernamespace refid="namespaceserif_1_1geometry">serif::geometry</innernamespace>
    <innernamespace refid="namespaceserif_1_1geometry_1_1deformation">serif::geometry::deformation</innernamespace>
    <innernamespace refid="namespaceserif_1_1geometry_1_1deformation_1_1backend">serif::geometry::deformation::backend</innernamespace>
    <innernamespace refid="namespaceserif_1_1geometry_1_1deformation_1_1backend_1_1mfem">serif::geometry::deformation::backend::mfem</innernamespace>
    <briefdescription>
    </briefdescription>
    <detaileddescription>
    </detaileddescription>
    <programlisting>
<codeline lineno="1"><highlight class="preprocessor">#pragma<sp/>once</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;memory&gt;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="4"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&lt;span&gt;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="5"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&lt;utility&gt;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="6"><highlight class="normal"></highlight></codeline>
<codeline lineno="7"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&quot;mfem/fem/fespace.hpp&quot;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="8"><highlight class="normal"></highlight><highlight class="preprocessor">#ifdef<sp/>MFEM_USE_MPI</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="9"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&quot;mfem/fem/pfespace.hpp&quot;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="10"><highlight class="normal"></highlight><highlight class="preprocessor">#endif</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="11"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&quot;mfem/linalg/sparsemat.hpp&quot;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="12"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&quot;<ref refid="radial__extension_8hpp" kindref="compound">serif/geometry/deformation/radial_extension.hpp</ref>&quot;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="13"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&quot;<ref refid="action__storage_8hpp" kindref="compound">serif/geometry/deformation/detail/action_storage.hpp</ref>&quot;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="14"><highlight class="normal"></highlight></codeline>
<codeline lineno="15"><highlight class="normal"></highlight><highlight class="keyword">namespace<sp/></highlight><highlight class="normal"><ref refid="namespaceserif_1_1geometry_1_1deformation_1_1backend_1_1mfem" kindref="compound">serif::geometry::deformation::backend::mfem</ref><sp/>{</highlight></codeline>
<codeline lineno="21" refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension" refkind="compound"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">class<sp/></highlight><highlight class="normal"><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1ae80591257dccff5948ea1bf361f5361d" kindref="member">ConstrainedRadialExtension</ref><sp/>{</highlight></codeline>
<codeline lineno="22"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">public</highlight><highlight class="normal">:</highlight></codeline>
<codeline lineno="23"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/></highlight><highlight class="normal"><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a78f1793cc8a36ca7036bb915820eacbb" kindref="member">InputSpace</ref><sp/>=<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1_radial_extension_1abdb8ed08c56c90bc7f195cec6ce518dc" kindref="member">RadialExtension::InputSpace</ref>;</highlight></codeline>
<codeline lineno="24"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/></highlight><highlight class="normal"><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a57726f110c4832e21b9573cc8ad71aff" kindref="member">OutputSpace</ref><sp/>=<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1_radial_extension_1a48b6a881cc9b709d3ee74e1fd73d96bc" kindref="member">RadialExtension::OutputSpace</ref>;</highlight></codeline>
<codeline lineno="25"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">static</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">constexpr</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">bool</highlight><highlight class="normal"><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a1d245d42e6ddae89d177e4df21a58931" kindref="member">is_linear</ref><sp/>=<sp/></highlight><highlight class="keyword">true</highlight><highlight class="normal">;</highlight></codeline>
<codeline lineno="26"><highlight class="normal"></highlight></codeline>
<codeline lineno="28" refid="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1_1_workspace" refkind="compound"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">struct<sp/></highlight><highlight class="normal"><ref refid="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1_1_workspace" kindref="compound">Workspace</ref><sp/>{</highlight></codeline>
<codeline lineno="29"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::Vector<sp/><ref refid="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1_1_workspace_1a0833ebad427c658beb06430db8107e92" kindref="member">local</ref>;</highlight></codeline>
<codeline lineno="30"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::Vector<sp/><ref refid="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1_1_workspace_1aa830b65d5cdb67345b61d839ce4f4909" kindref="member">independent</ref>;</highlight></codeline>
<codeline lineno="31"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="32"><highlight class="normal"></highlight></codeline>
<codeline lineno="34"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>std::shared_ptr&lt;const<sp/>RadialInterpolation&gt;&amp;<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a6151f50ac40bd03b6d4fd65fb376b68b" kindref="member">interpolation</ref>()<sp/>const<sp/>noexcept<sp/>{</highlight></codeline>
<codeline lineno="35"><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="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a53f1ed40b9c1baeb98bea41e8712d78d" kindref="member">m_radial</ref>.interpolation();</highlight></codeline>
<codeline lineno="36"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="37"><highlight class="normal"></highlight></codeline>
<codeline lineno="39"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/>std::span&lt;const<sp/>double&gt;<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a3133e300456b4b8be7ac644e84fb3108" kindref="member">radial_weights</ref>()<sp/>const<sp/>noexcept<sp/>{</highlight></codeline>
<codeline lineno="40"><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="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a53f1ed40b9c1baeb98bea41e8712d78d" kindref="member">m_radial</ref>.<ref refid="classserif_1_1geometry_1_1deformation_1_1_radial_extension_1acf26b497fe74d7dd8c28462bb64dabd7" kindref="member">radial_weights</ref>();</highlight></codeline>
<codeline lineno="41"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="42"><highlight class="normal"></highlight></codeline>
<codeline lineno="44"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/>std::size_t<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a3816438e98f5165b3bbf9060e4bb0ebf" kindref="member">input_size</ref>()<sp/>const<sp/>noexcept<sp/>{</highlight></codeline>
<codeline lineno="45"><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="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a53f1ed40b9c1baeb98bea41e8712d78d" kindref="member">m_radial</ref>.<ref refid="classserif_1_1geometry_1_1deformation_1_1_radial_extension_1ad7b113fe439b95bf7636b653b0047e86" kindref="member">input_size</ref>();</highlight></codeline>
<codeline lineno="46"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="47"><highlight class="normal"></highlight></codeline>
<codeline lineno="49"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/>std::size_t<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a6f7ab32c2cad0ab4904627705802a9d3" kindref="member">output_size</ref>()<sp/>const<sp/>noexcept<sp/>{</highlight></codeline>
<codeline lineno="50"><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="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a53f1ed40b9c1baeb98bea41e8712d78d" kindref="member">m_radial</ref>.<ref refid="classserif_1_1geometry_1_1deformation_1_1_radial_extension_1ac4a58eec559d81378ca020dfe156fa81" kindref="member">output_size</ref>();</highlight></codeline>
<codeline lineno="51"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="52"><highlight class="normal"></highlight></codeline>
<codeline lineno="54"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/>Workspace<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a79f4075553947ee93eb4266ea318e606" kindref="member">make_workspace</ref>()<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal">;</highlight></codeline>
<codeline lineno="55"><highlight class="normal"></highlight></codeline>
<codeline lineno="57"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a365909946a576f7dfcdc904300f7214c" kindref="member">apply</ref>(std::span&lt;const<sp/>double&gt;<sp/><ref refid="namespaceserif_1_1geometry_1_1surface" kindref="compound">surface</ref>,<sp/>std::span&lt;double&gt;<sp/>volume)<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal">;</highlight></codeline>
<codeline lineno="58"><highlight class="normal"></highlight></codeline>
<codeline lineno="60"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1aa7d29fe02f7649e3917ff2b7b69f33f4" kindref="member">apply_transpose</ref>(std::span&lt;const<sp/>double&gt;<sp/>volume_force,<sp/>std::span&lt;double&gt;<sp/>surface_force)<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal">;</highlight></codeline>
<codeline lineno="61"><highlight class="normal"></highlight></codeline>
<codeline lineno="63"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a365909946a576f7dfcdc904300f7214c" kindref="member">apply</ref>(std::span&lt;const<sp/>double&gt;<sp/><ref refid="namespaceserif_1_1geometry_1_1surface" kindref="compound">surface</ref>,<sp/>std::span&lt;double&gt;<sp/>volume,<sp/>Workspace&amp;<sp/>workspace)<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal">;</highlight></codeline>
<codeline lineno="64"><highlight class="normal"></highlight></codeline>
<codeline lineno="66"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1aa7d29fe02f7649e3917ff2b7b69f33f4" kindref="member">apply_transpose</ref>(std::span&lt;const<sp/>double&gt;<sp/>volume_force,<sp/>std::span&lt;double&gt;<sp/>surface_force,<sp/>Workspace&amp;<sp/>workspace)<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal">;</highlight></codeline>
<codeline lineno="67"><highlight class="normal"></highlight></codeline>
<codeline lineno="68"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">private</highlight><highlight class="normal">:</highlight></codeline>
<codeline lineno="70"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1ae80591257dccff5948ea1bf361f5361d" kindref="member">ConstrainedRadialExtension</ref>(RadialExtension<sp/>radial,<sp/>std::unique_ptr&lt;const<sp/>::mfem::SparseMatrix&gt;<sp/>prolongation,<sp/>std::unique_ptr&lt;const<sp/>::mfem::SparseMatrix&gt;<sp/>restriction);</highlight></codeline>
<codeline lineno="71"><highlight class="normal"></highlight></codeline>
<codeline lineno="73"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a5e89e83f6ce00f805f5e802c7d7676de" kindref="member">validate_workspace</ref>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>Workspace&amp;<sp/>workspace,<sp/>std::span&lt;const<sp/>double&gt;<sp/>input,<sp/>std::span&lt;double&gt;<sp/>output)<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal">;</highlight></codeline>
<codeline lineno="74"><highlight class="normal"></highlight></codeline>
<codeline lineno="75"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">friend</highlight><highlight class="normal"><sp/>std::expected&lt;ConstrainedRadialExtension,<sp/>DeformationFailure&gt;<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a27b36315f46f337ca2361833767b82c7" kindref="member">make_constrained_radial_extension</ref>(</highlight></codeline>
<codeline lineno="76"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::shared_ptr&lt;const<sp/>RadialInterpolation&gt;,<sp/>const<sp/>::mfem::FiniteElementSpace&amp;,<sp/>PowerLawRadialExtension,<sp/>FixedInfinityRadialExtension);</highlight></codeline>
<codeline lineno="77"><highlight class="normal"></highlight></codeline>
<codeline lineno="78"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>RadialExtension<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a53f1ed40b9c1baeb98bea41e8712d78d" kindref="member">m_radial</ref>;</highlight></codeline>
<codeline lineno="79"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::unique_ptr&lt;const<sp/>::mfem::SparseMatrix&gt;<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1ac4b00cd302ad31d545f4e5ebde67c8e5" kindref="member">m_prolongation</ref>;</highlight></codeline>
<codeline lineno="80"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::unique_ptr&lt;const<sp/>::mfem::SparseMatrix&gt;<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a74d566ff3095fc95ad230e0fa18f2c84" kindref="member">m_restriction</ref>;</highlight></codeline>
<codeline lineno="81"><highlight class="normal"><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="82"><highlight class="normal"></highlight></codeline>
<codeline lineno="84"><highlight class="normal"><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/>std::expected&lt;ConstrainedRadialExtension,<sp/>DeformationFailure&gt;<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1afd8573b009b4014cbb94493addfaf526" kindref="member">make_constrained_radial_extension</ref>(</highlight></codeline>
<codeline lineno="85"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::shared_ptr&lt;const<sp/>RadialInterpolation&gt;<sp/>interpolation,</highlight></codeline>
<codeline lineno="86"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>const<sp/>::mfem::FiniteElementSpace&amp;<sp/>displacement_space,</highlight></codeline>
<codeline lineno="87"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>PowerLawRadialExtension<sp/>interior<sp/>=<sp/>{},</highlight></codeline>
<codeline lineno="88"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>FixedInfinityRadialExtension<sp/>vacuum<sp/>=<sp/>{});</highlight></codeline>
<codeline lineno="89"><highlight class="normal">}<sp/></highlight><highlight class="comment">//<sp/>namespace<sp/>serif::geometry::deformation::backend::mfem</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="90"><highlight class="normal"></highlight></codeline>
<codeline lineno="91"><highlight class="normal"></highlight><highlight class="keyword">namespace<sp/></highlight><highlight class="normal"><ref refid="namespaceserif_1_1geometry_1_1deformation_1_1backend_1_1mfem" kindref="compound">serif::geometry::deformation::backend::mfem</ref><sp/>{</highlight></codeline>
<codeline lineno="92"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">inline</highlight><highlight class="normal"><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1ae80591257dccff5948ea1bf361f5361d" kindref="member">ConstrainedRadialExtension::ConstrainedRadialExtension</ref>(RadialExtension<sp/>radial,</highlight></codeline>
<codeline lineno="93"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::unique_ptr&lt;const<sp/>::mfem::SparseMatrix&gt;<sp/>prolongation,<sp/>std::unique_ptr&lt;const<sp/>::mfem::SparseMatrix&gt;<sp/>restriction)<sp/>:</highlight></codeline>
<codeline lineno="94"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>m_radial(std::move(radial)),<sp/>m_prolongation(std::move(prolongation)),<sp/>m_restriction(std::move(restriction))<sp/>{}</highlight></codeline>
<codeline lineno="95"><highlight class="normal"></highlight></codeline>
<codeline lineno="96"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">inline</highlight><highlight class="normal"><sp/>ConstrainedRadialExtension::Workspace<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a79f4075553947ee93eb4266ea318e606" kindref="member">ConstrainedRadialExtension::make_workspace</ref>()</highlight><highlight class="keyword"><sp/>const<sp/></highlight><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/>(!<ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1ac4b00cd302ad31d545f4e5ebde67c8e5" kindref="member">m_prolongation</ref>)<sp/>{</highlight></codeline>
<codeline lineno="98"><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></codeline>
<codeline lineno="99"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="100"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>{.local<sp/>=<sp/>::mfem::Vector(</highlight><highlight class="keyword">static_cast&lt;</highlight><highlight class="keywordtype">int</highlight><highlight class="keyword">&gt;</highlight><highlight class="normal">(<ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a6f7ab32c2cad0ab4904627705802a9d3" kindref="member">output_size</ref>())),<sp/>.independent<sp/>=<sp/>::mfem::Vector(<ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a74d566ff3095fc95ad230e0fa18f2c84" kindref="member">m_restriction</ref>-&gt;Height())};</highlight></codeline>
<codeline lineno="101"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="102"><highlight class="normal"></highlight></codeline>
<codeline lineno="103"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">inline</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a365909946a576f7dfcdc904300f7214c" kindref="member">ConstrainedRadialExtension::apply</ref>(std::span&lt;const<sp/>double&gt;<sp/>surface,<sp/>std::span&lt;double&gt;<sp/>volume)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="104"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1_1_workspace" kindref="compound">Workspace</ref><sp/>workspace<sp/>=<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a79f4075553947ee93eb4266ea318e606" kindref="member">make_workspace</ref>();</highlight></codeline>
<codeline lineno="105"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a365909946a576f7dfcdc904300f7214c" kindref="member">apply</ref>(surface,<sp/>volume,<sp/>workspace);</highlight></codeline>
<codeline lineno="106"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="107"><highlight class="normal"></highlight></codeline>
<codeline lineno="108"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">inline</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1aa7d29fe02f7649e3917ff2b7b69f33f4" kindref="member">ConstrainedRadialExtension::apply_transpose</ref>(std::span&lt;const<sp/>double&gt;<sp/>volume_force,<sp/>std::span&lt;double&gt;<sp/>surface_force)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="109"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1_1_workspace" kindref="compound">Workspace</ref><sp/>workspace<sp/>=<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a79f4075553947ee93eb4266ea318e606" kindref="member">make_workspace</ref>();</highlight></codeline>
<codeline lineno="110"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1aa7d29fe02f7649e3917ff2b7b69f33f4" kindref="member">apply_transpose</ref>(volume_force,<sp/>surface_force,<sp/>workspace);</highlight></codeline>
<codeline lineno="111"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="112"><highlight class="normal"></highlight></codeline>
<codeline lineno="113"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">inline</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a5e89e83f6ce00f805f5e802c7d7676de" kindref="member">ConstrainedRadialExtension::validate_workspace</ref>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1_1_workspace" kindref="compound">Workspace</ref>&amp;<sp/>workspace,</highlight></codeline>
<codeline lineno="114"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::span&lt;const<sp/>double&gt;<sp/>input,<sp/>std::span&lt;double&gt;<sp/>output)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="115"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(workspace.local.Size()<sp/>!=<sp/></highlight><highlight class="keyword">static_cast&lt;</highlight><highlight class="keywordtype">int</highlight><highlight class="keyword">&gt;</highlight><highlight class="normal">(<ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a6f7ab32c2cad0ab4904627705802a9d3" kindref="member">output_size</ref>())<sp/>||<sp/>workspace.independent.Size()<sp/>!=<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a74d566ff3095fc95ad230e0fa18f2c84" kindref="member">m_restriction</ref>-&gt;Height())<sp/>{</highlight></codeline>
<codeline lineno="116"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(DeformationFailure{.reason<sp/>=<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1abff36c0dd3d8b9c916a2e928ab69a871a0181b84c792eb178181854d09c8fe7b9" kindref="member">DeformationFailureReason::invalid_size</ref>,<sp/>.message<sp/>=<sp/></highlight><highlight class="stringliteral">&quot;Constrained<sp/>radial<sp/>action<sp/>requires<sp/>scratch<sp/>allocated<sp/>for<sp/>this<sp/>extension.&quot;</highlight><highlight class="normal">});</highlight></codeline>
<codeline lineno="117"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="118"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>std::span&lt;double&gt;<sp/>local(workspace.local.GetData(),<sp/>workspace.local.Size());</highlight></codeline>
<codeline lineno="119"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>std::span&lt;double&gt;<sp/><ref refid="namespaceserif_1_1solver_1_1preconditioner_1_1strategies_1_1coupling_1a1843f298d4e3fb5770085697b48ffa55ae2580777557122a236c803c1ccad7fa4" kindref="member">independent</ref>(workspace.independent.GetData(),<sp/>workspace.independent.Size());</highlight></codeline>
<codeline lineno="120"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/>status<sp/>=<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1_1detail_1a9bb57f09e703c8a56d6ce6bd658a752f" kindref="member">deformation::detail::validate_disjoint_output</ref>(local,<sp/>{input,<sp/>output,<sp/><ref refid="namespaceserif_1_1solver_1_1preconditioner_1_1strategies_1_1coupling_1a1843f298d4e3fb5770085697b48ffa55ae2580777557122a236c803c1ccad7fa4" kindref="member">independent</ref>});<sp/>!status)<sp/>{</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/>status;</highlight></codeline>
<codeline lineno="122"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="123"><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_1geometry_1_1deformation_1_1detail_1a9bb57f09e703c8a56d6ce6bd658a752f" kindref="member">deformation::detail::validate_disjoint_output</ref>(independent,<sp/>{input,<sp/>output});</highlight></codeline>
<codeline lineno="124"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="125"><highlight class="normal"></highlight></codeline>
<codeline lineno="126"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">inline</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a365909946a576f7dfcdc904300f7214c" kindref="member">ConstrainedRadialExtension::apply</ref>(std::span&lt;const<sp/>double&gt;<sp/>surface,</highlight></codeline>
<codeline lineno="127"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::span&lt;double&gt;<sp/>volume,<sp/><ref refid="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1_1_workspace" kindref="compound">Workspace</ref>&amp;<sp/>workspace)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="128"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1ac4b00cd302ad31d545f4e5ebde67c8e5" kindref="member">m_prolongation</ref>)<sp/>{</highlight></codeline>
<codeline lineno="129"><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="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a53f1ed40b9c1baeb98bea41e8712d78d" kindref="member">m_radial</ref>.apply(surface,<sp/>volume);</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/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/>status<sp/>=<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1_1detail_1a39ae230635f7915d0cb1942b333aac6f" kindref="member">deformation::detail::validate_action_arguments</ref>(volume,<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a6f7ab32c2cad0ab4904627705802a9d3" kindref="member">output_size</ref>(),<sp/>{{<ref refid="namespaceserif_1_1solver_1_1normalization_1ac75238439fd5d1610b177bf85e9376d9ad302e976c9dd527a9521a88c012437c5" kindref="member">surface</ref>,<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a3816438e98f5165b3bbf9060e4bb0ebf" kindref="member">input_size</ref>()}});<sp/>!status)<sp/>{</highlight></codeline>
<codeline lineno="132"><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/>status;</highlight></codeline>
<codeline lineno="133"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="134"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/>status<sp/>=<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a5e89e83f6ce00f805f5e802c7d7676de" kindref="member">validate_workspace</ref>(workspace,<sp/>surface,<sp/>volume);<sp/>!status)<sp/>{</highlight></codeline>
<codeline lineno="135"><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/>status;</highlight></codeline>
<codeline lineno="136"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="137"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/>status<sp/>=<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a53f1ed40b9c1baeb98bea41e8712d78d" kindref="member">m_radial</ref>.apply(surface,<sp/>{workspace.local.GetData(),<sp/>static_cast&lt;std::size_t&gt;(workspace.local.Size())});<sp/>!status)<sp/>{</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/>status;</highlight></codeline>
<codeline lineno="139"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="140"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a74d566ff3095fc95ad230e0fa18f2c84" kindref="member">m_restriction</ref>-&gt;Mult(workspace.local,<sp/>workspace.independent);</highlight></codeline>
<codeline lineno="141"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::Vector<sp/>output(<ref refid="namespaceserif_1_1solver_1_1normalization_1ac75238439fd5d1610b177bf85e9376d9a210ab9e731c9c36c2c38db15c28a8d1c" kindref="member">volume</ref>.data(),<sp/></highlight><highlight class="keyword">static_cast&lt;</highlight><highlight class="keywordtype">int</highlight><highlight class="keyword">&gt;</highlight><highlight class="normal">(<ref refid="namespaceserif_1_1solver_1_1normalization_1ac75238439fd5d1610b177bf85e9376d9a210ab9e731c9c36c2c38db15c28a8d1c" kindref="member">volume</ref>.size()));</highlight></codeline>
<codeline lineno="142"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1ac4b00cd302ad31d545f4e5ebde67c8e5" kindref="member">m_prolongation</ref>-&gt;Mult(workspace.independent,<sp/>output);</highlight></codeline>
<codeline lineno="143"><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_1geometry_1_1deformation_1a371fc6cc353f2df9f7752d802c2155da" kindref="member">validate_deformation_values</ref>(volume,<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a6f7ab32c2cad0ab4904627705802a9d3" kindref="member">output_size</ref>());</highlight></codeline>
<codeline lineno="144"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="145"><highlight class="normal"></highlight></codeline>
<codeline lineno="146"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">inline</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1aa7d29fe02f7649e3917ff2b7b69f33f4" kindref="member">ConstrainedRadialExtension::apply_transpose</ref>(std::span&lt;const<sp/>double&gt;<sp/>volume_force,</highlight></codeline>
<codeline lineno="147"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::span&lt;double&gt;<sp/>surface_force,<sp/><ref refid="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1_1_workspace" kindref="compound">Workspace</ref>&amp;<sp/>workspace)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="148"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1ac4b00cd302ad31d545f4e5ebde67c8e5" kindref="member">m_prolongation</ref>)<sp/>{</highlight></codeline>
<codeline lineno="149"><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="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a53f1ed40b9c1baeb98bea41e8712d78d" kindref="member">m_radial</ref>.apply_transpose(volume_force,<sp/>surface_force);</highlight></codeline>
<codeline lineno="150"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="151"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/>status<sp/>=<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1_1detail_1a39ae230635f7915d0cb1942b333aac6f" kindref="member">deformation::detail::validate_action_arguments</ref>(surface_force,<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a3816438e98f5165b3bbf9060e4bb0ebf" kindref="member">input_size</ref>(),<sp/>{{volume_force,<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a6f7ab32c2cad0ab4904627705802a9d3" kindref="member">output_size</ref>()}});<sp/>!status)<sp/>{</highlight></codeline>
<codeline lineno="152"><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/>status;</highlight></codeline>
<codeline lineno="153"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="154"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/>status<sp/>=<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a5e89e83f6ce00f805f5e802c7d7676de" kindref="member">validate_workspace</ref>(workspace,<sp/>volume_force,<sp/>surface_force);<sp/>!status)<sp/>{</highlight></codeline>
<codeline lineno="155"><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/>status;</highlight></codeline>
<codeline lineno="156"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="157"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>MFEM&apos;s<sp/>vector<sp/>view<sp/>requires<sp/>mutable<sp/>storage;<sp/>matrix<sp/>application<sp/>only<sp/>reads<sp/>this<sp/>force.</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="158"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>const<sp/>::mfem::Vector<sp/>input(</highlight><highlight class="keyword">const_cast&lt;</highlight><highlight class="keywordtype">double</highlight><highlight class="normal">*</highlight><highlight class="keyword">&gt;</highlight><highlight class="normal">(volume_force.data()),<sp/></highlight><highlight class="keyword">static_cast&lt;</highlight><highlight class="keywordtype">int</highlight><highlight class="keyword">&gt;</highlight><highlight class="normal">(volume_force.size()));</highlight></codeline>
<codeline lineno="159"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1ac4b00cd302ad31d545f4e5ebde67c8e5" kindref="member">m_prolongation</ref>-&gt;MultTranspose(input,<sp/>workspace.independent);</highlight></codeline>
<codeline lineno="160"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a74d566ff3095fc95ad230e0fa18f2c84" kindref="member">m_restriction</ref>-&gt;MultTranspose(workspace.independent,<sp/>workspace.local);</highlight></codeline>
<codeline lineno="161"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a53f1ed40b9c1baeb98bea41e8712d78d" kindref="member">m_radial</ref>.apply_transpose({workspace.local.GetData(),<sp/></highlight><highlight class="keyword">static_cast&lt;</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">&gt;</highlight><highlight class="normal">(workspace.local.Size())},<sp/>surface_force);</highlight></codeline>
<codeline lineno="162"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="163"><highlight class="normal"></highlight></codeline>
<codeline lineno="164"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">inline</highlight><highlight class="normal"><sp/>std::expected&lt;ConstrainedRadialExtension,<sp/>DeformationFailure&gt;<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a27b36315f46f337ca2361833767b82c7" kindref="member">make_constrained_radial_extension</ref>(</highlight></codeline>
<codeline lineno="165"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::shared_ptr&lt;const<sp/>RadialInterpolation&gt;<sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a6151f50ac40bd03b6d4fd65fb376b68b" kindref="member">interpolation</ref>,<sp/>const<sp/>::mfem::FiniteElementSpace&amp;<sp/>displacement_space,</highlight></codeline>
<codeline lineno="166"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>PowerLawRadialExtension<sp/>interior,<sp/>FixedInfinityRadialExtension<sp/>vacuum)<sp/>{</highlight></codeline>
<codeline lineno="167"><highlight class="normal"></highlight><highlight class="preprocessor">#ifdef<sp/>MFEM_USE_MPI</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="168"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keyword">dynamic_cast&lt;</highlight><highlight class="keyword">const<sp/></highlight><highlight class="normal">::mfem::ParFiniteElementSpace*</highlight><highlight class="keyword">&gt;</highlight><highlight class="normal">(&amp;displacement_space)<sp/>!=<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">)<sp/>{</highlight></codeline>
<codeline lineno="169"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(DeformationFailure{.reason<sp/>=<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1abff36c0dd3d8b9c916a2e928ab69a871ad73416fa06a45d748a7c8d4349054cbb" kindref="member">DeformationFailureReason::incompatible_discretization</ref>,<sp/>.message<sp/>=<sp/></highlight><highlight class="stringliteral">&quot;Constrained<sp/>radial<sp/>extension<sp/>requires<sp/>a<sp/>serial<sp/>displacement<sp/>space.&quot;</highlight><highlight class="normal">});</highlight></codeline>
<codeline lineno="170"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="171"><highlight class="normal"></highlight><highlight class="preprocessor">#endif</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="172"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>radial<sp/>=<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a9744ea82df8c839a005c29235aecdb03" kindref="member">make_radial_extension</ref>(std::move(<ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1a6151f50ac40bd03b6d4fd65fb376b68b" kindref="member">interpolation</ref>),<sp/>interior,<sp/>vacuum);</highlight></codeline>
<codeline lineno="173"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!radial)<sp/>{</highlight></codeline>
<codeline lineno="174"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(radial.error());</highlight></codeline>
<codeline lineno="175"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="176"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(displacement_space.GetVDim()<sp/>!=<sp/>3<sp/>||<sp/>displacement_space.GetOrdering()<sp/>!=<sp/>::mfem::Ordering::byNODES<sp/>||</highlight></codeline>
<codeline lineno="177"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>displacement_space.GetVSize()<sp/>!=<sp/></highlight><highlight class="keyword">static_cast&lt;</highlight><highlight class="keywordtype">int</highlight><highlight class="keyword">&gt;</highlight><highlight class="normal">(radial-&gt;output_size()))<sp/>{</highlight></codeline>
<codeline lineno="178"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(DeformationFailure{.reason<sp/>=<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1abff36c0dd3d8b9c916a2e928ab69a871ad73416fa06a45d748a7c8d4349054cbb" kindref="member">DeformationFailureReason::incompatible_discretization</ref>,<sp/>.message<sp/>=<sp/></highlight><highlight class="stringliteral">&quot;Constrained<sp/>radial<sp/>extension<sp/>requires<sp/>its<sp/>three-component<sp/>reference<sp/>space<sp/>in<sp/>component-block<sp/>ordering.&quot;</highlight><highlight class="normal">});</highlight></codeline>
<codeline lineno="179"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="180"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>const<sp/>::mfem::SparseMatrix*<sp/>prolongation<sp/>=<sp/>displacement_space.GetConformingProlongation();</highlight></codeline>
<codeline lineno="181"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(prolongation<sp/>==<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">)<sp/>{</highlight></codeline>
<codeline lineno="182"><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="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1ae80591257dccff5948ea1bf361f5361d" kindref="member">ConstrainedRadialExtension</ref>(std::move(*radial),<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">,<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">);</highlight></codeline>
<codeline lineno="183"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="184"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>const<sp/>::mfem::SparseMatrix*<sp/>restriction<sp/>=<sp/>displacement_space.GetConformingRestriction();</highlight></codeline>
<codeline lineno="185"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(displacement_space.IsVariableOrder()<sp/>||<sp/>restriction<sp/>==<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal"><sp/>||<sp/>prolongation-&gt;Height()<sp/>!=<sp/>displacement_space.GetVSize()<sp/>||</highlight></codeline>
<codeline lineno="186"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>prolongation-&gt;Width()<sp/>!=<sp/>displacement_space.GetTrueVSize()<sp/>||<sp/>restriction-&gt;Height()<sp/>!=<sp/>prolongation-&gt;Width()<sp/>||</highlight></codeline>
<codeline lineno="187"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>restriction-&gt;Width()<sp/>!=<sp/>prolongation-&gt;Height())<sp/>{</highlight></codeline>
<codeline lineno="188"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(DeformationFailure{.reason<sp/>=<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1abff36c0dd3d8b9c916a2e928ab69a871ad73416fa06a45d748a7c8d4349054cbb" kindref="member">DeformationFailureReason::incompatible_discretization</ref>,<sp/>.message<sp/>=<sp/></highlight><highlight class="stringliteral">&quot;Constrained<sp/>radial<sp/>extension<sp/>requires<sp/>fixed-order<sp/>serial<sp/>prolongation<sp/>and<sp/>matching<sp/>coefficient<sp/>restriction<sp/>matrices.&quot;</highlight><highlight class="normal">});</highlight></codeline>
<codeline lineno="189"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="190"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="classserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_constrained_radial_extension_1ae80591257dccff5948ea1bf361f5361d" kindref="member">ConstrainedRadialExtension</ref>(std::move(*radial),<sp/>std::make_unique&lt;const<sp/>::mfem::SparseMatrix&gt;(*prolongation),</highlight></codeline>
<codeline lineno="191"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::make_unique&lt;const<sp/>::mfem::SparseMatrix&gt;(*restriction));</highlight></codeline>
<codeline lineno="192"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="193"><highlight class="normal">}<sp/></highlight><highlight class="comment">//<sp/>namespace<sp/>serif::geometry::deformation::backend::mfem</highlight><highlight class="normal"></highlight></codeline>
    </programlisting>
    <location file="src/include/serif/backend/mfem/geometry/deformation/constrained_radial_extension.hpp"/>
  </compounddef>
</doxygen>