MeanField_SERiF/docs/xml/backend_2mfem_2serialize_2top_8cpp.xml

<?xml version='1.0' encoding='UTF-8' standalone='no'?>
<doxygen xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="compound.xsd" version="1.13.2" xml:lang="en-US">
  <compounddef id="backend_2mfem_2serialize_2top_8cpp" kind="file" language="C++">
    <compoundname>top.cpp</compoundname>
    <includes refid="backend_2mfem_2serialize_2top_8hpp" local="yes">serif/backend/mfem/serialize/top.hpp</includes>
    <includes local="no">algorithm</includes>
    <includes local="no">limits</includes>
    <includes local="no">map</includes>
    <includes local="no">optional</includes>
    <includes local="no">tuple</includes>
    <includes local="no">stdexcept</includes>
    <includes local="no">utility</includes>
    <includes local="no">vector</includes>
    <includes local="yes">mfem/fem/gridfunc.hpp</includes>
    <includes refid="field__values_8hpp" local="yes">serif/backend/mfem/discretization/field_values.hpp</includes>
    <incdepgraph>
      <node id="169">
        <label>serif/backend/mfem/discretization/basis_evaluation.hpp</label>
        <link refid="basis__evaluation_8hpp"/>
        <childnode refid="6" relation="include">
        </childnode>
        <childnode refid="15" relation="include">
        </childnode>
        <childnode refid="170" relation="include">
        </childnode>
        <childnode refid="171" relation="include">
        </childnode>
      </node>
      <node id="11">
        <label>serif/backend/mfem/discretization/dof_map.hpp</label>
        <link refid="dof__map_8hpp"/>
        <childnode refid="4" relation="include">
        </childnode>
        <childnode refid="12" relation="include">
        </childnode>
        <childnode refid="6" relation="include">
        </childnode>
        <childnode refid="13" relation="include">
        </childnode>
        <childnode refid="55" relation="include">
        </childnode>
        <childnode refid="56" relation="include">
        </childnode>
        <childnode refid="57" relation="include">
        </childnode>
        <childnode refid="41" relation="include">
        </childnode>
        <childnode refid="58" relation="include">
        </childnode>
      </node>
      <node id="162">
        <label>serif/backend/mfem/discretization/field_values.hpp</label>
        <link refid="field__values_8hpp"/>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="163" relation="include">
        </childnode>
        <childnode refid="165" relation="include">
        </childnode>
      </node>
      <node id="13">
        <label>serif/backend/mfem/discretization/reconstruction.hpp</label>
        <link refid="reconstruction_8hpp"/>
        <childnode refid="6" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
        <childnode refid="15" relation="include">
        </childnode>
        <childnode refid="16" relation="include">
        </childnode>
        <childnode refid="40" relation="include">
        </childnode>
      </node>
      <node id="171">
        <label>serif/backend/mfem/discretization/tensor_basis.hpp</label>
        <link refid="tensor__basis_8hpp"/>
        <childnode refid="44" relation="include">
        </childnode>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="4" relation="include">
        </childnode>
        <childnode refid="142" relation="include">
        </childnode>
        <childnode refid="6" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="172" relation="include">
        </childnode>
        <childnode refid="173" relation="include">
        </childnode>
        <childnode refid="41" relation="include">
        </childnode>
        <childnode refid="13" relation="include">
        </childnode>
      </node>
      <node id="164">
        <label>serif/backend/mfem/discretization/term_binding.hpp</label>
        <link refid="term__binding_8hpp"/>
        <childnode refid="44" relation="include">
        </childnode>
        <childnode refid="132" relation="include">
        </childnode>
        <childnode refid="94" relation="include">
        </childnode>
        <childnode refid="11" relation="include">
        </childnode>
        <childnode refid="41" relation="include">
        </childnode>
      </node>
      <node id="163">
        <label>serif/backend/mfem/discretization/term_coefficients.hpp</label>
        <link refid="term__coefficients_8hpp"/>
        <childnode refid="108" relation="include">
        </childnode>
        <childnode refid="164" relation="include">
        </childnode>
      </node>
      <node id="59">
        <label>serif/backend/mfem/geometry/deformation/computations.hpp</label>
        <link refid="backend_2mfem_2geometry_2deformation_2computations_8hpp"/>
        <childnode refid="44" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="60" relation="include">
        </childnode>
        <childnode refid="72" relation="include">
        </childnode>
        <childnode refid="112" relation="include">
        </childnode>
        <childnode refid="116" relation="include">
        </childnode>
        <childnode refid="115" relation="include">
        </childnode>
      </node>
      <node id="60">
        <label>serif/backend/mfem/geometry/deformation/constrained_radial_extension.hpp</label>
        <link refid="constrained__radial__extension_8hpp"/>
        <childnode refid="6" relation="include">
        </childnode>
        <childnode refid="8" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="56" relation="include">
        </childnode>
        <childnode refid="61" relation="include">
        </childnode>
        <childnode refid="62" relation="include">
        </childnode>
        <childnode refid="69" relation="include">
        </childnode>
      </node>
      <node id="72">
        <label>serif/backend/mfem/geometry/deformation/inputs.hpp</label>
        <link refid="geometry_2deformation_2inputs_8hpp"/>
        <childnode refid="44" relation="include">
        </childnode>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="73" relation="include">
        </childnode>
        <childnode refid="8" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
        <childnode refid="56" relation="include">
        </childnode>
        <childnode refid="74" relation="include">
        </childnode>
        <childnode refid="94" relation="include">
        </childnode>
        <childnode refid="66" relation="include">
        </childnode>
        <childnode refid="96" relation="include">
        </childnode>
        <childnode refid="97" relation="include">
        </childnode>
        <childnode refid="108" relation="include">
        </childnode>
      </node>
      <node id="112">
        <label>serif/backend/mfem/geometry/deformation/radial_interpolation.hpp</label>
        <link refid="backend_2mfem_2geometry_2deformation_2radial__interpolation_8hpp"/>
        <childnode refid="8" relation="include">
        </childnode>
        <childnode refid="62" relation="include">
        </childnode>
        <childnode refid="113" relation="include">
        </childnode>
      </node>
      <node id="113">
        <label>serif/backend/mfem/geometry/deformation/reference_geometry.hpp</label>
        <link refid="reference__geometry_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
        <childnode refid="64" relation="include">
        </childnode>
        <childnode refid="114" relation="include">
        </childnode>
        <childnode refid="115" relation="include">
        </childnode>
      </node>
      <node id="116">
        <label>serif/backend/mfem/geometry/deformation/surface_metric.hpp</label>
        <link refid="surface__metric_8hpp"/>
        <childnode refid="65" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
        <childnode refid="64" relation="include">
        </childnode>
        <childnode refid="114" relation="include">
        </childnode>
        <childnode refid="115" relation="include">
        </childnode>
      </node>
      <node id="178">
        <label>serif/backend/mfem/mapping/context.hpp</label>
        <link refid="context_8hpp"/>
        <childnode refid="176" relation="include">
        </childnode>
        <childnode refid="177" relation="include">
        </childnode>
      </node>
      <node id="174">
        <label>serif/backend/mfem/mapping/transformations.hpp</label>
        <link refid="transformations_8hpp"/>
        <childnode refid="175" relation="include">
        </childnode>
        <childnode refid="178" relation="include">
        </childnode>
      </node>
      <node id="175">
        <label>serif/backend/mfem/mapping/variation.hpp</label>
        <link refid="variation_8hpp"/>
        <childnode refid="176" relation="include">
        </childnode>
        <childnode refid="177" relation="include">
        </childnode>
      </node>
      <node id="165">
        <label>serif/backend/mfem/operators/operand_evaluation.hpp</label>
        <link refid="operand__evaluation_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="55" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="15" relation="include">
        </childnode>
        <childnode refid="166" relation="include">
        </childnode>
        <childnode refid="167" relation="include">
        </childnode>
        <childnode refid="169" relation="include">
        </childnode>
        <childnode refid="11" relation="include">
        </childnode>
        <childnode refid="174" relation="include">
        </childnode>
        <childnode refid="41" relation="include">
        </childnode>
        <childnode refid="179" relation="include">
        </childnode>
      </node>
      <node id="2">
        <label>serif/backend/mfem/serialize/top.hpp</label>
        <link refid="backend_2mfem_2serialize_2top_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="4" relation="include">
        </childnode>
        <childnode refid="5" relation="include">
        </childnode>
        <childnode refid="6" relation="include">
        </childnode>
        <childnode refid="7" relation="include">
        </childnode>
        <childnode refid="8" relation="include">
        </childnode>
        <childnode refid="9" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="11" relation="include">
        </childnode>
        <childnode refid="59" relation="include">
        </childnode>
        <childnode refid="117" relation="include">
        </childnode>
        <childnode refid="122" relation="include">
        </childnode>
        <childnode refid="137" relation="include">
        </childnode>
        <childnode refid="138" relation="include">
        </childnode>
        <childnode refid="139" relation="include">
        </childnode>
        <childnode refid="149" relation="include">
        </childnode>
        <childnode refid="157" relation="include">
        </childnode>
      </node>
      <node id="96">
        <label>serif/backend/mfem/shared/reference_tracking.hpp</label>
        <link refid="reference__tracking_8hpp"/>
        <childnode refid="14" relation="include">
        </childnode>
      </node>
      <node id="18">
        <label>serif/dimensions/dimensions.hpp</label>
        <link refid="dimensions_8hpp"/>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="19" relation="include">
        </childnode>
      </node>
      <node id="121">
        <label>serif/dimensions/partial.hpp</label>
        <link refid="partial_8hpp"/>
        <childnode refid="19" relation="include">
        </childnode>
        <childnode refid="24" relation="include">
        </childnode>
      </node>
      <node id="19">
        <label>serif/dimensions/quantities.hpp</label>
        <link refid="quantities_8hpp"/>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="20" relation="include">
        </childnode>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="23" relation="include">
        </childnode>
        <childnode refid="24" relation="include">
        </childnode>
      </node>
      <node id="120">
        <label>serif/dimensions/type_alias.hpp</label>
        <link refid="type__alias_8hpp"/>
        <childnode refid="19" relation="include">
        </childnode>
        <childnode refid="18" relation="include">
        </childnode>
      </node>
      <node id="25">
        <label>serif/discretization/blocks/base.hpp</label>
        <link refid="base_8hpp"/>
      </node>
      <node id="133">
        <label>serif/discretization/blocks/concepts.hpp</label>
        <link refid="discretization_2blocks_2concepts_8hpp"/>
        <childnode refid="16" relation="include">
        </childnode>
      </node>
      <node id="16">
        <label>serif/discretization/blocks/fields.hpp</label>
        <link refid="fields_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="18" relation="include">
        </childnode>
        <childnode refid="25" relation="include">
        </childnode>
        <childnode refid="26" relation="include">
        </childnode>
        <childnode refid="31" relation="include">
        </childnode>
        <childnode refid="36" relation="include">
        </childnode>
        <childnode refid="37" relation="include">
        </childnode>
        <childnode refid="38" relation="include">
        </childnode>
      </node>
      <node id="132">
        <label>serif/discretization/blocks/traits.hpp</label>
        <link refid="discretization_2blocks_2traits_8hpp"/>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="133" relation="include">
        </childnode>
      </node>
      <node id="31">
        <label>serif/discretization/domain/concepts.hpp</label>
        <link refid="discretization_2domain_2concepts_8hpp"/>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="8" relation="include">
        </childnode>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="27" relation="include">
        </childnode>
        <childnode refid="33" relation="include">
        </childnode>
        <childnode refid="34" relation="include">
        </childnode>
        <childnode refid="23" relation="include">
        </childnode>
        <childnode refid="35" relation="include">
        </childnode>
      </node>
      <node id="75">
        <label>serif/discretization/domain/ids/boundary.hpp</label>
        <link refid="ids_2boundary_8hpp"/>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="31" relation="include">
        </childnode>
      </node>
      <node id="79">
        <label>serif/discretization/domain/ids/concepts.hpp</label>
        <link refid="discretization_2domain_2ids_2concepts_8hpp"/>
        <childnode refid="75" relation="include">
        </childnode>
        <childnode refid="76" relation="include">
        </childnode>
      </node>
      <node id="76">
        <label>serif/discretization/domain/ids/domain.hpp</label>
        <link refid="domain_8hpp"/>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="31" relation="include">
        </childnode>
      </node>
      <node id="87">
        <label>serif/discretization/domain/ids/lists/concepts.hpp</label>
        <link refid="discretization_2domain_2ids_2lists_2concepts_8hpp"/>
        <childnode refid="88" relation="include">
        </childnode>
      </node>
      <node id="77">
        <label>serif/discretization/domain/ids/lists/lists.hpp</label>
        <link refid="lists_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="78" relation="include">
        </childnode>
        <childnode refid="35" relation="include">
        </childnode>
      </node>
      <node id="88">
        <label>serif/discretization/domain/ids/lists/traits.hpp</label>
        <link refid="discretization_2domain_2ids_2lists_2traits_8hpp"/>
        <childnode refid="77" relation="include">
        </childnode>
      </node>
      <node id="78">
        <label>serif/discretization/domain/ids/validation/concepts.hpp</label>
        <link refid="discretization_2domain_2ids_2validation_2concepts_8hpp"/>
        <childnode refid="79" relation="include">
        </childnode>
        <childnode refid="80" relation="include">
        </childnode>
        <childnode refid="81" relation="include">
        </childnode>
      </node>
      <node id="80">
        <label>serif/discretization/domain/ids/validation/unique_ids.hpp</label>
        <link refid="unique__ids_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="35" relation="include">
        </childnode>
      </node>
      <node id="36">
        <label>serif/discretization/domain/physical_domains.hpp</label>
        <link refid="physical__domains_8hpp"/>
        <childnode refid="34" relation="include">
        </childnode>
        <childnode refid="31" relation="include">
        </childnode>
      </node>
      <node id="84">
        <label>serif/discretization/domain/relation/concepts.hpp</label>
        <link refid="discretization_2domain_2relation_2concepts_8hpp"/>
        <childnode refid="85" relation="include">
        </childnode>
        <childnode refid="23" relation="include">
        </childnode>
      </node>
      <node id="89">
        <label>serif/discretization/domain/relation/lists/concepts.hpp</label>
        <link refid="discretization_2domain_2relation_2lists_2concepts_8hpp"/>
        <childnode refid="90" relation="include">
        </childnode>
      </node>
      <node id="83">
        <label>serif/discretization/domain/relation/lists/relation_list.hpp</label>
        <link refid="relation__list_8hpp"/>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="84" relation="include">
        </childnode>
        <childnode refid="35" relation="include">
        </childnode>
      </node>
      <node id="90">
        <label>serif/discretization/domain/relation/lists/traits.hpp</label>
        <link refid="discretization_2domain_2relation_2lists_2traits_8hpp"/>
        <childnode refid="83" relation="include">
        </childnode>
      </node>
      <node id="85">
        <label>serif/discretization/domain/relation/relations.hpp</label>
        <link refid="discretization_2domain_2relation_2relations_8hpp"/>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="31" relation="include">
        </childnode>
        <childnode refid="24" relation="include">
        </childnode>
        <childnode refid="35" relation="include">
        </childnode>
      </node>
      <node id="91">
        <label>serif/discretization/domain/relation/validation/concepts.hpp</label>
        <link refid="discretization_2domain_2relation_2validation_2concepts_8hpp"/>
        <childnode refid="92" relation="include">
        </childnode>
      </node>
      <node id="92">
        <label>serif/discretization/domain/relation/validation/validation.hpp</label>
        <link refid="discretization_2domain_2relation_2validation_2validation_8hpp"/>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="31" relation="include">
        </childnode>
        <childnode refid="87" relation="include">
        </childnode>
        <childnode refid="84" relation="include">
        </childnode>
        <childnode refid="89" relation="include">
        </childnode>
        <childnode refid="83" relation="include">
        </childnode>
        <childnode refid="93" relation="include">
        </childnode>
      </node>
      <node id="93">
        <label>serif/discretization/domain/resolver.hpp</label>
        <link refid="resolver_8hpp"/>
        <childnode refid="31" relation="include">
        </childnode>
        <childnode refid="87" relation="include">
        </childnode>
        <childnode refid="77" relation="include">
        </childnode>
      </node>
      <node id="94">
        <label>serif/discretization/domain/schema/concepts.hpp</label>
        <link refid="discretization_2domain_2schema_2concepts_8hpp"/>
        <childnode refid="95" relation="include">
        </childnode>
      </node>
      <node id="86">
        <label>serif/discretization/domain/schema/domain_schema.hpp</label>
        <link refid="domain__schema_8hpp"/>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="31" relation="include">
        </childnode>
        <childnode refid="87" relation="include">
        </childnode>
        <childnode refid="89" relation="include">
        </childnode>
        <childnode refid="91" relation="include">
        </childnode>
        <childnode refid="93" relation="include">
        </childnode>
      </node>
      <node id="74">
        <label>serif/discretization/domain/schema/schemas.hpp</label>
        <link refid="schemas_8hpp"/>
        <childnode refid="75" relation="include">
        </childnode>
        <childnode refid="76" relation="include">
        </childnode>
        <childnode refid="77" relation="include">
        </childnode>
        <childnode refid="36" relation="include">
        </childnode>
        <childnode refid="83" relation="include">
        </childnode>
        <childnode refid="85" relation="include">
        </childnode>
        <childnode refid="86" relation="include">
        </childnode>
        <childnode refid="34" relation="include">
        </childnode>
      </node>
      <node id="95">
        <label>serif/discretization/domain/schema/traits.hpp</label>
        <link refid="discretization_2domain_2schema_2traits_8hpp"/>
        <childnode refid="86" relation="include">
        </childnode>
      </node>
      <node id="34">
        <label>serif/discretization/domain/types.hpp</label>
        <link refid="discretization_2domain_2types_8hpp"/>
        <childnode refid="17" relation="include">
        </childnode>
      </node>
      <node id="134">
        <label>serif/discretization/forms/operands.hpp</label>
        <link refid="operands_8hpp"/>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="133" relation="include">
        </childnode>
      </node>
      <node id="131">
        <label>serif/discretization/forms/value.hpp</label>
        <link refid="value_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="132" relation="include">
        </childnode>
        <childnode refid="134" relation="include">
        </childnode>
      </node>
      <node id="26">
        <label>serif/discretization/functional_space/space.hpp</label>
        <link refid="space_8hpp"/>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="27" relation="include">
        </childnode>
        <childnode refid="28" relation="include">
        </childnode>
      </node>
      <node id="40">
        <label>serif/discretization/reconstruction/surface_power.hpp</label>
        <link refid="surface__power_8hpp"/>
        <childnode refid="4" relation="include">
        </childnode>
        <childnode refid="8" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="36" relation="include">
        </childnode>
        <childnode refid="41" relation="include">
        </childnode>
      </node>
      <node id="123">
        <label>serif/eos/evaluation.hpp</label>
        <link refid="eos_2evaluation_8hpp"/>
        <childnode refid="124" relation="include">
        </childnode>
      </node>
      <node id="128">
        <label>serif/eos/exceptions.hpp</label>
        <link refid="eos_2exceptions_8hpp"/>
        <childnode refid="43" relation="include">
        </childnode>
        <childnode refid="55" relation="include">
        </childnode>
        <childnode refid="9" relation="include">
        </childnode>
      </node>
      <node id="124">
        <label>serif/eos/models/concepts.hpp</label>
        <link refid="eos_2models_2concepts_8hpp"/>
        <childnode refid="125" relation="include">
        </childnode>
        <childnode refid="118" relation="include">
        </childnode>
      </node>
      <node id="117">
        <label>serif/eos/models/polytropic.hpp</label>
        <link refid="eos_2models_2polytropic_8hpp"/>
        <childnode refid="118" relation="include">
        </childnode>
        <childnode refid="120" relation="include">
        </childnode>
        <childnode refid="121" relation="include">
        </childnode>
      </node>
      <node id="125">
        <label>serif/eos/models/traits.hpp</label>
        <link refid="eos_2models_2traits_8hpp"/>
        <childnode refid="19" relation="include">
        </childnode>
        <childnode refid="121" relation="include">
        </childnode>
        <childnode refid="126" relation="include">
        </childnode>
      </node>
      <node id="119">
        <label>serif/eos/relations/catalog.hpp</label>
        <link refid="catalog_8hpp"/>
        <childnode refid="81" relation="include">
        </childnode>
        <childnode refid="35" relation="include">
        </childnode>
        <childnode refid="19" relation="include">
        </childnode>
      </node>
      <node id="126">
        <label>serif/eos/relations/concepts.hpp</label>
        <link refid="eos_2relations_2concepts_8hpp"/>
        <childnode refid="127" relation="include">
        </childnode>
      </node>
      <node id="118">
        <label>serif/eos/relations/relations.hpp</label>
        <link refid="eos_2relations_2relations_8hpp"/>
        <childnode refid="119" relation="include">
        </childnode>
        <childnode refid="18" relation="include">
        </childnode>
      </node>
      <node id="127">
        <label>serif/eos/relations/traits.hpp</label>
        <link refid="eos_2relations_2traits_8hpp"/>
        <childnode refid="119" relation="include">
        </childnode>
      </node>
      <node id="122">
        <label>serif/equations/barotropic_density.hpp</label>
        <link refid="barotropic__density_8hpp"/>
        <childnode refid="4" relation="include">
        </childnode>
        <childnode refid="12" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="120" relation="include">
        </childnode>
        <childnode refid="16" relation="include">
        </childnode>
        <childnode refid="123" relation="include">
        </childnode>
        <childnode refid="128" relation="include">
        </childnode>
        <childnode refid="118" relation="include">
        </childnode>
        <childnode refid="129" relation="include">
        </childnode>
        <childnode refid="135" relation="include">
        </childnode>
        <childnode refid="136" relation="include">
        </childnode>
        <childnode refid="130" relation="include">
        </childnode>
        <childnode refid="41" relation="include">
        </childnode>
      </node>
      <node id="129">
        <label>serif/equations/companion.hpp</label>
        <link refid="companion_8hpp"/>
        <childnode refid="4" relation="include">
        </childnode>
        <childnode refid="130" relation="include">
        </childnode>
        <childnode refid="41" relation="include">
        </childnode>
      </node>
      <node id="168">
        <label>serif/equations/concepts.hpp</label>
        <link refid="equations_2concepts_8hpp"/>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="135" relation="include">
        </childnode>
        <childnode refid="136" relation="include">
        </childnode>
      </node>
      <node id="135">
        <label>serif/equations/contributions.hpp</label>
        <link refid="equations_2contributions_8hpp"/>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="132" relation="include">
        </childnode>
        <childnode refid="134" relation="include">
        </childnode>
      </node>
      <node id="136">
        <label>serif/equations/point.hpp</label>
        <link refid="point_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="131" relation="include">
        </childnode>
        <childnode refid="37" relation="include">
        </childnode>
      </node>
      <node id="137">
        <label>serif/equations/primal_gravity.hpp</label>
        <link refid="primal__gravity_8hpp"/>
        <childnode refid="4" relation="include">
        </childnode>
        <childnode refid="55" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="16" relation="include">
        </childnode>
        <childnode refid="135" relation="include">
        </childnode>
        <childnode refid="136" relation="include">
        </childnode>
        <childnode refid="41" relation="include">
        </childnode>
      </node>
      <node id="130">
        <label>serif/equations/rigid_rotation.hpp</label>
        <link refid="rigid__rotation_8hpp"/>
        <childnode refid="4" relation="include">
        </childnode>
        <childnode refid="55" relation="include">
        </childnode>
        <childnode refid="131" relation="include">
        </childnode>
        <childnode refid="41" relation="include">
        </childnode>
      </node>
      <node id="99">
        <label>serif/geometry/deformation/actions.hpp</label>
        <link refid="geometry_2deformation_2actions_8hpp"/>
        <childnode refid="44" relation="include">
        </childnode>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="8" relation="include">
        </childnode>
        <childnode refid="64" relation="include">
        </childnode>
        <childnode refid="69" relation="include">
        </childnode>
        <childnode refid="67" relation="include">
        </childnode>
        <childnode refid="68" relation="include">
        </childnode>
      </node>
      <node id="98">
        <label>serif/geometry/deformation/composition.hpp</label>
        <link refid="composition_8hpp"/>
        <childnode refid="6" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
        <childnode refid="99" relation="include">
        </childnode>
      </node>
      <node id="97">
        <label>serif/geometry/deformation/computations.hpp</label>
        <link refid="geometry_2deformation_2computations_8hpp"/>
        <childnode refid="8" relation="include">
        </childnode>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="98" relation="include">
        </childnode>
        <childnode refid="100" relation="include">
        </childnode>
        <childnode refid="101" relation="include">
        </childnode>
      </node>
      <node id="69">
        <label>serif/geometry/deformation/detail/action_storage.hpp</label>
        <link refid="action__storage_8hpp"/>
        <childnode refid="44" relation="include">
        </childnode>
        <childnode refid="4" relation="include">
        </childnode>
        <childnode refid="70" relation="include">
        </childnode>
        <childnode refid="71" relation="include">
        </childnode>
        <childnode refid="8" relation="include">
        </childnode>
        <childnode refid="64" relation="include">
        </childnode>
      </node>
      <node id="64">
        <label>serif/geometry/deformation/failure.hpp</label>
        <link refid="geometry_2deformation_2failure_8hpp"/>
        <childnode refid="65" relation="include">
        </childnode>
        <childnode refid="9" relation="include">
        </childnode>
      </node>
      <node id="100">
        <label>serif/geometry/deformation/linearization.hpp</label>
        <link refid="deformation_2linearization_8hpp"/>
        <childnode refid="6" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
        <childnode refid="99" relation="include">
        </childnode>
      </node>
      <node id="66">
        <label>serif/geometry/deformation/models.hpp</label>
        <link refid="models_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="12" relation="include">
        </childnode>
      </node>
      <node id="62">
        <label>serif/geometry/deformation/radial_extension.hpp</label>
        <link refid="radial__extension_8hpp"/>
        <childnode refid="6" relation="include">
        </childnode>
        <childnode refid="8" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
        <childnode refid="63" relation="include">
        </childnode>
        <childnode refid="66" relation="include">
        </childnode>
        <childnode refid="67" relation="include">
        </childnode>
        <childnode refid="68" relation="include">
        </childnode>
      </node>
      <node id="63">
        <label>serif/geometry/deformation/radial_interpolation.hpp</label>
        <link refid="geometry_2deformation_2radial__interpolation_8hpp"/>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
        <childnode refid="64" relation="include">
        </childnode>
      </node>
      <node id="114">
        <label>serif/geometry/deformation/reference_surface.hpp</label>
        <link refid="reference__surface_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
      </node>
      <node id="67">
        <label>serif/geometry/deformation/spaces.hpp</label>
        <link refid="spaces_8hpp"/>
      </node>
      <node id="115">
        <label>serif/geometry/deformation/surface_embedding.hpp</label>
        <link refid="surface__embedding_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="65" relation="include">
        </childnode>
        <childnode refid="8" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
        <childnode refid="64" relation="include">
        </childnode>
        <childnode refid="66" relation="include">
        </childnode>
        <childnode refid="67" relation="include">
        </childnode>
        <childnode refid="68" relation="include">
        </childnode>
      </node>
      <node id="68">
        <label>serif/geometry/deformation/workspace.hpp</label>
        <link refid="geometry_2deformation_2workspace_8hpp"/>
      </node>
      <node id="138">
        <label>serif/geometry/surface/isobaric.hpp</label>
        <link refid="isobaric_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="18" relation="include">
        </childnode>
        <childnode refid="120" relation="include">
        </childnode>
      </node>
      <node id="28">
        <label>serif/meta/enums.hpp</label>
        <link refid="enums_8hpp"/>
        <childnode refid="29" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="30" relation="include">
        </childnode>
      </node>
      <node id="38">
        <label>serif/meta/names.hpp</label>
        <link refid="names_8hpp"/>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="39" relation="include">
        </childnode>
      </node>
      <node id="139">
        <label>serif/model/description.hpp</label>
        <link refid="description_8hpp"/>
        <childnode refid="140" relation="include">
        </childnode>
        <childnode refid="141" relation="include">
        </childnode>
        <childnode refid="143" relation="include">
        </childnode>
        <childnode refid="145" relation="include">
        </childnode>
        <childnode refid="146" relation="include">
        </childnode>
      </node>
      <node id="147">
        <label>serif/model/detail/specification_list.hpp</label>
        <link refid="specification__list_8hpp"/>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="141" relation="include">
        </childnode>
        <childnode refid="148" relation="include">
        </childnode>
      </node>
      <node id="140">
        <label>serif/model/equilibrium_equations.hpp</label>
        <link refid="equilibrium__equations_8hpp"/>
        <childnode refid="70" relation="include">
        </childnode>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
      </node>
      <node id="141">
        <label>serif/model/specification.hpp</label>
        <link refid="specification_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="142" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="19" relation="include">
        </childnode>
        <childnode refid="53" relation="include">
        </childnode>
        <childnode refid="37" relation="include">
        </childnode>
      </node>
      <node id="143">
        <label>serif/model/specifications/fixed_angular_momentum.hpp</label>
        <link refid="fixed__angular__momentum_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="120" relation="include">
        </childnode>
        <childnode refid="141" relation="include">
        </childnode>
        <childnode refid="144" relation="include">
        </childnode>
      </node>
      <node id="145">
        <label>serif/model/specifications/fixed_total_mass.hpp</label>
        <link refid="fixed__total__mass_8hpp"/>
        <childnode refid="120" relation="include">
        </childnode>
        <childnode refid="141" relation="include">
        </childnode>
        <childnode refid="144" relation="include">
        </childnode>
      </node>
      <node id="144">
        <label>serif/model/specifications/physical_names.hpp</label>
        <link refid="physical__names_8hpp"/>
      </node>
      <node id="146">
        <label>serif/model/stellar_model.hpp</label>
        <link refid="stellar__model_8hpp"/>
        <childnode refid="70" relation="include">
        </childnode>
        <childnode refid="8" relation="include">
        </childnode>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="147" relation="include">
        </childnode>
        <childnode refid="140" relation="include">
        </childnode>
      </node>
      <node id="167">
        <label>serif/operators/concepts.hpp</label>
        <link refid="operators_2concepts_8hpp"/>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="132" relation="include">
        </childnode>
        <childnode refid="26" relation="include">
        </childnode>
        <childnode refid="168" relation="include">
        </childnode>
      </node>
      <node id="104">
        <label>serif/utils/dependencies/binding.hpp</label>
        <link refid="utils_2dependencies_2binding_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="82" relation="include">
        </childnode>
      </node>
      <node id="102">
        <label>serif/utils/dependencies/computation.hpp</label>
        <link refid="computation_8hpp"/>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="70" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="103" relation="include">
        </childnode>
      </node>
      <node id="149">
        <label>serif/utils/dependencies/computation_cache.hpp</label>
        <link refid="computation__cache_8hpp"/>
        <childnode refid="41" relation="include">
        </childnode>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="43" relation="include">
        </childnode>
        <childnode refid="65" relation="include">
        </childnode>
        <childnode refid="6" relation="include">
        </childnode>
        <childnode refid="55" relation="include">
        </childnode>
        <childnode refid="9" relation="include">
        </childnode>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="150" relation="include">
        </childnode>
        <childnode refid="151" relation="include">
        </childnode>
        <childnode refid="154" relation="include">
        </childnode>
      </node>
      <node id="101">
        <label>serif/utils/dependencies/computation_graph.hpp</label>
        <link refid="computation__graph_8hpp"/>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="102" relation="include">
        </childnode>
        <childnode refid="105" relation="include">
        </childnode>
        <childnode refid="107" relation="include">
        </childnode>
      </node>
      <node id="108">
        <label>serif/utils/dependencies/dependency_context.hpp</label>
        <link refid="dependency__context_8hpp"/>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="109" relation="include">
        </childnode>
        <childnode refid="106" relation="include">
        </childnode>
        <childnode refid="111" relation="include">
        </childnode>
      </node>
      <node id="154">
        <label>serif/utils/dependencies/detail/cache_storage.hpp</label>
        <link refid="cache__storage_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="43" relation="include">
        </childnode>
        <childnode refid="6" relation="include">
        </childnode>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="155" relation="include">
        </childnode>
      </node>
      <node id="155">
        <label>serif/utils/dependencies/detail/cache_traits.hpp</label>
        <link refid="cache__traits_8hpp"/>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="156" relation="include">
        </childnode>
        <childnode refid="152" relation="include">
        </childnode>
        <childnode refid="148" relation="include">
        </childnode>
      </node>
      <node id="107">
        <label>serif/utils/dependencies/detail/graph.hpp</label>
        <link refid="utils_2dependencies_2detail_2graph_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="106" relation="include">
        </childnode>
      </node>
      <node id="153">
        <label>serif/utils/dependencies/detail/result.hpp</label>
        <link refid="utils_2dependencies_2detail_2result_8hpp"/>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="65" relation="include">
        </childnode>
        <childnode refid="6" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="37" relation="include">
        </childnode>
      </node>
      <node id="152">
        <label>serif/utils/dependencies/detail/schema.hpp</label>
        <link refid="utils_2dependencies_2detail_2schema_8hpp"/>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="108" relation="include">
        </childnode>
        <childnode refid="101" relation="include">
        </childnode>
        <childnode refid="153" relation="include">
        </childnode>
      </node>
      <node id="106">
        <label>serif/utils/dependencies/detail/type_lists.hpp</label>
        <link refid="type__lists_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="103" relation="include">
        </childnode>
      </node>
      <node id="103">
        <label>serif/utils/dependencies/key_list.hpp</label>
        <link refid="key__list_8hpp"/>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="104" relation="include">
        </childnode>
        <childnode refid="82" relation="include">
        </childnode>
        <childnode refid="37" relation="include">
        </childnode>
      </node>
      <node id="105">
        <label>serif/utils/dependencies/key_list_operations.hpp</label>
        <link refid="key__list__operations_8hpp"/>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="106" relation="include">
        </childnode>
      </node>
      <node id="151">
        <label>serif/utils/dependencies/observers.hpp</label>
        <link refid="utils_2dependencies_2observers_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="51" relation="include">
        </childnode>
        <childnode refid="70" relation="include">
        </childnode>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="27" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="152" relation="include">
        </childnode>
      </node>
      <node id="150">
        <label>serif/utils/dependencies/preparation.hpp</label>
        <link refid="utils_2dependencies_2preparation_8hpp"/>
        <childnode refid="44" relation="include">
        </childnode>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="29" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="106" relation="include">
        </childnode>
        <childnode refid="28" relation="include">
        </childnode>
      </node>
      <node id="109">
        <label>serif/utils/dependencies/resource.hpp</label>
        <link refid="resource_8hpp"/>
        <childnode refid="70" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="110" relation="include">
        </childnode>
      </node>
      <node id="110">
        <label>serif/utils/dependencies/snapshot.hpp</label>
        <link refid="snapshot_8hpp"/>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="70" relation="include">
        </childnode>
        <childnode refid="6" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
      </node>
      <node id="52">
        <label>serif/utils/error/codes.hpp</label>
        <link refid="codes_8hpp"/>
        <childnode refid="43" relation="include">
        </childnode>
        <childnode refid="28" relation="include">
        </childnode>
      </node>
      <node id="41">
        <label>serif/utils/error/errors.hpp</label>
        <link refid="errors_8hpp"/>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="27" relation="include">
        </childnode>
        <childnode refid="42" relation="include">
        </childnode>
        <childnode refid="9" relation="include">
        </childnode>
        <childnode refid="43" relation="include">
        </childnode>
        <childnode refid="44" relation="include">
        </childnode>
        <childnode refid="45" relation="include">
        </childnode>
        <childnode refid="46" relation="include">
        </childnode>
        <childnode refid="50" relation="include">
        </childnode>
        <childnode refid="52" relation="include">
        </childnode>
        <childnode refid="53" relation="include">
        </childnode>
      </node>
      <node id="58">
        <label>serif/utils/error/nonfinite.hpp</label>
        <link refid="nonfinite_8hpp"/>
        <childnode refid="55" relation="include">
        </childnode>
      </node>
      <node id="179">
        <label>serif/utils/misc/compensated_sum.hpp</label>
        <link refid="compensated__sum_8hpp"/>
        <childnode refid="4" relation="include">
        </childnode>
      </node>
      <node id="111">
        <label>serif/utils/misc/concepts/concat.hpp</label>
        <link refid="concat_8hpp"/>
        <childnode refid="37" relation="include">
        </childnode>
      </node>
      <node id="50">
        <label>serif/utils/misc/concepts/exceptions.hpp</label>
        <link refid="utils_2misc_2concepts_2exceptions_8hpp"/>
        <childnode refid="51" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="24" relation="include">
        </childnode>
      </node>
      <node id="148">
        <label>serif/utils/misc/concepts/filter.hpp</label>
        <link refid="filter_8hpp"/>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="37" relation="include">
        </childnode>
      </node>
      <node id="24">
        <label>serif/utils/misc/concepts/numeric.hpp</label>
        <link refid="numeric_8hpp"/>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
      </node>
      <node id="81">
        <label>serif/utils/misc/concepts/type_uniqueness.hpp</label>
        <link refid="concepts_2type__uniqueness_8hpp"/>
        <childnode refid="82" relation="include">
        </childnode>
      </node>
      <node id="23">
        <label>serif/utils/misc/std_helper/cleaning.hpp</label>
        <link refid="cleaning_8hpp"/>
        <childnode refid="10" relation="include">
        </childnode>
      </node>
      <node id="35">
        <label>serif/utils/misc/std_helper/variadic_cleaning.hpp</label>
        <link refid="variadic__cleaning_8hpp"/>
        <childnode refid="32" relation="include">
        </childnode>
      </node>
      <node id="46">
        <label>serif/utils/misc/terminal/colors.hpp</label>
        <link refid="colors_8hpp"/>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="42" relation="include">
        </childnode>
        <childnode refid="47" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="49" relation="include">
        </childnode>
      </node>
      <node id="82">
        <label>serif/utils/misc/traits/type_uniqueness.hpp</label>
        <link refid="traits_2type__uniqueness_8hpp"/>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="23" relation="include">
        </childnode>
      </node>
      <node id="158">
        <label>serif/utils/serialize/record.hpp</label>
        <link refid="record_8hpp"/>
        <childnode refid="44" relation="include">
        </childnode>
        <childnode refid="3" relation="include">
        </childnode>
        <childnode refid="159" relation="include">
        </childnode>
        <childnode refid="21" relation="include">
        </childnode>
        <childnode refid="43" relation="include">
        </childnode>
        <childnode refid="51" relation="include">
        </childnode>
        <childnode refid="65" relation="include">
        </childnode>
        <childnode refid="70" relation="include">
        </childnode>
        <childnode refid="142" relation="include">
        </childnode>
        <childnode refid="160" relation="include">
        </childnode>
        <childnode refid="12" relation="include">
        </childnode>
        <childnode refid="9" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="10" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="156" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
      </node>
      <node id="157">
        <label>serif/utils/serialize/top.hpp</label>
        <link refid="utils_2serialize_2top_8hpp"/>
        <childnode refid="32" relation="include">
        </childnode>
        <childnode refid="5" relation="include">
        </childnode>
        <childnode refid="70" relation="include">
        </childnode>
        <childnode refid="158" relation="include">
        </childnode>
      </node>
      <node id="53">
        <label>serif/utils/types/strings.hpp</label>
        <link refid="types_2strings_8hpp"/>
        <childnode refid="44" relation="include">
        </childnode>
        <childnode refid="54" relation="include">
        </childnode>
        <childnode refid="17" relation="include">
        </childnode>
      </node>
      <node id="37">
        <label>serif/utils/types/type_list.hpp</label>
        <link refid="type__list_8hpp"/>
      </node>
      <node id="1">
        <label>src/lib/serif/backend/mfem/serialize/top.cpp</label>
        <link refid="backend_2mfem_2serialize_2top_8cpp"/>
        <childnode refid="2" relation="include">
        </childnode>
        <childnode refid="44" relation="include">
        </childnode>
        <childnode refid="142" relation="include">
        </childnode>
        <childnode refid="160" relation="include">
        </childnode>
        <childnode refid="12" relation="include">
        </childnode>
        <childnode refid="22" relation="include">
        </childnode>
        <childnode refid="55" relation="include">
        </childnode>
        <childnode refid="48" relation="include">
        </childnode>
        <childnode refid="14" relation="include">
        </childnode>
        <childnode refid="161" relation="include">
        </childnode>
        <childnode refid="162" relation="include">
        </childnode>
      </node>
      <node id="44">
        <label>algorithm</label>
      </node>
      <node id="3">
        <label>array</label>
      </node>
      <node id="159">
        <label>bit</label>
      </node>
      <node id="4">
        <label>cmath</label>
      </node>
      <node id="20">
        <label>compare</label>
      </node>
      <node id="21">
        <label>concepts</label>
      </node>
      <node id="32">
        <label>cstddef</label>
      </node>
      <node id="43">
        <label>cstdint</label>
      </node>
      <node id="42">
        <label>cstdio</label>
      </node>
      <node id="47">
        <label>cstdlib</label>
      </node>
      <node id="51">
        <label>exception</label>
      </node>
      <node id="65">
        <label>expected</label>
      </node>
      <node id="5">
        <label>filesystem</label>
      </node>
      <node id="27">
        <label>format</label>
      </node>
      <node id="70">
        <label>functional</label>
      </node>
      <node id="71">
        <label>initializer_list</label>
      </node>
      <node id="54">
        <label>iostream</label>
      </node>
      <node id="142">
        <label>limits</label>
      </node>
      <node id="30">
        <label>magic_enum/magic_enum_all.hpp</label>
      </node>
      <node id="160">
        <label>map</label>
      </node>
      <node id="6">
        <label>memory</label>
      </node>
      <node id="166">
        <label>mfem/fem/eltrans.hpp</label>
      </node>
      <node id="15">
        <label>mfem/fem/fe/fe_base.hpp</label>
      </node>
      <node id="172">
        <label>mfem/fem/fe/fe_h1.hpp</label>
      </node>
      <node id="173">
        <label>mfem/fem/fe/fe_l2.hpp</label>
      </node>
      <node id="56">
        <label>mfem/fem/fespace.hpp</label>
      </node>
      <node id="161">
        <label>mfem/fem/gridfunc.hpp</label>
      </node>
      <node id="170">
        <label>mfem/fem/intrules.hpp</label>
      </node>
      <node id="57">
        <label>mfem/fem/pfespace.hpp</label>
      </node>
      <node id="176">
        <label>mfem/linalg/densemat.hpp</label>
      </node>
      <node id="61">
        <label>mfem/linalg/sparsemat.hpp</label>
      </node>
      <node id="177">
        <label>mfem/linalg/vector.hpp</label>
      </node>
      <node id="39">
        <label>nameof.hpp</label>
      </node>
      <node id="7">
        <label>numbers</label>
      </node>
      <node id="33">
        <label>numeric</label>
      </node>
      <node id="12">
        <label>optional</label>
      </node>
      <node id="73">
        <label>ranges</label>
      </node>
      <node id="8">
        <label>span</label>
      </node>
      <node id="55">
        <label>stdexcept</label>
      </node>
      <node id="9">
        <label>string</label>
      </node>
      <node id="17">
        <label>string_view</label>
      </node>
      <node id="22">
        <label>tuple</label>
      </node>
      <node id="10">
        <label>type_traits</label>
      </node>
      <node id="49">
        <label>unistd.h</label>
      </node>
      <node id="29">
        <label>unordered_map</label>
      </node>
      <node id="48">
        <label>utility</label>
      </node>
      <node id="156">
        <label>variant</label>
      </node>
      <node id="14">
        <label>vector</label>
      </node>
      <node id="45">
        <label>version</label>
      </node>
    </incdepgraph>
    <innerclass refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl" prot="public">serif::backend::mfem::serialize::detail::TopModelSampler::Impl</innerclass>
    <innerclass refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_error" prot="public">serif::backend::mfem::serialize::detail::Error</innerclass>
    <innernamespace refid="namespaceserif">serif</innernamespace>
    <innernamespace refid="namespaceserif_1_1backend">serif::backend</innernamespace>
    <innernamespace refid="namespaceserif_1_1backend_1_1mfem">serif::backend::mfem</innernamespace>
    <innernamespace refid="namespaceserif_1_1backend_1_1mfem_1_1serialize">serif::backend::mfem::serialize</innernamespace>
    <innernamespace refid="namespaceserif_1_1backend_1_1mfem_1_1serialize_1_1detail">serif::backend::mfem::serialize::detail</innernamespace>
    <briefdescription>
    </briefdescription>
    <detaileddescription>
    </detaileddescription>
    <programlisting>
