field_values.hpp concepts serif/backend/mfem/discretization/term_coefficients.hpp serif/backend/mfem/operators/operand_evaluation.hpp src/lib/serif/backend/mfem/serialize/top.cpp serif/backend/mfem/discretization/basis_evaluation.hpp serif/backend/mfem/discretization/dof_map.hpp src/include/serif/backend/mfem/discretization/field_values.hpp serif/backend/mfem/discretization/reconstruction.hpp serif/backend/mfem/discretization/tensor_basis.hpp serif/backend/mfem/discretization/term_binding.hpp serif/backend/mfem/discretization/term_coefficients.hpp serif/backend/mfem/mapping/context.hpp serif/backend/mfem/mapping/transformations.hpp serif/backend/mfem/mapping/variation.hpp serif/backend/mfem/operators/operand_evaluation.hpp serif/dimensions/dimensions.hpp serif/dimensions/quantities.hpp serif/discretization/blocks/base.hpp serif/discretization/blocks/concepts.hpp serif/discretization/blocks/fields.hpp serif/discretization/blocks/traits.hpp serif/discretization/domain/concepts.hpp serif/discretization/domain/ids/boundary.hpp serif/discretization/domain/ids/concepts.hpp serif/discretization/domain/ids/domain.hpp serif/discretization/domain/ids/lists/concepts.hpp serif/discretization/domain/ids/lists/lists.hpp serif/discretization/domain/ids/lists/traits.hpp serif/discretization/domain/ids/validation/concepts.hpp serif/discretization/domain/ids/validation/unique_ids.hpp serif/discretization/domain/physical_domains.hpp serif/discretization/domain/relation/concepts.hpp serif/discretization/domain/relation/lists/concepts.hpp serif/discretization/domain/relation/lists/relation_list.hpp serif/discretization/domain/relation/lists/traits.hpp serif/discretization/domain/relation/relations.hpp serif/discretization/domain/relation/validation/concepts.hpp serif/discretization/domain/relation/validation/validation.hpp serif/discretization/domain/resolver.hpp serif/discretization/domain/schema/concepts.hpp serif/discretization/domain/schema/domain_schema.hpp serif/discretization/domain/schema/traits.hpp serif/discretization/domain/types.hpp serif/discretization/forms/operands.hpp serif/discretization/forms/value.hpp serif/discretization/functional_space/space.hpp serif/discretization/reconstruction/surface_power.hpp serif/equations/concepts.hpp serif/equations/contributions.hpp serif/equations/point.hpp serif/meta/enums.hpp serif/meta/names.hpp serif/operators/concepts.hpp serif/utils/dependencies/binding.hpp serif/utils/dependencies/dependency_context.hpp serif/utils/dependencies/detail/type_lists.hpp serif/utils/dependencies/key_list.hpp serif/utils/dependencies/resource.hpp serif/utils/dependencies/snapshot.hpp serif/utils/error/codes.hpp serif/utils/error/errors.hpp serif/utils/error/nonfinite.hpp serif/utils/misc/compensated_sum.hpp serif/utils/misc/concepts/concat.hpp serif/utils/misc/concepts/exceptions.hpp serif/utils/misc/concepts/numeric.hpp serif/utils/misc/concepts/type_uniqueness.hpp serif/utils/misc/std_helper/cleaning.hpp serif/utils/misc/std_helper/variadic_cleaning.hpp serif/utils/misc/terminal/colors.hpp serif/utils/misc/traits/type_uniqueness.hpp serif/utils/types/strings.hpp serif/utils/types/type_list.hpp algorithm array cmath compare concepts cstddef cstdint cstdio cstdlib exception format functional iostream limits magic_enum/magic_enum_all.hpp memory mfem/fem/eltrans.hpp mfem/fem/fe/fe_base.hpp mfem/fem/fe/fe_h1.hpp mfem/fem/fe/fe_l2.hpp mfem/fem/fespace.hpp mfem/fem/intrules.hpp mfem/fem/pfespace.hpp mfem/linalg/densemat.hpp mfem/linalg/vector.hpp nameof.hpp numeric optional span stdexcept string string_view tuple type_traits unistd.h unordered_map utility vector version src/include/serif/backend/mfem/discretization/field_values.hpp src/lib/serif/backend/mfem/serialize/top.cpp serif::discretization::backend::mfem::FieldValues serif serif::discretization serif::discretization::backend serif::discretization::backend::mfem #pragmaonce #include<concepts> #include"serif/backend/mfem/discretization/term_coefficients.hpp" #include"serif/backend/mfem/operators/operand_evaluation.hpp" namespaceserif::discretization::backend::mfem{ template<blocks::SpatialTermTerm> classFieldValuesfinal{ public: usingValueOperand=forms::Operand<Term>; usingValue=forms::OperandValue<ValueOperand>; FieldValues(constDofMap&binding,const::mfem::Vector&coefficients): m_binding(binding),m_coefficients(coefficients),m_elements(binding.validated_elements()),m_rule(1){ usingnamespaceutils::errors; usingTraits=blocks::SpatialTermTraits<Term>; binding.validate_values(coefficients); constexprintcopies=Traits::rank==0?1:3; constexprintcontinuity=Traits::family==functional_space::Family::h1?::mfem::FiniteElementCollection::CONTINUOUS:::mfem::FiniteElementCollection::DISCONTINUOUS; serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(Traits::rank<=1&& (Traits::family==functional_space::Family::h1||Traits::family==functional_space::Family::l2), "FieldValuessupportsscalar/vectorH1andL2fields."); serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(binding.space().GetVDim()==copies&& binding.space().FEColl()->GetContType()==continuity, "FieldValuestermfamilyorrankdisagreeswiththesuppliedbinding;passthematchingtypedcoefficientview."); for(intelement=0;element<binding.space().GetNE();++element){ constauto&finite_element=*binding.space().GetFE(element); serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(finite_element.GetOrder()==Traits::order&& finite_element.GetRangeType()==::mfem::FiniteElement::SCALAR&& finite_element.GetMapType()==::mfem::FiniteElement::VALUE, "FieldValuestermorderormappingdisagreeswiththesuppliedbinding;passthematchingtypedcoefficientview."); } } template<typenameBinding>requiresstd::same_as<typenameBinding::Term,Term> explicitFieldValues(constoperators::backend::mfem::CoefficientView<Binding>&field): FieldValues(field.binding,field.values){} FieldValues(constDofMap&,::mfem::Vector&&)=delete; FieldValues(constDofMap&,const::mfem::Vector&&)=delete; [[nodiscard]]Valuevalue(intelement,const::mfem::IntegrationPoint&point){ returnsample<forms::Operation::Value>(element,point); } voidvalues(intelement,const::mfem::IntegrationRule&points,::mfem::Vector&result){ usingEvaluator=operators::backend::mfem::OperandEvaluator<ValueOperand>; m_elements.element_values(element,m_coefficients,m_local); auto&transformation=*m_binding.mesh().GetElementTransformation(element); Evaluator::sample_basis(m_binding,element,points,m_basis,m_cache, ContractionMode::automatic,transformation,m_workspace); Evaluator::evaluate(m_basis,m_local,result,m_workspace); } [[nodiscard]]Valuevalue(intelement,const::mfem::IntegrationPoint&point,constmapping::MappingPointContext&context){ usingEvaluator=operators::backend::mfem::OperandEvaluator<ValueOperand>; returnEvaluator::physical_value(sample<forms::Operation::Value>(element,point),context); } [[nodiscard]]std::array<double,3>gradient(intelement,const::mfem::IntegrationPoint&point) requires(blocks::SpatialTermTraits<Term>::rank==0&&blocks::SpatialTermTraits<Term>::family==functional_space::Family::h1){ returnsample<forms::Operation::Gradient>(element,point); } [[nodiscard]]std::array<double,3>gradient(intelement,const::mfem::IntegrationPoint&point,constmapping::MappingPointContext&context) requires(blocks::SpatialTermTraits<Term>::rank==0&&blocks::SpatialTermTraits<Term>::family==functional_space::Family::h1){ usingOperand=forms::Operand<Term, forms::Operation::Gradient>; usingEvaluator=operators::backend::mfem::OperandEvaluator<Operand>; returnEvaluator::physical_value(sample<forms::Operation::Gradient>(element,point),context); } private: template<typenameOperation> [[nodiscard]]autosample(intelement,const::mfem::IntegrationPoint&point){ usingOperand=forms::Operand<Term, Operation>; usingEvaluator=operators::backend::mfem::OperandEvaluator<Operand>; m_elements.element_values(element,m_coefficients,m_local); m_rule.IntPoint(0)=point; auto&transformation=*m_binding.mesh().GetElementTransformation(element); //Arbitraryone-pointqueriesmustnotretainaseparatetensorplanforeveryquery. Evaluator::sample_basis(m_binding,element,m_rule,m_basis,m_cache, ContractionMode::dense,transformation,m_workspace); Evaluator::evaluate(m_basis,m_local,m_values,m_workspace); returnEvaluator::reference_value(m_values,0); } constDofMap&m_binding; const::mfem::Vector&m_coefficients; DofMap::Elementsm_elements; ::mfem::IntegrationRulem_rule; ::mfem::Vectorm_local,m_values; BasisEvaluationm_basis; TensorBasisCachem_cache; BasisWorkspacem_workspace; }; }