261 lines
41 KiB
XML
261 lines
41 KiB
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="harmonics_8cpp" kind="file" language="C++">
|
||
|
|
<compoundname>harmonics.cpp</compoundname>
|
||
|
|
<includes refid="harmonics_8hpp" local="yes">serif/geometry/deformation/detail/harmonics.hpp</includes>
|
||
|
|
<includes local="no">algorithm</includes>
|
||
|
|
<includes local="no">cmath</includes>
|
||
|
|
<includes local="no">cstdint</includes>
|
||
|
|
<includes local="no">limits</includes>
|
||
|
|
<includes local="no">numbers</includes>
|
||
|
|
<incdepgraph>
|
||
|
|
<node id="2">
|
||
|
|
<label>serif/geometry/deformation/detail/harmonics.hpp</label>
|
||
|
|
<link refid="harmonics_8hpp"/>
|
||
|
|
<childnode refid="3" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="4" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="5" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="6" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="7" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="8" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="9">
|
||
|
|
<label>serif/geometry/deformation/failure.hpp</label>
|
||
|
|
<link refid="geometry_2deformation_2failure_8hpp"/>
|
||
|
|
<childnode refid="5" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="10" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="11">
|
||
|
|
<label>serif/geometry/deformation/models.hpp</label>
|
||
|
|
<link refid="models_8hpp"/>
|
||
|
|
<childnode refid="3" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="12" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="13">
|
||
|
|
<label>serif/geometry/deformation/spaces.hpp</label>
|
||
|
|
<link refid="spaces_8hpp"/>
|
||
|
|
</node>
|
||
|
|
<node id="8">
|
||
|
|
<label>serif/geometry/deformation/surface_embedding.hpp</label>
|
||
|
|
<link refid="surface__embedding_8hpp"/>
|
||
|
|
<childnode refid="3" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="4" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="5" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="6" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="7" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="9" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="11" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="14" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="14">
|
||
|
|
<label>serif/geometry/deformation/workspace.hpp</label>
|
||
|
|
<link refid="geometry_2deformation_2workspace_8hpp"/>
|
||
|
|
</node>
|
||
|
|
<node id="1">
|
||
|
|
<label>src/lib/serif/geometry/deformation/harmonics.cpp</label>
|
||
|
|
<link refid="harmonics_8cpp"/>
|
||
|
|
<childnode refid="2" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="15" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="16" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="17" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="18" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="19" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="15">
|
||
|
|
<label>algorithm</label>
|
||
|
|
</node>
|
||
|
|
<node id="3">
|
||
|
|
<label>array</label>
|
||
|
|
</node>
|
||
|
|
<node id="16">
|
||
|
|
<label>cmath</label>
|
||
|
|
</node>
|
||
|
|
<node id="4">
|
||
|
|
<label>cstddef</label>
|
||
|
|
</node>
|
||
|
|
<node id="17">
|
||
|
|
<label>cstdint</label>
|
||
|
|
</node>
|
||
|
|
<node id="5">
|
||
|
|
<label>expected</label>
|
||
|
|
</node>
|
||
|
|
<node id="18">
|
||
|
|
<label>limits</label>
|
||
|
|
</node>
|
||
|
|
<node id="19">
|
||
|
|
<label>numbers</label>
|
||
|
|
</node>
|
||
|
|
<node id="12">
|
||
|
|
<label>optional</label>
|
||
|
|
</node>
|
||
|
|
<node id="6">
|
||
|
|
<label>span</label>
|
||
|
|
</node>
|
||
|
|
<node id="10">
|
||
|
|
<label>string</label>
|
||
|
|
</node>
|
||
|
|
<node id="7">
|
||
|
|
<label>vector</label>
|
||
|
|
</node>
|
||
|
|
</incdepgraph>
|
||
|
|
<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_1detail">serif::geometry::deformation::detail</innernamespace>
|
||
|
|
<briefdescription>
|
||
|
|
</briefdescription>
|
||
|
|
<detaileddescription>
|
||
|
|
</detaileddescription>
|
||
|
|
<programlisting>
|
||
|
|
<codeline lineno="1"><highlight class="preprocessor">#include<sp/>"<ref refid="harmonics_8hpp" kindref="compound">serif/geometry/deformation/detail/harmonics.hpp</ref>"</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/><algorithm></highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="4"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/><cmath></highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="5"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/><cstdint></highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="6"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/><limits></highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="7"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/><numbers></highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="8"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="9"><highlight class="normal"></highlight><highlight class="keyword">namespace<sp/></highlight><highlight class="normal"><ref refid="namespaceserif_1_1geometry_1_1deformation_1_1detail" kindref="compound">serif::geometry::deformation::detail</ref><sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="10"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">namespace<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
||
|
|
<codeline lineno="12"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>column_product(std::span<const<sp/>double><sp/>basis,<sp/>std::size_t<sp/>rows,<sp/>std::size_t<sp/>columns,<sp/>std::size_t<sp/>first,<sp/>std::size_t<sp/>second)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="13"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">long</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>result<sp/>=<sp/>0.0L;</highlight></codeline>
|
||
|
|
<codeline lineno="14"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(std::size_t<sp/>row<sp/>=<sp/>0;<sp/>row<sp/><<sp/>rows;<sp/>++row)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="15"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>result<sp/>+=<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">long</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(basis[row<sp/>*<sp/>columns<sp/>+<sp/>first])<sp/>*<sp/>basis[row<sp/>*<sp/>columns<sp/>+<sp/>second];</highlight></codeline>
|
||
|
|
<codeline lineno="16"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="17"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">double</highlight><highlight class="keyword">></highlight><highlight class="normal">(result);</highlight></codeline>
|
||
|
|
<codeline lineno="18"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="19"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="21"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">void</highlight><highlight class="normal"><sp/>evaluate_legendre(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>std::array<double,<sp/>3>&<sp/>direction,<sp/></highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>maximum_degree,<sp/></highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>maximum_order,<sp/>std::span<double><sp/>legendre)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="22"><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/>std::size_t<sp/>stride<sp/>=<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(maximum_order)<sp/>+<sp/>1;</highlight></codeline>
|
||
|
|
<codeline lineno="23"><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/>cylindrical_radius<sp/>=<sp/>std::hypot(direction[0],<sp/>direction[1]);</highlight></codeline>
|
||
|
|
<codeline lineno="25"><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="keyword">auto</highlight><highlight class="normal"><sp/>polynomial<sp/>=<sp/>[&](</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>degree,<sp/></highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>order)<sp/>-><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal">&<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="26"><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/>legendre[</highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(degree)<sp/>*<sp/>stride<sp/>+<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(order)];</highlight></codeline>
|
||
|
|
<codeline lineno="27"><highlight class="normal"><sp/><sp/><sp/><sp/><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/>polynomial(0,<sp/>0)<sp/>=<sp/>1.0<sp/>/<sp/>std::sqrt(4.0<sp/>*<sp/>std::numbers::pi);</highlight></codeline>
|
||
|
|
<codeline lineno="29"><highlight class="normal"><sp/><sp/><sp/><sp/><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/>order<sp/>=<sp/>0;<sp/>order<sp/><=<sp/>maximum_order;<sp/>++order)<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">if</highlight><highlight class="normal"><sp/>(order<sp/>><sp/>0)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="31"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>polynomial(order,<sp/>order)<sp/>=<sp/>-std::sqrt((2.0<sp/>*<sp/>order<sp/>+<sp/>1.0)<sp/>/<sp/>(2.0<sp/>*<sp/>order))<sp/>*<sp/>cylindrical_radius<sp/>*<sp/>polynomial(order<sp/>-<sp/>1,<sp/>order<sp/>-<sp/>1);</highlight></codeline>
|
||
|
|
<codeline lineno="32"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="33"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(order<sp/><<sp/>maximum_degree)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="34"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>polynomial(order<sp/>+<sp/>1,<sp/>order)<sp/>=<sp/>std::sqrt(2.0<sp/>*<sp/>order<sp/>+<sp/>3.0)<sp/>*<sp/>direction[2]<sp/>*<sp/>polynomial(order,<sp/>order);</highlight></codeline>
|
||
|
|
<codeline lineno="35"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="36"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><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/>degree<sp/>=<sp/>order<sp/>+<sp/>2;<sp/>degree<sp/><=<sp/>maximum_degree;<sp/>++degree)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="37"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><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/>degree_value<sp/>=<sp/>degree;</highlight></codeline>
|
||
|
|
<codeline lineno="38"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><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/>order_value<sp/>=<sp/>order;</highlight></codeline>
|
||
|
|
<codeline lineno="39"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><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/>leading<sp/>=<sp/>std::sqrt((4.0<sp/>*<sp/>degree_value<sp/>*<sp/>degree_value<sp/>-<sp/>1.0)<sp/>/<sp/>(degree_value<sp/>*<sp/>degree_value<sp/>-<sp/>order_value<sp/>*<sp/>order_value));</highlight></codeline>
|
||
|
|
<codeline lineno="40"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><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/>trailing<sp/>=<sp/>std::sqrt(((2.0<sp/>*<sp/>degree_value<sp/>+<sp/>1.0)<sp/>*<sp/>((degree_value<sp/>-<sp/>1.0)<sp/>*<sp/>(degree_value<sp/>-<sp/>1.0)<sp/>-<sp/>order_value<sp/>*<sp/>order_value))<sp/>/<sp/>((2.0<sp/>*<sp/>degree_value<sp/>-<sp/>3.0)<sp/>*<sp/>(degree_value<sp/>*<sp/>degree_value<sp/>-<sp/>order_value<sp/>*<sp/>order_value)));</highlight></codeline>
|
||
|
|
<codeline lineno="41"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>polynomial(degree,<sp/>order)<sp/>=<sp/>leading<sp/>*<sp/>direction[2]<sp/>*<sp/>polynomial(degree<sp/>-<sp/>1,<sp/>order)<sp/>-<sp/>trailing<sp/>*<sp/>polynomial(degree<sp/>-<sp/>2,<sp/>order);</highlight></codeline>
|
||
|
|
<codeline lineno="42"><highlight class="normal"><sp/><sp/><sp/><sp/><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/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="44"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="45"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="46"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="47"><highlight class="normal"><sp/><sp/><sp/><sp/>std::expected<std::vector<SphericalHarmonicMode>,<sp/>DeformationFailure><sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1_1detail_1a2b90b4a4ff8a4b85379079d55b2eb3b1" kindref="member">make_harmonic_modes</ref>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>SphericalHarmonicSurface&<sp/>model,<sp/>std::size_t<sp/>node_count)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="48"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>maximum_order<sp/>=<sp/>model.maximum_azimuthal_order.value_or(model.maximum_degree);</highlight></codeline>
|
||
|
|
<codeline lineno="49"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(model.maximum_degree<sp/><<sp/>0<sp/>||<sp/>maximum_order<sp/><<sp/>0<sp/>||<sp/>maximum_order<sp/>><sp/>model.maximum_degree)<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/>std::unexpected(DeformationFailure{.reason<sp/>=<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1abff36c0dd3d8b9c916a2e928ab69a871a9cc57114b319994e90322314a00521d4" kindref="member">DeformationFailureReason::invalid_configuration</ref>,<sp/>.message<sp/>=<sp/></highlight><highlight class="stringliteral">"The<sp/>maximum<sp/>harmonic<sp/>degree<sp/>must<sp/>be<sp/>nonnegative<sp/>and<sp/>the<sp/>azimuthal<sp/>cutoff<sp/>must<sp/>lie<sp/>between<sp/>zero<sp/>and<sp/>that<sp/>degree."</highlight><highlight class="normal">});</highlight></codeline>
|
||
|
|
<codeline lineno="51"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="52"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>std::int64_t<sp/>degree<sp/>=<sp/>model.maximum_degree;</highlight></codeline>
|
||
|
|
<codeline lineno="53"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>std::int64_t<sp/>order<sp/>=<sp/>maximum_order;</highlight></codeline>
|
||
|
|
<codeline lineno="54"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>std::int64_t<sp/>count<sp/>=<sp/>(order<sp/>+<sp/>1)<sp/>*<sp/>(order<sp/>+<sp/>1)<sp/>+<sp/>(degree<sp/>-<sp/>order)<sp/>*<sp/>(2<sp/>*<sp/>order<sp/>+<sp/>1);</highlight></codeline>
|
||
|
|
<codeline lineno="55"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(count<sp/>><sp/>std::numeric_limits<int>::max())<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="56"><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_1abff36c0dd3d8b9c916a2e928ab69a871a9cc57114b319994e90322314a00521d4" kindref="member">DeformationFailureReason::invalid_configuration</ref>,<sp/>.message<sp/>=<sp/></highlight><highlight class="stringliteral">"The<sp/>requested<sp/>harmonic<sp/>mode<sp/>count<sp/>exceeds<sp/>supported<sp/>vector<sp/>indexing."</highlight><highlight class="normal">});</highlight></codeline>
|
||
|
|
<codeline lineno="57"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="58"><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">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(count)<sp/>><sp/>node_count)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="59"><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_1abff36c0dd3d8b9c916a2e928ab69a871a9a8a2ae86e6fce17b9be87ec37249370" kindref="member">DeformationFailureReason::unresolved_surface</ref>,<sp/>.message<sp/>=<sp/></highlight><highlight class="stringliteral">"The<sp/>requested<sp/>harmonic<sp/>mode<sp/>count<sp/>exceeds<sp/>the<sp/>number<sp/>of<sp/>surface<sp/>nodes;<sp/>reduce<sp/>the<sp/>degree<sp/>or<sp/>refine<sp/>the<sp/>surface."</highlight><highlight class="normal">});</highlight></codeline>
|
||
|
|
<codeline lineno="60"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="61"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::vector<SphericalHarmonicMode><sp/>modes;</highlight></codeline>
|
||
|
|
<codeline lineno="62"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>modes.reserve(</highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(count));</highlight></codeline>
|
||
|
|
<codeline lineno="63"><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/>degree_index<sp/>=<sp/>0;<sp/>degree_index<sp/><=<sp/>model.maximum_degree;<sp/>++degree_index)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="64"><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">int</highlight><highlight class="normal"><sp/>cutoff<sp/>=<sp/>std::min(degree_index,<sp/>maximum_order);</highlight></codeline>
|
||
|
|
<codeline lineno="65"><highlight class="normal"><sp/><sp/><sp/><sp/><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/>order_index<sp/>=<sp/>-cutoff;<sp/>order_index<sp/><=<sp/>cutoff;<sp/>++order_index)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="66"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>modes.push_back(SphericalHarmonicMode{.degree<sp/>=<sp/>degree_index,<sp/>.order<sp/>=<sp/>order_index});</highlight></codeline>
|
||
|
|
<codeline lineno="67"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="68"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="69"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>modes;</highlight></codeline>
|
||
|
|
<codeline lineno="70"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="71"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="72"><highlight class="normal"><sp/><sp/><sp/><sp/>std::expected<std::vector<double>,<sp/>DeformationFailure><sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1_1detail_1a02c8c6fceb651f31e86585b8a794562a" kindref="member">sample_harmonic_basis</ref>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>SphericalHarmonicSurface&<sp/>model,<sp/>std::span<const<sp/>SphericalHarmonicMode><sp/>modes,<sp/>std::span<</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>std::array<double,<sp/>3>><sp/>directions)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="73"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(modes.empty()<sp/>||<sp/>directions.size()<sp/>><sp/>std::vector<double>{}.max_size()<sp/>/<sp/>modes.size())<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="74"><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">"The<sp/>requested<sp/>sampled<sp/>harmonic<sp/>basis<sp/>exceeds<sp/>supported<sp/>storage<sp/>indexing."</highlight><highlight class="normal">});</highlight></codeline>
|
||
|
|
<codeline lineno="75"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="76"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>maximum_order<sp/>=<sp/>model.maximum_azimuthal_order.value_or(model.maximum_degree);</highlight></codeline>
|
||
|
|
<codeline lineno="77"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>std::size_t<sp/>stride<sp/>=<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(maximum_order)<sp/>+<sp/>1;</highlight></codeline>
|
||
|
|
<codeline lineno="78"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::vector<double><sp/>legendre((</highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(model.maximum_degree)<sp/>+<sp/>1)<sp/>*<sp/>stride,<sp/>0.0);</highlight></codeline>
|
||
|
|
<codeline lineno="79"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::vector<double><sp/>basis(directions.size()<sp/>*<sp/>modes.size());</highlight></codeline>
|
||
|
|
<codeline lineno="80"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(std::size_t<sp/>node<sp/>=<sp/>0;<sp/>node<sp/><<sp/>directions.size();<sp/>++node)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="81"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>evaluate_legendre(directions[node],<sp/>model.maximum_degree,<sp/>maximum_order,<sp/>legendre);</highlight></codeline>
|
||
|
|
<codeline lineno="82"><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/>azimuth<sp/>=<sp/>std::atan2(directions[node][1],<sp/>directions[node][0]);</highlight></codeline>
|
||
|
|
<codeline lineno="83"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(std::size_t<sp/>column<sp/>=<sp/>0;<sp/>column<sp/><<sp/>modes.size();<sp/>++column)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="84"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>SphericalHarmonicMode<sp/>mode<sp/>=<sp/>modes[column];</highlight></codeline>
|
||
|
|
<codeline lineno="85"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1discretization_1_1backend_1_1mfem_1af75f98db90618153ea4ca169945733f7a2063c1608d6e0baf80249c42e2be5804" kindref="member">value</ref><sp/>=<sp/>legendre[</highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(mode.degree)<sp/>*<sp/>stride<sp/>+<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(std::abs(mode.order))];</highlight></codeline>
|
||
|
|
<codeline lineno="86"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(mode.order<sp/><<sp/>0)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="87"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1discretization_1_1backend_1_1mfem_1af75f98db90618153ea4ca169945733f7a2063c1608d6e0baf80249c42e2be5804" kindref="member">value</ref><sp/>*=<sp/>std::sqrt(2.0)<sp/>*<sp/>std::sin(-mode.order<sp/>*<sp/>azimuth);</highlight></codeline>
|
||
|
|
<codeline lineno="88"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="89"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(mode.order<sp/>><sp/>0)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="90"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1discretization_1_1backend_1_1mfem_1af75f98db90618153ea4ca169945733f7a2063c1608d6e0baf80249c42e2be5804" kindref="member">value</ref><sp/>*=<sp/>std::sqrt(2.0)<sp/>*<sp/>std::cos(mode.order<sp/>*<sp/>azimuth);</highlight></codeline>
|
||
|
|
<codeline lineno="91"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="92"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!std::isfinite(value))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="93"><highlight class="normal"><sp/><sp/><sp/><sp/><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/>std::unexpected(DeformationFailure{.reason<sp/>=<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1abff36c0dd3d8b9c916a2e928ab69a871a4013d0a2d04acba5257678a15922a094" kindref="member">DeformationFailureReason::non_finite</ref>,<sp/>.message<sp/>=<sp/></highlight><highlight class="stringliteral">"Sampling<sp/>the<sp/>spherical<sp/>harmonic<sp/>basis<sp/>produced<sp/>a<sp/>non-finite<sp/>value."</highlight><highlight class="normal">});</highlight></codeline>
|
||
|
|
<codeline lineno="94"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="95"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>basis[node<sp/>*<sp/>modes.size()<sp/>+<sp/>column]<sp/>=<sp/><ref refid="namespaceserif_1_1discretization_1_1backend_1_1mfem_1af75f98db90618153ea4ca169945733f7a2063c1608d6e0baf80249c42e2be5804" kindref="member">value</ref>;</highlight></codeline>
|
||
|
|
<codeline lineno="96"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="97"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="98"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>basis;</highlight></codeline>
|
||
|
|
<codeline lineno="99"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="100"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="101"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1a08306ee593c172f4877facc061e98f97" kindref="member">DeformationStatus</ref><sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1_1detail_1a7628b540a8b43925a269627be0e1f8a7" kindref="member">check_independent_harmonics</ref>(std::span<const<sp/>double><sp/>basis,<sp/>std::size_t<sp/>node_count,<sp/>std::size_t<sp/>mode_count)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="102"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::vector<double><sp/>orthogonal(basis.begin(),<sp/>basis.end());</highlight></codeline>
|
||
|
|
<codeline lineno="103"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">constexpr</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>minimum_relative_pivot<sp/>=<sp/>1e-12;</highlight></codeline>
|
||
|
|
<codeline lineno="104"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(std::size_t<sp/>column<sp/>=<sp/>0;<sp/>column<sp/><<sp/>mode_count;<sp/>++column)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="105"><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/>original<sp/>=<sp/>std::sqrt(column_product(basis,<sp/>node_count,<sp/>mode_count,<sp/>column,<sp/>column));</highlight></codeline>
|
||
|
|
<codeline lineno="106"><highlight class="normal"><sp/><sp/><sp/><sp/><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/>pass<sp/>=<sp/>0;<sp/>pass<sp/><<sp/>2;<sp/>++pass)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="107"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(std::size_t<sp/>previous<sp/>=<sp/>0;<sp/>previous<sp/><<sp/>column;<sp/>++previous)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="108"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><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/><ref refid="namespaceserif_1_1discretization_1_1quadrature_1a2170ba7c65dd986e68afb2cc85f2ad88adfa555786d19f339a13e6719241e9d6f" kindref="member">projection</ref><sp/>=<sp/>column_product(orthogonal,<sp/>node_count,<sp/>mode_count,<sp/>previous,<sp/>column);</highlight></codeline>
|
||
|
|
<codeline lineno="109"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(std::size_t<sp/>row<sp/>=<sp/>0;<sp/>row<sp/><<sp/>node_count;<sp/>++row)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="110"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>orthogonal[row<sp/>*<sp/>mode_count<sp/>+<sp/>column]<sp/>-=<sp/><ref refid="namespaceserif_1_1discretization_1_1quadrature_1a2170ba7c65dd986e68afb2cc85f2ad88adfa555786d19f339a13e6719241e9d6f" kindref="member">projection</ref><sp/>*<sp/>orthogonal[row<sp/>*<sp/>mode_count<sp/>+<sp/>previous];</highlight></codeline>
|
||
|
|
<codeline lineno="111"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="112"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="113"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="114"><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/>pivot<sp/>=<sp/>std::sqrt(column_product(orthogonal,<sp/>node_count,<sp/>mode_count,<sp/>column,<sp/>column));</highlight></codeline>
|
||
|
|
<codeline lineno="115"><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(original)<sp/>||<sp/>!(original<sp/>><sp/>0.0)<sp/>||<sp/>!std::isfinite(pivot)<sp/>||<sp/>!(pivot<sp/>><sp/>minimum_relative_pivot<sp/>*<sp/>original))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="116"><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/>std::unexpected(DeformationFailure{.reason<sp/>=<sp/><ref refid="namespaceserif_1_1geometry_1_1deformation_1abff36c0dd3d8b9c916a2e928ab69a871a9a8a2ae86e6fce17b9be87ec37249370" kindref="member">DeformationFailureReason::unresolved_surface</ref>,<sp/>.message<sp/>=<sp/></highlight><highlight class="stringliteral">"The<sp/>selected<sp/>harmonics<sp/>are<sp/>numerically<sp/>dependent<sp/>on<sp/>the<sp/>surface<sp/>nodes;<sp/>reduce<sp/>their<sp/>degree<sp/>or<sp/>order,<sp/>or<sp/>refine<sp/>the<sp/>surface."</highlight><highlight class="normal">});</highlight></codeline>
|
||
|
|
<codeline lineno="117"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="118"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(std::size_t<sp/>row<sp/>=<sp/>0;<sp/>row<sp/><<sp/>node_count;<sp/>++row)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="119"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>orthogonal[row<sp/>*<sp/>mode_count<sp/>+<sp/>column]<sp/>/=<sp/>pivot;</highlight></codeline>
|
||
|
|
<codeline lineno="120"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="121"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="122"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>{};</highlight></codeline>
|
||
|
|
<codeline lineno="123"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="124"><highlight class="normal">}</highlight></codeline>
|
||
|
|
</programlisting>
|
||
|
|
<location file="src/lib/serif/geometry/deformation/harmonics.cpp"/>
|
||
|
|
</compounddef>
|
||
|
|
</doxygen>
|