<codeline lineno="1"><highlight class="preprocessor">#include<sp/>&quot;<ref refid="backend_2mfem_2serialize_2top_8hpp" kindref="compound">serif/backend/mfem/serialize/top.hpp</ref>&quot;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="2"><highlight class="normal"></highlight></codeline>
<codeline lineno="3"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&lt;algorithm&gt;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="4"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&lt;limits&gt;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="5"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&lt;map&gt;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="6"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&lt;optional&gt;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="7"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&lt;tuple&gt;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="8"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&lt;stdexcept&gt;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="9"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&lt;utility&gt;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="10"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&lt;vector&gt;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="11"><highlight class="normal"></highlight></codeline>
<codeline lineno="12"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&quot;mfem/fem/gridfunc.hpp&quot;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="13"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>&quot;<ref refid="field__values_8hpp" kindref="compound">serif/backend/mfem/discretization/field_values.hpp</ref>&quot;</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="14"><highlight class="normal"></highlight></codeline>
<codeline lineno="15"><highlight class="normal"></highlight><highlight class="keyword">namespace<sp/></highlight><highlight class="normal"><ref refid="namespaceserif_1_1backend_1_1mfem_1_1serialize_1_1detail" kindref="compound">serif::backend::mfem::serialize::detail</ref><sp/>{</highlight></codeline>
<codeline lineno="16"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">namespace<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="17"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/></highlight><highlight class="normal">Vec<sp/>=<sp/>std::array&lt;double,<sp/>3&gt;;</highlight></codeline>
<codeline lineno="18"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/></highlight><highlight class="normal">utils::serialize::Error;</highlight></codeline>
<codeline lineno="19"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/></highlight><highlight class="normal"><ref refid="namespaceserif_1_1utils_1_1serialize_1ace9e0af29889a767ef1c95a09a4b7265" kindref="member">utils::serialize::Result</ref>;</highlight></codeline>
<codeline lineno="20"><highlight class="normal"></highlight></codeline>
<codeline lineno="22"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">class<sp/></highlight><highlight class="normal">DensityBinding<sp/>final<sp/>:<sp/></highlight><highlight class="keyword">public</highlight><highlight class="normal"><sp/>discretization::backend::mfem::DofMap<sp/>{</highlight></codeline>
<codeline lineno="23"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">public</highlight><highlight class="normal">:</highlight></codeline>
<codeline lineno="25"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">explicit</highlight><highlight class="normal"><sp/><ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1abde2910012b0374e8c072cb30981371b" kindref="member">DensityBinding</ref>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>DofMap&amp;<sp/>source)<sp/>:</highlight></codeline>
<codeline lineno="26"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>DofMap(source.space(),<sp/>support(source),<sp/>source.reconstruction())<sp/>{}</highlight></codeline>
<codeline lineno="27"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">private</highlight><highlight class="normal">:</highlight></codeline>
<codeline lineno="29"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>static<sp/>::mfem::Array&lt;int&gt;<sp/>support(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>DofMap&amp;<sp/>source)<sp/>{</highlight></codeline>
<codeline lineno="30"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>source.validate();</highlight></codeline>
<codeline lineno="31"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::Array&lt;int&gt;<sp/>result(source.mesh().GetNE());</highlight></codeline>
<codeline lineno="32"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>element<sp/>=<sp/>0;<sp/>element<sp/>&lt;<sp/>result.Size();<sp/>++element)<sp/>result[element]<sp/>=<sp/>source.supports(element);</highlight></codeline>
<codeline lineno="33"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>result;</highlight></codeline>
<codeline lineno="34"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="35"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="36"><highlight class="normal"></highlight></codeline>
<codeline lineno="38"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1backend_1_1mfem_1afdd4d9966e5040b2de0684e018dce66e" kindref="member">::mfem::Vector</ref><sp/>copy_density(std::span&lt;const<sp/>double&gt;<sp/>values)<sp/>{</highlight></codeline>
<codeline lineno="39"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1backend_1_1mfem_1afdd4d9966e5040b2de0684e018dce66e" kindref="member">::mfem::Vector</ref><sp/>result(</highlight><highlight class="keyword">static_cast&lt;</highlight><highlight class="keywordtype">int</highlight><highlight class="keyword">&gt;</highlight><highlight class="normal">(<ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1adf1d7a29193c8b8a6cd17745686a04d8" kindref="member">values</ref>.size()));</highlight></codeline>
<codeline lineno="40"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::copy(<ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1adf1d7a29193c8b8a6cd17745686a04d8" kindref="member">values</ref>.begin(),<sp/><ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1adf1d7a29193c8b8a6cd17745686a04d8" kindref="member">values</ref>.end(),<sp/>result.GetData());</highlight></codeline>
<codeline lineno="41"><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/>result;</highlight></codeline>
<codeline lineno="42"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="43"><highlight class="normal"></highlight></codeline>
<codeline lineno="44"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>dot(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>Vec&amp;<sp/>a,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>Vec&amp;<sp/>b)<sp/>{</highlight></codeline>
<codeline lineno="45"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>a[0]<sp/>*<sp/>b[0]<sp/>+<sp/>a[1]<sp/>*<sp/>b[1]<sp/>+<sp/>a[2]<sp/>*<sp/>b[2];</highlight></codeline>
<codeline lineno="46"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="47"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>Vec<sp/>cross(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>Vec&amp;<sp/>a,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>Vec&amp;<sp/>b)<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/>{a[1]<sp/>*<sp/>b[2]<sp/>-<sp/>a[2]<sp/>*<sp/>b[1],<sp/>a[2]<sp/>*<sp/>b[0]<sp/>-<sp/>a[0]<sp/>*<sp/>b[2],<sp/>a[0]<sp/>*<sp/>b[1]<sp/>-<sp/>a[1]<sp/>*<sp/>b[0]};</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/>Vec<sp/>normalized(Vec<sp/>value)<sp/>{</highlight></codeline>
<codeline lineno="51"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>norm<sp/>=<sp/>std::sqrt(dot(value,<sp/>value));</highlight></codeline>
<codeline lineno="52"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!std::isfinite(norm)<sp/>||<sp/>norm<sp/>&lt;=<sp/>0.0)<sp/>{<sp/></highlight><highlight class="keywordflow">throw</highlight><highlight class="normal"><sp/>std::invalid_argument(</highlight><highlight class="stringliteral">&quot;TOP<sp/>sampling<sp/>requires<sp/>a<sp/>finite<sp/>nonzero<sp/>direction.&quot;</highlight><highlight class="normal">);<sp/>}</highlight></codeline>
<codeline lineno="53"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">double</highlight><highlight class="normal">&amp;<sp/>component<sp/>:<sp/>value)<sp/>{<sp/><ref refid="namespaceserif_1_1discretization_1_1forms_1aa4daa912e2157d80dc731a6ce2349dfe" kindref="member">component</ref><sp/>/=<sp/>norm;<sp/>}</highlight></codeline>
<codeline lineno="54"><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_1discretization_1_1backend_1_1mfem_1af75f98db90618153ea4ca169945733f7a2063c1608d6e0baf80249c42e2be5804" kindref="member">value</ref>;</highlight></codeline>
<codeline lineno="55"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="56"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::pair&lt;Vec,<sp/>Vec&gt;<sp/>transverse_basis(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>Vec&amp;<sp/>axis)<sp/>{</highlight></codeline>
<codeline lineno="57"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>helper<sp/>=<sp/>std::abs(axis[2])<sp/>&lt;<sp/>0.8<sp/>?<sp/>Vec{0.0,<sp/>0.0,<sp/>1.0}<sp/>:<sp/>Vec{0.0,<sp/>1.0,<sp/>0.0};</highlight></codeline>
<codeline lineno="58"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>Vec<sp/>first<sp/>=<sp/>normalized(cross(helper,<sp/>axis));</highlight></codeline>
<codeline lineno="59"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>{first,<sp/>cross(axis,<sp/>first)};</highlight></codeline>
<codeline lineno="60"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="61"><highlight class="normal"></highlight></codeline>
<codeline lineno="62"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">struct<sp/></highlight><highlight class="normal">Bounds<sp/>{</highlight></codeline>
<codeline lineno="63"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>Vec<sp/>lower,<sp/>upper;</highlight></codeline>
<codeline lineno="64"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">bool</highlight><highlight class="normal"><sp/>contains(const<sp/>::mfem::Vector&amp;<sp/>point)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="65"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>k<sp/>=<sp/>0;<sp/>k<sp/>&lt;<sp/>3;<sp/>++k)<sp/>{</highlight></codeline>
<codeline lineno="66"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(point[k]<sp/>&lt;<sp/>lower[k]<sp/>||<sp/>point[k]<sp/>&gt;<sp/>upper[k])<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">false</highlight><highlight class="normal">;<sp/>}</highlight></codeline>
<codeline lineno="67"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="68"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">true</highlight><highlight class="normal">;</highlight></codeline>
<codeline lineno="69"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="70"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">bool</highlight><highlight class="normal"><sp/>intersects_ray(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>Vec&amp;<sp/>center,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>Vec&amp;<sp/>direction)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="71"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>near<sp/>=<sp/>0.0,<sp/>far<sp/>=<sp/>std::numeric_limits&lt;double&gt;::infinity();</highlight></codeline>
<codeline lineno="72"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>k<sp/>=<sp/>0;<sp/>k<sp/>&lt;<sp/>3;<sp/>++k)<sp/>{</highlight></codeline>
<codeline lineno="73"><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">if</highlight><highlight class="normal"><sp/>(std::abs(direction[k])<sp/>&lt;<sp/>1.0e-15)<sp/>{</highlight></codeline>
<codeline lineno="74"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(center[k]<sp/>&lt;<sp/>lower[k]<sp/>||<sp/>center[k]<sp/>&gt;<sp/>upper[k])<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">false</highlight><highlight class="normal">;<sp/>}</highlight></codeline>
<codeline lineno="75"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}<sp/></highlight><highlight class="keywordflow">else</highlight><highlight class="normal"><sp/>{</highlight></codeline>
<codeline lineno="76"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>a<sp/>=<sp/>(lower[k]<sp/>-<sp/>center[k])<sp/>/<sp/>direction[k];</highlight></codeline>
<codeline lineno="77"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>b<sp/>=<sp/>(upper[k]<sp/>-<sp/>center[k])<sp/>/<sp/>direction[k];</highlight></codeline>
<codeline lineno="78"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(a<sp/>&gt;<sp/>b)<sp/>{<sp/>std::swap(a,<sp/>b);<sp/>}</highlight></codeline>
<codeline lineno="79"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>near<sp/>=<sp/>std::max(near,<sp/>a);</highlight></codeline>
<codeline lineno="80"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>far<sp/>=<sp/>std::min(far,<sp/>b);</highlight></codeline>
<codeline lineno="81"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(near<sp/>&gt;<sp/>far)<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">false</highlight><highlight class="normal">;<sp/>}</highlight></codeline>
<codeline lineno="82"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="83"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="84"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>far<sp/>&gt;<sp/>0.0;</highlight></codeline>
<codeline lineno="85"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="86"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="87"><highlight class="normal"></highlight></codeline>
<codeline lineno="88"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>These<sp/>boxes<sp/>accelerate<sp/>the<sp/>search,<sp/>but<sp/>are<sp/>not<sp/>treated<sp/>as<sp/>rigorous<sp/>polynomial<sp/>bounds.</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="89"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Every<sp/>failed<sp/>accelerated<sp/>search<sp/>retries<sp/>all<sp/>elements,<sp/>so<sp/>curved<sp/>overshoot<sp/>cannot<sp/>cause<sp/>a<sp/>false<sp/>exclusion.</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="90"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>Bounds<sp/>bounds(::mfem::ElementTransformation&amp;<sp/>transformation)<sp/>{</highlight></codeline>
<codeline lineno="91"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>Bounds<sp/>result;</highlight></codeline>
<codeline lineno="92"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>result.lower.fill(std::numeric_limits&lt;double&gt;::infinity());</highlight></codeline>
<codeline lineno="93"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>result.upper.fill(-std::numeric_limits&lt;double&gt;::infinity());</highlight></codeline>
<codeline lineno="94"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal">&amp;<sp/>rule<sp/>=<sp/>::mfem::IntRules.Get(transformation.GetGeometryType(),<sp/>2<sp/>*<sp/>transformation.Order()<sp/>+<sp/>4);</highlight></codeline>
<codeline lineno="95"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal">*<sp/>vertices<sp/>=<sp/>::mfem::Geometries.GetVertices(transformation.GetGeometryType());</highlight></codeline>
<codeline lineno="96"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1backend_1_1mfem_1afdd4d9966e5040b2de0684e018dce66e" kindref="member">::mfem::Vector</ref><sp/>point(3);</highlight></codeline>
<codeline lineno="97"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>visit<sp/>=<sp/>[&amp;](const<sp/>::mfem::IntegrationPoint&amp;<sp/>ip)<sp/>{</highlight></codeline>
<codeline lineno="98"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>transformation.Transform(ip,<sp/>point);</highlight></codeline>
<codeline lineno="99"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>k<sp/>=<sp/>0;<sp/>k<sp/>&lt;<sp/>3;<sp/>++k)<sp/>{</highlight></codeline>
<codeline lineno="100"><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">if</highlight><highlight class="normal"><sp/>(!std::isfinite(point[k]))<sp/>{<sp/></highlight><highlight class="keywordflow">throw</highlight><highlight class="normal"><sp/>std::invalid_argument(</highlight><highlight class="stringliteral">&quot;TOP<sp/>sampling<sp/>found<sp/>nonfinite<sp/>deformed<sp/>geometry.&quot;</highlight><highlight class="normal">);<sp/>}</highlight></codeline>
<codeline lineno="101"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>result.lower[k]<sp/>=<sp/>std::min(result.lower[k],<sp/>point[k]);</highlight></codeline>
<codeline lineno="102"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>result.upper[k]<sp/>=<sp/>std::max(result.upper[k],<sp/>point[k]);</highlight></codeline>
<codeline lineno="103"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="104"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="105"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>q<sp/>=<sp/>0;<sp/>q<sp/>&lt;<sp/>rule.GetNPoints();<sp/>++q)<sp/>{<sp/>visit(rule.IntPoint(q));<sp/>}</highlight></codeline>
<codeline lineno="106"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>q<sp/>=<sp/>0;<sp/>q<sp/>&lt;<sp/>vertices-&gt;GetNPoints();<sp/>++q)<sp/>{<sp/>visit(vertices-&gt;IntPoint(q));<sp/>}</highlight></codeline>
<codeline lineno="107"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>diameter<sp/>=<sp/>0.0;</highlight></codeline>
<codeline lineno="108"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>k<sp/>=<sp/>0;<sp/>k<sp/>&lt;<sp/>3;<sp/>++k)<sp/>{<sp/>diameter<sp/>=<sp/>std::max(diameter,<sp/>result.upper[k]<sp/>-<sp/>result.lower[k]);<sp/>}</highlight></codeline>
<codeline lineno="109"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>k<sp/>=<sp/>0;<sp/>k<sp/>&lt;<sp/>3;<sp/>++k)<sp/>{</highlight></codeline>
<codeline lineno="110"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>result.lower[k]<sp/>-=<sp/>0.05<sp/>*<sp/>diameter;</highlight></codeline>
<codeline lineno="111"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>result.upper[k]<sp/>+=<sp/>0.05<sp/>*<sp/>diameter;</highlight></codeline>
<codeline lineno="112"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="113"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>result;</highlight></codeline>
<codeline lineno="114"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="115"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="116"><highlight class="normal"></highlight></codeline>
<codeline lineno="117" refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl" refkind="compound"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">struct<sp/></highlight><highlight class="normal"><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl" kindref="compound">TopModelSampler::Impl</ref><sp/>{</highlight></codeline>
<codeline lineno="118"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::Mesh<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>;</highlight></codeline>
<codeline lineno="119"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>DensityBinding<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a83c029bd3792b0a85315bc1905f5e7e6" kindref="member">density_binding</ref>;</highlight></codeline>
<codeline lineno="120"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1backend_1_1mfem_1afdd4d9966e5040b2de0684e018dce66e" kindref="member">::mfem::Vector</ref><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a22e7e90fd0232dcd294dc20a83bf2f14" kindref="member">density_coefficients</ref>;</highlight></codeline>
<codeline lineno="121"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="classserif_1_1discretization_1_1backend_1_1mfem_1_1_field_values" kindref="compound">discretization::backend::mfem::FieldValues&lt;discretization::blocks::density::mass&gt;</ref><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a8b3cd740d35e3dcc8e87881ed16f3994" kindref="member">density_values</ref>;</highlight></codeline>
<codeline lineno="122"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>Vec<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a9289593e42872bfc91200ebcea927304" kindref="member">axis</ref>,<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad695bf1897581c9ae922380b2351d3ca" kindref="member">center</ref>,<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a043b20664ac109ac3db0d771b35757ab" kindref="member">first</ref>,<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a613797a28ac2d08bef5e0fe8c386a025" kindref="member">second</ref>;</highlight></codeline>
<codeline lineno="123"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::vector&lt;int&gt;<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1acc2d88650d1b9c9e274445b56c30dde8" kindref="member">material</ref>,<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a9112c1297929fc007097ad2124a382f6" kindref="member">boundary</ref>;</highlight></codeline>
<codeline lineno="124"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::vector&lt;Bounds&gt;<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a2de173cda192a08dda21c73b48732e5f" kindref="member">material_bounds</ref>,<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a9b9dee8b7f254c922764664f12a72749" kindref="member">boundary_bounds</ref>;</highlight></codeline>
<codeline lineno="125"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::map&lt;std::pair&lt;double,<sp/>double&gt;,<sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal">&gt;<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad729491cb2e4686a5a876c136f18f77e" kindref="member">radii</ref>;</highlight></codeline>
<codeline lineno="126"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a5826e7f5b36ddfcaef9db42e1cb59189" kindref="member">density_scale</ref>{0.0},<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ab84dd7daeca474df1e393ffc0f3e694f" kindref="member">total_mass</ref>{0.0},<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a6f0d903bf13dd4c76ed6eb3cf4fafdbc" kindref="member">length_scale</ref>{0.0};</highlight></codeline>
<codeline lineno="127"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::size_t<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ac6fa8d17bbc3747c6296a220de74b0f1" kindref="member">clamps</ref>{0};</highlight></codeline>
<codeline lineno="128"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad98c4d933c14ea6f809b03aa38da8e11" kindref="member">last_element</ref>{-1};</highlight></codeline>
<codeline lineno="129"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::IntegrationPoint<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1afc485b408dc4164446231108e3a05713" kindref="member">last_point</ref>;</highlight></codeline>
<codeline lineno="130"><highlight class="normal"></highlight></codeline>
<codeline lineno="131"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a09ffe52d54e5978bb64abc8d3007afc8" kindref="member">Impl</ref>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="classserif_1_1discretization_1_1backend_1_1mfem_1_1_dof_map" kindref="compound">discretization::backend::mfem::DofMap</ref>&amp;<sp/>binding,<sp/>std::span&lt;const<sp/>double&gt;<sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1a6803a4e0f33dc359c6ad55ae43d22a15" kindref="member">density</ref>,</highlight></codeline>
<codeline lineno="132"><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="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_reference_discretization" kindref="compound">geometry::deformation::backend::mfem::ReferenceDiscretization</ref>&amp;<sp/>reference,</highlight></codeline>
<codeline lineno="133"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::span&lt;const<sp/>double&gt;<sp/>displacement,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>Vec&amp;<sp/>supplied_axis,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>Vec&amp;<sp/>supplied_center)<sp/>:</highlight></codeline>
<codeline lineno="134"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>(binding.<ref refid="classserif_1_1discretization_1_1backend_1_1mfem_1_1_dof_map_1a43309a9a969cda1821eca7c08fc3167d" kindref="member">mesh</ref>(),<sp/></highlight><highlight class="keyword">true</highlight><highlight class="normal">),<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a83c029bd3792b0a85315bc1905f5e7e6" kindref="member">density_binding</ref>(binding),<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a22e7e90fd0232dcd294dc20a83bf2f14" kindref="member">density_coefficients</ref>(copy_density(<ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1a6803a4e0f33dc359c6ad55ae43d22a15" kindref="member">density</ref>)),</highlight></codeline>
<codeline lineno="135"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a8b3cd740d35e3dcc8e87881ed16f3994" kindref="member">density_values</ref>(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a83c029bd3792b0a85315bc1905f5e7e6" kindref="member">density_binding</ref>,<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a22e7e90fd0232dcd294dc20a83bf2f14" kindref="member">density_coefficients</ref>),<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a9289593e42872bfc91200ebcea927304" kindref="member">axis</ref>(normalized(supplied_axis)),<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad695bf1897581c9ae922380b2351d3ca" kindref="member">center</ref>(supplied_center)<sp/>{</highlight></codeline>
<codeline lineno="136"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>binding.<ref refid="classserif_1_1discretization_1_1backend_1_1mfem_1_1_dof_map_1a6b3f1ddbb971cfb23f779930a1417785" kindref="member">validate</ref>();</highlight></codeline>
<codeline lineno="137"><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/>(binding.<ref refid="classserif_1_1discretization_1_1backend_1_1mfem_1_1_dof_map_1ad7309ffab5c37db444ecd53450f3cfdc" kindref="member">space</ref>().GetVDim()<sp/>!=<sp/>1<sp/>||<sp/>binding.<ref refid="classserif_1_1discretization_1_1backend_1_1mfem_1_1_dof_map_1ad7309ffab5c37db444ecd53450f3cfdc" kindref="member">space</ref>().IsVariableOrder()<sp/>||</highlight></codeline>
<codeline lineno="138"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>reference.<ref refid="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_reference_discretization_1ae51ef5f73a4def1a7866cb664ace3985" kindref="member">displacement_space</ref>.IsVariableOrder()<sp/>||<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>.Dimension()<sp/>!=<sp/>3<sp/>||<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>.SpaceDimension()<sp/>!=<sp/>3<sp/>||</highlight></codeline>
<codeline lineno="139"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>binding.<ref refid="classserif_1_1discretization_1_1backend_1_1mfem_1_1_dof_map_1ad7309ffab5c37db444ecd53450f3cfdc" kindref="member">space</ref>().GetMesh()<sp/>!=<sp/>reference.<ref refid="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_reference_discretization_1ae51ef5f73a4def1a7866cb664ace3985" kindref="member">displacement_space</ref>.GetMesh()<sp/>||</highlight></codeline>
<codeline lineno="140"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1a6803a4e0f33dc359c6ad55ae43d22a15" kindref="member">density</ref>.size()<sp/>!=<sp/></highlight><highlight class="keyword">static_cast&lt;</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">&gt;</highlight><highlight class="normal">(binding.<ref refid="classserif_1_1discretization_1_1backend_1_1mfem_1_1_dof_map_1a8a6ce2a8958a84dfc99ca0484c64f302" kindref="member">size</ref>()))<sp/>{</highlight></codeline>
<codeline lineno="141"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">throw</highlight><highlight class="normal"><sp/>std::invalid_argument(</highlight><highlight class="stringliteral">&quot;TOP<sp/>sampling<sp/>requires<sp/>matching<sp/>fixed-order<sp/>serial<sp/>scalar<sp/>density<sp/>and<sp/>three-dimensional<sp/>displacement<sp/>spaces.&quot;</highlight><highlight class="normal">);</highlight></codeline>
<codeline lineno="142"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><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">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>value<sp/>:<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad695bf1897581c9ae922380b2351d3ca" kindref="member">center</ref>)<sp/>{</highlight></codeline>
<codeline lineno="144"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!std::isfinite(value))<sp/>{<sp/></highlight><highlight class="keywordflow">throw</highlight><highlight class="normal"><sp/>std::invalid_argument(</highlight><highlight class="stringliteral">&quot;TOP<sp/>sampling<sp/>requires<sp/>a<sp/>finite<sp/>rotation<sp/>center.&quot;</highlight><highlight class="normal">);<sp/>}</highlight></codeline>
<codeline lineno="145"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="146"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::tie(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a043b20664ac109ac3db0d771b35757ab" kindref="member">first</ref>,<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a613797a28ac2d08bef5e0fe8c386a025" kindref="member">second</ref>)<sp/>=<sp/>transverse_basis(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a9289593e42872bfc91200ebcea927304" kindref="member">axis</ref>);</highlight></codeline>
<codeline lineno="147"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::GridFunction<sp/>motion(</highlight><highlight class="keyword">const_cast&lt;</highlight><highlight class="normal">::mfem::FiniteElementSpace*</highlight><highlight class="keyword">&gt;</highlight><highlight class="normal">(&amp;reference.<ref refid="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_reference_discretization_1ae51ef5f73a4def1a7866cb664ace3985" kindref="member">displacement_space</ref>));</highlight></codeline>
<codeline lineno="148"><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/>(motion.Size()<sp/>!=<sp/></highlight><highlight class="keyword">static_cast&lt;</highlight><highlight class="keywordtype">int</highlight><highlight class="keyword">&gt;</highlight><highlight class="normal">(displacement.size())<sp/>||<sp/>motion.FESpace()-&gt;GetVDim()<sp/>!=<sp/>3)<sp/>{</highlight></codeline>
<codeline lineno="149"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">throw</highlight><highlight class="normal"><sp/>std::invalid_argument(</highlight><highlight class="stringliteral">&quot;TOP<sp/>displacement<sp/>coefficients<sp/>do<sp/>not<sp/>match<sp/>the<sp/>deformation<sp/>reference<sp/>space.&quot;</highlight><highlight class="normal">);</highlight></codeline>
<codeline lineno="150"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="151"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>nodes<sp/>=<sp/>motion.FESpace()-&gt;GetNDofs();</highlight></codeline>
<codeline lineno="152"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>node<sp/>=<sp/>0;<sp/>node<sp/>&lt;<sp/>nodes;<sp/>++node)<sp/>{</highlight></codeline>
<codeline lineno="153"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>component<sp/>=<sp/>0;<sp/>component<sp/>&lt;<sp/>3;<sp/>++component)<sp/>{</highlight></codeline>
<codeline lineno="154"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>value<sp/>=<sp/>displacement[node<sp/>+<sp/>component<sp/>*<sp/>nodes];</highlight></codeline>
<codeline lineno="155"><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">if</highlight><highlight class="normal"><sp/>(!std::isfinite(value))<sp/>{<sp/></highlight><highlight class="keywordflow">throw</highlight><highlight class="normal"><sp/>std::invalid_argument(</highlight><highlight class="stringliteral">&quot;TOP<sp/>displacement<sp/>contains<sp/>a<sp/>nonfinite<sp/>coefficient.&quot;</highlight><highlight class="normal">);<sp/>}</highlight></codeline>
<codeline lineno="156"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>motion[motion.FESpace()-&gt;DofToVDof(node,<sp/>component)]<sp/>=<sp/>value;</highlight></codeline>
<codeline lineno="157"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="158"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="159"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Preserve<sp/>all<sp/>reference-geometry<sp/>modes<sp/>before<sp/>adding<sp/>the<sp/>lower-<sp/>or<sp/>higher-order<sp/>displacement.</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="160"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>order<sp/>=<sp/>reference.<ref refid="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_reference_discretization_1ae51ef5f73a4def1a7866cb664ace3985" kindref="member">displacement_space</ref>.GetMaxElementOrder();</highlight></codeline>
<codeline lineno="161"><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="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>.GetNodalFESpace())<sp/>{<sp/>order<sp/>=<sp/>std::max(order,<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>.GetNodalFESpace()-&gt;GetMaxElementOrder());<sp/>}</highlight></codeline>
<codeline lineno="162"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>.SetCurvature(order,<sp/></highlight><highlight class="keyword">false</highlight><highlight class="normal">,<sp/>3,<sp/>::mfem::Ordering::byNODES);</highlight></codeline>
<codeline lineno="163"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::GridFunction<sp/>nodal_motion(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>.GetNodes()-&gt;FESpace());</highlight></codeline>
<codeline lineno="164"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>nodal_motion.ProjectGridFunction(motion);</highlight></codeline>
<codeline lineno="165"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>*<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>.GetNodes()<sp/>+=<sp/>nodal_motion;</highlight></codeline>
<codeline lineno="166"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>.NodesUpdated();</highlight></codeline>
<codeline lineno="167"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>element<sp/>=<sp/>0;<sp/>element<sp/>&lt;<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>.GetNE();<sp/>++element)<sp/>{</highlight></codeline>
<codeline lineno="168"><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/>(binding.<ref refid="classserif_1_1discretization_1_1backend_1_1mfem_1_1_dof_map_1aef238cf43b66e8499399b520d190f80a" kindref="member">supports</ref>(element))<sp/>{</highlight></codeline>
<codeline lineno="169"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1acc2d88650d1b9c9e274445b56c30dde8" kindref="member">material</ref>.push_back(element);</highlight></codeline>
<codeline lineno="170"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a2de173cda192a08dda21c73b48732e5f" kindref="member">material_bounds</ref>.push_back(bounds(*<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>.GetElementTransformation(element)));</highlight></codeline>
<codeline lineno="171"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="172"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="173"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>face<sp/>=<sp/>0;<sp/>face<sp/>&lt;<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>.GetNBE();<sp/>++face)<sp/>{</highlight></codeline>
<codeline lineno="174"><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/>(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>.GetBdrAttribute(face)<sp/>==<sp/>reference.<ref refid="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_reference_discretization_1a31f0502da71188954ce9150d835c2956" kindref="member">surface_attribute</ref>)<sp/>{</highlight></codeline>
<codeline lineno="175"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a9112c1297929fc007097ad2124a382f6" kindref="member">boundary</ref>.push_back(face);</highlight></codeline>
<codeline lineno="176"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a9b9dee8b7f254c922764664f12a72749" kindref="member">boundary_bounds</ref>.push_back(bounds(*<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>.GetBdrElementTransformation(face)));</highlight></codeline>
<codeline lineno="177"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="178"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="179"><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="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1acc2d88650d1b9c9e274445b56c30dde8" kindref="member">material</ref>.empty()<sp/>||<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a9112c1297929fc007097ad2124a382f6" kindref="member">boundary</ref>.empty())<sp/>{<sp/></highlight><highlight class="keywordflow">throw</highlight><highlight class="normal"><sp/>std::invalid_argument(</highlight><highlight class="stringliteral">&quot;TOP<sp/>sampling<sp/>found<sp/>no<sp/>material<sp/>or<sp/>stellar<sp/>surface.&quot;</highlight><highlight class="normal">);<sp/>}</highlight></codeline>
<codeline lineno="180"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal">&amp;<sp/>box<sp/>:<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a9b9dee8b7f254c922764664f12a72749" kindref="member">boundary_bounds</ref>)<sp/>{</highlight></codeline>
<codeline lineno="181"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>k<sp/>=<sp/>0;<sp/>k<sp/>&lt;<sp/>3;<sp/>++k)<sp/>{</highlight></codeline>
<codeline lineno="182"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a6f0d903bf13dd4c76ed6eb3cf4fafdbc" kindref="member">length_scale</ref><sp/>=<sp/>std::max({<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a6f0d903bf13dd4c76ed6eb3cf4fafdbc" kindref="member">length_scale</ref>,<sp/>std::abs(box.lower[k]<sp/>-<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad695bf1897581c9ae922380b2351d3ca" kindref="member">center</ref>[k]),<sp/>std::abs(box.upper[k]<sp/>-<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad695bf1897581c9ae922380b2351d3ca" kindref="member">center</ref>[k])});</highlight></codeline>
<codeline lineno="183"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="184"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="185"><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="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a6f0d903bf13dd4c76ed6eb3cf4fafdbc" kindref="member">length_scale</ref><sp/>&gt;<sp/>0.0)<sp/>||<sp/>!std::isfinite(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a6f0d903bf13dd4c76ed6eb3cf4fafdbc" kindref="member">length_scale</ref>))<sp/>{<sp/></highlight><highlight class="keywordflow">throw</highlight><highlight class="normal"><sp/>std::invalid_argument(</highlight><highlight class="stringliteral">&quot;Invalid<sp/>stellar<sp/>length<sp/>scale<sp/>for<sp/>TOP<sp/>sampling.&quot;</highlight><highlight class="normal">);<sp/>}</highlight></codeline>
<codeline lineno="186"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1backend_1_1mfem_1afdd4d9966e5040b2de0684e018dce66e" kindref="member">::mfem::Vector</ref><sp/>origin(3);</highlight></codeline>
<codeline lineno="187"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>k<sp/>=<sp/>0;<sp/>k<sp/>&lt;<sp/>3;<sp/>++k)<sp/>{<sp/>origin[k]<sp/>=<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad695bf1897581c9ae922380b2351d3ca" kindref="member">center</ref>[k];<sp/>}</highlight></codeline>
<codeline lineno="188"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>central<sp/>=<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a7c43ff965f8652cbe4f640749f5eaed1" kindref="member">sample_raw</ref>(origin);</highlight></codeline>
<codeline lineno="189"><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/>(!central<sp/>||<sp/>!std::isfinite(*central)<sp/>||<sp/>*central<sp/>&lt;=<sp/>0.0)<sp/>{</highlight></codeline>
<codeline lineno="190"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">throw</highlight><highlight class="normal"><sp/>std::invalid_argument(</highlight><highlight class="stringliteral">&quot;TOP<sp/>requires<sp/>positive<sp/>density<sp/>at<sp/>a<sp/>rotation<sp/>center<sp/>inside<sp/>the<sp/>stellar<sp/>material.&quot;</highlight><highlight class="normal">);</highlight></codeline>
<codeline lineno="191"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="192"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a5826e7f5b36ddfcaef9db42e1cb59189" kindref="member">density_scale</ref><sp/>=<sp/>*central;</highlight></codeline>
<codeline lineno="193"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>correction<sp/>=<sp/>0.0;</highlight></codeline>
<codeline lineno="194"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1backend_1_1mfem_1afdd4d9966e5040b2de0684e018dce66e" kindref="member">::mfem::Vector</ref><sp/>density_samples;</highlight></codeline>
<codeline lineno="195"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>element<sp/>:<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1acc2d88650d1b9c9e274445b56c30dde8" kindref="member">material</ref>)<sp/>{</highlight></codeline>
<codeline lineno="196"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal">&amp;<sp/>transformation<sp/>=<sp/>*<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>.GetElementTransformation(element);</highlight></codeline>
<codeline lineno="197"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Fractional<sp/>radial<sp/>powers<sp/>are<sp/>not<sp/>polynomials;<sp/>overintegrate<sp/>only<sp/>the<sp/>reconstructed<sp/>collar.</highlight><highlight class="normal"></highlight></codeline>
<codeline lineno="198"><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/></highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>additional_order<sp/>=<sp/>binding.<ref refid="classserif_1_1discretization_1_1backend_1_1mfem_1_1_dof_map_1ab9bbd3372044a045d75ce1c9b2bb3097" kindref="member">reconstruction_element</ref>(element)<sp/>==<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal"><sp/>?<sp/>0<sp/>:<sp/>binding.<ref refid="classserif_1_1discretization_1_1backend_1_1mfem_1_1_dof_map_1a575fa7bf1d77adeb2e06fababd825f72" kindref="member">reconstruction</ref>()-&gt;integration_order_increment();</highlight></codeline>
<codeline lineno="199"><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/></highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>quadrature_order<sp/>=<sp/>binding.<ref refid="classserif_1_1discretization_1_1backend_1_1mfem_1_1_dof_map_1ad7309ffab5c37db444ecd53450f3cfdc" kindref="member">space</ref>().GetFE(element)-&gt;GetOrder()<sp/>+<sp/>3<sp/>*<sp/>order<sp/>+<sp/>3<sp/>+<sp/>additional_order;</highlight></codeline>
<codeline lineno="200"><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/></highlight><highlight class="keyword">auto</highlight><highlight class="normal">&amp;<sp/>rule<sp/>=<sp/>::mfem::IntRules.Get(transformation.GetGeometryType(),<sp/>quadrature_order);</highlight></codeline>
<codeline lineno="201"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a8b3cd740d35e3dcc8e87881ed16f3994" kindref="member">density_values</ref>.values(element,<sp/>rule,<sp/>density_samples);</highlight></codeline>
<codeline lineno="202"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>q<sp/>=<sp/>0;<sp/>q<sp/>&lt;<sp/>rule.GetNPoints();<sp/>++q)<sp/>{</highlight></codeline>
<codeline lineno="203"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal">&amp;<sp/>ip<sp/>=<sp/>rule.IntPoint(q);</highlight></codeline>
<codeline lineno="204"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>transformation.SetIntPoint(&amp;ip);</highlight></codeline>
<codeline lineno="205"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>determinant<sp/>=<sp/>transformation.Jacobian().Det();</highlight></codeline>
<codeline lineno="206"><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">if</highlight><highlight class="normal"><sp/>(!std::isfinite(determinant)<sp/>||<sp/>determinant<sp/>&lt;=<sp/>0.0)<sp/>{</highlight></codeline>
<codeline lineno="207"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">throw</highlight><highlight class="normal"><sp/>std::invalid_argument(</highlight><highlight class="stringliteral">&quot;TOP<sp/>sampling<sp/>rejected<sp/>a<sp/>nonpositive<sp/>deformed<sp/>material<sp/>Jacobian.&quot;</highlight><highlight class="normal">);</highlight></codeline>
<codeline lineno="208"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="209"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>rho<sp/>=<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a14439c331108199bdfd79c5cc677f6b6" kindref="member">checked_density</ref>(density_samples[q]);</highlight></codeline>
<codeline lineno="210"><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">if</highlight><highlight class="normal"><sp/>(!rho)<sp/>{<sp/></highlight><highlight class="keywordflow">throw</highlight><highlight class="normal"><sp/>std::invalid_argument(rho.error().message);<sp/>}</highlight></codeline>
<codeline lineno="211"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>term<sp/>=<sp/>*rho<sp/>*<sp/>determinant<sp/>*<sp/>ip.weight<sp/>-<sp/>correction;</highlight></codeline>
<codeline lineno="212"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>next<sp/>=<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ab84dd7daeca474df1e393ffc0f3e694f" kindref="member">total_mass</ref><sp/>+<sp/>term;</highlight></codeline>
<codeline lineno="213"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>correction<sp/>=<sp/>(next<sp/>-<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ab84dd7daeca474df1e393ffc0f3e694f" kindref="member">total_mass</ref>)<sp/>-<sp/>term;</highlight></codeline>
<codeline lineno="214"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ab84dd7daeca474df1e393ffc0f3e694f" kindref="member">total_mass</ref><sp/>=<sp/>next;</highlight></codeline>
<codeline lineno="215"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="216"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="217"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!std::isfinite(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ab84dd7daeca474df1e393ffc0f3e694f" kindref="member">total_mass</ref>)<sp/>||<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ab84dd7daeca474df1e393ffc0f3e694f" kindref="member">total_mass</ref><sp/>&lt;=<sp/>0.0)<sp/>{<sp/></highlight><highlight class="keywordflow">throw</highlight><highlight class="normal"><sp/>std::invalid_argument(</highlight><highlight class="stringliteral">&quot;TOP<sp/>requires<sp/>positive<sp/>finite<sp/>integrated<sp/>stellar<sp/>mass.&quot;</highlight><highlight class="normal">);<sp/>}</highlight></codeline>
<codeline lineno="218"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="219"><highlight class="normal"></highlight></codeline>
<codeline lineno="220"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>Vec<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ae32afa65a346f5f33469021c75efe46d" kindref="member">direction</ref>(</highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>mu,<sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>phi)</highlight><highlight class="keyword"><sp/>const<sp/></highlight><highlight class="normal">{</highlight></codeline>
<codeline lineno="221"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!std::isfinite(mu)<sp/>||<sp/>!std::isfinite(phi)<sp/>||<sp/>mu<sp/>&lt;<sp/>-1.0<sp/>||<sp/>mu<sp/>&gt;<sp/>1.0)<sp/>{</highlight></codeline>
<codeline lineno="222"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">throw</highlight><highlight class="normal"><sp/>std::invalid_argument(</highlight><highlight class="stringliteral">&quot;TOP<sp/>sampling<sp/>requires<sp/>finite<sp/>angles<sp/>with<sp/>-1<sp/>&lt;=<sp/>mu<sp/>&lt;=<sp/>1.&quot;</highlight><highlight class="normal">);</highlight></codeline>
<codeline lineno="223"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="224"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>sine<sp/>=<sp/>std::sqrt(std::max(0.0,<sp/>1.0<sp/>-<sp/>mu<sp/>*<sp/>mu));</highlight></codeline>
<codeline lineno="225"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>Vec<sp/>result;</highlight></codeline>
<codeline lineno="226"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>k<sp/>=<sp/>0;<sp/>k<sp/>&lt;<sp/>3;<sp/>++k)<sp/>{<sp/>result[k]<sp/>=<sp/>mu<sp/>*<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a9289593e42872bfc91200ebcea927304" kindref="member">axis</ref>[k]<sp/>+<sp/>sine<sp/>*<sp/>(std::cos(phi)<sp/>*<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a043b20664ac109ac3db0d771b35757ab" kindref="member">first</ref>[k]<sp/>+<sp/>std::sin(phi)<sp/>*<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a613797a28ac2d08bef5e0fe8c386a025" kindref="member">second</ref>[k]);<sp/>}</highlight></codeline>
<codeline lineno="227"><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/>result;</highlight></codeline>
<codeline lineno="228"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="229"><highlight class="normal"></highlight></codeline>
<codeline lineno="230"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1ace9e0af29889a767ef1c95a09a4b7265" kindref="member">Result&lt;double&gt;</ref><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a14439c331108199bdfd79c5cc677f6b6" kindref="member">checked_density</ref>(</highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>value)<sp/>{</highlight></codeline>
<codeline lineno="231"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(!std::isfinite(value))<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_error" kindref="compound">Error</ref>{</highlight><highlight class="stringliteral">&quot;top/state/density&quot;</highlight><highlight class="normal">,<sp/></highlight><highlight class="stringliteral">&quot;The<sp/>sampled<sp/>FE<sp/>density<sp/>is<sp/>nonfinite.&quot;</highlight><highlight class="normal">});<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">if</highlight><highlight class="normal"><sp/>(value<sp/>&lt;<sp/>0.0)<sp/>{</highlight></codeline>
<codeline lineno="233"><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/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>tolerance<sp/>=<sp/>128.0<sp/>*<sp/>std::numeric_limits&lt;double&gt;::epsilon()<sp/>*<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a5826e7f5b36ddfcaef9db42e1cb59189" kindref="member">density_scale</ref>;</highlight></codeline>
<codeline lineno="234"><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/>(-value<sp/>&gt;<sp/>tolerance)<sp/>{</highlight></codeline>
<codeline lineno="235"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_error" kindref="compound">Error</ref>{</highlight><highlight class="stringliteral">&quot;top/state/density&quot;</highlight><highlight class="normal">,<sp/></highlight><highlight class="stringliteral">&quot;The<sp/>sampled<sp/>FE<sp/>density<sp/>is<sp/>negative<sp/>beyond<sp/>128<sp/>machine<sp/>epsilons<sp/>times<sp/>the<sp/>central<sp/>density;<sp/>TOP<sp/>export<sp/>does<sp/>not<sp/>repair<sp/>the<sp/>stellar<sp/>solution.&quot;</highlight><highlight class="normal">});</highlight></codeline>
<codeline lineno="236"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="237"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>++<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ac6fa8d17bbc3747c6296a220de74b0f1" kindref="member">clamps</ref>;</highlight></codeline>
<codeline lineno="238"><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/>0.0;</highlight></codeline>
<codeline lineno="239"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="240"><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/>value;</highlight></codeline>
<codeline lineno="241"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="242"><highlight class="normal"></highlight></codeline>
<codeline lineno="243"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1ace9e0af29889a767ef1c95a09a4b7265" kindref="member">Result&lt;double&gt;</ref><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a7c43ff965f8652cbe4f640749f5eaed1" kindref="member">sample_raw</ref>(const<sp/>::mfem::Vector&amp;<sp/>point)<sp/>{</highlight></codeline>
<codeline lineno="244"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>attempt<sp/>=<sp/>[&amp;](</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>element,<sp/></highlight><highlight class="keywordtype">bool</highlight><highlight class="normal"><sp/>previous,<sp/></highlight><highlight class="keywordtype">bool</highlight><highlight class="normal"><sp/>robust)<sp/>-&gt;<sp/></highlight><highlight class="keywordtype">bool</highlight><highlight class="normal"><sp/>{</highlight></codeline>
<codeline lineno="245"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal">&amp;<sp/>transformation<sp/>=<sp/>*<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>.GetElementTransformation(element);</highlight></codeline>
<codeline lineno="246"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::InverseElementTransformation<sp/>inverse(&amp;transformation);</highlight></codeline>
<codeline lineno="247"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>inverse.SetMaxIter(40);</highlight></codeline>
<codeline lineno="248"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>inverse.SetPhysicalRelTol(2.0e-13);</highlight></codeline>
<codeline lineno="249"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>inverse.SetReferenceTol(2.0e-13);</highlight></codeline>
<codeline lineno="250"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>inverse.SetPrintLevel(-1);</highlight></codeline>
<codeline lineno="251"><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/>(previous)<sp/>{<sp/>inverse.SetInitialGuess(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1afc485b408dc4164446231108e3a05713" kindref="member">last_point</ref>);<sp/>}</highlight></codeline>
<codeline lineno="252"><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/>(robust)<sp/>{<sp/>inverse.SetInitialGuessType(::mfem::InverseElementTransformation::ClosestPhysNode);<sp/>}</highlight></codeline>
<codeline lineno="253"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::IntegrationPoint<sp/>ip;</highlight></codeline>
<codeline lineno="254"><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/>(inverse.Transform(point,<sp/>ip)<sp/>!=<sp/>::mfem::InverseElementTransformation::Inside)<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">false</highlight><highlight class="normal">;<sp/>}</highlight></codeline>
<codeline lineno="255"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1backend_1_1mfem_1afdd4d9966e5040b2de0684e018dce66e" kindref="member">::mfem::Vector</ref><sp/>reconstructed(3);</highlight></codeline>
<codeline lineno="256"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>transformation.Transform(ip,<sp/>reconstructed);</highlight></codeline>
<codeline lineno="257"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>reconstructed<sp/>-=<sp/>point;</highlight></codeline>
<codeline lineno="258"><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/>(reconstructed.Norml2()<sp/>&gt;<sp/>2.0e-11<sp/>*<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a6f0d903bf13dd4c76ed6eb3cf4fafdbc" kindref="member">length_scale</ref>)<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">false</highlight><highlight class="normal">;<sp/>}</highlight></codeline>
<codeline lineno="259"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad98c4d933c14ea6f809b03aa38da8e11" kindref="member">last_element</ref><sp/>=<sp/>element;</highlight></codeline>
<codeline lineno="260"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1afc485b408dc4164446231108e3a05713" kindref="member">last_point</ref><sp/>=<sp/>ip;</highlight></codeline>
<codeline lineno="261"><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/></highlight><highlight class="keyword">true</highlight><highlight class="normal">;</highlight></codeline>
<codeline lineno="262"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="263"><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="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad98c4d933c14ea6f809b03aa38da8e11" kindref="member">last_element</ref><sp/>&gt;=<sp/>0<sp/>&amp;&amp;<sp/>attempt(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad98c4d933c14ea6f809b03aa38da8e11" kindref="member">last_element</ref>,<sp/></highlight><highlight class="keyword">true</highlight><highlight class="normal">,<sp/></highlight><highlight class="keyword">false</highlight><highlight class="normal">))<sp/>{</highlight></codeline>
<codeline lineno="264"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a8b3cd740d35e3dcc8e87881ed16f3994" kindref="member">density_values</ref>.value(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad98c4d933c14ea6f809b03aa38da8e11" kindref="member">last_element</ref>,<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1afc485b408dc4164446231108e3a05713" kindref="member">last_point</ref>);</highlight></codeline>
<codeline lineno="265"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="266"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(std::size_t<sp/>i<sp/>=<sp/>0;<sp/>i<sp/>&lt;<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1acc2d88650d1b9c9e274445b56c30dde8" kindref="member">material</ref>.size();<sp/>++i)<sp/>{</highlight></codeline>
<codeline lineno="267"><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/>(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a2de173cda192a08dda21c73b48732e5f" kindref="member">material_bounds</ref>[i].contains(point)<sp/>&amp;&amp;<sp/>attempt(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1acc2d88650d1b9c9e274445b56c30dde8" kindref="member">material</ref>[i],<sp/></highlight><highlight class="keyword">false</highlight><highlight class="normal">,<sp/></highlight><highlight class="keyword">false</highlight><highlight class="normal">))<sp/>{</highlight></codeline>
<codeline lineno="268"><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/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a8b3cd740d35e3dcc8e87881ed16f3994" kindref="member">density_values</ref>.value(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad98c4d933c14ea6f809b03aa38da8e11" kindref="member">last_element</ref>,<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1afc485b408dc4164446231108e3a05713" kindref="member">last_point</ref>);</highlight></codeline>
<codeline lineno="269"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="270"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="271"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>element<sp/>:<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1acc2d88650d1b9c9e274445b56c30dde8" kindref="member">material</ref>)<sp/>{</highlight></codeline>
<codeline lineno="272"><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/>(attempt(element,<sp/></highlight><highlight class="keyword">false</highlight><highlight class="normal">,<sp/></highlight><highlight class="keyword">true</highlight><highlight class="normal">))<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a8b3cd740d35e3dcc8e87881ed16f3994" kindref="member">density_values</ref>.value(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad98c4d933c14ea6f809b03aa38da8e11" kindref="member">last_element</ref>,<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1afc485b408dc4164446231108e3a05713" kindref="member">last_point</ref>);<sp/>}</highlight></codeline>
<codeline lineno="273"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="274"><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(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_error" kindref="compound">Error</ref>{</highlight><highlight class="stringliteral">&quot;top/sampling&quot;</highlight><highlight class="normal">,<sp/></highlight><highlight class="stringliteral">&quot;A<sp/>requested<sp/>physical<sp/>point<sp/>could<sp/>not<sp/>be<sp/>located<sp/>inside<sp/>the<sp/>deformed<sp/>stellar<sp/>material;<sp/>no<sp/>field<sp/>extrapolation<sp/>was<sp/>performed.&quot;</highlight><highlight class="normal">});</highlight></codeline>
<codeline lineno="275"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="276"><highlight class="normal"></highlight></codeline>
<codeline lineno="277"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::optional&lt;double&gt;<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ae647a22f6d2c0fef01cf039c91474fbe" kindref="member">intersect</ref>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>face,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>Vec&amp;<sp/>ray,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>Vec&amp;<sp/>a,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>Vec&amp;<sp/>b,<sp/></highlight><highlight class="keywordtype">bool</highlight><highlight class="normal"><sp/>robust)<sp/>{</highlight></codeline>
<codeline lineno="278"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal">&amp;<sp/>transformation<sp/>=<sp/>*<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a47bc85e37dbfae7b1bd20a0e8ba8a04c" kindref="member">mesh</ref>.GetBdrElementTransformation(face);</highlight></codeline>
<codeline lineno="279"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>solve<sp/>=<sp/>[&amp;](::mfem::IntegrationPoint<sp/>ip)<sp/>-&gt;<sp/>std::optional&lt;double&gt;<sp/>{</highlight></codeline>
<codeline lineno="280"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1backend_1_1mfem_1afdd4d9966e5040b2de0684e018dce66e" kindref="member">::mfem::Vector</ref><sp/>point(3);</highlight></codeline>
<codeline lineno="281"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>iteration<sp/>=<sp/>0;<sp/>iteration<sp/>&lt;<sp/>35;<sp/>++iteration)<sp/>{</highlight></codeline>
<codeline lineno="282"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>transformation.Transform(ip,<sp/>point);</highlight></codeline>
<codeline lineno="283"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>Vec<sp/>relative{point[0]<sp/>-<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad695bf1897581c9ae922380b2351d3ca" kindref="member">center</ref>[0],<sp/>point[1]<sp/>-<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad695bf1897581c9ae922380b2351d3ca" kindref="member">center</ref>[1],<sp/>point[2]<sp/>-<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad695bf1897581c9ae922380b2351d3ca" kindref="member">center</ref>[2]};</highlight></codeline>
<codeline lineno="284"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>f<sp/>=<sp/>dot(relative,<sp/>a),<sp/>g<sp/>=<sp/>dot(relative,<sp/>b);</highlight></codeline>
<codeline lineno="285"><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">if</highlight><highlight class="normal"><sp/>(std::hypot(f,<sp/>g)<sp/>&lt;=<sp/>5.0e-13<sp/>*<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a6f0d903bf13dd4c76ed6eb3cf4fafdbc" kindref="member">length_scale</ref>)<sp/>{</highlight></codeline>
<codeline lineno="286"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>radius<sp/>=<sp/>dot(relative,<sp/>ray);</highlight></codeline>
<codeline lineno="287"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(radius<sp/>&gt;<sp/>0.0<sp/>&amp;&amp;<sp/>::mfem::Geometry::CheckPoint(transformation.GetGeometryType(),<sp/>ip,<sp/>2.0e-10))<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>radius;<sp/>}</highlight></codeline>
<codeline lineno="288"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>{};</highlight></codeline>
<codeline lineno="289"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="290"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>transformation.SetIntPoint(&amp;ip);</highlight></codeline>
<codeline lineno="291"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal">&amp;<sp/>jacobian<sp/>=<sp/>transformation.Jacobian();</highlight></codeline>
<codeline lineno="292"><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="keywordtype">double</highlight><highlight class="normal"><sp/>j00<sp/>=<sp/>0.0,<sp/>j01<sp/>=<sp/>0.0,<sp/>j10<sp/>=<sp/>0.0,<sp/>j11<sp/>=<sp/>0.0;</highlight></codeline>
<codeline lineno="293"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>k<sp/>=<sp/>0;<sp/>k<sp/>&lt;<sp/>3;<sp/>++k)<sp/>{</highlight></codeline>
<codeline lineno="294"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>j00<sp/>+=<sp/>a[k]<sp/>*<sp/>jacobian(k,<sp/>0);<sp/>j01<sp/>+=<sp/>a[k]<sp/>*<sp/>jacobian(k,<sp/>1);</highlight></codeline>
<codeline lineno="295"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>j10<sp/>+=<sp/>b[k]<sp/>*<sp/>jacobian(k,<sp/>0);<sp/>j11<sp/>+=<sp/>b[k]<sp/>*<sp/>jacobian(k,<sp/>1);</highlight></codeline>
<codeline lineno="296"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="297"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>determinant<sp/>=<sp/>j00<sp/>*<sp/>j11<sp/>-<sp/>j01<sp/>*<sp/>j10;</highlight></codeline>
<codeline lineno="298"><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">if</highlight><highlight class="normal"><sp/>(!std::isfinite(determinant)<sp/>||<sp/>std::abs(determinant)<sp/>&lt;<sp/>1.0e-16<sp/>*<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a6f0d903bf13dd4c76ed6eb3cf4fafdbc" kindref="member">length_scale</ref><sp/>*<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a6f0d903bf13dd4c76ed6eb3cf4fafdbc" kindref="member">length_scale</ref>)<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>{};<sp/>}</highlight></codeline>
<codeline lineno="299"><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="keywordtype">double</highlight><highlight class="normal"><sp/>du<sp/>=<sp/>(j11<sp/>*<sp/>f<sp/>-<sp/>j01<sp/>*<sp/>g)<sp/>/<sp/>determinant;</highlight></codeline>
<codeline lineno="300"><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="keywordtype">double</highlight><highlight class="normal"><sp/>dv<sp/>=<sp/>(j00<sp/>*<sp/>g<sp/>-<sp/>j10<sp/>*<sp/>f)<sp/>/<sp/>determinant;</highlight></codeline>
<codeline lineno="301"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>scale<sp/>=<sp/>std::max({1.0,<sp/>std::abs(du),<sp/>std::abs(dv)});</highlight></codeline>
<codeline lineno="302"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>ip.x<sp/>-=<sp/>du<sp/>/<sp/>scale;</highlight></codeline>
<codeline lineno="303"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>ip.y<sp/>-=<sp/>dv<sp/>/<sp/>scale;</highlight></codeline>
<codeline lineno="304"><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">if</highlight><highlight class="normal"><sp/>(!std::isfinite(ip.x)<sp/>||<sp/>!std::isfinite(ip.y)<sp/>||<sp/>ip.x<sp/>&lt;<sp/>-0.5<sp/>||<sp/>ip.x<sp/>&gt;<sp/>1.5<sp/>||<sp/>ip.y<sp/>&lt;<sp/>-0.5<sp/>||<sp/>ip.y<sp/>&gt;<sp/>1.5)<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>{};<sp/>}</highlight></codeline>
<codeline lineno="305"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="306"><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/>{};</highlight></codeline>
<codeline lineno="307"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="308"><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/>(</highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>result<sp/>=<sp/>solve(::mfem::Geometries.GetCenter(transformation.GetGeometryType())))<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>result;<sp/>}</highlight></codeline>
<codeline lineno="309"><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/>(robust)<sp/>{</highlight></codeline>
<codeline lineno="310"><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/></highlight><highlight class="keyword">auto</highlight><highlight class="normal">&amp;<sp/>rule<sp/>=<sp/>::mfem::IntRules.Get(transformation.GetGeometryType(),<sp/>4);</highlight></codeline>
<codeline lineno="311"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>q<sp/>=<sp/>0;<sp/>q<sp/>&lt;<sp/>rule.GetNPoints();<sp/>++q)<sp/>{<sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>result<sp/>=<sp/>solve(rule.IntPoint(q)))<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>result;<sp/>}<sp/>}</highlight></codeline>
<codeline lineno="312"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="313"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>{};</highlight></codeline>
<codeline lineno="314"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="315"><highlight class="normal"></highlight></codeline>
<codeline lineno="316"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1ace9e0af29889a767ef1c95a09a4b7265" kindref="member">Result&lt;double&gt;</ref><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ab1dc30f4ac2fa2c196bbe28e1afff3ca" kindref="member">surface</ref>(</highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>mu,<sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>phi)<sp/>{</highlight></codeline>
<codeline lineno="317"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>key<sp/>=<sp/>std::pair{mu,<sp/>phi};</highlight></codeline>
<codeline lineno="318"><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/>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>found<sp/>=<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad729491cb2e4686a5a876c136f18f77e" kindref="member">radii</ref>.find(key);<sp/>found<sp/>!=<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad729491cb2e4686a5a876c136f18f77e" kindref="member">radii</ref>.end())<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>found-&gt;second;<sp/>}</highlight></codeline>
<codeline lineno="319"><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/>Vec<sp/>ray<sp/>=<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ae32afa65a346f5f33469021c75efe46d" kindref="member">direction</ref>(mu,<sp/>phi);</highlight></codeline>
<codeline lineno="320"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>[a,<sp/>b]<sp/>=<sp/>transverse_basis(ray);</highlight></codeline>
<codeline lineno="321"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::optional&lt;double&gt;<sp/>result;</highlight></codeline>
<codeline lineno="322"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(std::size_t<sp/>i<sp/>=<sp/>0;<sp/>i<sp/>&lt;<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a9112c1297929fc007097ad2124a382f6" kindref="member">boundary</ref>.size();<sp/>++i)<sp/>{</highlight></codeline>
<codeline lineno="323"><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/>(!<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a9b9dee8b7f254c922764664f12a72749" kindref="member">boundary_bounds</ref>[i].intersects_ray(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad695bf1897581c9ae922380b2351d3ca" kindref="member">center</ref>,<sp/>ray))<sp/>{<sp/></highlight><highlight class="keywordflow">continue</highlight><highlight class="normal">;<sp/>}</highlight></codeline>
<codeline lineno="324"><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/>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>radius<sp/>=<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ae647a22f6d2c0fef01cf039c91474fbe" kindref="member">intersect</ref>(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a9112c1297929fc007097ad2124a382f6" kindref="member">boundary</ref>[i],<sp/>ray,<sp/>a,<sp/>b,<sp/></highlight><highlight class="keyword">true</highlight><highlight class="normal">))<sp/>{</highlight></codeline>
<codeline lineno="325"><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">if</highlight><highlight class="normal"><sp/>(result<sp/>&amp;&amp;<sp/>std::abs(*result<sp/>-<sp/>*radius)<sp/>&gt;<sp/>2.0e-9<sp/>*<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a6f0d903bf13dd4c76ed6eb3cf4fafdbc" kindref="member">length_scale</ref>)<sp/>{</highlight></codeline>
<codeline lineno="326"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_error" kindref="compound">Error</ref>{</highlight><highlight class="stringliteral">&quot;top/geometry/surface&quot;</highlight><highlight class="normal">,<sp/></highlight><highlight class="stringliteral">&quot;The<sp/>stellar<sp/>boundary<sp/>has<sp/>multiple<sp/>radial<sp/>intersections<sp/>and<sp/>cannot<sp/>be<sp/>represented<sp/>as<sp/>a<sp/>TOP<sp/>star-shaped<sp/>surface.&quot;</highlight><highlight class="normal">});</highlight></codeline>
<codeline lineno="327"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="328"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>result<sp/>=<sp/>radius;</highlight></codeline>
<codeline lineno="329"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="330"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="331"><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/>(!result)<sp/>{</highlight></codeline>
<codeline lineno="332"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>face<sp/>:<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a9112c1297929fc007097ad2124a382f6" kindref="member">boundary</ref>)<sp/>{</highlight></codeline>
<codeline lineno="333"><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">if</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>radius<sp/>=<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ae647a22f6d2c0fef01cf039c91474fbe" kindref="member">intersect</ref>(face,<sp/>ray,<sp/>a,<sp/>b,<sp/></highlight><highlight class="keyword">true</highlight><highlight class="normal">))<sp/>{</highlight></codeline>
<codeline lineno="334"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(result<sp/>&amp;&amp;<sp/>std::abs(*result<sp/>-<sp/>*radius)<sp/>&gt;<sp/>2.0e-9<sp/>*<sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1a6f0d903bf13dd4c76ed6eb3cf4fafdbc" kindref="member">length_scale</ref>)<sp/>{</highlight></codeline>
<codeline lineno="335"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_error" kindref="compound">Error</ref>{</highlight><highlight class="stringliteral">&quot;top/geometry/surface&quot;</highlight><highlight class="normal">,<sp/></highlight><highlight class="stringliteral">&quot;The<sp/>stellar<sp/>boundary<sp/>is<sp/>not<sp/>single-valued<sp/>in<sp/>radius.&quot;</highlight><highlight class="normal">});</highlight></codeline>
<codeline lineno="336"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="337"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>result<sp/>=<sp/>radius;</highlight></codeline>
<codeline lineno="338"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="339"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="340"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="341"><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/>(!result)<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_error" kindref="compound">Error</ref>{</highlight><highlight class="stringliteral">&quot;top/geometry/surface&quot;</highlight><highlight class="normal">,<sp/></highlight><highlight class="stringliteral">&quot;Could<sp/>not<sp/>intersect<sp/>the<sp/>requested<sp/>ray<sp/>with<sp/>the<sp/>actual<sp/>deformed<sp/>FE<sp/>stellar<sp/>surface.&quot;</highlight><highlight class="normal">});<sp/>}</highlight></codeline>
<codeline lineno="342"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl_1ad729491cb2e4686a5a876c136f18f77e" kindref="member">radii</ref>.emplace(key,<sp/>*result);</highlight></codeline>
<codeline lineno="343"><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/>*result;</highlight></codeline>
<codeline lineno="344"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="345"><highlight class="normal"><sp/><sp/><sp/><sp/>};</highlight></codeline>
<codeline lineno="346"><highlight class="normal"></highlight></codeline>
<codeline lineno="347"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1a638d2d9807f0d652dadeeeef02194b02" kindref="member">TopModelSampler::TopModelSampler</ref>(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="classserif_1_1discretization_1_1backend_1_1mfem_1_1_dof_map" kindref="compound">discretization::backend::mfem::DofMap</ref>&amp;<sp/>binding,</highlight></codeline>
<codeline lineno="348"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::span&lt;const<sp/>double&gt;<sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1a6803a4e0f33dc359c6ad55ae43d22a15" kindref="member">density</ref>,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/><ref refid="structserif_1_1geometry_1_1deformation_1_1backend_1_1mfem_1_1_reference_discretization" kindref="compound">geometry::deformation::backend::mfem::ReferenceDiscretization</ref>&amp;<sp/>reference,</highlight></codeline>
<codeline lineno="349"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::span&lt;const<sp/>double&gt;<sp/>displacement,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>Vec&amp;<sp/>axis,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>Vec&amp;<sp/>center)<sp/>:</highlight></codeline>
<codeline lineno="350"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1abbd87041d750fb2fdd34bd1621192dd5" kindref="member">m_impl</ref>(std::make_unique&lt;<ref refid="structserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1_1_impl" kindref="compound">Impl</ref>&gt;(binding,<sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1a6803a4e0f33dc359c6ad55ae43d22a15" kindref="member">density</ref>,<sp/>reference,<sp/>displacement,<sp/>axis,<sp/>center))<sp/>{}</highlight></codeline>
<codeline lineno="351"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1a4b9ba12ef9f1d1c1fe91e4b95f26c9f0" kindref="member">TopModelSampler::~TopModelSampler</ref>()<sp/>=<sp/></highlight><highlight class="keywordflow">default</highlight><highlight class="normal">;</highlight></codeline>
<codeline lineno="352"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1ace9e0af29889a767ef1c95a09a4b7265" kindref="member">Result&lt;double&gt;</ref><sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1a877f31e990f85c7013eac540dfa3b080" kindref="member">TopModelSampler::surface_radius</ref>(</highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>mu,<sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>phi)<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1abbd87041d750fb2fdd34bd1621192dd5" kindref="member">m_impl</ref>-&gt;surface(mu,<sp/>phi);<sp/>}</highlight></codeline>
<codeline lineno="353"><highlight class="normal"><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1serialize_1ace9e0af29889a767ef1c95a09a4b7265" kindref="member">Result&lt;double&gt;</ref><sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1a6803a4e0f33dc359c6ad55ae43d22a15" kindref="member">TopModelSampler::density</ref>(</highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>radius,<sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>mu,<sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>phi)<sp/>{</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/>(!std::isfinite(radius)<sp/>||<sp/>radius<sp/>&lt;<sp/>0.0)<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(<ref refid="structserif_1_1utils_1_1serialize_1_1_error" kindref="compound">Error</ref>{</highlight><highlight class="stringliteral">&quot;top/sampling/radius&quot;</highlight><highlight class="normal">,<sp/></highlight><highlight class="stringliteral">&quot;TOP<sp/>sampling<sp/>requires<sp/>a<sp/>finite<sp/>nonnegative<sp/>physical<sp/>radius.&quot;</highlight><highlight class="normal">});<sp/>}</highlight></codeline>
<codeline lineno="355"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(radius<sp/>==<sp/>0.0)<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1abbd87041d750fb2fdd34bd1621192dd5" kindref="member">m_impl</ref>-&gt;density_scale;<sp/>}</highlight></codeline>
<codeline lineno="356"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>Vec<sp/>direction<sp/>=<sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1abbd87041d750fb2fdd34bd1621192dd5" kindref="member">m_impl</ref>-&gt;direction(mu,<sp/>phi);</highlight></codeline>
<codeline lineno="357"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1backend_1_1mfem_1afdd4d9966e5040b2de0684e018dce66e" kindref="member">::mfem::Vector</ref><sp/>point(3);</highlight></codeline>
<codeline lineno="358"><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/>k<sp/>=<sp/>0;<sp/>k<sp/>&lt;<sp/>3;<sp/>++k)<sp/>{<sp/>point[k]<sp/>=<sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1abbd87041d750fb2fdd34bd1621192dd5" kindref="member">m_impl</ref>-&gt;center[k]<sp/>+<sp/>radius<sp/>*<sp/>direction[k];<sp/>}</highlight></codeline>
<codeline lineno="359"><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/>value<sp/>=<sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1abbd87041d750fb2fdd34bd1621192dd5" kindref="member">m_impl</ref>-&gt;sample_raw(point);</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/>(!value)<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(value.error());<sp/>}</highlight></codeline>
<codeline lineno="361"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1abbd87041d750fb2fdd34bd1621192dd5" kindref="member">m_impl</ref>-&gt;checked_density(*value);</highlight></codeline>
<codeline lineno="362"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
<codeline lineno="363"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1a4a3c41443e79e3056900878169478dd0" kindref="member">TopModelSampler::mass</ref>()<sp/>const<sp/>noexcept<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1abbd87041d750fb2fdd34bd1621192dd5" kindref="member">m_impl</ref>-&gt;total_mass;<sp/>}</highlight></codeline>
<codeline lineno="364"><highlight class="normal"><sp/><sp/><sp/><sp/>std::size_t<sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1aa3f8c1628586cb0aefbdf9d708abaf94" kindref="member">TopModelSampler::roundoff_density_clamps</ref>()<sp/>const<sp/>noexcept<sp/>{<sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/><ref refid="classserif_1_1backend_1_1mfem_1_1serialize_1_1detail_1_1_top_model_sampler_1abbd87041d750fb2fdd34bd1621192dd5" kindref="member">m_impl</ref>-&gt;clamps;<sp/>}</highlight></codeline>
<codeline lineno="365"><highlight class="normal">}</highlight></codeline>
    </programlisting>
    <location file="src/lib/serif/backend/mfem/serialize/top.cpp"/>
  </compounddef>
</doxygen>