831 lines
63 KiB
XML
831 lines
63 KiB
XML
|
<?xml version='1.0' encoding='UTF-8' standalone='no'?>
|
||
|
|
<doxygen xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="compound.xsd" version="1.13.2" xml:lang="en-US">
|
||
|
|
<compounddef id="mass__diagonal_8cpp" kind="file" language="C++">
|
||
|
|
<compoundname>mass_diagonal.cpp</compoundname>
|
||
|
|
<includes refid="mass__diagonal_8hpp" local="yes">serif/backend/mfem/discretization/mass_diagonal.hpp</includes>
|
||
|
|
<includes local="no">cmath</includes>
|
||
|
|
<includes local="no">memory</includes>
|
||
|
|
<includes local="no">algorithm</includes>
|
||
|
|
<includes local="no">utility</includes>
|
||
|
|
<includes local="no">vector</includes>
|
||
|
|
<includes local="yes">mfem/fem/bilinearform.hpp</includes>
|
||
|
|
<includes local="yes">mfem/fem/bilininteg.hpp</includes>
|
||
|
|
<includes local="yes">mfem/general/forall.hpp</includes>
|
||
|
|
<includes refid="errors_8hpp" local="yes">serif/utils/error/errors.hpp</includes>
|
||
|
|
<includes refid="element__coordinates_8hpp" local="yes">serif/backend/mfem/discretization/element_coordinates.hpp</includes>
|
||
|
|
<includes refid="parallel_8hpp" local="yes">serif/backend/mfem/shared/parallel.hpp</includes>
|
||
|
|
<incdepgraph>
|
||
|
|
<node id="58">
|
||
|
|
<label>serif/backend/mfem/discretization/dof_map.hpp</label>
|
||
|
|
<link refid="dof__map_8hpp"/>
|
||
|
|
<childnode refid="59" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="60" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="61" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="62" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="78" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="79" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="80" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="65" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="81" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="85">
|
||
|
|
<label>serif/backend/mfem/discretization/element_coordinates.hpp</label>
|
||
|
|
<link refid="element__coordinates_8hpp"/>
|
||
|
|
<childnode refid="71" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="86" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="87" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="58" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="2">
|
||
|
|
<label>serif/backend/mfem/discretization/mass_diagonal.hpp</label>
|
||
|
|
<link refid="mass__diagonal_8hpp"/>
|
||
|
|
<childnode refid="3" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="4" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="5" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="6" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="62">
|
||
|
|
<label>serif/backend/mfem/discretization/reconstruction.hpp</label>
|
||
|
|
<link refid="reconstruction_8hpp"/>
|
||
|
|
<childnode refid="61" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="5" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="63" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="11" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="64" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="6">
|
||
|
|
<label>serif/backend/mfem/discretization/term_binding.hpp</label>
|
||
|
|
<link refid="term__binding_8hpp"/>
|
||
|
|
<childnode refid="7" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="8" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="37" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="58" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="65" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="88">
|
||
|
|
<label>serif/backend/mfem/shared/parallel.hpp</label>
|
||
|
|
<link refid="parallel_8hpp"/>
|
||
|
|
<childnode refid="74" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="5" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="14">
|
||
|
|
<label>serif/dimensions/dimensions.hpp</label>
|
||
|
|
<link refid="dimensions_8hpp"/>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="15" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="15">
|
||
|
|
<label>serif/dimensions/quantities.hpp</label>
|
||
|
|
<link refid="quantities_8hpp"/>
|
||
|
|
<childnode refid="9" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="16" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="17" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="18" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="19" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="20" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="21">
|
||
|
|
<label>serif/discretization/blocks/base.hpp</label>
|
||
|
|
<link refid="base_8hpp"/>
|
||
|
|
</node>
|
||
|
|
<node id="10">
|
||
|
|
<label>serif/discretization/blocks/concepts.hpp</label>
|
||
|
|
<link refid="discretization_2blocks_2concepts_8hpp"/>
|
||
|
|
<childnode refid="11" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="11">
|
||
|
|
<label>serif/discretization/blocks/fields.hpp</label>
|
||
|
|
<link refid="fields_8hpp"/>
|
||
|
|
<childnode refid="12" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="14" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="21" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="22" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="27" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="33" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="34" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="35" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="8">
|
||
|
|
<label>serif/discretization/blocks/traits.hpp</label>
|
||
|
|
<link refid="discretization_2blocks_2traits_8hpp"/>
|
||
|
|
<childnode refid="9" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="10" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="27">
|
||
|
|
<label>serif/discretization/domain/concepts.hpp</label>
|
||
|
|
<link refid="discretization_2domain_2concepts_8hpp"/>
|
||
|
|
<childnode refid="28" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="29" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="12" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="23" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="30" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="31" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="19" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="32" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="45">
|
||
|
|
<label>serif/discretization/domain/ids/boundary.hpp</label>
|
||
|
|
<link refid="ids_2boundary_8hpp"/>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="27" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="44">
|
||
|
|
<label>serif/discretization/domain/ids/concepts.hpp</label>
|
||
|
|
<link refid="discretization_2domain_2ids_2concepts_8hpp"/>
|
||
|
|
<childnode refid="45" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="46" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="46">
|
||
|
|
<label>serif/discretization/domain/ids/domain.hpp</label>
|
||
|
|
<link refid="domain_8hpp"/>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="27" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="40">
|
||
|
|
<label>serif/discretization/domain/ids/lists/concepts.hpp</label>
|
||
|
|
<link refid="discretization_2domain_2ids_2lists_2concepts_8hpp"/>
|
||
|
|
<childnode refid="41" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="42">
|
||
|
|
<label>serif/discretization/domain/ids/lists/lists.hpp</label>
|
||
|
|
<link refid="lists_8hpp"/>
|
||
|
|
<childnode refid="12" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="28" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="43" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="32" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="41">
|
||
|
|
<label>serif/discretization/domain/ids/lists/traits.hpp</label>
|
||
|
|
<link refid="discretization_2domain_2ids_2lists_2traits_8hpp"/>
|
||
|
|
<childnode refid="42" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="43">
|
||
|
|
<label>serif/discretization/domain/ids/validation/concepts.hpp</label>
|
||
|
|
<link refid="discretization_2domain_2ids_2validation_2concepts_8hpp"/>
|
||
|
|
<childnode refid="44" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="47" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="48" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="47">
|
||
|
|
<label>serif/discretization/domain/ids/validation/unique_ids.hpp</label>
|
||
|
|
<link refid="unique__ids_8hpp"/>
|
||
|
|
<childnode refid="12" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="17" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="9" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="32" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="33">
|
||
|
|
<label>serif/discretization/domain/physical_domains.hpp</label>
|
||
|
|
<link refid="physical__domains_8hpp"/>
|
||
|
|
<childnode refid="31" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="27" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="53">
|
||
|
|
<label>serif/discretization/domain/relation/concepts.hpp</label>
|
||
|
|
<link refid="discretization_2domain_2relation_2concepts_8hpp"/>
|
||
|
|
<childnode refid="54" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="19" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="50">
|
||
|
|
<label>serif/discretization/domain/relation/lists/concepts.hpp</label>
|
||
|
|
<link refid="discretization_2domain_2relation_2lists_2concepts_8hpp"/>
|
||
|
|
<childnode refid="51" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="52">
|
||
|
|
<label>serif/discretization/domain/relation/lists/relation_list.hpp</label>
|
||
|
|
<link refid="relation__list_8hpp"/>
|
||
|
|
<childnode refid="28" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="53" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="32" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="51">
|
||
|
|
<label>serif/discretization/domain/relation/lists/traits.hpp</label>
|
||
|
|
<link refid="discretization_2domain_2relation_2lists_2traits_8hpp"/>
|
||
|
|
<childnode refid="52" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="54">
|
||
|
|
<label>serif/discretization/domain/relation/relations.hpp</label>
|
||
|
|
<link refid="discretization_2domain_2relation_2relations_8hpp"/>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="27" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="20" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="32" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="55">
|
||
|
|
<label>serif/discretization/domain/relation/validation/concepts.hpp</label>
|
||
|
|
<link refid="discretization_2domain_2relation_2validation_2concepts_8hpp"/>
|
||
|
|
<childnode refid="56" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="56">
|
||
|
|
<label>serif/discretization/domain/relation/validation/validation.hpp</label>
|
||
|
|
<link refid="discretization_2domain_2relation_2validation_2validation_8hpp"/>
|
||
|
|
<childnode refid="9" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="27" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="40" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="53" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="50" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="52" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="57" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="57">
|
||
|
|
<label>serif/discretization/domain/resolver.hpp</label>
|
||
|
|
<link refid="resolver_8hpp"/>
|
||
|
|
<childnode refid="27" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="40" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="42" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="37">
|
||
|
|
<label>serif/discretization/domain/schema/concepts.hpp</label>
|
||
|
|
<link refid="discretization_2domain_2schema_2concepts_8hpp"/>
|
||
|
|
<childnode refid="38" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="39">
|
||
|
|
<label>serif/discretization/domain/schema/domain_schema.hpp</label>
|
||
|
|
<link refid="domain__schema_8hpp"/>
|
||
|
|
<childnode refid="28" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="27" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="40" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="50" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="55" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="57" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="38">
|
||
|
|
<label>serif/discretization/domain/schema/traits.hpp</label>
|
||
|
|
<link refid="discretization_2domain_2schema_2traits_8hpp"/>
|
||
|
|
<childnode refid="39" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="31">
|
||
|
|
<label>serif/discretization/domain/types.hpp</label>
|
||
|
|
<link refid="discretization_2domain_2types_8hpp"/>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="22">
|
||
|
|
<label>serif/discretization/functional_space/space.hpp</label>
|
||
|
|
<link refid="space_8hpp"/>
|
||
|
|
<childnode refid="9" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="23" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="24" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="64">
|
||
|
|
<label>serif/discretization/reconstruction/surface_power.hpp</label>
|
||
|
|
<link refid="surface__power_8hpp"/>
|
||
|
|
<childnode refid="59" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="29" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="33" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="65" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="24">
|
||
|
|
<label>serif/meta/enums.hpp</label>
|
||
|
|
<link refid="enums_8hpp"/>
|
||
|
|
<childnode refid="25" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="9" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="26" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="35">
|
||
|
|
<label>serif/meta/names.hpp</label>
|
||
|
|
<link refid="names_8hpp"/>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="36" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="75">
|
||
|
|
<label>serif/utils/error/codes.hpp</label>
|
||
|
|
<link refid="codes_8hpp"/>
|
||
|
|
<childnode refid="67" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="24" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="65">
|
||
|
|
<label>serif/utils/error/errors.hpp</label>
|
||
|
|
<link refid="errors_8hpp"/>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="23" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="66" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="4" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="67" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="7" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="68" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="69" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="73" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="75" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="76" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="81">
|
||
|
|
<label>serif/utils/error/nonfinite.hpp</label>
|
||
|
|
<link refid="nonfinite_8hpp"/>
|
||
|
|
<childnode refid="78" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="73">
|
||
|
|
<label>serif/utils/misc/concepts/exceptions.hpp</label>
|
||
|
|
<link refid="utils_2misc_2concepts_2exceptions_8hpp"/>
|
||
|
|
<childnode refid="74" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="9" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="20" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="20">
|
||
|
|
<label>serif/utils/misc/concepts/numeric.hpp</label>
|
||
|
|
<link refid="numeric_8hpp"/>
|
||
|
|
<childnode refid="17" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="9" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="48">
|
||
|
|
<label>serif/utils/misc/concepts/type_uniqueness.hpp</label>
|
||
|
|
<link refid="concepts_2type__uniqueness_8hpp"/>
|
||
|
|
<childnode refid="49" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="19">
|
||
|
|
<label>serif/utils/misc/std_helper/cleaning.hpp</label>
|
||
|
|
<link refid="cleaning_8hpp"/>
|
||
|
|
<childnode refid="9" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="32">
|
||
|
|
<label>serif/utils/misc/std_helper/variadic_cleaning.hpp</label>
|
||
|
|
<link refid="variadic__cleaning_8hpp"/>
|
||
|
|
<childnode refid="28" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="69">
|
||
|
|
<label>serif/utils/misc/terminal/colors.hpp</label>
|
||
|
|
<link refid="colors_8hpp"/>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="66" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="70" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="71" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="72" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="49">
|
||
|
|
<label>serif/utils/misc/traits/type_uniqueness.hpp</label>
|
||
|
|
<link refid="traits_2type__uniqueness_8hpp"/>
|
||
|
|
<childnode refid="12" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="19" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="76">
|
||
|
|
<label>serif/utils/types/strings.hpp</label>
|
||
|
|
<link refid="types_2strings_8hpp"/>
|
||
|
|
<childnode refid="7" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="77" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="13" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="34">
|
||
|
|
<label>serif/utils/types/type_list.hpp</label>
|
||
|
|
<link refid="type__list_8hpp"/>
|
||
|
|
</node>
|
||
|
|
<node id="1">
|
||
|
|
<label>src/lib/serif/backend/mfem/discretization/mass_diagonal.cpp</label>
|
||
|
|
<link refid="mass__diagonal_8cpp"/>
|
||
|
|
<childnode refid="2" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="59" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="61" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="7" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="71" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="5" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="82" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="83" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="84" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="65" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="85" relation="include">
|
||
|
|
</childnode>
|
||
|
|
<childnode refid="88" relation="include">
|
||
|
|
</childnode>
|
||
|
|
</node>
|
||
|
|
<node id="7">
|
||
|
|
<label>algorithm</label>
|
||
|
|
</node>
|
||
|
|
<node id="12">
|
||
|
|
<label>array</label>
|
||
|
|
</node>
|
||
|
|
<node id="59">
|
||
|
|
<label>cmath</label>
|
||
|
|
</node>
|
||
|
|
<node id="16">
|
||
|
|
<label>compare</label>
|
||
|
|
</node>
|
||
|
|
<node id="17">
|
||
|
|
<label>concepts</label>
|
||
|
|
</node>
|
||
|
|
<node id="28">
|
||
|
|
<label>cstddef</label>
|
||
|
|
</node>
|
||
|
|
<node id="67">
|
||
|
|
<label>cstdint</label>
|
||
|
|
</node>
|
||
|
|
<node id="66">
|
||
|
|
<label>cstdio</label>
|
||
|
|
</node>
|
||
|
|
<node id="70">
|
||
|
|
<label>cstdlib</label>
|
||
|
|
</node>
|
||
|
|
<node id="74">
|
||
|
|
<label>exception</label>
|
||
|
|
</node>
|
||
|
|
<node id="3">
|
||
|
|
<label>expected</label>
|
||
|
|
</node>
|
||
|
|
<node id="23">
|
||
|
|
<label>format</label>
|
||
|
|
</node>
|
||
|
|
<node id="77">
|
||
|
|
<label>iostream</label>
|
||
|
|
</node>
|
||
|
|
<node id="26">
|
||
|
|
<label>magic_enum/magic_enum_all.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="61">
|
||
|
|
<label>memory</label>
|
||
|
|
</node>
|
||
|
|
<node id="82">
|
||
|
|
<label>mfem/fem/bilinearform.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="83">
|
||
|
|
<label>mfem/fem/bilininteg.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="87">
|
||
|
|
<label>mfem/fem/doftrans.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="63">
|
||
|
|
<label>mfem/fem/fe/fe_base.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="79">
|
||
|
|
<label>mfem/fem/fespace.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="86">
|
||
|
|
<label>mfem/fem/geom.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="80">
|
||
|
|
<label>mfem/fem/pfespace.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="84">
|
||
|
|
<label>mfem/general/forall.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="36">
|
||
|
|
<label>nameof.hpp</label>
|
||
|
|
</node>
|
||
|
|
<node id="30">
|
||
|
|
<label>numeric</label>
|
||
|
|
</node>
|
||
|
|
<node id="60">
|
||
|
|
<label>optional</label>
|
||
|
|
</node>
|
||
|
|
<node id="29">
|
||
|
|
<label>span</label>
|
||
|
|
</node>
|
||
|
|
<node id="78">
|
||
|
|
<label>stdexcept</label>
|
||
|
|
</node>
|
||
|
|
<node id="4">
|
||
|
|
<label>string</label>
|
||
|
|
</node>
|
||
|
|
<node id="13">
|
||
|
|
<label>string_view</label>
|
||
|
|
</node>
|
||
|
|
<node id="18">
|
||
|
|
<label>tuple</label>
|
||
|
|
</node>
|
||
|
|
<node id="9">
|
||
|
|
<label>type_traits</label>
|
||
|
|
</node>
|
||
|
|
<node id="72">
|
||
|
|
<label>unistd.h</label>
|
||
|
|
</node>
|
||
|
|
<node id="25">
|
||
|
|
<label>unordered_map</label>
|
||
|
|
</node>
|
||
|
|
<node id="71">
|
||
|
|
<label>utility</label>
|
||
|
|
</node>
|
||
|
|
<node id="5">
|
||
|
|
<label>vector</label>
|
||
|
|
</node>
|
||
|
|
<node id="68">
|
||
|
|
<label>version</label>
|
||
|
|
</node>
|
||
|
|
</incdepgraph>
|
||
|
|
<innernamespace refid="namespaceserif">serif</innernamespace>
|
||
|
|
<innernamespace refid="namespaceserif_1_1discretization">serif::discretization</innernamespace>
|
||
|
|
<innernamespace refid="namespaceserif_1_1discretization_1_1backend">serif::discretization::backend</innernamespace>
|
||
|
|
<innernamespace refid="namespaceserif_1_1discretization_1_1backend_1_1mfem">serif::discretization::backend::mfem</innernamespace>
|
||
|
|
<briefdescription>
|
||
|
|
</briefdescription>
|
||
|
|
<detaileddescription>
|
||
|
|
</detaileddescription>
|
||
|
|
<programlisting>
|
||
|
|
<codeline lineno="1"><highlight class="preprocessor">#include<sp/>"<ref refid="mass__diagonal_8hpp" kindref="compound">serif/backend/mfem/discretization/mass_diagonal.hpp</ref>"</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="2"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="3"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/><cmath></highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="4"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/><memory></highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="5"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/><algorithm></highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="6"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/><utility></highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="7"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/><vector></highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="8"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="9"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>"mfem/fem/bilinearform.hpp"</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="10"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>"mfem/fem/bilininteg.hpp"</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="11"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>"mfem/general/forall.hpp"</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="12"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>"<ref refid="errors_8hpp" kindref="compound">serif/utils/error/errors.hpp</ref>"</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="13"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>"<ref refid="element__coordinates_8hpp" kindref="compound">serif/backend/mfem/discretization/element_coordinates.hpp</ref>"</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="14"><highlight class="normal"></highlight><highlight class="preprocessor">#include<sp/>"<ref refid="parallel_8hpp" kindref="compound">serif/backend/mfem/shared/parallel.hpp</ref>"</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="15"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="16"><highlight class="normal"></highlight><highlight class="keyword">namespace<sp/></highlight><highlight class="normal"><ref refid="namespaceserif_1_1discretization_1_1backend_1_1mfem" kindref="compound">serif::discretization::backend::mfem</ref><sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="17"><highlight class="normal"><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">namespace<sp/></highlight><highlight class="normal">{</highlight></codeline>
|
||
|
|
<codeline lineno="25"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>[[nodiscard]]<sp/>::mfem::Vector<sp/>scalar_mass_diagonal(</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">DofMap</ref>&<sp/>binding,<sp/>const<sp/>::mfem::Array<int>&<sp/>marker)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="26"><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/>elements<sp/>=<sp/>binding.validated_elements();</highlight></codeline>
|
||
|
|
<codeline lineno="27"><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/>mesh<sp/>=<sp/>binding.mesh();</highlight></codeline>
|
||
|
|
<codeline lineno="28"><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/>plan<sp/>=<sp/>binding.reconstruction().get();</highlight></codeline>
|
||
|
|
<codeline lineno="29"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">struct<sp/></highlight><highlight class="normal">ElementDiagonal<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/>std::vector<int><sp/>indices;</highlight></codeline>
|
||
|
|
<codeline lineno="31"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::vector<double><sp/><ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1adf1d7a29193c8b8a6cd17745686a04d8" kindref="member">values</ref>;</highlight></codeline>
|
||
|
|
<codeline lineno="32"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
|
||
|
|
<codeline lineno="33"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">struct<sp/></highlight><highlight class="normal">Workspace<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="34"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::IsoparametricTransformation<sp/>transformation;</highlight></codeline>
|
||
|
|
<codeline lineno="35"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::Vector<sp/>shape,<sp/>weighted;</highlight></codeline>
|
||
|
|
<codeline lineno="36"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::DenseMatrix<sp/>derivatives;</highlight></codeline>
|
||
|
|
<codeline lineno="37"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::vector<int><sp/>local_of_reduced;</highlight></codeline>
|
||
|
|
<codeline lineno="38"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>};</highlight></codeline>
|
||
|
|
<codeline lineno="39"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::vector<ElementDiagonal><sp/>results(</highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(mesh.GetNE()));</highlight></codeline>
|
||
|
|
<codeline lineno="40"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::vector<Workspace><sp/>workspaces(</highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(<ref refid="namespaceserif_1_1utils_1_1parallel_1a845c655da35a46d0be87880f0e4fa2d9" kindref="member">utils::parallel::worker_count</ref>()));</highlight></codeline>
|
||
|
|
<codeline lineno="41"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1parallel_1a453320d288c9e42300bdc59ee8dccd2d" kindref="member">utils::parallel::for_each_index</ref>(mesh.GetNE(),<sp/>[&](</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>element,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>worker)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="42"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>if<sp/>(marker[mesh.GetAttribute(element)<sp/>-<sp/>1]<sp/>==<sp/>0)<sp/>{<sp/>return;<sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="43"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>Workspace&<sp/>work<sp/>=<sp/>workspaces[</highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(worker)];</highlight></codeline>
|
||
|
|
<codeline lineno="44"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>work.local_of_reduced.resize(</highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(binding.size()),<sp/>-1);</highlight></codeline>
|
||
|
|
<codeline lineno="45"><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">&<sp/>finite_element<sp/>=<sp/>*binding.space().GetFE(element);</highlight></codeline>
|
||
|
|
<codeline lineno="46"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><ref refid="namespaceserif_1_1utils_1_1parallel_1a8e0e78d02ab6c49ae588644e579bed79" kindref="member">utils::parallel::element_transformation</ref>(mesh,<sp/>element,<sp/>work.transformation);</highlight></codeline>
|
||
|
|
<codeline lineno="47"><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"><sp/>coordinates<sp/>=<sp/>elements.coordinates(element);</highlight></codeline>
|
||
|
|
<codeline lineno="48"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>ElementDiagonal&<sp/>result<sp/>=<sp/>results[</highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(element)];</highlight></codeline>
|
||
|
|
<codeline lineno="49"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::vector<std::vector<std::pair<int,<sp/>double>>><sp/>constraint_rows(</highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(finite_element.GetDof()));</highlight></codeline>
|
||
|
|
<codeline lineno="50"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>for<sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>node<sp/>=<sp/>0;<sp/>node<sp/><<sp/>finite_element.GetDof();<sp/>++node)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="51"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>coordinates.for_each_coefficient(node,<sp/>[&](</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>reduced,<sp/></highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">double</highlight><highlight class="normal"><sp/>weight)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="52"><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">int</highlight><highlight class="normal">&<sp/>local<sp/>=<sp/>work.local_of_reduced[</highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(reduced)];</highlight></codeline>
|
||
|
|
<codeline lineno="53"><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/>(local<sp/><<sp/>0)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="54"><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/>local<sp/>=<sp/></highlight><highlight class="keyword">static_cast<</highlight><highlight class="keywordtype">int</highlight><highlight class="keyword">></highlight><highlight class="normal">(result.indices.size());</highlight></codeline>
|
||
|
|
<codeline lineno="55"><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/>result.indices.push_back(reduced);</highlight></codeline>
|
||
|
|
<codeline lineno="56"><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="57"><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/>constraint_rows[</highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(node)].emplace_back(local,<sp/>weight);</highlight></codeline>
|
||
|
|
<codeline lineno="58"><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="59"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="60"><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="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>reduced<sp/>:<sp/>result.indices)<sp/>{<sp/>work.local_of_reduced[static_cast<std::size_t>(reduced)]<sp/>=<sp/>-1;<sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="61"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>result.values.assign(result.indices.size(),<sp/>0.0);</highlight></codeline>
|
||
|
|
<codeline lineno="62"><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">*<sp/>reconstruction<sp/>=<sp/>plan<sp/>!=<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal"><sp/>?<sp/>binding.reconstruction_element(element)<sp/>:<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">;</highlight></codeline>
|
||
|
|
<codeline lineno="63"><highlight class="normal"><sp/><sp/><sp/><sp/><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/>::mfem::MassIntegrator::GetRule(finite_element,<sp/>finite_element,<sp/>work.transformation).GetOrder();</highlight></codeline>
|
||
|
|
<codeline lineno="64"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>if<sp/>(reconstruction<sp/>!=<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">)<sp/>{<sp/>order<sp/>+=<sp/>plan->integration_order_increment();<sp/>}</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="keyword">const</highlight><highlight class="normal"><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal">&<sp/>rule<sp/>=<sp/>::mfem::IntRules.Get(finite_element.GetGeomType(),<sp/>order);</highlight></codeline>
|
||
|
|
<codeline lineno="66"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::vector<double><sp/>combined(result.indices.size());</highlight></codeline>
|
||
|
|
<codeline lineno="67"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>work.shape.SetSize(finite_element.GetDof());</highlight></codeline>
|
||
|
|
<codeline lineno="68"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>for<sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>point<sp/>=<sp/>0;<sp/>point<sp/><<sp/>rule.GetNPoints();<sp/>++point)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="69"><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">&<sp/>location<sp/>=<sp/>rule.IntPoint(point);</highlight></codeline>
|
||
|
|
<codeline lineno="70"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>const<sp/>::mfem::Vector*<sp/>shape<sp/>=<sp/>&work.shape;</highlight></codeline>
|
||
|
|
<codeline lineno="71"><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/>(reconstruction<sp/>!=<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="72"><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/>reconstruction->sample(finite_element,<sp/>location,<sp/>work.weighted,<sp/>work.derivatives);</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/><sp/><sp/><sp/><sp/>shape<sp/>=<sp/>&work.weighted;</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/></highlight><highlight class="keywordflow">else</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/><sp/><sp/><sp/>finite_element.CalcShape(location,<sp/>work.shape);</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/>}</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/>std::fill(combined.begin(),<sp/>combined.end(),<sp/>0.0);</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/></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/><<sp/>finite_element.GetDof();<sp/>++node)<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/></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">&<sp/>[local,<sp/>coefficient]<sp/>:<sp/>constraint_rows[</highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(node)])<sp/>{</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/><sp/><sp/><sp/><sp/>combined[</highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(local)]<sp/>+=<sp/>coefficient<sp/>*<sp/>(*shape)[node];</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></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/><sp/><sp/><sp/><sp/>work.transformation.SetIntPoint(&location);</highlight></codeline>
|
||
|
|
<codeline lineno="84"><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/>weight<sp/>=<sp/>location.weight<sp/>*<sp/>work.transformation.Weight();</highlight></codeline>
|
||
|
|
<codeline lineno="85"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(std::size_t<sp/>local<sp/>=<sp/>0;<sp/>local<sp/><<sp/>combined.size();<sp/>++local)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="86"><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.values[local]<sp/>+=<sp/>weight<sp/>*<sp/>combined[local]<sp/>*<sp/>combined[local];</highlight></codeline>
|
||
|
|
<codeline lineno="87"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="88"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="89"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>});</highlight></codeline>
|
||
|
|
<codeline lineno="90"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::Vector<sp/>diagonal(binding.size());</highlight></codeline>
|
||
|
|
<codeline lineno="91"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>diagonal<sp/>=<sp/>0.0;</highlight></codeline>
|
||
|
|
<codeline lineno="92"><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">&<sp/>result<sp/>:<sp/>results)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="93"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(std::size_t<sp/>k<sp/>=<sp/>0;<sp/>k<sp/><<sp/>result.indices.size();<sp/>++k)<sp/>{<sp/>diagonal[result.indices[k]]<sp/>+=<sp/>result.values[k];<sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="94"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="95"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>diagonal;</highlight></codeline>
|
||
|
|
<codeline lineno="96"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="97"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="98"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="99"><highlight class="normal"><sp/><sp/><sp/><sp/>std::expected<std::vector<double>,<sp/>std::string><sp/>detail::reference_mass_diagonal(</highlight><highlight class="keyword">const</highlight><highlight class="normal"><sp/>DofMap&<sp/>binding)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="100"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">using<sp/>namespace<sp/></highlight><highlight class="normal">utils::errors;</highlight></codeline>
|
||
|
|
<codeline lineno="101"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="102"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>binding.validate();</highlight></codeline>
|
||
|
|
<codeline lineno="103"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(binding.size()<sp/>==<sp/>0)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="104"><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::vector<double>{};</highlight></codeline>
|
||
|
|
<codeline lineno="105"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="106"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="107"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::FiniteElementSpace&<sp/>space<sp/>=<sp/>binding.space();</highlight></codeline>
|
||
|
|
<codeline lineno="108"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::Array<int><sp/>marker(binding.mesh().attributes.Max());</highlight></codeline>
|
||
|
|
<codeline lineno="109"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>marker<sp/>=<sp/>0;</highlight></codeline>
|
||
|
|
<codeline lineno="110"><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/>element<sp/>=<sp/>0;<sp/>element<sp/><<sp/>space.GetNE();<sp/>++element)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="111"><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.supports(element))<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="112"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>marker[binding.mesh().GetAttribute(element)<sp/>-<sp/>1]<sp/>=<sp/>1;</highlight></codeline>
|
||
|
|
<codeline lineno="113"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="114"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="115"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="116"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(space.GetVDim()<sp/>==<sp/>1)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="117"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>The<sp/>scalar<sp/>fields<sp/>(density,<sp/>potential)<sp/>need<sp/>only<sp/>diag(P^T<sp/>M<sp/>P);<sp/>assembling<sp/>M<sp/>and<sp/>P^T<sp/>M<sp/>P<sp/>to<sp/>read<sp/>it</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="118"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>cost<sp/>6<sp/>s<sp/>of<sp/>a<sp/>19<sp/>s<sp/>benchmark<sp/>solve.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="119"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>const<sp/>::mfem::Vector<sp/><ref refid="namespaceserif_1_1geometry_1_1surface_1a739898ab5e92990aaeb57582842b1d1ca42bd4c0819c498d1c1ab622de74d6f82" kindref="member">supported</ref><sp/>=<sp/>scalar_mass_diagonal(binding,<sp/>marker);</highlight></codeline>
|
||
|
|
<codeline lineno="120"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::vector<double><sp/>result(</highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(<ref refid="namespaceserif_1_1geometry_1_1surface_1a739898ab5e92990aaeb57582842b1d1ca42bd4c0819c498d1c1ab622de74d6f82" kindref="member">supported</ref>.Size()));</highlight></codeline>
|
||
|
|
<codeline lineno="121"><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/>index<sp/>=<sp/>0;<sp/><ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1ab297c6eedd8df878b6af2135b05e4616" kindref="member">index</ref><sp/><<sp/><ref refid="namespaceserif_1_1geometry_1_1surface_1a739898ab5e92990aaeb57582842b1d1ca42bd4c0819c498d1c1ab622de74d6f82" kindref="member">supported</ref>.Size();<sp/>++<ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1ab297c6eedd8df878b6af2135b05e4616" kindref="member">index</ref>)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="122"><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(supported[index])<sp/>||<sp/>supported[index]<sp/><=<sp/>0.0)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="123"><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(format_error_message<SERiFErrorCode::NORMALIZATION_ERROR>(</highlight><highlight class="stringliteral">"Reference<sp/>mass<sp/>diagonal<sp/>entry<sp/>{}<sp/>is<sp/>{},<sp/>but<sp/>every<sp/>supported<sp/>coefficient<sp/>requires<sp/>positive<sp/>finite<sp/>entries.<sp/>Check<sp/>the<sp/>reference<sp/>mesh<sp/>geometry<sp/>and<sp/>finite-element<sp/>support."</highlight><highlight class="normal">,<sp/>index,<sp/>supported[index]));</highlight></codeline>
|
||
|
|
<codeline lineno="124"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="125"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>result[</highlight><highlight class="keyword">static_cast<</highlight><highlight class="normal">std::size_t</highlight><highlight class="keyword">></highlight><highlight class="normal">(<ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1ab297c6eedd8df878b6af2135b05e4616" kindref="member">index</ref>)]<sp/>=<sp/>supported[index];</highlight></codeline>
|
||
|
|
<codeline lineno="126"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="127"><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="128"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="129"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="130"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>serif_ensure<SERiFErrorCode::NORMALIZATION_ERROR>(binding.reconstruction()<sp/>==<sp/></highlight><highlight class="keyword">nullptr</highlight><highlight class="normal">,</highlight></codeline>
|
||
|
|
<codeline lineno="131"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="stringliteral">"Reconstructed<sp/>reference<sp/>metrics<sp/>currently<sp/>require<sp/>a<sp/>scalar<sp/>H1<sp/>density<sp/>field."</highlight><highlight class="normal">);</highlight></codeline>
|
||
|
|
<codeline lineno="132"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="comment">//<sp/>Vector<sp/>H1<sp/>fields<sp/>assemble<sp/>their<sp/>componentwise<sp/>physical<sp/>inner<sp/>product.</highlight><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="133"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::BilinearForm<sp/>mass(&space);</highlight></codeline>
|
||
|
|
<codeline lineno="134"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keyword">auto</highlight><highlight class="normal"><sp/>integrator<sp/>=<sp/>std::make_unique<::mfem::VectorMassIntegrator>();</highlight></codeline>
|
||
|
|
<codeline lineno="135"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>integrator->SetVDim(space.GetVDim());</highlight></codeline>
|
||
|
|
<codeline lineno="136"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mass.AddDomainIntegrator(integrator.release(),<sp/>marker);</highlight></codeline>
|
||
|
|
<codeline lineno="137"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mass.Assemble();</highlight></codeline>
|
||
|
|
<codeline lineno="138"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mass.Finalize();</highlight></codeline>
|
||
|
|
<codeline lineno="139"><highlight class="normal"></highlight></codeline>
|
||
|
|
<codeline lineno="140"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">if</highlight><highlight class="normal"><sp/>(binding.constrained())<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="141"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>const<sp/>::mfem::Array<int><sp/>essential_dofs;</highlight></codeline>
|
||
|
|
<codeline lineno="142"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::OperatorHandle<sp/>system;</highlight></codeline>
|
||
|
|
<codeline lineno="143"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mass.FormSystemMatrix(essential_dofs,<sp/>system);</highlight></codeline>
|
||
|
|
<codeline lineno="144"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="145"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::Vector<sp/>diagonal(binding.constrained()<sp/>?<sp/>space.GetTrueVSize()<sp/>:<sp/>space.GetVSize());</highlight></codeline>
|
||
|
|
<codeline lineno="146"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>mass.AssembleDiagonal(diagonal);</highlight></codeline>
|
||
|
|
<codeline lineno="147"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>::mfem::Vector<sp/><ref refid="namespaceserif_1_1geometry_1_1surface_1a739898ab5e92990aaeb57582842b1d1ca42bd4c0819c498d1c1ab622de74d6f82" kindref="member">supported</ref><sp/>=<sp/>binding.constrained()<sp/>?<sp/>binding.gather_true(diagonal)<sp/>:<sp/>binding.gather(diagonal);</highlight></codeline>
|
||
|
|
<codeline lineno="148"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>std::vector<double><sp/>result(<ref refid="namespaceserif_1_1geometry_1_1surface_1a739898ab5e92990aaeb57582842b1d1ca42bd4c0819c498d1c1ab622de74d6f82" kindref="member">supported</ref>.Size());</highlight></codeline>
|
||
|
|
<codeline lineno="149"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">for</highlight><highlight class="normal"><sp/>(</highlight><highlight class="keywordtype">int</highlight><highlight class="normal"><sp/>index<sp/>=<sp/>0;<sp/><ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1ab297c6eedd8df878b6af2135b05e4616" kindref="member">index</ref><sp/><<sp/><ref refid="namespaceserif_1_1geometry_1_1surface_1a739898ab5e92990aaeb57582842b1d1ca42bd4c0819c498d1c1ab622de74d6f82" kindref="member">supported</ref>.Size();<sp/>++<ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1ab297c6eedd8df878b6af2135b05e4616" kindref="member">index</ref>)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="150"><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(supported[index])<sp/>||<sp/>supported[index]<sp/><=<sp/>0.0)<sp/>{</highlight></codeline>
|
||
|
|
<codeline lineno="151"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>std::unexpected(format_error_message<SERiFErrorCode::NORMALIZATION_ERROR>(</highlight><highlight class="stringliteral">"Reference<sp/>mass<sp/>diagonal<sp/>entry<sp/>{}<sp/>is<sp/>{},<sp/>but<sp/>every<sp/>supported<sp/>coefficient<sp/>requires<sp/>positive<sp/>finite<sp/>entries.<sp/>Check<sp/>the<sp/>reference<sp/>mesh<sp/>geometry<sp/>and<sp/>finite-element<sp/>support."</highlight><highlight class="normal">,<sp/>index,<sp/>supported[index]));</highlight></codeline>
|
||
|
|
<codeline lineno="152"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="153"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>result[<ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1ab297c6eedd8df878b6af2135b05e4616" kindref="member">index</ref>]<sp/>=<sp/><ref refid="namespaceserif_1_1geometry_1_1surface_1a739898ab5e92990aaeb57582842b1d1ca42bd4c0819c498d1c1ab622de74d6f82" kindref="member">supported</ref>[<ref refid="namespaceserif_1_1operators_1_1backend_1_1mfem_1ab297c6eedd8df878b6af2135b05e4616" kindref="member">index</ref>];</highlight></codeline>
|
||
|
|
<codeline lineno="154"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="155"><highlight class="normal"><sp/><sp/><sp/><sp/><sp/><sp/><sp/><sp/></highlight><highlight class="keywordflow">return</highlight><highlight class="normal"><sp/>result;</highlight></codeline>
|
||
|
|
<codeline lineno="156"><highlight class="normal"><sp/><sp/><sp/><sp/>}</highlight></codeline>
|
||
|
|
<codeline lineno="157"><highlight class="normal">}</highlight></codeline>
|
||
|
|
</programlisting>
|
||
|
|
<location file="src/lib/serif/backend/mfem/discretization/mass_diagonal.cpp"/>
|
||
|
|
</compounddef>
|
||
|
|
</doxygen>
|