basis_evaluation.hpp memory mfem/fem/fe/fe_base.hpp mfem/fem/intrules.hpp serif/backend/mfem/discretization/tensor_basis.hpp src/include/serif/backend/mfem/operators/element_data.hpp src/include/serif/backend/mfem/operators/operand_evaluation.hpp src/include/serif/backend/mfem/operators/workspace.hpp src/lib/serif/backend/mfem/discretization/basis_evaluation.cpp src/include/serif/backend/mfem/discretization/basis_evaluation.hpp serif/backend/mfem/discretization/reconstruction.hpp serif/backend/mfem/discretization/tensor_basis.hpp serif/dimensions/dimensions.hpp serif/dimensions/quantities.hpp serif/discretization/blocks/base.hpp serif/discretization/blocks/fields.hpp serif/discretization/domain/concepts.hpp serif/discretization/domain/physical_domains.hpp serif/discretization/domain/types.hpp serif/discretization/functional_space/space.hpp serif/discretization/reconstruction/surface_power.hpp serif/meta/enums.hpp serif/meta/names.hpp serif/utils/error/codes.hpp serif/utils/error/errors.hpp serif/utils/misc/concepts/exceptions.hpp serif/utils/misc/concepts/numeric.hpp serif/utils/misc/std_helper/cleaning.hpp serif/utils/misc/std_helper/variadic_cleaning.hpp serif/utils/misc/terminal/colors.hpp serif/utils/types/strings.hpp serif/utils/types/type_list.hpp algorithm array cmath compare concepts cstddef cstdint cstdio cstdlib exception format iostream limits magic_enum/magic_enum_all.hpp memory mfem/fem/fe/fe_base.hpp mfem/fem/fe/fe_h1.hpp mfem/fem/fe/fe_l2.hpp mfem/fem/intrules.hpp nameof.hpp numeric span string string_view tuple type_traits unistd.h unordered_map utility vector version src/compile_time_checks/serif/model/finite_elements.cpp src/compile_time_checks/serif/model/scalar_row.cpp src/compile_time_checks/serif/model/state.cpp src/compile_time_checks/serif/operators/computation_owners.cpp src/compile_time_checks/serif/operators/plan.cpp src/include/serif/backend/mfem/discretization/basis_evaluation.hpp src/include/serif/backend/mfem/discretization/field_values.hpp src/include/serif/backend/mfem/model/binding.hpp src/include/serif/backend/mfem/model/computations/interior.hpp src/include/serif/backend/mfem/model/computations/potential_infinity.hpp src/include/serif/backend/mfem/model/computations/pressure_boundary.hpp src/include/serif/backend/mfem/model/computations/scalar_row.hpp src/include/serif/backend/mfem/model/computations/surface_potential.hpp src/include/serif/backend/mfem/model/discretization.hpp src/include/serif/backend/mfem/model/equilibrium.hpp src/include/serif/backend/mfem/model/equilibrium_evaluator.hpp src/include/serif/backend/mfem/model/finite_elements.hpp src/include/serif/backend/mfem/model/graph.hpp src/include/serif/backend/mfem/model/inputs.hpp src/include/serif/backend/mfem/model/volume.hpp src/include/serif/backend/mfem/model/volume_actions.hpp src/include/serif/backend/mfem/model/volume_displacement.hpp src/include/serif/backend/mfem/operators/element_data.hpp src/include/serif/backend/mfem/operators/element_evaluation.hpp src/include/serif/backend/mfem/operators/element_preparation.hpp src/include/serif/backend/mfem/operators/element_transpose.hpp src/include/serif/backend/mfem/operators/geometry_derivatives.hpp src/include/serif/backend/mfem/operators/jacobian_operator.hpp src/include/serif/backend/mfem/operators/operand_evaluation.hpp src/include/serif/backend/mfem/operators/volume_application.hpp src/include/serif/backend/mfem/operators/volume_computations.hpp src/include/serif/backend/mfem/operators/volume_data.hpp src/include/serif/backend/mfem/operators/volume_field_block.hpp src/include/serif/backend/mfem/operators/volume_inputs.hpp src/include/serif/backend/mfem/operators/volume_stages.hpp src/include/serif/backend/mfem/operators/workspace.hpp src/include/serif/backend/mfem/serialize/configuration.hpp src/include/serif/backend/mfem/serialize/discretization.hpp src/include/serif/backend/mfem/serialize/options.hpp src/include/serif/backend/mfem/solver/admissibility/geometry.hpp src/include/serif/backend/mfem/solver/equilibrium.hpp src/include/serif/backend/mfem/solver/nonlinear/stellar_problem.hpp src/include/serif/backend/mfem/solver/seed/state.hpp src/include/serif/model/finite_elements.hpp src/include/serif/solver/continuation/equilibrium.hpp src/include/serif/solver/equilibrium.hpp src/include/serif/solver/seed/state.hpp src/include/serif/solver/solver.hpp src/lib/serif/backend/mfem/discretization/basis_evaluation.cpp src/lib/serif/backend/mfem/model/finite_elements.cpp src/lib/serif/backend/mfem/model/system_selection.cpp src/lib/serif/backend/mfem/serialize/top.cpp serif::discretization::backend::mfem::BasisWorkspace serif::discretization::backend::mfem::BasisEvaluation serif serif::discretization serif::discretization::backend serif::discretization::backend::mfem Basis evaluation, pointwise mapping, and transpose application using dense or tensor-product representations. Emily M. Boudreaux #pragmaonce #include<memory> #include"mfem/fem/fe/fe_base.hpp" #include"mfem/fem/intrules.hpp" #include"serif/backend/mfem/discretization/tensor_basis.hpp" namespaceserif::discretization::backend::mfem{ structBasisWorkspace{ TensorWorkspacetensor; ::mfem::Vectorreference_values; ::mfem::Vectormapping_input; ::mfem::Vectormapping_output; }; classBasisEvaluationfinal{ public: voidsample( const::mfem::FiniteElement&element, const::mfem::IntegrationRule&rule, ReferenceOperationoperation, intcopies, TensorBasisCache&cache, ContractionModemode, constElementReconstruction*reconstruction=nullptr ); voidset_constant( intpoints ); voidreset_from( constBasisEvaluation&reference, intcomponents ); template<typenameMapT> voidmap_point( intpoint, constBasisEvaluation&reference, MapT&&map, BasisWorkspace&workspace ); voidevaluate( const::mfem::Vector&coefficients, ::mfem::Vector&values, BasisWorkspace&workspace )const; voidevaluate_many( const::mfem::Vector&coefficients, ::mfem::Vector&values, intlanes, BasisWorkspace&workspace )const; voidintegrate( const::mfem::Vector&values, ::mfem::Vector&coefficients, BasisWorkspace&workspace )const; [[nodiscard]]intheight()constnoexcept; [[nodiscard]]intwidth()constnoexcept; [[nodiscard]]booluses_tensor()constnoexcept; private: intm_points{0}; intm_components{0}; intm_width{0}; boolm_identity_mapping{false}; std::shared_ptr<constTensorBasis>m_tensor; ::mfem::DenseMatrixm_dense; ::mfem::DenseMatrixm_mapping; }; }//namespaceserif::discretization::backend::mfem template<typenameMapT> voidserif::discretization::backend::mfem::BasisEvaluation::map_point(constintpoint,constBasisEvaluation&reference,MapT&&map,BasisWorkspace&workspace){ constboolcanonical_reference=reference.m_identity_mapping; if(m_identity_mapping){ //Directlymappingacanonicaldestinationmustpreservetheuntouchedpoints. m_mapping.SetSize(m_points*m_components,m_components); m_mapping=0.0; for(intq=0;q<m_points;++q){ for(intcomponent=0;component<m_components;++component){ m_mapping(q*m_components+component,component)=1.0; } } } m_identity_mapping=false; const::mfem::DenseMatrix&source=reference.m_tensor?reference.m_mapping:reference.m_dense; ::mfem::DenseMatrix&destination=m_tensor?m_mapping:m_dense; constintcolumns=canonical_reference?reference.m_components:source.Width(); ::mfem::Vector&input=workspace.mapping_input; ::mfem::Vector&output=workspace.mapping_output; input.SetSize(reference.m_components); output.SetSize(m_components); for(intcolumn=0;column<columns;++column){ for(intaxis=0;axis<reference.m_components;++axis){ input[axis]=canonical_reference?static_cast<double>(axis==column):source(point*reference.m_components+axis,column); } map(input,output); for(intaxis=0;axis<m_components;++axis){ destination(point*m_components+axis,column)=output[axis]; } } }