dof_selection.hpp serif/backend/mfem/discretization/term_binding.hpp serif/backend/mfem/discretization/dof_map.hpp src/include/serif/backend/mfem/discretization/boundary_trace.hpp src/include/serif/backend/mfem/model/computations/potential_infinity.hpp src/include/serif/backend/mfem/model/normalization/computations.hpp src/include/serif/backend/mfem/solver/seed/state.hpp serif/backend/mfem/discretization/dof_map.hpp src/include/serif/backend/mfem/discretization/dof_selection.hpp serif/backend/mfem/discretization/reconstruction.hpp serif/backend/mfem/discretization/term_binding.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/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/error/nonfinite.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 memory mfem/fem/fe/fe_base.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 vector version src/compile_time_checks/serif/geometry/surface_boundary.cpp src/compile_time_checks/serif/model/state.cpp src/include/serif/backend/mfem/discretization/boundary_trace.hpp src/include/serif/backend/mfem/discretization/dof_selection.hpp src/include/serif/backend/mfem/geometry/surface/boundary_computations.hpp src/include/serif/backend/mfem/geometry/surface/boundary_data.hpp src/include/serif/backend/mfem/geometry/surface/boundary_inputs.hpp src/include/serif/backend/mfem/geometry/surface/boundary_rows.hpp src/include/serif/backend/mfem/geometry/surface/pressure_boundary.hpp src/include/serif/backend/mfem/geometry/surface/pressure_field_block.hpp src/include/serif/backend/mfem/model/binding.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/equilibrium.hpp src/include/serif/backend/mfem/model/equilibrium_evaluator.hpp src/include/serif/backend/mfem/model/graph.hpp src/include/serif/backend/mfem/model/normalization/computations.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/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/boundary_trace.cpp serif::discretization::backend::mfem::DofSelection serif serif::discretization serif::discretization::backend serif::discretization::backend::mfem This file contains the DofSelection which allows for the selection of a subset of DOFs from a term binding, with boundary attribute information. Emily M. Boudreaux #pragmaonce #include"serif/backend/mfem/discretization/term_binding.hpp" #include"serif/backend/mfem/discretization/dof_map.hpp" namespaceserif::discretization::backend::mfem{ classDofSelectionfinal{ public: DofSelection( constDofMap&binding, const::mfem::Array<int>&reduced_dofs, constbooltrack_coordinates=false ): m_binding(binding), m_dofs(reduced_dofs), m_boundary_attributes(binding.mesh().GetNBE()), m_nodes_sequence(track_coordinates?binding.mesh().GetNodesSequence():-1){ usingnamespaceutils::errors; binding.validate(); intprevious=-1; for(constintdof:m_dofs){ serif_ensure<SERiFErrorCode::INVALID_DOF_COUNT,std::invalid_argument>(dof>previous&&dof<binding.size(),"SelectedDOFsmustbesorted,unique,andinsidethetermbinding."); previous=dof; } for(intboundary=0;boundary<m_boundary_attributes.Size();++boundary){ m_boundary_attributes[boundary]=binding.mesh().GetBdrAttribute(boundary); } } [[nodiscard]]constDofMap&binding()constnoexcept{returnm_binding;} [[nodiscard]]const::mfem::Array<int>&reduced_dofs()constnoexcept{returnm_dofs;} [[nodiscard]]intsize()constnoexcept{returnm_dofs.Size();} [[nodiscard]]const::mfem::Array<int>&boundary_attributes()constnoexcept{returnm_boundary_attributes;} [[nodiscard]]booltracks_coordinates()constnoexcept{returnm_nodes_sequence>=0;} voidvalidate()const{ usingnamespaceutils::errors; m_binding.validate(); constauto&mesh=m_binding.mesh(); serif_ensure<SERiFErrorCode::CACHE_FE_BINDING_ERROR,std::logic_error>(m_nodes_sequence<0||mesh.GetNodesSequence()==m_nodes_sequence,"Rebuildpointselectionsafterfixed-meshcoordinatechanges."); serif_ensure<SERiFErrorCode::CACHE_FE_BINDING_ERROR,std::logic_error>(mesh.GetNBE()==m_boundary_attributes.Size(),"RebuildDOFselectionsafterboundarychanges."); for(intboundary=0;boundary<m_boundary_attributes.Size();++boundary){ serif_ensure<SERiFErrorCode::CACHE_FE_BINDING_ERROR,std::logic_error>(mesh.GetBdrAttribute(boundary)==m_boundary_attributes[boundary],"RebuildDOFselectionsafterboundaryattributeschange."); } } private: constDofMap&m_binding; ::mfem::Array<int>m_dofs; ::mfem::Array<int>m_boundary_attributes; longm_nodes_sequence; }; template<domain::IsBoundaryBoundaryT,blocks::SpatialTermTermT,domain::schema::IsSchemaSchemaT> [[nodiscard]]DofSelectionselect_boundary_dofs( constTermBinding<TermT, SchemaT>&binding ){ usingnamespaceutils::errors; static_assert(SchemaT::templatecontains_boundary<BoundaryT>(),"Theboundarymustberegisteredinthedomainschema."); binding.validate(); ::mfem::Array<int>marker(binding.size()); marker=0; ::mfem::Array<int>dofs; ::mfem::Array<int>selected; auto&space=binding.space(); constauto&mesh=binding.mesh(); for(intboundary=0;boundary<mesh.GetNBE();++boundary){ if(!SchemaT::templateboundary_id_matches<BoundaryT>(mesh.GetBdrAttribute(boundary))){ continue; } space.GetBdrElementVDofs(boundary,dofs); for(constintencoded:dofs){ boolsupported=false; binding.for_each_coefficient(::mfem::FiniteElementSpace::DecodeDof(encoded),[&](constintindex,constdouble){ marker[index]=1; supported=true; }); serif_ensure<SERiFErrorCode::INVALID_DOMAIN_ERROR,std::invalid_argument>(supported,"AselectedboundaryDOFliesoutsidetheterm'ssupport."); } } ::mfem::FiniteElementSpace::MarkerToList(marker,selected); serif_ensure<SERiFErrorCode::INVALID_DOMAIN_ERROR,std::invalid_argument>(selected.Size()>0,"TherequestedboundaryhasnoDOFsonthisserialmesh."); returnDofSelection(binding,selected); } }//namespaceserif::discretization::backend::mfem