pressure_field_block.hpp cmath memory utility vector serif/backend/mfem/geometry/surface/boundary_rows.hpp serif/backend/mfem/operators/replaced_rows.hpp serif/utils/dependencies/key_list_operations.hpp src/include/serif/backend/mfem/model/computations/pressure_boundary.hpp serif/backend/mfem/discretization/boundary_trace.hpp serif/backend/mfem/discretization/dof_map.hpp serif/backend/mfem/discretization/dof_selection.hpp serif/backend/mfem/discretization/reconstruction.hpp serif/backend/mfem/discretization/term_binding.hpp serif/backend/mfem/discretization/term_coefficients.hpp serif/backend/mfem/geometry/surface/boundary_data.hpp serif/backend/mfem/geometry/surface/boundary_rows.hpp src/include/serif/backend/mfem/geometry/surface/pressure_field_block.hpp serif/backend/mfem/operators/linearized_block.hpp serif/backend/mfem/operators/replaced_rows.hpp serif/backend/mfem/shared/finite.hpp serif/dimensions/dimensions.hpp serif/dimensions/partial.hpp serif/dimensions/quantities.hpp serif/dimensions/type_alias.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/eos/models/concepts.hpp serif/eos/models/traits.hpp serif/eos/relations/catalog.hpp serif/eos/relations/concepts.hpp serif/eos/relations/relations.hpp serif/eos/relations/traits.hpp serif/equations/point.hpp serif/geometry/surface/compilation.hpp serif/geometry/surface/detail/relations.hpp serif/geometry/surface/isobaric.hpp serif/geometry/surface/linearization.hpp serif/meta/enums.hpp serif/meta/names.hpp serif/utils/checks/finite.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/key_list_operations.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/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/fe/fe_base.hpp mfem/fem/fespace.hpp mfem/fem/pfespace.hpp mfem/linalg/densemat.hpp mfem/linalg/operator.hpp mfem/linalg/sparsemat.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/compile_time_checks/serif/model/state.cpp src/include/serif/backend/mfem/geometry/surface/pressure_field_block.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/graph.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/solver/continuation/equilibrium.hpp src/include/serif/solver/equilibrium.hpp src/include/serif/solver/solver.hpp serif serif::geometry serif::geometry::surface serif::geometry::surface::backend serif::geometry::surface::backend::mfem #pragmaonce #include<cmath> #include<memory> #include<utility> #include<vector> #include"serif/backend/mfem/geometry/surface/boundary_rows.hpp" #include"serif/backend/mfem/operators/replaced_rows.hpp" #include"serif/utils/dependencies/key_list_operations.hpp" namespaceserif::geometry::surface::backend::mfem{ template<typenameField,PressureBoundaryFormulationFormulation,typenameTraces> requiresutils::dependencies::ContainsKey<Field,typenameFormulation::StateTerms> [[nodiscard]]std::shared_ptr<const::mfem::SparseMatrix>pressure_field_trace(constPressureBoundary<Formulation, Traces>&boundary,conststd::vector<SurfaceLinearization<Formulation>>&linearizations){ usingnamespaceutils::errors; boundary.validate(); detail::validate_boundary_samples(boundary,linearizations); autoresult=std::make_shared<::mfem::SparseMatrix>(boundary.traces.templateget<Field>().interpolation()); ::mfem::Vectorderivatives(boundary.selected_rows->size()); for(introw=0;row<derivatives.Size();++row){ derivatives[row]=linearizations[row].field_derivatives.templateget<Field>(); serif_ensure<SERiFErrorCode::NON_FINITE_VALUE_ERROR>(std::isfinite(derivatives[row]),"Pressure-boundaryderivative{}isnonfinite.Supplyanacceptedsurfacelinearization.",row); } result->ScaleRows(derivatives); returnresult; } template<typenameField,PressureBoundaryFormulationFormulation,typenameTraces> [[nodiscard]]operators::backend::mfem::LinearizedBlockpressure_field_block(operators::backend::mfem::LinearizedBlockvolume,constPressureBoundary<Formulation, Traces>&boundary,conststd::vector<SurfaceLinearization<Formulation>>&linearizations){ usingnamespaceutils::errors; boundary.validate(); detail::validate_boundary_samples(boundary,linearizations); serif_ensure<SERiFErrorCode::MATRIX_SHAPE_ERROR>(volume.action()->Height()==boundary.selected_rows->binding().size(),"Theselectedphysicalblockandpressureboundarymustusethesamefullresidualcoordinatecount."); std::shared_ptr<const::mfem::SparseMatrix>replacement; ifconstexpr(utils::dependencies::ContainsKey<Field,typenameFormulation::StateTerms>){ replacement=pressure_field_trace<Field>(boundary,linearizations); serif_ensure<SERiFErrorCode::MATRIX_SHAPE_ERROR>(replacement->Width()==volume.action()->Width(),"Theselectedphysicalfieldanditspressuretracehavedifferentcoefficientcounts."); } else{ autozero=std::make_shared<::mfem::SparseMatrix>(boundary.selected_rows->size(),volume.action()->Width()); zero->Finalize(); replacement=std::move(zero); } const::mfem::Array<int>&selected=boundary.selected_rows->reduced_dofs(); returnoperators::backend::mfem::replace_block_rows(std::move(volume),std::vector<int>(selected.begin(),selected.end()),std::move(replacement)); } }//namespaceserif::geometry::surface::backend::mfem