1096 lines
172 KiB
XML
1096 lines
172 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="domain__mapper_8cpp" kind="file" language="C++">
|
||
|
|
<compoundname>domain_mapper.cpp</compoundname>
|
||
|
|
<includes local="yes">mfem/linalg/densemat.hpp</includes>
|
||
|
|
<includes refid="backend_2mfem_2mapping_2domain__mapper_8hpp" local="yes">serif/backend/mfem/mapping/domain_mapper.hpp</includes>
|
||
|
|
<includes refid="errors_8hpp" local="yes">serif/utils/error/errors.hpp</includes>
|
||
|
|
<includes refid="backend_2mfem_2shared_2finite_8hpp" local="yes">serif/backend/mfem/shared/finite.hpp</includes>
|
||
|
|
<includes refid="invertible_8hpp" local="yes">serif/backend/mfem/shared/invertible.hpp</includes>
|
||
|
|
<incdepgraph>
|
||
|
|
<node id="19">
|
||
|
|
<label>serif/backend/mfem/mapping/compactification/compactification.hpp</label>
|
||
|
|
<link refid="compactification_8hpp"/>
|
||
|
|
<childnode refid="10" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="5" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="11" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="20" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="21" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="20">
|
||
|
|
<label>serif/backend/mfem/mapping/compactification/input_result.hpp</label>
|
||
|
|
<link refid="input__result_8hpp"/>
|
||
|
|
<childnode refid="6" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="2" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="9">
|
||
|
|
<label>serif/backend/mfem/mapping/compactification/kelvin.hpp</label>
|
||
|
|
<link refid="backend_2mfem_2mapping_2compactification_2kelvin_8hpp"/>
|
||
|
|
<childnode refid="10" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="6" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="11" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="19" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="22" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="21">
|
||
|
|
<label>serif/backend/mfem/mapping/compactification/variation.hpp</label>
|
||
|
|
<link refid="compactification_2variation_8hpp"/>
|
||
|
|
<childnode refid="6" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="2" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="23">
|
||
|
|
<label>serif/backend/mfem/mapping/context.hpp</label>
|
||
|
|
<link refid="context_8hpp"/>
|
||
|
|
<childnode refid="2" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="6" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="24">
|
||
|
|
<label>serif/backend/mfem/mapping/data.hpp</label>
|
||
|
|
<link refid="data_8hpp"/>
|
||
|
|
<childnode refid="25" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="26" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="6" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="2" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="27" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="9" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="23" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="24" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="3">
|
||
|
|
<label>serif/backend/mfem/mapping/domain_mapper.hpp</label>
|
||
|
|
<link refid="backend_2mfem_2mapping_2domain__mapper_8hpp"/>
|
||
|
|
<childnode refid="4" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="5" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="6" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="2" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="7" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="8" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="9" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="23" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="24" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="28" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="29" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="11">
|
||
|
|
<label>serif/backend/mfem/mapping/status.hpp</label>
|
||
|
|
<link refid="status_8hpp"/>
|
||
|
|
<childnode refid="12" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="14" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="15" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="16" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="28">
|
||
|
|
<label>serif/backend/mfem/mapping/variation.hpp</label>
|
||
|
|
<link refid="variation_8hpp"/>
|
||
|
|
<childnode refid="2" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="6" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="45">
|
||
|
|
<label>serif/backend/mfem/shared/finite.hpp</label>
|
||
|
|
<link refid="backend_2mfem_2shared_2finite_8hpp"/>
|
||
|
|
<childnode refid="46" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="6" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="2" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="48">
|
||
|
|
<label>serif/backend/mfem/shared/invertible.hpp</label>
|
||
|
|
<link refid="invertible_8hpp"/>
|
||
|
|
<childnode refid="45" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="16">
|
||
|
|
<label>serif/meta/enums.hpp</label>
|
||
|
|
<link refid="enums_8hpp"/>
|
||
|
|
<childnode refid="12" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="17" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="18" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="46">
|
||
|
|
<label>serif/utils/checks/finite.hpp</label>
|
||
|
|
<link refid="utils_2checks_2finite_8hpp"/>
|
||
|
|
<childnode refid="47" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="41" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="15">
|
||
|
|
<label>serif/utils/error/codes.hpp</label>
|
||
|
|
<link refid="codes_8hpp"/>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="16" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="30">
|
||
|
|
<label>serif/utils/error/errors.hpp</label>
|
||
|
|
<link refid="errors_8hpp"/>
|
||
|
|
<childnode refid="10" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="31" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="32" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="33" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="34" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="35" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="36" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="39" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="15" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="43" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="39">
|
||
|
|
<label>serif/utils/misc/concepts/exceptions.hpp</label>
|
||
|
|
<link refid="utils_2misc_2concepts_2exceptions_8hpp"/>
|
||
|
|
<childnode refid="40" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="17" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="41" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="41">
|
||
|
|
<label>serif/utils/misc/concepts/numeric.hpp</label>
|
||
|
|
<link refid="numeric_8hpp"/>
|
||
|
|
<childnode refid="42" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="17" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="36">
|
||
|
|
<label>serif/utils/misc/terminal/colors.hpp</label>
|
||
|
|
<link refid="colors_8hpp"/>
|
||
|
|
<childnode refid="10" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="32" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="37" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="14" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="38" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="43">
|
||
|
|
<label>serif/utils/types/strings.hpp</label>
|
||
|
|
<link refid="types_2strings_8hpp"/>
|
||
|
|
<childnode refid="34" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="44" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="10" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="22">
|
||
|
|
<label>serif/utils/user_options/mapping/compactification/kelvin.hpp</label>
|
||
|
|
<link refid="utils_2user__options_2mapping_2compactification_2kelvin_8hpp"/>
|
||
|
|
</node>
|
||
|
|
<node id="29">
|
||
|
|
<label>serif/utils/user_options/mapping/domain_mapper.hpp</label>
|
||
|
|
<link refid="utils_2user__options_2mapping_2domain__mapper_8hpp"/>
|
||
|
|
</node>
|
||
|
|
<node id="1">
|
||
|
|
<label>src/lib/serif/backend/mfem/mapping/domain_mapper.cpp</label>
|
||
|
|
<link refid="domain__mapper_8cpp"/>
|
||
|
|
<childnode refid="2" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="3" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="30" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="45" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="48" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="34">
|
||
|
|
<label>algorithm</label>
|
||
|
|
</node>
|
||
|
|
<node id="47">
|
||
|
|
<label>cmath</label>
|
||
|
|
</node>
|
||
|
|
<node id="42">
|
||
|
|
<label>concepts</label>
|
||
|
|
</node>
|
||
|
|
<node id="13">
|
||
|
|
<label>cstdint</label>
|
||
|
|
</node>
|
||
|
|
<node id="32">
|
||
|
|
<label>cstdio</label>
|
||
|
|
</node>
|
||
|
|
<node id="37">
|
||
|
|
<label>cstdlib</label>
|
||
|
|
</node>
|
||
|
|
<node id="40">
|
||
|
|
<label>exception</label>
|
||
|
|
</node>
|
||
|
|
<node id="4">
|
||
|
|
<label>expected</label>
|
||
|
|
</node>
|
||
|
|
<node id="31">
|
||
|
|
<label>format</label>
|
||
|
|
</node>
|
||
|
|
<node id="44">
|
||
|
|
<label>iostream</label>
|
||
|
|
</node>
|
||
|
|
<node id="18">
|
||
|
|
<label>magic_enum/magic_enum_all.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="5">
|
||
|
|
<label>memory</label>
|
||
|
|
</node>
|
||
|
|
<node id="7">
|
||
|
|
<label>mfem/fem/eltrans.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="25">
|
||
|
|
<label>mfem/fem/fe/fe_base.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="26">
|
||
|
|
<label>mfem/fem/fespace.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="8">
|
||
|
|
<label>mfem/fem/intrules.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="2">
|
||
|
|
<label>mfem/linalg/densemat.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="6">
|
||
|
|
<label>mfem/linalg/vector.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="27">
|
||
|
|
<label>mfem/mesh/mesh.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="33">
|
||
|
|
<label>string</label>
|
||
|
|
</node>
|
||
|
|
<node id="10">
|
||
|
|
<label>string_view</label>
|
||
|
|
</node>
|
||
|
|
<node id="17">
|
||
|
|
<label>type_traits</label>
|
||
|
|
</node>
|
||
|
|
<node id="38">
|
||
|
|
<label>unistd.h</label>
|
||
|
|
</node>
|
||
|
|
<node id="12">
|
||
|
|
<label>unordered_map</label>
|
||
|
|
</node>
|
||
|
|
<node id="14">
|
||
|
|
<label>utility</label>
|
||
|
|
</node>
|
||
|
|
<node id="35">
|
||
|
|
<label>version</label>
|
||
|
|
</node>
|
||
|
|
</incdepgraph>
|
||
|
|
<innernamespace refid="namespaceserif">serif</innernamespace>
|
||
|
|
<innernamespace refid="namespaceserif_1_1mapping">serif::mapping</innernamespace>
|
||
|
|
<briefdescription>
|
||
|
|
</briefdescription>
|
||
|
|
<detaileddescription>
|
||
|
|
</detaileddescription>
|
||
|
|
<programlisting>
|
||
|
|
<codeline lineno="1"><highlight class="preprocessor">#include<sp/>"mfem/linalg/densemat.hpp"</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/>"<ref refid="backend_2mfem_2mapping_2domain__mapper_8hpp" kindref="compound">serif/backend/mfem/mapping/domain_mapper.hpp</ref>"</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="4"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>"<ref refid="errors_8hpp" kindref="compound">serif/utils/error/errors.hpp</ref>"</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="5"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="6"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>"<ref refid="backend_2mfem_2shared_2finite_8hpp" kindref="compound">serif/backend/mfem/shared/finite.hpp</ref>"</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="7"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>"<ref refid="invertible_8hpp" kindref="compound">serif/backend/mfem/shared/invertible.hpp</ref>"</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="8"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="9"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="10"><highlight class="normal"></highlight><highlight class="comment">//<sp/>Anonymous<sp/>namespace<sp/>for<sp/>translation<sp/>unit-local<sp/>helper<sp/>functions<sp/>and<sp/>types</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="11"><highlight class="normal"></highlight><highlight class="keyword">namespace<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
||
|
|
<codeline lineno="12"><highlight class="normal"><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>TODO:<sp/>Put<sp/>something<sp/>in<sp/>here<sp/>using<sp/>variadic<sp/>args<sp/>to<sp/>validate<sp/>multiple<sp/>things<sp/>satisfy<sp/>some<sp/>predicate</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="13"><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_1mapping" kindref="compound">serif::mapping</ref><sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="16"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/>namespace<sp/></highlight><highlight class="normal">serif::utils::errors;</highlight></codeline>
|
||
|
|
<codeline lineno="17"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/>namespace<sp/></highlight><highlight class="normal">serif::utils::checks;</highlight></codeline>
|
||
|
|
<codeline lineno="18"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="19"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="classserif_1_1mapping_1_1_domain_mapper_1abe3f3d0dafe5425f6a8d2d74b5c20176" kindref="member">DomainMapper::DomainMapper</ref>(</highlight></codeline>
|
||
|
|
<codeline lineno="20"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>utils::user_options::DomainMapperOptions<sp/>options,</highlight></codeline>
|
||
|
|
<codeline lineno="21"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::unique_ptr<const<sp/>compactification::ExteriorDomainMap><sp/>exterior_map</highlight></codeline>
|
||
|
|
<codeline lineno="22"><highlight class="normal"><sp/><sp/><sp/><sp/>)<sp/>:</highlight></codeline>
|
||
|
|
<codeline lineno="23"><highlight class="normal"><sp/><sp/><sp/><sp/>m_options(options),</highlight></codeline>
|
||
|
|
<codeline lineno="24"><highlight class="normal"><sp/><sp/><sp/><sp/>m_exterior_map(std::move(exterior_map))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="25"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(m_options.dimension<sp/><=<sp/>0)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="26"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>serif_error<std::invalid_argument,<sp/>SERiFErrorCode::INVALID_DIMENSION>(</highlight><highlight class="stringliteral">"Domain<sp/>Mapper<sp/>requires<sp/>a<sp/>positive,<sp/>non-zero,<sp/>dimension.<sp/>Received<sp/>dimension:<sp/>{}"</highlight><highlight class="normal">,<sp/>m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="27"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="28"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="29"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>if<sp/>(m_options.vacuum_element_attribute_ID<sp/><=<sp/>0)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="30"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>serif_error<std::invalid_argument,<sp/>SERiFErrorCode::INVALID_VACUUM_ELEMENT_ATTRIBUTE_ID>(</highlight><highlight class="stringliteral">"Domain<sp/>Mapper<sp/>requires<sp/>a<sp/>positive,<sp/>non-zero,<sp/>vacuum<sp/>element<sp/>attribute<sp/>ID.<sp/>Received<sp/>attribute<sp/>ID:<sp/>{}"</highlight><highlight class="normal">,<sp/>m_options.vacuum_element_attribute_ID);</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="33"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>if<sp/>(!m_exterior_map)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="34"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a2447eb80947534cdf0fbd205251cfd9c" kindref="member">serif_error<std::invalid_argument, SERiFErrorCode::MISSING_EXTERIOR_MAP></ref>(</highlight><highlight class="stringliteral">"Domain<sp/>Mapper<sp/>requires<sp/>a<sp/>valid<sp/>exterior<sp/>map.<sp/>One<sp/>was<sp/>not<sp/>provided<sp/>or<sp/>a<sp/>nullptr<sp/>was<sp/>provided."</highlight><highlight class="normal">);</highlight></codeline>
|
||
|
|
<codeline lineno="35"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="36"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="37"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="38"><highlight class="normal"><sp/><sp/><sp/><sp/>std::expected<DMFixedMeshGeometry,<sp/>MappingStatus><sp/><ref refid="namespaceserif_1_1mapping_1a8935271a2117c2232be8306651d5ce9f" kindref="member">sample_fixed_mesh_geometry</ref>(</highlight></codeline>
|
||
|
|
<codeline lineno="39"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::ElementTransformation&<sp/>transformation,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint&<sp/>point</highlight></codeline>
|
||
|
|
<codeline lineno="40"><highlight class="normal"><sp/><sp/><sp/><sp/>)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="41"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>transformation.SetIntPoint(&point);</highlight></codeline>
|
||
|
|
<codeline lineno="42"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>DMFixedMeshGeometry<sp/><ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>;</highlight></codeline>
|
||
|
|
<codeline lineno="43"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>transformation.Transform(point,<sp/><ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.position);</highlight></codeline>
|
||
|
|
<codeline lineno="44"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.jacobian<sp/>=<sp/>transformation.Jacobian();</highlight></codeline>
|
||
|
|
<codeline lineno="45"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.weight<sp/>=<sp/>transformation.Weight();</highlight></codeline>
|
||
|
|
<codeline lineno="46"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Check<sp/>before<sp/>MFEM<sp/>attempts<sp/>an<sp/>inverse:<sp/>singular<sp/>reference<sp/>elements<sp/>are<sp/>numerical<sp/>failures.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="47"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(<ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.position)<sp/>||<sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(<ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.jacobian)<sp/>||<sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(<ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.weight))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="48"><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(MappingStatus::NON_FINITE_MAPPING_INPUT);</highlight></codeline>
|
||
|
|
<codeline lineno="49"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="50"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(<ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.weight<sp/><=<sp/>0.0<sp/>||</highlight></codeline>
|
||
|
|
<codeline lineno="51"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>(<ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.jacobian.Height()<sp/>==<sp/><ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.jacobian.Width()<sp/>&&<sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a3471cc276b4a977cff023828197af29e" kindref="member">has_stable_inverse</ref>(<ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.jacobian)))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="52"><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(MappingStatus::NON_POSITIVE_DETERMINATE);</highlight></codeline>
|
||
|
|
<codeline lineno="53"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="54"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.inverse_jacobian<sp/>=<sp/>transformation.InverseJacobian();</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/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(<ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.inverse_jacobian))<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(MappingStatus::NON_FINITE_MAPPING_INPUT);</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">return</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>;</highlight></codeline>
|
||
|
|
<codeline lineno="59"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="60"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="61"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>DomainMapper::evaluate_point(</highlight></codeline>
|
||
|
|
<codeline lineno="62"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementMappingData<sp/>&element_data,</highlight></codeline>
|
||
|
|
<codeline lineno="63"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::ElementTransformation<sp/>&transformation,</highlight></codeline>
|
||
|
|
<codeline lineno="64"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint<sp/>&ip,</highlight></codeline>
|
||
|
|
<codeline lineno="65"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>DMWorkspace<sp/>&workspace,</highlight></codeline>
|
||
|
|
<codeline lineno="66"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>MappingPointContext<sp/>&context</highlight></codeline>
|
||
|
|
<codeline lineno="67"><highlight class="normal"><sp/><sp/><sp/><sp/>)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
||
|
|
<codeline lineno="68"><highlight class="normal"><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/><ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref><sp/>=<sp/><ref refid="namespaceserif_1_1mapping_1a8935271a2117c2232be8306651d5ce9f" kindref="member">sample_fixed_mesh_geometry</ref>(transformation,<sp/>ip);</highlight></codeline>
|
||
|
|
<codeline lineno="69"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="70"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.error();</highlight></codeline>
|
||
|
|
<codeline lineno="71"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="72"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>evaluate_point_impl(element_data,<sp/>transformation,<sp/>ip,<sp/>workspace,<sp/>context,<sp/>*<ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>);</highlight></codeline>
|
||
|
|
<codeline lineno="73"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="74"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="75"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>DomainMapper::evaluate_point(</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/>ElementMappingData&<sp/>element_data,</highlight></codeline>
|
||
|
|
<codeline lineno="77"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::ElementTransformation&<sp/>transformation,</highlight></codeline>
|
||
|
|
<codeline lineno="78"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint&<sp/>ip,</highlight></codeline>
|
||
|
|
<codeline lineno="79"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>DMWorkspace&<sp/>workspace,</highlight></codeline>
|
||
|
|
<codeline lineno="80"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>MappingPointContext&<sp/>context,</highlight></codeline>
|
||
|
|
<codeline lineno="81"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>DMFixedMeshGeometry&<sp/>fixed_geometry</highlight></codeline>
|
||
|
|
<codeline lineno="82"><highlight class="normal"><sp/><sp/><sp/><sp/>)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
||
|
|
<codeline lineno="83"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(fixed_geometry.position.Size()<sp/>!=<sp/>m_options.dimension<sp/>||</highlight></codeline>
|
||
|
|
<codeline lineno="84"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>fixed_geometry.jacobian.Height()<sp/>!=<sp/>m_options.dimension<sp/>||<sp/>fixed_geometry.jacobian.Width()<sp/>!=<sp/>m_options.dimension<sp/>||</highlight></codeline>
|
||
|
|
<codeline lineno="85"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>fixed_geometry.inverse_jacobian.Height()<sp/>!=<sp/>m_options.dimension<sp/>||<sp/>fixed_geometry.inverse_jacobian.Width()<sp/>!=<sp/>m_options.dimension)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="86"><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/>MappingStatus::INVALID_DIMENSION;</highlight></codeline>
|
||
|
|
<codeline lineno="87"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="88"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(fixed_geometry.position)<sp/>||<sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(fixed_geometry.jacobian)<sp/>||</highlight></codeline>
|
||
|
|
<codeline lineno="89"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(fixed_geometry.inverse_jacobian)<sp/>||<sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(fixed_geometry.weight))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="90"><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/>MappingStatus::NON_FINITE_MAPPING_INPUT;</highlight></codeline>
|
||
|
|
<codeline lineno="91"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="92"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(fixed_geometry.weight<sp/><=<sp/>0.0<sp/>||<sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a3471cc276b4a977cff023828197af29e" kindref="member">has_stable_inverse</ref>(fixed_geometry.jacobian))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="93"><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/>MappingStatus::NON_POSITIVE_DETERMINATE;</highlight></codeline>
|
||
|
|
<codeline lineno="94"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="95"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>transformation.SetIntPoint(&ip);</highlight></codeline>
|
||
|
|
<codeline lineno="96"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>evaluate_point_impl(element_data,<sp/>transformation,<sp/>ip,<sp/>workspace,<sp/>context,<sp/>fixed_geometry);</highlight></codeline>
|
||
|
|
<codeline lineno="97"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="98"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="99"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>DomainMapper::evaluate_volume(</highlight></codeline>
|
||
|
|
<codeline lineno="100"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementMappingData<sp/>&element_data,</highlight></codeline>
|
||
|
|
<codeline lineno="101"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::ElementTransformation<sp/>&transformation,</highlight></codeline>
|
||
|
|
<codeline lineno="102"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint<sp/>&ip,</highlight></codeline>
|
||
|
|
<codeline lineno="103"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>DMWorkspace<sp/>&workspace,</highlight></codeline>
|
||
|
|
<codeline lineno="104"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>VolumeMappingContext<sp/>&context,</highlight></codeline>
|
||
|
|
<codeline lineno="105"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>DMReferenceDisplacementSample<sp/>*displacement_sample,</highlight></codeline>
|
||
|
|
<codeline lineno="106"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>DMFixedMeshGeometry*<sp/>fixed_geometry</highlight></codeline>
|
||
|
|
<codeline lineno="107"><highlight class="normal"><sp/><sp/><sp/><sp/>)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
||
|
|
<codeline lineno="108"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>transformation.SetIntPoint(&ip);</highlight></codeline>
|
||
|
|
<codeline lineno="109"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::expected<DMFixedMeshGeometry,<sp/>MappingStatus><sp/>sampled_geometry;</highlight></codeline>
|
||
|
|
<codeline lineno="110"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(fixed_geometry<sp/>==<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="111"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>sampled_geometry<sp/>=<sp/><ref refid="namespaceserif_1_1mapping_1a8935271a2117c2232be8306651d5ce9f" kindref="member">sample_fixed_mesh_geometry</ref>(transformation,<sp/>ip);</highlight></codeline>
|
||
|
|
<codeline lineno="112"><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/>(!sampled_geometry)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="113"><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/>sampled_geometry.error();</highlight></codeline>
|
||
|
|
<codeline lineno="114"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="115"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>fixed_geometry<sp/>=<sp/>&*sampled_geometry;</highlight></codeline>
|
||
|
|
<codeline lineno="116"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="117"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>point_status<sp/>=<sp/>evaluate_point_impl(</highlight></codeline>
|
||
|
|
<codeline lineno="118"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>element_data,</highlight></codeline>
|
||
|
|
<codeline lineno="119"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>transformation,</highlight></codeline>
|
||
|
|
<codeline lineno="120"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>ip,</highlight></codeline>
|
||
|
|
<codeline lineno="121"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace,</highlight></codeline>
|
||
|
|
<codeline lineno="122"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.mapping,</highlight></codeline>
|
||
|
|
<codeline lineno="123"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>*fixed_geometry,</highlight></codeline>
|
||
|
|
<codeline lineno="124"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>displacement_sample</highlight></codeline>
|
||
|
|
<codeline lineno="125"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>);</highlight></codeline>
|
||
|
|
<codeline lineno="126"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="127"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(point_status<sp/>!=<sp/>MappingStatus::VALID)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="128"><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/>point_status;</highlight></codeline>
|
||
|
|
<codeline lineno="129"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="130"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="131"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::Mult(context.mapping.mapping_jacobian,fixed_geometry->jacobian,<sp/>workspace.m_full_element_jacobian);</highlight></codeline>
|
||
|
|
<codeline lineno="132"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="133"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.quadrature.J_element_to_physical_inverse.SetSize(m_options.dimension,<sp/>m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="134"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::CalcInverse(workspace.m_full_element_jacobian,<sp/>context.quadrature.J_element_to_physical_inverse);</highlight></codeline>
|
||
|
|
<codeline lineno="135"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="136"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.quadrature.mapping_determinate<sp/>=<sp/>context.mapping.mapping_determinate;</highlight></codeline>
|
||
|
|
<codeline lineno="137"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.quadrature.physical_integration_weight<sp/>=<sp/>ip.weight<sp/>*<sp/>fixed_geometry->weight<sp/>*<sp/>context.mapping.mapping_determinate;</highlight></codeline>
|
||
|
|
<codeline lineno="138"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="139"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(context.quadrature.J_element_to_physical_inverse)<sp/>||<sp/>!std::isfinite(context.quadrature.physical_integration_weight))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="140"><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/>MappingStatus::NON_FINITE_MAPPING_RESULT;</highlight></codeline>
|
||
|
|
<codeline lineno="141"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="142"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(context.quadrature.physical_integration_weight<sp/><=<sp/>0.0)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="143"><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/>MappingStatus::NON_POSITIVE_DETERMINATE;</highlight></codeline>
|
||
|
|
<codeline lineno="144"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="145"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="146"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>MappingStatus::VALID;</highlight></codeline>
|
||
|
|
<codeline lineno="147"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="148"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="149"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>DomainMapper::evaluate_face(</highlight></codeline>
|
||
|
|
<codeline lineno="150"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementMappingData<sp/>&element_data,</highlight></codeline>
|
||
|
|
<codeline lineno="151"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::FaceElementTransformations<sp/>&transformation,</highlight></codeline>
|
||
|
|
<codeline lineno="152"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>FaceElementSide<sp/>side,</highlight></codeline>
|
||
|
|
<codeline lineno="153"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint<sp/>&ip,</highlight></codeline>
|
||
|
|
<codeline lineno="154"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>DMWorkspace<sp/>&workspace,</highlight></codeline>
|
||
|
|
<codeline lineno="155"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>FaceMappingContext<sp/>&context</highlight></codeline>
|
||
|
|
<codeline lineno="156"><highlight class="normal"><sp/><sp/><sp/><sp/>)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
||
|
|
<codeline lineno="157"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>transformation.SetAllIntPoints(&ip);</highlight></codeline>
|
||
|
|
<codeline lineno="158"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::ElementTransformation&<sp/><ref refid="namespaceserif_1_1utils_1_1parallel_1a8e0e78d02ab6c49ae588644e579bed79" kindref="member">element_transformation</ref><sp/>=<sp/>select_face_element_transformation(transformation,<sp/>side);</highlight></codeline>
|
||
|
|
<codeline lineno="159"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint&<sp/>element_integration_point<sp/>=<sp/>select_face_element_integration_point(transformation,<sp/>side);</highlight></codeline>
|
||
|
|
<codeline lineno="160"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="161"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>point_status<sp/>=<sp/>evaluate_point(element_data,<sp/>element_transformation,<sp/>element_integration_point,<sp/>workspace,<sp/>context.mapping);</highlight></codeline>
|
||
|
|
<codeline lineno="162"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(point_status<sp/>!=<sp/>MappingStatus::VALID)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="163"><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/>point_status;</highlight></codeline>
|
||
|
|
<codeline lineno="164"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="165"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="166"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Get<sp/>the<sp/>normal<sp/>vector<sp/>from<sp/>the<sp/>face<sp/>in<sp/>reference<sp/>space</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="167"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_reference_normal.SetSize(m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="168"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::CalcOrtho(transformation.Jacobian(),<sp/>workspace.m_reference_normal);</highlight></codeline>
|
||
|
|
<codeline lineno="169"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(side<sp/>==<sp/>FaceElementSide::ELEMENT_2)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="170"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_reference_normal<sp/>*=<sp/>-1.0;<sp/></highlight><highlight class="comment">//<sp/>Flip<sp/>the<sp/>normal<sp/>for<sp/>the<sp/>second<sp/>element<sp/>(e.g.<sp/>the<sp/>"other"<sp/>side<sp/>of<sp/>the<sp/>face)</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="171"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="172"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="173"><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">double</highlight><highlight class="normal"><sp/>reference_normal_magnitude<sp/>=<sp/>workspace.m_reference_normal.Norml2();</highlight></codeline>
|
||
|
|
<codeline lineno="174"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(reference_normal_magnitude)<sp/>||<sp/>reference_normal_magnitude<sp/><=<sp/>0.0)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="175"><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/>MappingStatus::NON_FINITE_MAPPING_RESULT;</highlight></codeline>
|
||
|
|
<codeline lineno="176"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="177"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.reference_normal.SetSize(m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="178"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.reference_normal<sp/>=<sp/>workspace.m_reference_normal;</highlight></codeline>
|
||
|
|
<codeline lineno="179"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.reference_normal<sp/>/=<sp/>reference_normal_magnitude;</highlight></codeline>
|
||
|
|
<codeline lineno="180"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="181"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Get<sp/>the<sp/>normal<sp/>vector<sp/>from<sp/>the<sp/>face<sp/>in<sp/>physical<sp/>space</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="182"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.mapping.inverse_mapping_jacobian.MultTranspose(workspace.m_reference_normal,<sp/>workspace.m_mapped_normal);</highlight></codeline>
|
||
|
|
<codeline lineno="183"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_mapped_normal<sp/>*=<sp/>context.mapping.mapping_determinate;</highlight></codeline>
|
||
|
|
<codeline lineno="184"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="185"><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">double</highlight><highlight class="normal"><sp/>mapped_normal_magnitude<sp/>=<sp/>workspace.m_mapped_normal.Norml2();</highlight></codeline>
|
||
|
|
<codeline lineno="186"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(mapped_normal_magnitude)<sp/>||<sp/>mapped_normal_magnitude<sp/><=<sp/>0.0)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="187"><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/>MappingStatus::NON_FINITE_MAPPING_RESULT;</highlight></codeline>
|
||
|
|
<codeline lineno="188"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="189"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.quadrature.normal.SetSize(m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="190"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.quadrature.normal<sp/>=<sp/>workspace.m_mapped_normal;</highlight></codeline>
|
||
|
|
<codeline lineno="191"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.quadrature.normal<sp/>/=<sp/>mapped_normal_magnitude;</highlight></codeline>
|
||
|
|
<codeline lineno="192"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="193"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.reference_surface_weight<sp/>=<sp/>ip.weight<sp/>*<sp/>reference_normal_magnitude;</highlight></codeline>
|
||
|
|
<codeline lineno="194"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.physical_surface_weight<sp/><sp/>=<sp/>ip.weight<sp/>*<sp/>mapped_normal_magnitude;</highlight></codeline>
|
||
|
|
<codeline lineno="195"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="196"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.quadrature.ds<sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>=<sp/>context.reference_surface_weight;</highlight></codeline>
|
||
|
|
<codeline lineno="197"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.quadrature.v_dot_n_scale<sp/>=<sp/>mapped_normal_magnitude<sp/>/<sp/>reference_normal_magnitude;</highlight></codeline>
|
||
|
|
<codeline lineno="198"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="199"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(context.quadrature.normal)<sp/>||<sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(context.reference_surface_weight)<sp/>||<sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(context.physical_surface_weight)<sp/>||<sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(context.quadrature.v_dot_n_scale))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="200"><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/>MappingStatus::NON_FINITE_MAPPING_RESULT;</highlight></codeline>
|
||
|
|
<codeline lineno="201"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="202"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="203"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>MappingStatus::VALID;</highlight></codeline>
|
||
|
|
<codeline lineno="204"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="205"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="206"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>DomainMapper::evaluate_point_variation(</highlight></codeline>
|
||
|
|
<codeline lineno="207"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementMappingData<sp/>&element_data,</highlight></codeline>
|
||
|
|
<codeline lineno="208"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementDisplacementData<sp/>&direction,</highlight></codeline>
|
||
|
|
<codeline lineno="209"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::ElementTransformation<sp/>&transformation,</highlight></codeline>
|
||
|
|
<codeline lineno="210"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint<sp/>&ip,</highlight></codeline>
|
||
|
|
<codeline lineno="211"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>MappingPointContext<sp/>&context,</highlight></codeline>
|
||
|
|
<codeline lineno="212"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>DMWorkspace<sp/>&workspace,</highlight></codeline>
|
||
|
|
<codeline lineno="213"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>MappingPointVariation<sp/>&variation</highlight></codeline>
|
||
|
|
<codeline lineno="214"><highlight class="normal"><sp/><sp/><sp/><sp/>)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
||
|
|
<codeline lineno="215"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>evaluate_point_variation_impl(element_data,<sp/>direction,<sp/>transformation,<sp/>ip,<sp/>context,<sp/>workspace,<sp/>variation,<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">);</highlight></codeline>
|
||
|
|
<codeline lineno="216"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="217"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="218"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>DomainMapper::evaluate_point_variation(</highlight></codeline>
|
||
|
|
<codeline lineno="219"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementMappingData&<sp/>element_data,</highlight></codeline>
|
||
|
|
<codeline lineno="220"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementDisplacementData&<sp/>direction,</highlight></codeline>
|
||
|
|
<codeline lineno="221"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::ElementTransformation&<sp/>transformation,</highlight></codeline>
|
||
|
|
<codeline lineno="222"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint&<sp/>ip,</highlight></codeline>
|
||
|
|
<codeline lineno="223"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>MappingPointContext&<sp/>context,</highlight></codeline>
|
||
|
|
<codeline lineno="224"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>DMWorkspace&<sp/>workspace,</highlight></codeline>
|
||
|
|
<codeline lineno="225"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>MappingPointVariation&<sp/>variation,</highlight></codeline>
|
||
|
|
<codeline lineno="226"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::DenseMatrix&<sp/>inverse_mesh_jacobian</highlight></codeline>
|
||
|
|
<codeline lineno="227"><highlight class="normal"><sp/><sp/><sp/><sp/>)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
||
|
|
<codeline lineno="228"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(inverse_mesh_jacobian.Height()<sp/>!=<sp/>m_options.dimension<sp/>||<sp/>inverse_mesh_jacobian.Width()<sp/>!=<sp/>m_options.dimension)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="229"><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/>MappingStatus::INVALID_DIMENSION;</highlight></codeline>
|
||
|
|
<codeline lineno="230"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="231"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(inverse_mesh_jacobian))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="232"><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/>MappingStatus::NON_FINITE_MAPPING_INPUT;</highlight></codeline>
|
||
|
|
<codeline lineno="233"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="234"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(direction.get_element().GetMapType()<sp/>!=<sp/>mfem::FiniteElement::VALUE)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="235"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Only<sp/>VALUE<sp/>fields<sp/>admit<sp/>the<sp/>reference-gradient<sp/>pullback<sp/>used<sp/>by<sp/>the<sp/>reuse<sp/>path.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="236"><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/>evaluate_point_variation(element_data,<sp/>direction,<sp/>transformation,<sp/>ip,<sp/>context,<sp/>workspace,<sp/>variation);</highlight></codeline>
|
||
|
|
<codeline lineno="237"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="238"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>evaluate_point_variation_impl(element_data,<sp/>direction,<sp/>transformation,<sp/>ip,<sp/>context,<sp/>workspace,<sp/>variation,<sp/>&inverse_mesh_jacobian);</highlight></codeline>
|
||
|
|
<codeline lineno="239"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="240"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="241"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>DomainMapper::evaluate_volume_direction_variation(</highlight></codeline>
|
||
|
|
<codeline lineno="242"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementMappingData<sp/>&element_data,</highlight></codeline>
|
||
|
|
<codeline lineno="243"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementDisplacementData<sp/>&direction,</highlight></codeline>
|
||
|
|
<codeline lineno="244"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::ElementTransformation<sp/>&transformation,</highlight></codeline>
|
||
|
|
<codeline lineno="245"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint<sp/>&ip,</highlight></codeline>
|
||
|
|
<codeline lineno="246"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>VolumeMappingContext<sp/>&context,</highlight></codeline>
|
||
|
|
<codeline lineno="247"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>DMWorkspace<sp/>&workspace,</highlight></codeline>
|
||
|
|
<codeline lineno="248"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>VolumeMappingVariation<sp/>&variation,</highlight></codeline>
|
||
|
|
<codeline lineno="249"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>DMReferenceDisplacementSample<sp/>&direction_sample,</highlight></codeline>
|
||
|
|
<codeline lineno="250"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>CompactificationPointData*<sp/>compactification_sample,</highlight></codeline>
|
||
|
|
<codeline lineno="251"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::DenseMatrix<sp/>&inverse_mesh_jacobian</highlight></codeline>
|
||
|
|
<codeline lineno="252"><highlight class="normal"><sp/><sp/><sp/><sp/>)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
||
|
|
<codeline lineno="253"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(inverse_mesh_jacobian.Height()<sp/>!=<sp/>m_options.dimension<sp/>||<sp/>inverse_mesh_jacobian.Width()<sp/>!=<sp/>m_options.dimension)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="254"><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/>MappingStatus::INVALID_DIMENSION;</highlight></codeline>
|
||
|
|
<codeline lineno="255"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="256"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>point_status<sp/>=<sp/>evaluate_point_variation_impl(element_data,<sp/>direction,<sp/>transformation,<sp/>ip,<sp/>context.mapping,<sp/>workspace,<sp/>variation.mapping,<sp/>&inverse_mesh_jacobian,<sp/>&direction_sample,<sp/>compactification_sample);</highlight></codeline>
|
||
|
|
<codeline lineno="257"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(point_status<sp/>!=<sp/>MappingStatus::VALID)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="258"><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/>point_status;</highlight></codeline>
|
||
|
|
<codeline lineno="259"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="260"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="261"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.weight_variation<sp/>=<sp/>(context.quadrature.physical_integration_weight<sp/>/<sp/>context.mapping.mapping_determinate)<sp/>*<sp/>variation.mapping.mapping_determinate_variation;</highlight></codeline>
|
||
|
|
<codeline lineno="262"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="263"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(variation.weight_variation))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="264"><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/>MappingStatus::NON_FINITE_MAPPING_RESULT;</highlight></codeline>
|
||
|
|
<codeline lineno="265"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="266"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="267"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>MappingStatus::VALID;</highlight></codeline>
|
||
|
|
<codeline lineno="268"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="269"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="270"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>DomainMapper::evaluate_face_variation(</highlight></codeline>
|
||
|
|
<codeline lineno="271"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementMappingData<sp/>&element_data,</highlight></codeline>
|
||
|
|
<codeline lineno="272"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementDisplacementData<sp/>&direction,</highlight></codeline>
|
||
|
|
<codeline lineno="273"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::FaceElementTransformations<sp/>&transformation,</highlight></codeline>
|
||
|
|
<codeline lineno="274"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>FaceElementSide<sp/>side,</highlight></codeline>
|
||
|
|
<codeline lineno="275"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint<sp/>&ip,</highlight></codeline>
|
||
|
|
<codeline lineno="276"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>FaceMappingContext<sp/>&context,</highlight></codeline>
|
||
|
|
<codeline lineno="277"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>DMWorkspace<sp/>&workspace,</highlight></codeline>
|
||
|
|
<codeline lineno="278"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>FaceMappingVariation<sp/>&variation</highlight></codeline>
|
||
|
|
<codeline lineno="279"><highlight class="normal"><sp/><sp/><sp/><sp/>)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
||
|
|
<codeline lineno="280"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>transformation.SetAllIntPoints(&ip);</highlight></codeline>
|
||
|
|
<codeline lineno="281"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::ElementTransformation<sp/>&<ref refid="namespaceserif_1_1utils_1_1parallel_1a8e0e78d02ab6c49ae588644e579bed79" kindref="member">element_transformation</ref><sp/>=<sp/>select_face_element_transformation(transformation,<sp/>side);</highlight></codeline>
|
||
|
|
<codeline lineno="282"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint<sp/>&element_integration_point<sp/>=<sp/>select_face_element_integration_point(transformation,<sp/>side);</highlight></codeline>
|
||
|
|
<codeline lineno="283"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="284"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>point_status<sp/>=<sp/>evaluate_point_variation(element_data,<sp/>direction,<sp/>element_transformation,<sp/>element_integration_point,<sp/>context.mapping,<sp/>workspace,<sp/>variation.mapping);</highlight></codeline>
|
||
|
|
<codeline lineno="285"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(point_status<sp/>!=<sp/>MappingStatus::VALID)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="286"><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/>point_status;</highlight></codeline>
|
||
|
|
<codeline lineno="287"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="288"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="289"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_reference_normal.SetSize(m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="290"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::CalcOrtho(transformation.Jacobian(),<sp/>workspace.m_reference_normal);</highlight></codeline>
|
||
|
|
<codeline lineno="291"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(side<sp/>==<sp/>FaceElementSide::ELEMENT_2)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="292"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_reference_normal<sp/>*=<sp/>-1.0;<sp/></highlight><highlight class="comment">//<sp/>Flip<sp/>the<sp/>normal<sp/>for<sp/>the<sp/>second<sp/>element<sp/>(e.g.<sp/>the<sp/>"other"<sp/>side<sp/>of<sp/>the<sp/>face)</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="293"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="294"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="295"><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">double</highlight><highlight class="normal"><sp/>reference_normal_magnitude<sp/>=<sp/>workspace.m_reference_normal.Norml2();</highlight></codeline>
|
||
|
|
<codeline lineno="296"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(reference_normal_magnitude)<sp/>||<sp/>reference_normal_magnitude<sp/><=<sp/>0)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="297"><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/>MappingStatus::NON_FINITE_MAPPING_RESULT;</highlight></codeline>
|
||
|
|
<codeline lineno="298"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="299"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="300"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.mapping.inverse_mapping_jacobian.MultTranspose(workspace.m_reference_normal,<sp/>workspace.m_vector_temp);</highlight></codeline>
|
||
|
|
<codeline lineno="301"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_mapped_normal<sp/>=<sp/>workspace.m_vector_temp;</highlight></codeline>
|
||
|
|
<codeline lineno="302"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_mapped_normal<sp/>*=<sp/>context.mapping.mapping_determinate;</highlight></codeline>
|
||
|
|
<codeline lineno="303"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="304"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.physical_normal_variation.SetSize(m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="305"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.mapping.inverse_mapping_jacobian_variation.MultTranspose(workspace.m_reference_normal,<sp/>variation.physical_normal_variation);</highlight></codeline>
|
||
|
|
<codeline lineno="306"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.physical_normal_variation<sp/>*=<sp/>context.mapping.mapping_determinate;</highlight></codeline>
|
||
|
|
<codeline lineno="307"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.physical_normal_variation.Add(variation.mapping.mapping_determinate_variation,<sp/>workspace.m_vector_temp);</highlight></codeline>
|
||
|
|
<codeline lineno="308"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="309"><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">double</highlight><highlight class="normal"><sp/>mapped_normal_magnitude<sp/>=<sp/>workspace.m_mapped_normal.Norml2();</highlight></codeline>
|
||
|
|
<codeline lineno="310"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(mapped_normal_magnitude)<sp/>||<sp/>mapped_normal_magnitude<sp/><=<sp/>0)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="311"><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/>MappingStatus::NON_FINITE_MAPPING_RESULT;</highlight></codeline>
|
||
|
|
<codeline lineno="312"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="313"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="314"><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">double</highlight><highlight class="normal"><sp/>mapped_normal_magnitude_variation<sp/>=<sp/>context.quadrature.normal<sp/>*<sp/>variation.physical_normal_variation;</highlight></codeline>
|
||
|
|
<codeline lineno="315"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="316"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.physical_normal_variation.Add(-mapped_normal_magnitude_variation,<sp/>context.quadrature.normal);</highlight></codeline>
|
||
|
|
<codeline lineno="317"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.physical_normal_variation<sp/>/=<sp/>mapped_normal_magnitude;</highlight></codeline>
|
||
|
|
<codeline lineno="318"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="319"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.physical_surface_weight_variation<sp/>=<sp/>ip.weight<sp/>*<sp/>mapped_normal_magnitude_variation;</highlight></codeline>
|
||
|
|
<codeline lineno="320"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.normal_flux_scale_variation<sp/>=<sp/>mapped_normal_magnitude_variation<sp/>/<sp/>reference_normal_magnitude;</highlight></codeline>
|
||
|
|
<codeline lineno="321"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="322"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(variation.physical_normal_variation)<sp/>||<sp/>!std::isfinite(variation.physical_surface_weight_variation)<sp/>||<sp/>!std::isfinite(variation.normal_flux_scale_variation))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="323"><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/>MappingStatus::NON_FINITE_MAPPING_RESULT;</highlight></codeline>
|
||
|
|
<codeline lineno="324"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="325"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="326"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>MappingStatus::VALID;</highlight></codeline>
|
||
|
|
<codeline lineno="327"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="328"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="329"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">bool</highlight><highlight class="normal"><sp/>DomainMapper::is_compactified_element(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::ElementTransformation<sp/>&transformation)<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">noexcept</highlight><highlight class="normal"><sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="330"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>transformation.Attribute<sp/>==<sp/>m_options.vacuum_element_attribute_ID;</highlight></codeline>
|
||
|
|
<codeline lineno="331"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="332"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="333"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>DomainMapper::get_dimension()<sp/>const<sp/>noexcept<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="334"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>m_options.dimension;</highlight></codeline>
|
||
|
|
<codeline lineno="335"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="336"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="337"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>compactification::ExteriorDomainMap&<sp/>DomainMapper::get_exterior_map()<sp/>const<sp/>noexcept<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="338"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>*m_exterior_map;</highlight></codeline>
|
||
|
|
<codeline lineno="339"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="340"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="341"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>DomainMapper::evaluate_point_impl(</highlight></codeline>
|
||
|
|
<codeline lineno="342"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementMappingData<sp/>&element_data,</highlight></codeline>
|
||
|
|
<codeline lineno="343"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::ElementTransformation<sp/>&transformation,</highlight></codeline>
|
||
|
|
<codeline lineno="344"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint<sp/>&ip,</highlight></codeline>
|
||
|
|
<codeline lineno="345"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>DMWorkspace<sp/>&workspace,</highlight></codeline>
|
||
|
|
<codeline lineno="346"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>MappingPointContext<sp/>&context,</highlight></codeline>
|
||
|
|
<codeline lineno="347"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>DMFixedMeshGeometry&<sp/><ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>,</highlight></codeline>
|
||
|
|
<codeline lineno="348"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>DMReferenceDisplacementSample<sp/>*displacement_sample</highlight></codeline>
|
||
|
|
<codeline lineno="349"><highlight class="normal"><sp/><sp/><sp/><sp/>)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
||
|
|
<codeline lineno="350"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>validate_element_data(element_data);</highlight></codeline>
|
||
|
|
<codeline lineno="351"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="352"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//TODO:<sp/>Factor<sp/>into<sp/>its<sp/>own<sp/>function</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="353"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//TODO:<sp/>Consider<sp/>if<sp/>these<sp/>should<sp/>return<sp/>or<sp/>throw.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="354"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(workspace.get_dimension()<sp/>!=<sp/>m_options.dimension)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="355"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a2447eb80947534cdf0fbd205251cfd9c" kindref="member">serif_error<std::invalid_argument, SERiFErrorCode::INVALID_DIMENSION></ref>(</highlight><highlight class="stringliteral">"Workspace<sp/>dimension<sp/>does<sp/>not<sp/>match<sp/>Domain<sp/>Mapper<sp/>dimension.<sp/>Workspace<sp/>dimension:<sp/>{},<sp/>Domain<sp/>Mapper<sp/>dimension:<sp/>{}"</highlight><highlight class="normal">,<sp/>workspace.get_dimension(),<sp/>m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="356"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="357"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(transformation.GetSpaceDim()<sp/>!=<sp/>m_options.dimension)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="358"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a2447eb80947534cdf0fbd205251cfd9c" kindref="member">serif_error<std::invalid_argument, SERiFErrorCode::INVALID_DIMENSION></ref>(</highlight><highlight class="stringliteral">"Element<sp/>transformation<sp/>space<sp/>dimension<sp/>does<sp/>not<sp/>match<sp/>Domain<sp/>Mapper<sp/>dimension.<sp/>Element<sp/>transformation<sp/>space<sp/>dimension:<sp/>{},<sp/>Domain<sp/>Mapper<sp/>dimension:<sp/>{}"</highlight><highlight class="normal">,<sp/>transformation.GetSpaceDim(),<sp/>m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="359"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="360"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(transformation.GetGeometryType()<sp/>!=<sp/>element_data.displacement.get_element().GetGeomType())<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="361"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a2447eb80947534cdf0fbd205251cfd9c" kindref="member">serif_error<std::invalid_argument, SERiFErrorCode::INVALID_ELEMENT_GEOMETRY></ref>(</highlight><highlight class="stringliteral">"Element<sp/>transformation<sp/>geometry<sp/>type<sp/>does<sp/>not<sp/>match<sp/>element<sp/>geometry<sp/>type.<sp/>Element<sp/>transformation<sp/>geometry<sp/>type:<sp/>{},<sp/>Element<sp/>geometry<sp/>type:<sp/>{}"</highlight><highlight class="normal">,<sp/>transformation.GetGeometryType(),<sp/>element_data.displacement.get_element().GetGeomType());</highlight></codeline>
|
||
|
|
<codeline lineno="362"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="363"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="364"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="365"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.reference_position<sp/>=<sp/><ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.position;</highlight></codeline>
|
||
|
|
<codeline lineno="366"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="367"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Get<sp/>the<sp/>displacement<sp/>field<sp/>value<sp/>and<sp/>its<sp/>Jacobian<sp/>at<sp/>the<sp/>integration<sp/>point.<sp/>Note<sp/>these<sp/>may<sp/>be<sp/>in<sp/>current<sp/>workspace<sp/>to<sp/>avoid<sp/>repeated<sp/>allocations<sp/>and<sp/>Multiplies.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="368"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(displacement_sample)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="369"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>TODO:<sp/>Generally<sp/>go<sp/>through<sp/>and<sp/>refactor<sp/>these<sp/>large<sp/>checks<sp/>into<sp/>their<sp/>own<sp/>functions,<sp/>at<sp/>the<sp/>very<sp/>least<sp/>it<sp/>will<sp/>be<sp/>easier<sp/>to<sp/>read</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="370"><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/>(element_data.displacement.get_element().GetMapType()<sp/>!=<sp/>mfem::FiniteElement::VALUE<sp/>||<sp/>displacement_sample->value.Size()<sp/>!=<sp/>m_options.dimension<sp/>||<sp/>displacement_sample->gradient.Height()<sp/>!=<sp/>m_options.dimension<sp/>||<sp/>displacement_sample->gradient.Width()<sp/>!=<sp/>m_options.dimension)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="371"><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/>MappingStatus::INVALID_DIMENSION;</highlight></codeline>
|
||
|
|
<codeline lineno="372"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="373"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="374"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_field_value<sp/>=<sp/>displacement_sample->value;</highlight></codeline>
|
||
|
|
<codeline lineno="375"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_field_jacobian.SetSize(m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="376"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::Mult(displacement_sample->gradient,<sp/><ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.inverse_jacobian,<sp/>workspace.m_field_jacobian);</highlight></codeline>
|
||
|
|
<codeline lineno="377"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}<sp/></highlight><highlight class="keywordflow">else</highlight><highlight class="normal"><sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="378"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>evaluate_displacement_field(element_data.displacement,<sp/>transformation,<sp/>ip,<sp/>workspace,<sp/>workspace.m_field_value,<sp/>workspace.m_field_jacobian,<sp/>&<ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.inverse_jacobian);</highlight></codeline>
|
||
|
|
<codeline lineno="379"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="380"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="381"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(context.reference_position)<sp/>||<sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(workspace.m_field_value)<sp/>||<sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(workspace.m_field_jacobian))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="382"><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/>MappingStatus::NON_FINITE_MAPPING_INPUT;</highlight></codeline>
|
||
|
|
<codeline lineno="383"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="384"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="385"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.displaced_position.SetSize(m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="386"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="387"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Get<sp/>the<sp/>position<sp/>of<sp/>the<sp/>point<sp/>in<sp/>physical<sp/>space<sp/>by<sp/>adding<sp/>the<sp/>displacement<sp/>to<sp/>the<sp/>reference<sp/>position.<sp/>Note<sp/>MFEM<sp/>really<sp/>dislikes<sp/>raw<sp/>arithmetic<sp/>operators</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="388"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>so<sp/>we<sp/>need<sp/>to<sp/>first<sp/>assign<sp/>the<sp/>reference<sp/>position<sp/>then<sp/>use<sp/>the<sp/>in<sp/>place<sp/>+=<sp/>operator.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="389"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.displaced_position<sp/>=<sp/>context.reference_position;</highlight></codeline>
|
||
|
|
<codeline lineno="390"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.displaced_position<sp/>+=<sp/>workspace.m_field_value;</highlight></codeline>
|
||
|
|
<codeline lineno="391"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="392"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.displacement_jacobian.SetSize(m_options.dimension,<sp/>m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="393"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.displacement_jacobian<sp/>=<sp/>workspace.m_field_jacobian;</highlight></codeline>
|
||
|
|
<codeline lineno="394"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="395"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Ensure<sp/>that<sp/>the<sp/>diagonal<sp/>of<sp/>the<sp/>displacement<sp/>Jacobian<sp/>is<sp/>incremented<sp/>by<sp/>1.0<sp/>to<sp/>account<sp/>for<sp/>the<sp/>identity<sp/>mapping<sp/>from<sp/>reference<sp/>to<sp/>physical<sp/>space.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="396"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>recall<sp/>that<sp/>r<sp/>=<sp/>x<sp/>+<sp/>d<sp/>(where<sp/>d<sp/>is<sp/>the<sp/>workspace.m_field_value<sp/>and<sp/>x<sp/>is<sp/>context.reference_position)<sp/>then<sp/>we<sp/>can<sp/>differentiate<sp/>this</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="397"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>component<sp/>wise<sp/>to<sp/>find<sp/>the<sp/>gradient<sp/>of<sp/>the<sp/>displaced<sp/>position<sp/>wrt.<sp/>the<sp/>mesh<sp/>coordinate<sp/>(reference<sp/>position).<sp/>E.g<sp/>as<sp/>you<sp/>move<sp/>along</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="398"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>the<sp/>mesh<sp/>coordinate<sp/>how<sp/>much<sp/>does<sp/>the<sp/>physical<sp/>coordinate<sp/>change<sp/>and<sp/>in<sp/>what<sp/>direction.<sp/>Lets<sp/>call<sp/>this<sp/>F</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="399"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>F<sp/>=<sp/>\frac{\partial<sp/>r_i}{\partial<sp/>x_j}<sp/>where<sp/>r<sp/>is<sp/>the<sp/>displaced<sp/>position<sp/>and<sp/>x<sp/>is<sp/>the<sp/>mesh<sp/>position.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="400"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>We<sp/>then<sp/>have<sp/>F<sp/>=<sp/>\frac{\partial<sp/>x_i}{\partial<sp/>x_j}<sp/>+<sp/>\frac{\partial<sp/>d_i}{x_j}<sp/>where<sp/>d<sp/>is<sp/>the<sp/>displacement<sp/>(recall<sp/>r<sp/>=<sp/>x<sp/>+<sp/>d)</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="401"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>By<sp/>definition<sp/>the<sp/>first<sp/>term<sp/>is<sp/>the<sp/>identity<sp/>matrix.<sp/>The<sp/>second<sp/>term<sp/>we<sp/>get<sp/>out<sp/>of<sp/>EvaluateField.<sp/>Thus<sp/>why<sp/>we<sp/>need<sp/>to<sp/>add<sp/>the<sp/>identity<sp/>matrix<sp/>here</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="402"><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/>i<sp/>=<sp/>0;<sp/>i<sp/><<sp/>m_options.dimension;<sp/>++i)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="403"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.displacement_jacobian(i,<sp/>i)<sp/>+=<sp/>1.0;</highlight></codeline>
|
||
|
|
<codeline lineno="404"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="405"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="406"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.compactified<sp/>=<sp/>is_compactified_element(transformation);</highlight></codeline>
|
||
|
|
<codeline lineno="407"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="408"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>This<sp/>branch<sp/>only<sp/>runs<sp/>for<sp/>vacuum<sp/>elements</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="409"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(context.compactified)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="410"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>There<sp/>are<sp/>two<sp/>things<sp/>that<sp/>we<sp/>need<sp/>to<sp/>the<sp/>mapping.<sp/>First<sp/>is<sp/>a<sp/>reference<sp/>coordinate<sp/>which<sp/>stroid<sp/>embeds<sp/>into<sp/>the<sp/>mesh<sp/>at<sp/>mesh<sp/>generation<sp/>time,<sp/>this<sp/>is</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="411"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>parameterized<sp/>from<sp/>0<sp/>-<sp/>1<sp/>where<sp/>0<sp/>is<sp/>the<sp/>model<sp/>surface<sp/>and<sp/>1<sp/>is<sp/>the<sp/>mesh<sp/>exterior<sp/>(what<sp/>will<sp/>becomes<sp/>the<sp/>compactified<sp/>infinity,<sp/>note<sp/>also<sp/>we<sp/>never<sp/>actually<sp/>evaluate<sp/>at<sp/>s=1;<sp/>rather<sp/>we<sp/>define<sp/>some<sp/>arbitrary<sp/>small<sp/>tolerance<sp/>to<sp/>approach<sp/>s=1).<sp/>Lets<sp/>call<sp/>this<sp/>s.<sp/>We<sp/>also<sp/>need</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="412"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>the<sp/>gradient<sp/>of<sp/>s<sp/>as<sp/>we<sp/>move<sp/>along<sp/>the<sp/>mesh<sp/>coordinates.<sp/>All<sp/>of<sp/>this<sp/>is<sp/>stashes<sp/>within<sp/>workspace.m_compactification_point.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="413"><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/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>coordinate_status<sp/>=<sp/>evaluate_compactification_coordinate(element_data.compactification,<sp/>transformation,<sp/>ip,<sp/>workspace,<sp/>workspace.m_compactification_point,<sp/>&<ref refid="namespaceserif_1_1geometry" kindref="compound">geometry</ref>.inverse_jacobian);</highlight></codeline>
|
||
|
|
<codeline lineno="414"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="415"><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/>(coordinate_status<sp/>!=<sp/>MappingStatus::VALID)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="416"><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/>coordinate_status;</highlight></codeline>
|
||
|
|
<codeline lineno="417"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="418"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="419"><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/>compactification::ExteriorMapInput<sp/>exterior_input{</highlight></codeline>
|
||
|
|
<codeline lineno="420"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.reference_position<sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>=<sp/>context.reference_position,</highlight></codeline>
|
||
|
|
<codeline lineno="421"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.displaced_position<sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>=<sp/>context.displaced_position,</highlight></codeline>
|
||
|
|
<codeline lineno="422"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.displacement_jacobian<sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>=<sp/>context.displacement_jacobian,</highlight></codeline>
|
||
|
|
<codeline lineno="423"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.compactification_coordinate_gradient<sp/>=<sp/>workspace.m_compactification_point.coordinate_gradient,</highlight></codeline>
|
||
|
|
<codeline lineno="424"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.compactification_coordinate<sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>=<sp/>workspace.m_compactification_point.compactification_coordinate,</highlight></codeline>
|
||
|
|
<codeline lineno="425"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
|
||
|
|
<codeline lineno="426"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="427"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>This<sp/>apply<sp/>whatever<sp/>the<sp/>exterior<sp/>map<sp/>is<sp/>to<sp/>generate<sp/>the<sp/>new<sp/>physical<sp/>exterior<sp/>coordinate<sp/>and<sp/>jacobian<sp/>between<sp/>physical<sp/>and<sp/>reference<sp/>space.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="428"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>In<sp/>general<sp/>we<sp/>have<sp/>only<sp/>implemented<sp/>a<sp/>kelvin<sp/>mapping;<sp/>however,<sp/>in<sp/>future<sp/>additional<sp/>mappings<sp/>may<sp/>be<sp/>implemented.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="429"><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/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>exterior_status<sp/>=<sp/>m_exterior_map->evaluate(exterior_input,<sp/>workspace.m_exterior_result);</highlight></codeline>
|
||
|
|
<codeline lineno="430"><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/>(exterior_status<sp/>!=<sp/>MappingStatus::VALID)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="431"><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/>exterior_status;</highlight></codeline>
|
||
|
|
<codeline lineno="432"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="433"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="434"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.physical_position<sp/>=<sp/>workspace.m_exterior_result.physical_position;</highlight></codeline>
|
||
|
|
<codeline lineno="435"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.mapping_jacobian<sp/><sp/>=<sp/>workspace.m_exterior_result.mapping_jacobian;</highlight></codeline>
|
||
|
|
<codeline lineno="436"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}<sp/></highlight><highlight class="keywordflow">else</highlight><highlight class="normal"><sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="437"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.physical_position<sp/>=<sp/>context.displaced_position;</highlight></codeline>
|
||
|
|
<codeline lineno="438"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.mapping_jacobian<sp/><sp/>=<sp/>context.displacement_jacobian;</highlight></codeline>
|
||
|
|
<codeline lineno="439"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="440"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="441"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Validation<sp/>work</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="442"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(context.physical_position)<sp/>||<sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(context.mapping_jacobian))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="443"><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/>MappingStatus::NON_FINITE_MAPPING_RESULT;</highlight></codeline>
|
||
|
|
<codeline lineno="444"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="445"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="446"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.mapping_determinate<sp/>=<sp/>context.mapping_jacobian.Det();</highlight></codeline>
|
||
|
|
<codeline lineno="447"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!std::isfinite(context.mapping_determinate))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="448"><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/>MappingStatus::NON_FINITE_MAPPING_RESULT;</highlight></codeline>
|
||
|
|
<codeline lineno="449"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="450"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(context.mapping_determinate<sp/><=<sp/>0.0)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="451"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>This<sp/>is<sp/>the<sp/>most<sp/>common<sp/>error<sp/>we<sp/>see<sp/>come<sp/>out<sp/>of<sp/>this<sp/>function,<sp/>specifically<sp/>it<sp/>is<sp/>common<sp/>when<sp/>we<sp/>try<sp/>to<sp/>deform<sp/>the<sp/>mesh<sp/>too<sp/>much<sp/>in<sp/>one<sp/>step.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="452"><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/>MappingStatus::NON_POSITIVE_DETERMINATE;</highlight></codeline>
|
||
|
|
<codeline lineno="453"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="454"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="455"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>context.inverse_mapping_jacobian.SetSize(m_options.dimension,<sp/>m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="456"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="457"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Finally<sp/>we<sp/>get<sp/>the<sp/>inverse<sp/>jacobian.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="458"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::CalcInverse(context.mapping_jacobian,<sp/>context.inverse_mapping_jacobian);</highlight></codeline>
|
||
|
|
<codeline lineno="459"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="460"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(context.inverse_mapping_jacobian))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="461"><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/>MappingStatus::NON_FINITE_MAPPING_RESULT;</highlight></codeline>
|
||
|
|
<codeline lineno="462"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="463"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="464"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>MappingStatus::VALID;</highlight></codeline>
|
||
|
|
<codeline lineno="465"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="466"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="467"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">void</highlight><highlight class="normal"><sp/>DomainMapper::validate_element_data(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementMappingData<sp/>&element_data)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
||
|
|
<codeline lineno="468"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementDisplacementData<sp/>&displacement<sp/>=<sp/>element_data.displacement;</highlight></codeline>
|
||
|
|
<codeline lineno="469"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementCompactificationData<sp/>&compactification<sp/>=<sp/>element_data.compactification;</highlight></codeline>
|
||
|
|
<codeline lineno="470"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="471"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(displacement.get_dimension()<sp/>!=<sp/>m_options.dimension)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="472"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a2447eb80947534cdf0fbd205251cfd9c" kindref="member">serif_error<std::invalid_argument, SERiFErrorCode::INVALID_DIMENSION></ref>(</highlight><highlight class="stringliteral">"Element<sp/>displacement<sp/>dimension<sp/>does<sp/>not<sp/>match<sp/>Domain<sp/>Mapper<sp/>dimension.<sp/>Element<sp/>displacement<sp/>dimension:<sp/>{},<sp/>Domain<sp/>Mapper<sp/>dimension:<sp/>{}"</highlight><highlight class="normal">,<sp/>displacement.get_dimension(),<sp/>m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="473"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="474"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="475"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(displacement.get_element().GetDim()<sp/>!=<sp/>m_options.dimension)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="476"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a2447eb80947534cdf0fbd205251cfd9c" kindref="member">serif_error<std::invalid_argument, SERiFErrorCode::INVALID_DIMENSION></ref>(</highlight><highlight class="stringliteral">"Element<sp/>displacement<sp/>finite<sp/>element<sp/>dimension<sp/>does<sp/>not<sp/>match<sp/>Domain<sp/>Mapper<sp/>dimension.<sp/>Element<sp/>displacement<sp/>finite<sp/>element<sp/>dimension:<sp/>{},<sp/>Domain<sp/>Mapper<sp/>dimension:<sp/>{}"</highlight><highlight class="normal">,<sp/>displacement.get_element().GetDim(),<sp/>m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="477"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="478"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="479"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(compactification.get_element().GetDim()<sp/>!=<sp/>m_options.dimension)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="480"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a2447eb80947534cdf0fbd205251cfd9c" kindref="member">serif_error<std::invalid_argument, SERiFErrorCode::INVALID_DIMENSION></ref>(</highlight><highlight class="stringliteral">"Element<sp/>compactification<sp/>finite<sp/>element<sp/>dimension<sp/>does<sp/>not<sp/>match<sp/>Domain<sp/>Mapper<sp/>dimension.<sp/>Element<sp/>compactification<sp/>finite<sp/>element<sp/>dimension:<sp/>{},<sp/>Domain<sp/>Mapper<sp/>dimension:<sp/>{}"</highlight><highlight class="normal">,<sp/>compactification.get_element().GetDim(),<sp/>m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="481"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="482"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="483"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(displacement.get_element().GetGeomType()<sp/>!=<sp/>compactification.get_element().GetGeomType())<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="484"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a2447eb80947534cdf0fbd205251cfd9c" kindref="member">serif_error<std::invalid_argument, SERiFErrorCode::INVALID_ELEMENT_GEOMETRY></ref>(</highlight><highlight class="stringliteral">"Element<sp/>displacement<sp/>finite<sp/>element<sp/>geometry<sp/>type<sp/>does<sp/>not<sp/>match<sp/>element<sp/>compactification<sp/>finite<sp/>element<sp/>geometry<sp/>type.<sp/>Element<sp/>displacement<sp/>finite<sp/>element<sp/>geometry<sp/>type:<sp/>{},<sp/>Element<sp/>compactification<sp/>finite<sp/>element<sp/>geometry<sp/>type:<sp/>{}"</highlight><highlight class="normal">,<sp/>displacement.get_element().GetGeomType(),<sp/>compactification.get_element().GetGeomType());</highlight></codeline>
|
||
|
|
<codeline lineno="485"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="486"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="487"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(compactification.get_element().GetRangeType()<sp/>!=<sp/>mfem::FiniteElement::SCALAR)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="488"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a2447eb80947534cdf0fbd205251cfd9c" kindref="member">serif_error<std::invalid_argument, SERiFErrorCode::INVALID_COMPACTIFICATION_RANK></ref>(</highlight><highlight class="stringliteral">"Element<sp/>compactification<sp/>finite<sp/>element<sp/>range<sp/>type<sp/>is<sp/>not<sp/>scalar.<sp/>Element<sp/>compactification<sp/>finite<sp/>element<sp/>range<sp/>type:<sp/>{}"</highlight><highlight class="normal">,<sp/>compactification.get_element().GetRangeType());</highlight></codeline>
|
||
|
|
<codeline lineno="489"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="490"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="491"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(compactification.get_element().GetMapType()<sp/>!=<sp/>mfem::FiniteElement::VALUE)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="492"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a2447eb80947534cdf0fbd205251cfd9c" kindref="member">serif_error<std::invalid_argument, SERiFErrorCode::INVALID_COMPACTIFICATION_TYPE></ref>(</highlight><highlight class="stringliteral">"Element<sp/>compactification<sp/>finite<sp/>element<sp/>map<sp/>type<sp/>is<sp/>not<sp/>value.<sp/>Element<sp/>compactification<sp/>finite<sp/>element<sp/>map<sp/>type:<sp/>{}"</highlight><highlight class="normal">,<sp/>compactification.get_element().GetMapType());</highlight></codeline>
|
||
|
|
<codeline lineno="493"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="494"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="495"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(compactification.get_element().GetDerivType()<sp/>!=<sp/>mfem::FiniteElement::GRAD)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="496"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a2447eb80947534cdf0fbd205251cfd9c" kindref="member">serif_error<std::invalid_argument, SERiFErrorCode::INVALID_COMPACTIFICATION_DERIVATIVE_TYPE></ref>(</highlight><highlight class="stringliteral">"Element<sp/>compactification<sp/>finite<sp/>element<sp/>derivative<sp/>type<sp/>is<sp/>not<sp/>gradient.<sp/>Element<sp/>compactification<sp/>finite<sp/>element<sp/>derivative<sp/>type:<sp/>{}"</highlight><highlight class="normal">,<sp/>compactification.get_element().GetDerivType());</highlight></codeline>
|
||
|
|
<codeline lineno="497"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="498"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="499"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(compactification.get_dof_count()<sp/>!=<sp/>compactification.get_element().GetDof())<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="500"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a2447eb80947534cdf0fbd205251cfd9c" kindref="member">serif_error<std::invalid_argument, SERiFErrorCode::INVALID_DOF_COUNT></ref>(</highlight><highlight class="stringliteral">"Element<sp/>compactification<sp/>finite<sp/>element<sp/>DOF<sp/>count<sp/>does<sp/>not<sp/>match<sp/>the<sp/>number<sp/>of<sp/>DOFs<sp/>in<sp/>the<sp/>finite<sp/>element.<sp/>Element<sp/>compactification<sp/>finite<sp/>element<sp/>DOF<sp/>count:<sp/>{},<sp/>Finite<sp/>element<sp/>DOF<sp/>count:<sp/>{}"</highlight><highlight class="normal">,<sp/>compactification.get_dof_count(),<sp/>compactification.get_element().GetDof());</highlight></codeline>
|
||
|
|
<codeline lineno="501"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="502"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="503"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="504"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Note<sp/>that<sp/>there<sp/>is<sp/>nothing<sp/>preventing<sp/>this<sp/>function<sp/>from<sp/>being<sp/>used<sp/>with<sp/>some<sp/>field<sp/>other<sp/>than<sp/>displacement;<sp/>perhaps<sp/>we<sp/>should<sp/>tighten<sp/>this?</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="505"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">void</highlight><highlight class="normal"><sp/>DomainMapper::evaluate_displacement_field(</highlight></codeline>
|
||
|
|
<codeline lineno="506"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementDisplacementData<sp/>&field,</highlight></codeline>
|
||
|
|
<codeline lineno="507"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::ElementTransformation<sp/>&transformation,</highlight></codeline>
|
||
|
|
<codeline lineno="508"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint<sp/>&ip,</highlight></codeline>
|
||
|
|
<codeline lineno="509"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>DMWorkspace<sp/>&workspace,</highlight></codeline>
|
||
|
|
<codeline lineno="510"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::Vector<sp/>&value,</highlight></codeline>
|
||
|
|
<codeline lineno="511"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::DenseMatrix<sp/>&jacobian,</highlight></codeline>
|
||
|
|
<codeline lineno="512"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::DenseMatrix<sp/>*inverse_mesh_jacobian</highlight></codeline>
|
||
|
|
<codeline lineno="513"><highlight class="normal"><sp/><sp/><sp/><sp/>)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
||
|
|
<codeline lineno="514"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::FiniteElement<sp/>&element<sp/><sp/>=<sp/><ref refid="namespaceserif_1_1utils_1_1serialize_1a73c789fd95f83d2f547facf07bc3ea79" kindref="member">field</ref>.get_element();</highlight></codeline>
|
||
|
|
<codeline lineno="515"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::DenseMatrix<sp/>&dof_matrix<sp/>=<sp/><ref refid="namespaceserif_1_1utils_1_1serialize_1a73c789fd95f83d2f547facf07bc3ea79" kindref="member">field</ref>.get_dof_matrix();</highlight></codeline>
|
||
|
|
<codeline lineno="516"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="517"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_shape.SetSize(element.GetDof());</highlight></codeline>
|
||
|
|
<codeline lineno="518"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="519"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>element.CalcShape(ip,<sp/>workspace.m_shape);</highlight></codeline>
|
||
|
|
<codeline lineno="520"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1discretization_1_1backend_1_1mfem_1af75f98db90618153ea4ca169945733f7a2063c1608d6e0baf80249c42e2be5804" kindref="member">value</ref>.SetSize(m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="521"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>dof_matrix.MultTranspose(workspace.m_shape,<sp/>value);</highlight></codeline>
|
||
|
|
<codeline lineno="522"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>jacobian.SetSize(m_options.dimension,<sp/>m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="523"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="524"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(inverse_mesh_jacobian<sp/>!=<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal"><sp/>||<sp/>element.GetMapType()<sp/>==<sp/>mfem::FiniteElement::VALUE)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="525"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_reference_dshape.SetSize(element.GetDof(),<sp/>m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="526"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>element.CalcDShape(ip,<sp/>workspace.m_reference_dshape);</highlight></codeline>
|
||
|
|
<codeline lineno="527"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Contract<sp/>DOFs<sp/>before<sp/>applying<sp/>fixed-mesh<sp/>geometry.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="528"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::MultAtB(dof_matrix,<sp/>workspace.m_reference_dshape,<sp/>workspace.m_reference_field_jacobian);</highlight></codeline>
|
||
|
|
<codeline lineno="529"><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/>mfem::DenseMatrix<sp/>&inverseMeshJacobian<sp/>=<sp/>inverse_mesh_jacobian<sp/>!=<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal"><sp/>?<sp/>*inverse_mesh_jacobian<sp/>:<sp/>transformation.InverseJacobian();</highlight></codeline>
|
||
|
|
<codeline lineno="530"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::Mult(workspace.m_reference_field_jacobian,<sp/>inverseMeshJacobian,<sp/>jacobian);</highlight></codeline>
|
||
|
|
<codeline lineno="531"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}<sp/></highlight><highlight class="keywordflow">else</highlight><highlight class="normal"><sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="532"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Retain<sp/>the<sp/>original<sp/>finite-element-specific<sp/>physical-gradient</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="533"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>path<sp/>for<sp/>mapping<sp/>types<sp/>without<sp/>the<sp/>ordinary<sp/>VALUE<sp/>pullback.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="534"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_mesh_dshape.SetSize(element.GetDof(),<sp/>m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="535"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>element.CalcPhysDShape(transformation,<sp/>workspace.m_mesh_dshape);</highlight></codeline>
|
||
|
|
<codeline lineno="536"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::MultAtB(dof_matrix,<sp/>workspace.m_mesh_dshape,<sp/>jacobian);</highlight></codeline>
|
||
|
|
<codeline lineno="537"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="538"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="539"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="540"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>DomainMapper::evaluate_compactification_coordinate(</highlight></codeline>
|
||
|
|
<codeline lineno="541"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementCompactificationData<sp/>&compactification,</highlight></codeline>
|
||
|
|
<codeline lineno="542"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::ElementTransformation<sp/>&transformation,</highlight></codeline>
|
||
|
|
<codeline lineno="543"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint<sp/>&ip,</highlight></codeline>
|
||
|
|
<codeline lineno="544"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>DMWorkspace<sp/>&workspace,</highlight></codeline>
|
||
|
|
<codeline lineno="545"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>CompactificationPointData<sp/>&point_data,</highlight></codeline>
|
||
|
|
<codeline lineno="546"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::DenseMatrix<sp/>*inverse_mesh_jacobian</highlight></codeline>
|
||
|
|
<codeline lineno="547"><highlight class="normal"><sp/><sp/><sp/><sp/>)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
||
|
|
<codeline lineno="548"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::FiniteElement&<sp/>element<sp/>=<sp/>compactification.get_element();</highlight></codeline>
|
||
|
|
<codeline lineno="549"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::Vector&<sp/>dofs<sp/>=<sp/>compactification.get_dofs();</highlight></codeline>
|
||
|
|
<codeline lineno="550"><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/>dof_count<sp/>=<sp/>element.GetDof();</highlight></codeline>
|
||
|
|
<codeline lineno="551"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="552"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(workspace.get_dimension()<sp/>!=<sp/>m_options.dimension<sp/>||<sp/>transformation.GetSpaceDim()<sp/>!=<sp/>m_options.dimension<sp/>||<sp/>element.GetDim()<sp/>!=<sp/>m_options.dimension)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="553"><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/>MappingStatus::INVALID_DIMENSION;</highlight></codeline>
|
||
|
|
<codeline lineno="554"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="555"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="556"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(dofs.Size()<sp/>!=<sp/>dof_count)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="557"><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/>MappingStatus::INVALID_DIMENSION;</highlight></codeline>
|
||
|
|
<codeline lineno="558"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="559"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="560"><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/>i<sp/>=<sp/>0;<sp/>i<sp/><<sp/>dof_count;<sp/>++i)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="561"><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/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(dofs(i)))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="562"><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/>MappingStatus::NON_FINITE_MAPPING_INPUT;</highlight></codeline>
|
||
|
|
<codeline lineno="563"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="564"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="565"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="566"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(inverse_mesh_jacobian<sp/>!=<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="567"><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/>(inverse_mesh_jacobian->Height()<sp/>!=<sp/>m_options.dimension<sp/>||<sp/>inverse_mesh_jacobian->Width()<sp/>!=<sp/>m_options.dimension)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="568"><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/>MappingStatus::INVALID_DIMENSION;</highlight></codeline>
|
||
|
|
<codeline lineno="569"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="570"><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/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(*inverse_mesh_jacobian))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="571"><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/>MappingStatus::NON_FINITE_MAPPING_INPUT;</highlight></codeline>
|
||
|
|
<codeline lineno="572"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="573"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="574"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>transformation.SetIntPoint(&ip);</highlight></codeline>
|
||
|
|
<codeline lineno="575"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_compactification_shape.SetSize(dof_count);</highlight></codeline>
|
||
|
|
<codeline lineno="576"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_compactification_dshape.SetSize(dof_count,<sp/>m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="577"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="578"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>element.CalcShape(ip,<sp/>workspace.m_compactification_shape);</highlight></codeline>
|
||
|
|
<codeline lineno="579"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(inverse_mesh_jacobian<sp/>!=<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="580"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_reference_dshape.SetSize(dof_count,<sp/>m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="581"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>element.CalcDShape(ip,<sp/>workspace.m_reference_dshape);</highlight></codeline>
|
||
|
|
<codeline lineno="582"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::Mult(workspace.m_reference_dshape,<sp/>*inverse_mesh_jacobian,<sp/>workspace.m_compactification_dshape);</highlight></codeline>
|
||
|
|
<codeline lineno="583"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}<sp/></highlight><highlight class="keywordflow">else</highlight><highlight class="normal"><sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="584"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>element.CalcPhysDShape(transformation,<sp/>workspace.m_compactification_dshape);</highlight></codeline>
|
||
|
|
<codeline lineno="585"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="586"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="587"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>point_data.compactification_coordinate<sp/>=<sp/>dofs<sp/>*<sp/>workspace.m_compactification_shape;</highlight></codeline>
|
||
|
|
<codeline lineno="588"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>point_data.coordinate_gradient.SetSize(m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="589"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_compactification_dshape.MultTranspose(dofs,<sp/>point_data.coordinate_gradient);</highlight></codeline>
|
||
|
|
<codeline lineno="590"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="591"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(point_data.compactification_coordinate))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="592"><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/>MappingStatus::NON_FINITE_MAPPING_RESULT;</highlight></codeline>
|
||
|
|
<codeline lineno="593"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="594"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="595"><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/>d<sp/>=<sp/>0;<sp/>d<sp/><<sp/>point_data.coordinate_gradient.Size();<sp/>++d)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="596"><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/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(point_data.coordinate_gradient(d)))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="597"><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/>MappingStatus::NON_FINITE_MAPPING_RESULT;</highlight></codeline>
|
||
|
|
<codeline lineno="598"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="599"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="600"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="601"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>MappingStatus::VALID;</highlight></codeline>
|
||
|
|
<codeline lineno="602"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="603"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="604"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>DomainMapper::evaluate_point_variation_impl(</highlight></codeline>
|
||
|
|
<codeline lineno="605"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementMappingData<sp/>&element_data,</highlight></codeline>
|
||
|
|
<codeline lineno="606"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>ElementDisplacementData<sp/>&direction,</highlight></codeline>
|
||
|
|
<codeline lineno="607"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::ElementTransformation<sp/>&transformation,</highlight></codeline>
|
||
|
|
<codeline lineno="608"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint<sp/>&ip,</highlight></codeline>
|
||
|
|
<codeline lineno="609"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>MappingPointContext<sp/>&context,</highlight></codeline>
|
||
|
|
<codeline lineno="610"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>DMWorkspace<sp/>&workspace,</highlight></codeline>
|
||
|
|
<codeline lineno="611"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>MappingPointVariation<sp/>&variation,</highlight></codeline>
|
||
|
|
<codeline lineno="612"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::DenseMatrix<sp/>*inverse_mesh_jacobian,</highlight></codeline>
|
||
|
|
<codeline lineno="613"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>DMReferenceDisplacementSample*<sp/>direction_sample,</highlight></codeline>
|
||
|
|
<codeline lineno="614"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>CompactificationPointData*<sp/>compactification_sample</highlight></codeline>
|
||
|
|
<codeline lineno="615"><highlight class="normal"><sp/><sp/><sp/><sp/>)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
||
|
|
<codeline lineno="616"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>validate_element_data(element_data);</highlight></codeline>
|
||
|
|
<codeline lineno="617"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>A<sp/>direction<sp/>that<sp/>is<sp/>the<sp/>element's<sp/>own<sp/>displacement<sp/>data<sp/>was<sp/>validated<sp/>just<sp/>above.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="618"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(&direction<sp/>!=<sp/>&element_data.displacement)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="619"><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/>ElementMappingData<sp/>direction_data<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="620"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.displacement<sp/>=<sp/>direction,</highlight></codeline>
|
||
|
|
<codeline lineno="621"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.compactification<sp/>=<sp/>element_data.compactification</highlight></codeline>
|
||
|
|
<codeline lineno="622"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
|
||
|
|
<codeline lineno="623"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="624"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>validate_element_data(direction_data);</highlight></codeline>
|
||
|
|
<codeline lineno="625"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="626"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="627"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(element_data.displacement.get_dof_count()<sp/>!=<sp/>direction.get_dof_count())<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="628"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a2447eb80947534cdf0fbd205251cfd9c" kindref="member">serif_error<std::invalid_argument, SERiFErrorCode::INVALID_DOF_COUNT></ref>(</highlight><highlight class="stringliteral">"Element<sp/>displacement<sp/>DOF<sp/>count<sp/>does<sp/>not<sp/>match<sp/>direction<sp/>DOF<sp/>count.<sp/>Element<sp/>displacement<sp/>DOF<sp/>count:<sp/>{},<sp/>Direction<sp/>DOF<sp/>count:<sp/>{}"</highlight><highlight class="normal">,<sp/>element_data.displacement.get_dof_count(),<sp/>direction.get_dof_count());</highlight></codeline>
|
||
|
|
<codeline lineno="629"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="630"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="631"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(workspace.get_dimension()<sp/>!=<sp/>m_options.dimension)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="632"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a2447eb80947534cdf0fbd205251cfd9c" kindref="member">serif_error<std::invalid_argument, SERiFErrorCode::INVALID_DIMENSION></ref>(</highlight><highlight class="stringliteral">"Workspace<sp/>dimension<sp/>does<sp/>not<sp/>match<sp/>Domain<sp/>Mapper<sp/>dimension.<sp/>Workspace<sp/>dimension:<sp/>{},<sp/>Domain<sp/>Mapper<sp/>dimension:<sp/>{}"</highlight><highlight class="normal">,<sp/>workspace.get_dimension(),<sp/>m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="633"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="634"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="635"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(context.compactified<sp/>!=<sp/>is_compactified_element(transformation))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="636"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a2447eb80947534cdf0fbd205251cfd9c" kindref="member">serif_error<std::invalid_argument, SERiFErrorCode::COMPACTIFICATION_MISMATCH></ref>(</highlight><highlight class="stringliteral">"Mapping<sp/>context<sp/>compactification<sp/>status<sp/>does<sp/>not<sp/>match<sp/>element<sp/>transformation<sp/>compactification<sp/>status.<sp/>Mapping<sp/>context<sp/>compactified:<sp/>{},<sp/>Element<sp/>transformation<sp/>compactified:<sp/>{}"</highlight><highlight class="normal">,<sp/>context.compactified,<sp/>is_compactified_element(transformation));</highlight></codeline>
|
||
|
|
<codeline lineno="637"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="638"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="639"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>transformation.SetIntPoint(&ip);</highlight></codeline>
|
||
|
|
<codeline lineno="640"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Tensor<sp/>contractions<sp/>and<sp/>transpose<sp/>actions<sp/>can<sp/>supply<sp/>canonical<sp/>value/gradient<sp/>samples<sp/>directly.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="641"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(direction_sample)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="642"><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/>(direction.get_element().GetMapType()<sp/>!=<sp/>mfem::FiniteElement::VALUE<sp/>||<sp/>direction_sample->value.Size()<sp/>!=<sp/>m_options.dimension<sp/>||<sp/>direction_sample->gradient.Height()<sp/>!=<sp/>m_options.dimension<sp/>||<sp/>direction_sample->gradient.Width()<sp/>!=<sp/>m_options.dimension)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="643"><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/>MappingStatus::INVALID_DIMENSION;</highlight></codeline>
|
||
|
|
<codeline lineno="644"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="645"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_field_value<sp/>=<sp/>direction_sample->value;</highlight></codeline>
|
||
|
|
<codeline lineno="646"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_field_jacobian.SetSize(m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="647"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::Mult(direction_sample->gradient,<sp/>inverse_mesh_jacobian<sp/>?<sp/>*inverse_mesh_jacobian<sp/>:<sp/>transformation.InverseJacobian(),<sp/>workspace.m_field_jacobian);</highlight></codeline>
|
||
|
|
<codeline lineno="648"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}<sp/></highlight><highlight class="keywordflow">else</highlight><highlight class="normal"><sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="649"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>evaluate_displacement_field(direction,<sp/>transformation,<sp/>ip,<sp/>workspace,<sp/>workspace.m_field_value,<sp/>workspace.m_field_jacobian,<sp/>inverse_mesh_jacobian);</highlight></codeline>
|
||
|
|
<codeline lineno="650"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="651"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="652"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(workspace.m_field_jacobian)<sp/>||<sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(workspace.m_field_value))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="653"><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/>MappingStatus::NON_FINITE_MAPPING_INPUT;</highlight></codeline>
|
||
|
|
<codeline lineno="654"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="655"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="656"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.displacement_variation<sp/>=<sp/>workspace.m_field_value;</highlight></codeline>
|
||
|
|
<codeline lineno="657"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.displacement_jacobian_variation<sp/>=<sp/>workspace.m_field_jacobian;</highlight></codeline>
|
||
|
|
<codeline lineno="658"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="659"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(context.compactified)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="660"><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/>(compactification_sample<sp/>==<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="661"><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/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>coordinate_status<sp/>=<sp/>evaluate_compactification_coordinate(element_data.compactification,<sp/>transformation,<sp/>ip,<sp/>workspace,<sp/>workspace.m_compactification_point,<sp/>inverse_mesh_jacobian);</highlight></codeline>
|
||
|
|
<codeline lineno="662"><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/>(coordinate_status<sp/>!=<sp/>MappingStatus::VALID)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="663"><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/>coordinate_status;</highlight></codeline>
|
||
|
|
<codeline lineno="664"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="665"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>compactification_sample<sp/>=<sp/>&workspace.m_compactification_point;</highlight></codeline>
|
||
|
|
<codeline lineno="666"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="667"><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/>(compactification_sample->coordinate_gradient.Size()<sp/>!=<sp/>m_options.dimension)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="668"><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/>MappingStatus::INVALID_DIMENSION;</highlight></codeline>
|
||
|
|
<codeline lineno="669"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="670"><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/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(compactification_sample->compactification_coordinate)<sp/>||<sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(compactification_sample->coordinate_gradient))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="671"><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/>MappingStatus::NON_FINITE_MAPPING_INPUT;</highlight></codeline>
|
||
|
|
<codeline lineno="672"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="673"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="674"><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/>compactification::ExteriorMapInput<sp/>exterior_input{</highlight></codeline>
|
||
|
|
<codeline lineno="675"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.reference_position<sp/>=<sp/>context.reference_position,</highlight></codeline>
|
||
|
|
<codeline lineno="676"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.displaced_position<sp/>=<sp/>context.displaced_position,</highlight></codeline>
|
||
|
|
<codeline lineno="677"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.displacement_jacobian<sp/>=<sp/>context.displacement_jacobian,</highlight></codeline>
|
||
|
|
<codeline lineno="678"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.compactification_coordinate_gradient<sp/>=<sp/>compactification_sample->coordinate_gradient,</highlight></codeline>
|
||
|
|
<codeline lineno="679"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.compactification_coordinate<sp/>=<sp/>compactification_sample->compactification_coordinate,</highlight></codeline>
|
||
|
|
<codeline lineno="680"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
|
||
|
|
<codeline lineno="681"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="682"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_exterior_result.physical_position<sp/>=<sp/>context.physical_position;</highlight></codeline>
|
||
|
|
<codeline lineno="683"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>workspace.m_exterior_result.mapping_jacobian<sp/>=<sp/>context.mapping_jacobian;</highlight></codeline>
|
||
|
|
<codeline lineno="684"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="685"><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/>compactification::ExteriorMapDirection<sp/>exterior_direction<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="686"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.displaced_position_variation<sp/>=<sp/>variation.displacement_variation,</highlight></codeline>
|
||
|
|
<codeline lineno="687"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>.displacement_jacobian_variation<sp/>=<sp/>variation.displacement_jacobian_variation</highlight></codeline>
|
||
|
|
<codeline lineno="688"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
|
||
|
|
<codeline lineno="689"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="690"><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/><ref refid="namespaceserif_1_1mapping_1a4d8ed5ebf76c79667efd6777f3ee3b97" kindref="member">MappingStatus</ref><sp/>exterior_status<sp/>=<sp/>m_exterior_map->evaluate_variation(exterior_input,<sp/>workspace.m_exterior_result,<sp/>exterior_direction,<sp/>workspace.m_exterior_variation);</highlight></codeline>
|
||
|
|
<codeline lineno="691"><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/>(exterior_status<sp/>!=<sp/>MappingStatus::VALID)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="692"><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/>exterior_status;</highlight></codeline>
|
||
|
|
<codeline lineno="693"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="694"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="695"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.physical_position_variation<sp/>=<sp/>workspace.m_exterior_variation.physical_position_variation;</highlight></codeline>
|
||
|
|
<codeline lineno="696"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.mapping_jacobian_variation<sp/>=<sp/>workspace.m_exterior_variation.mapping_jacobian_variation;</highlight></codeline>
|
||
|
|
<codeline lineno="697"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}<sp/></highlight><highlight class="keywordflow">else</highlight><highlight class="normal"><sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="698"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.physical_position_variation<sp/>=<sp/>variation.displacement_variation;</highlight></codeline>
|
||
|
|
<codeline lineno="699"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.mapping_jacobian_variation<sp/>=<sp/>variation.displacement_jacobian_variation;</highlight></codeline>
|
||
|
|
<codeline lineno="700"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="701"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="702"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::Mult(context.inverse_mapping_jacobian,<sp/>variation.mapping_jacobian_variation,<sp/>workspace.m_matrix_temp_1);</highlight></codeline>
|
||
|
|
<codeline lineno="703"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>trace<sp/>=<sp/>0.0;</highlight></codeline>
|
||
|
|
<codeline lineno="704"><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/>dimID<sp/>=<sp/>0;<sp/>dimID<sp/><<sp/>m_options.dimension;<sp/>++dimID)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="705"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>trace<sp/>+=<sp/>workspace.m_matrix_temp_1(dimID,<sp/>dimID);</highlight></codeline>
|
||
|
|
<codeline lineno="706"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="707"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.mapping_determinate_variation<sp/>=<sp/>context.mapping_determinate<sp/>*<sp/>trace;</highlight></codeline>
|
||
|
|
<codeline lineno="708"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="709"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.inverse_mapping_jacobian_variation.SetSize(m_options.dimension,<sp/>m_options.dimension);</highlight></codeline>
|
||
|
|
<codeline lineno="710"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::Mult(workspace.m_matrix_temp_1,<sp/>context.inverse_mapping_jacobian,<sp/>variation.inverse_mapping_jacobian_variation);</highlight></codeline>
|
||
|
|
<codeline lineno="711"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>variation.inverse_mapping_jacobian_variation<sp/>*=<sp/>-1.0;</highlight></codeline>
|
||
|
|
<codeline lineno="712"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="713"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(variation.physical_position_variation)<sp/>||<sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(variation.mapping_jacobian_variation)<sp/>||<sp/>!std::isfinite(variation.mapping_determinate_variation)<sp/>||<sp/>!<ref refid="namespaceserif_1_1utils_1_1checks_1a4eb8e1d2bdd4bdae0242814f9018537e" kindref="member">is_finite</ref>(variation.inverse_mapping_jacobian_variation))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="714"><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/>MappingStatus::NON_FINITE_MAPPING_RESULT;</highlight></codeline>
|
||
|
|
<codeline lineno="715"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="716"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="717"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>MappingStatus::VALID;</highlight></codeline>
|
||
|
|
<codeline lineno="718"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="719"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="720"><highlight class="normal"><sp/><sp/><sp/><sp/>mfem::ElementTransformation<sp/>&DomainMapper::select_face_element_transformation(</highlight></codeline>
|
||
|
|
<codeline lineno="721"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::FaceElementTransformations<sp/>&transformations,</highlight></codeline>
|
||
|
|
<codeline lineno="722"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>FaceElementSide<sp/>side</highlight></codeline>
|
||
|
|
<codeline lineno="723"><highlight class="normal"><sp/><sp/><sp/><sp/>)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="724"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(side<sp/>==<sp/>FaceElementSide::ELEMENT_1)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="725"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a980a92d4199f7d6b867ff79cdb9638b6" kindref="member">serif_assert<SERiFErrorCode::MISSING_FACE_ELEMENT></ref>(transformations.Elem1<sp/>!=<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">,<sp/></highlight><highlight class="stringliteral">"Face<sp/>element<sp/>transformation<sp/>for<sp/>ELEMENT_1<sp/>is<sp/>null."</highlight><highlight class="normal">);</highlight></codeline>
|
||
|
|
<codeline lineno="726"><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/>*transformations.Elem1;</highlight></codeline>
|
||
|
|
<codeline lineno="727"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="728"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="729"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1errors_1a980a92d4199f7d6b867ff79cdb9638b6" kindref="member">serif_assert<SERiFErrorCode::MISSING_FACE_ELEMENT></ref>(transformations.Elem2<sp/>!=<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">,<sp/></highlight><highlight class="stringliteral">"Face<sp/>element<sp/>transformation<sp/>for<sp/>ELEMENT_2<sp/>is<sp/>null."</highlight><highlight class="normal">);</highlight></codeline>
|
||
|
|
<codeline lineno="730"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>*transformations.Elem2;</highlight></codeline>
|
||
|
|
<codeline lineno="731"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="732"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="733"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::IntegrationPoint<sp/>&DomainMapper::select_face_element_integration_point(</highlight></codeline>
|
||
|
|
<codeline lineno="734"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>mfem::FaceElementTransformations<sp/>&transformations,</highlight></codeline>
|
||
|
|
<codeline lineno="735"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>FaceElementSide<sp/>side</highlight></codeline>
|
||
|
|
<codeline lineno="736"><highlight class="normal"><sp/><sp/><sp/><sp/>)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="737"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mfem::ElementTransformation&<sp/><ref refid="namespaceserif_1_1utils_1_1parallel_1a8e0e78d02ab6c49ae588644e579bed79" kindref="member">element_transformation</ref><sp/>=<sp/>select_face_element_transformation(transformations,<sp/>side);</highlight></codeline>
|
||
|
|
<codeline lineno="738"><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_1utils_1_1parallel_1a8e0e78d02ab6c49ae588644e579bed79" kindref="member">element_transformation</ref>.GetIntPoint();</highlight></codeline>
|
||
|
|
<codeline lineno="739"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="740"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="741"><highlight class="normal">}</highlight></codeline>
|
||
|
|
</programlisting>
|
||
|
|
<location file="src/lib/serif/backend/mfem/mapping/domain_mapper.cpp"/>
|
||
|
|
</compounddef>
|
||
|
|
</doxygen>
|