point_dofs.hpp array cmath serif/discretization/backend/mfem/dof_selection.hpp src/include/serif/discretization/backend/mfem/mfem_backend.hpp src/include/serif/operators/backend/mfem/nodal_inputs.hpp serif/dimensions/dimensions.hpp serif/dimensions/quantities.hpp serif/discretization/backend/mfem/dof_map.hpp serif/discretization/backend/mfem/dof_selection.hpp src/include/serif/discretization/backend/mfem/point_dofs.hpp serif/discretization/backend/mfem/term_binding.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/functional_space/space.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/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 iostream magic_enum/magic_enum_all.hpp mfem/fem/fespace.hpp mfem/fem/pfespace.hpp nameof.hpp numeric optional span stdexcept string string_view tuple type_traits unistd.h unordered_map utility version src/compile_time_checks/serif/operators/computation_owners.cpp src/compile_time_checks/serif/operators/serial_backend.cpp src/include/serif/discretization/backend/mfem/mfem_backend.hpp src/include/serif/discretization/backend/mfem/point_dofs.hpp src/include/serif/model/backend/mfem/binding.hpp src/include/serif/model/backend/mfem/equilibrium.hpp src/include/serif/model/computations/point_constraints.hpp src/include/serif/model/graph.hpp src/include/serif/operators/backend/mfem/nodal_computations.hpp src/include/serif/operators/backend/mfem/nodal_inputs.hpp serif serif::discretization serif::discretization::backend serif::discretization::backend::mfem This file contains the implementation of function select_point_dofs. This function is used to select the components from a single H1 vertex in the fixed computational mesh. This is used to implement point constraints in the discretization (e.g. fixed center of mass) Emily M. Boudreaux #pragmaonce #include<array> #include<cmath> #include"serif/discretization/backend/mfem/dof_selection.hpp" namespaceserif::discretization::backend::mfem{ template<blocks::SpatialTermTermT,domain::schema::IsSchemaSchemaT> requires(blocks::SpatialTermTraits<TermT>::family==functional_space::Family::h1) [[nodiscard]]DofSelectionselect_point_dofs( constTermBinding<TermT, SchemaT>&binding, conststd::array<double,3>&position={}, constdoubletolerance=1.0e-12 ){ usingnamespaceutils::errors; //Mostlythestartofthisisjustvalidationdownto... binding.validate(); serif_ensure<SERiFErrorCode::INVALID_ARGUMENT,std::invalid_argument>(std::isfinite(tolerance)&&tolerance>=0.0,"Pointselectionrequiresafinitenonnegativecoordinatetolerance."); for(constdoublecoordinate:position){ serif_ensure<SERiFErrorCode::NON_FINITE_VALUE_ERROR,std::invalid_argument>(std::isfinite(coordinate),"Pointselectionrequiresfinitereferencecoordinates."); } auto&space=binding.space(); for(intelement=0;element<space.GetNE();++element){ serif_ensure<SERiFErrorCode::INVALID_ARGUMENT,std::invalid_argument>(dynamic_cast<const::mfem::NodalFiniteElement*>(space.GetFE(element))!=nullptr,"PointvaluesrequireanodalH1basis;positive/Bernsteincoefficientsarenotnodalvalues."); } //downtohere. //Thenwecomethroughandwetrytofindthevertexthatmatchestherequestedpositionwithinthenumericaltolerance.e.g.wemightpass //apositionof(0,0,0)hereandatoleranceof1e-12andthiswillthenmatchanyvertexwithinaeuclideandistanceof1e-12from(0,0,0) //Therearebetterwaystodothisbutthisworksfornow. intselected_vertex=-1; for(intvertex=0;vertex<binding.mesh().GetNV();++vertex){ std::array<double,3>coordinates; binding.mesh().GetNode(vertex,coordinates.data()); boolmatches=true; for(intaxis=0;axis<3;++axis){ matches=matches&&std::abs(coordinates[axis]-position[axis])<=tolerance; } if(matches){ serif_ensure<SERiFErrorCode::INVALID_ARGUMENT,std::invalid_argument>(selected_vertex<0,"Therequestedpointmatchesmultiplemeshvertices;useauniquematerialpointorasmallertolerance."); selected_vertex=vertex; } } //Notethathereweenforcethattheselectedpositioniswithintoleranceofameshvertex.Afutureimplementationmightconsidersomemore //complexschemethatinvolvesfindingtheboundingelementsandenforcingsomesortofinterpolationscheme. serif_ensure<SERiFErrorCode::INVALID_ARGUMENT,std::invalid_argument>(selected_vertex>=0,"Therequestedpointisnotameshvertex;pointconstraintscurrentlyrequireauniquenodalH1vertex."); //Thenweloopthroughandselectthedofsassociatedwiththeselectedvertexinthereduceddofspaceoftheprovidedtermbinding. //Ifthatpositionisoutsidethesupportofthetermbindingthiswillthrow. ::mfem::Array<int>vertex_dofs,selected; space.GetVertexDofs(selected_vertex,vertex_dofs); serif_ensure<SERiFErrorCode::INVALID_DOF_COUNT,std::invalid_argument>(vertex_dofs.Size()==1,"ApointconstraintrequiresonescalarnodalDOFatitsvertex."); for(intcomponent=0;component<space.GetVDim();++component){ constautoreduced=binding.reduced_dof(space.DofToVDof(vertex_dofs[0],component)); serif_ensure<SERiFErrorCode::INVALID_DOMAIN_ERROR,std::invalid_argument>(reduced.has_value(),"Theselectedpointliesoutsidetheterm'ssupport."); selected.Append(*reduced); } returnDofSelection(binding,selected,true); } }//namespaceserif::discretization::backend::mfem