reconstruction.hpp memory vector mfem/fem/fe/fe_base.hpp serif/discretization/blocks/fields.hpp serif/discretization/reconstruction/surface_power.hpp src/include/serif/backend/mfem/discretization/dof_map.hpp src/include/serif/backend/mfem/discretization/tensor_basis.hpp src/lib/serif/backend/mfem/discretization/reconstruction.cpp src/include/serif/backend/mfem/discretization/reconstruction.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 magic_enum/magic_enum_all.hpp memory mfem/fem/fe/fe_base.hpp nameof.hpp numeric span 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/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/boundary_trace.hpp src/include/serif/backend/mfem/discretization/dof_map.hpp src/include/serif/backend/mfem/discretization/dof_selection.hpp src/include/serif/backend/mfem/discretization/element_coordinates.hpp src/include/serif/backend/mfem/discretization/field_values.hpp src/include/serif/backend/mfem/discretization/mass_diagonal.hpp src/include/serif/backend/mfem/discretization/reconstruction.hpp src/include/serif/backend/mfem/discretization/tensor_basis.hpp src/include/serif/backend/mfem/discretization/term_binding.hpp src/include/serif/backend/mfem/discretization/term_coefficients.hpp src/include/serif/backend/mfem/geometry/deformation/displacement_mapping.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/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/normalization/computations.hpp src/include/serif/backend/mfem/model/state.hpp src/include/serif/backend/mfem/model/surface_potential.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/serialize/resources.hpp src/include/serif/backend/mfem/serialize/top.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/normalization/metric_computations.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/discretization/boundary_trace.cpp src/lib/serif/backend/mfem/discretization/element_coordinates.cpp src/lib/serif/backend/mfem/discretization/mass_diagonal.cpp src/lib/serif/backend/mfem/discretization/reconstruction.cpp src/lib/serif/backend/mfem/model/deformation.cpp src/lib/serif/backend/mfem/model/finite_elements.cpp src/lib/serif/backend/mfem/model/surface_potential.cpp src/lib/serif/backend/mfem/model/system_selection.cpp src/lib/serif/backend/mfem/serialize/resources.cpp src/lib/serif/backend/mfem/serialize/top.cpp serif::discretization::backend::mfem::ElementReconstruction serif::discretization::backend::mfem::ReconstructionPlan mfem serif serif::discretization serif::discretization::backend serif::discretization::backend::mfem #pragmaonce #include<memory> #include<vector> #include"mfem/fem/fe/fe_base.hpp" #include"serif/discretization/blocks/fields.hpp" #include"serif/discretization/reconstruction/surface_power.hpp" namespacemfem{classFiniteElementSpace;} namespaceserif::discretization::backend::mfem{ classDofMap; structElementReconstruction{ constreconstruction::RadialProfile*profile{nullptr}; intinward_axis{0}; boolsurface_at_high_coordinate{false}; booloperator==(constElementReconstruction&)const=default; voidradial(const::mfem::FiniteElement&element,doublecoordinate,::mfem::Vector&values,::mfem::Vector&derivatives)const; voidsample(const::mfem::FiniteElement&element,const::mfem::IntegrationPoint&point, ::mfem::Vector&values,::mfem::DenseMatrix&derivatives)const; }; classReconstructionPlanfinal{ public: ReconstructionPlan(constDofMap&binding,intboundary_attribute,std::shared_ptr<constreconstruction::RadialProfile>profile); voidvalidate_binding(constDofMap&binding)const; [[nodiscard]]constElementReconstruction*element(intindex)const; [[nodiscard]]constreconstruction::RadialProfile&profile()constnoexcept{return*m_profile;} [[nodiscard]]intboundary_attribute()constnoexcept{returnm_boundary_attribute;} [[nodiscard]]intintegration_order_increment()constnoexcept{return8;} private: voidvalidate_interfaces(constDofMap&binding)const; const::mfem::FiniteElementSpace*m_space; longm_space_sequence; longm_mesh_sequence; std::vector<int>m_domain_attributes,m_boundary_attributes,m_element_support; intm_boundary_attribute; std::shared_ptr<constreconstruction::RadialProfile>m_profile; std::vector<ElementReconstruction>m_elements; }; } namespaceserif::discretization::backend::mfem{ template<typenameTerm,typenameBinding> [[nodiscard]]std::shared_ptr<constReconstructionPlan>make_reconstruction( constreconstruction::SurfacePower<Term>&specification,constBinding&binding){ usingnamespaceutils::errors; ifconstexpr(!std::same_as<Term,blocks::density::mass>){ serif_error<std::invalid_argument,SERiFErrorCode::INVALID_ARGUMENT>("SurfacePowersupportsonlyscalarH1massdensity;therequestedfieldhasnovalidatedreconstruction."); return{}; }else{ intboundary_attribute=0; for(constauto&boundary:Binding::Schema::boundary_descriptors()){ if(boundary.name==specification.region.boundary_name){boundary_attribute=boundary.ID;} } serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(boundary_attribute!=0, "SurfacePowerboundary'{}'isnotregisteredinthisdiscretization'sdomainschema.",specification.region.boundary_name); autoprofile=std::make_shared<constreconstruction::SurfacePowerProfile>(specification.density_exponent); if(specification.density_exponent==1.0){return{};} returnstd::make_shared<constReconstructionPlan>(binding,boundary_attribute,std::move(profile)); } } }