pressure_boundary.hpp cmath limits serif/discretization/blocks/system_layout.hpp serif/backend/mfem/geometry/surface/boundary_inputs.hpp serif/backend/mfem/model/computations/interior.hpp serif/backend/mfem/geometry/surface/pressure_field_block.hpp src/include/serif/backend/mfem/model/graph.hpp serif/backend/mfem/discretization/basis_evaluation.hpp serif/backend/mfem/discretization/boundary_trace.hpp serif/backend/mfem/discretization/coefficient_tracking.hpp serif/backend/mfem/discretization/dof_map.hpp serif/backend/mfem/discretization/dof_selection.hpp serif/backend/mfem/discretization/element_coordinates.hpp serif/backend/mfem/discretization/quadrature/concepts.hpp serif/backend/mfem/discretization/quadrature/configure.hpp serif/backend/mfem/discretization/quadrature/mfem_resolver.hpp serif/backend/mfem/discretization/quadrature/mfem_static_order.hpp serif/backend/mfem/discretization/quadrature/order.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/geometry/surface/boundary_data.hpp serif/backend/mfem/geometry/surface/boundary_inputs.hpp serif/backend/mfem/geometry/surface/boundary_rows.hpp serif/backend/mfem/geometry/surface/pressure_field_block.hpp serif/backend/mfem/mapping/compactification/compactification.hpp serif/backend/mfem/mapping/compactification/input_result.hpp serif/backend/mfem/mapping/compactification/kelvin.hpp serif/backend/mfem/mapping/compactification/variation.hpp serif/backend/mfem/mapping/context.hpp serif/backend/mfem/mapping/data.hpp serif/backend/mfem/mapping/domain_mapper.hpp serif/backend/mfem/mapping/status.hpp serif/backend/mfem/mapping/transformations.hpp serif/backend/mfem/mapping/variation.hpp serif/backend/mfem/mapping/volume_cache.hpp serif/backend/mfem/model/computations/interior.hpp src/include/serif/backend/mfem/model/computations/pressure_boundary.hpp serif/backend/mfem/model/computations/scalar_row.hpp serif/backend/mfem/model/computations/shape_derivatives.hpp serif/backend/mfem/model/deformation.hpp serif/backend/mfem/model/scalar_row.hpp serif/backend/mfem/model/system.hpp serif/backend/mfem/model/volume.hpp serif/backend/mfem/model/volume_actions.hpp serif/backend/mfem/operators/actions.hpp serif/backend/mfem/operators/element_data.hpp serif/backend/mfem/operators/element_evaluation.hpp serif/backend/mfem/operators/element_preparation.hpp serif/backend/mfem/operators/element_transpose.hpp serif/backend/mfem/operators/geometry_derivatives.hpp serif/backend/mfem/operators/linearized_block.hpp serif/backend/mfem/operators/operand_evaluation.hpp serif/backend/mfem/operators/preparation.hpp serif/backend/mfem/operators/quadrature.hpp serif/backend/mfem/operators/rejection.hpp serif/backend/mfem/operators/replaced_rows.hpp serif/backend/mfem/operators/volume_application.hpp serif/backend/mfem/operators/volume_computations.hpp serif/backend/mfem/operators/volume_data.hpp serif/backend/mfem/operators/volume_field_block.hpp serif/backend/mfem/operators/volume_inputs.hpp serif/backend/mfem/operators/volume_stages.hpp serif/backend/mfem/operators/workspace.hpp serif/backend/mfem/shared/finite.hpp serif/backend/mfem/shared/parallel.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/runtime.hpp serif/discretization/blocks/system_layout.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/concept.hpp serif/discretization/forms/form_traits.hpp serif/discretization/forms/forms.hpp serif/discretization/forms/operands.hpp serif/discretization/forms/runtime.hpp serif/discretization/forms/value.hpp serif/discretization/functional_space/space.hpp serif/discretization/quadrature/policy.hpp serif/discretization/quadrature/query.hpp serif/discretization/quadrature/rules.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/concepts.hpp serif/equations/contributions.hpp serif/equations/point.hpp serif/geometry/deformation/actions.hpp serif/geometry/deformation/composition.hpp serif/geometry/deformation/computations.hpp serif/geometry/deformation/detail/action_storage.hpp serif/geometry/deformation/failure.hpp serif/geometry/deformation/linearization.hpp serif/geometry/deformation/spaces.hpp serif/geometry/deformation/workspace.hpp serif/geometry/surface/compilation.hpp serif/geometry/surface/detail/relations.hpp serif/geometry/surface/isobaric.hpp serif/geometry/surface/linearization.hpp serif/geometry/surface/samples.hpp serif/meta/enums.hpp serif/meta/names.hpp serif/operators/concepts.hpp serif/operators/contributions.hpp serif/operators/plan.hpp serif/utils/checks/finite.hpp serif/utils/dependencies/binding.hpp serif/utils/dependencies/computation.hpp serif/utils/dependencies/computation_cache.hpp serif/utils/dependencies/computation_graph.hpp serif/utils/dependencies/dependencies.hpp serif/utils/dependencies/dependency_context.hpp serif/utils/dependencies/detail/cache_storage.hpp serif/utils/dependencies/detail/cache_traits.hpp serif/utils/dependencies/detail/graph.hpp serif/utils/dependencies/detail/result.hpp serif/utils/dependencies/detail/schema.hpp serif/utils/dependencies/detail/type_lists.hpp serif/utils/dependencies/key_list.hpp serif/utils/dependencies/key_list_operations.hpp serif/utils/dependencies/keys.hpp serif/utils/dependencies/observers.hpp serif/utils/dependencies/preparation.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/format/strings.hpp serif/utils/misc/compensated_sum.hpp serif/utils/misc/concepts/concat.hpp serif/utils/misc/concepts/exceptions.hpp serif/utils/misc/concepts/filter.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 serif/utils/user_options/mapping/compactification/kelvin.hpp serif/utils/user_options/mapping/domain_mapper.hpp algorithm array atomic cmath compare concepts cstddef cstdint cstdio cstdlib exception expected format functional initializer_list iostream limits magic_enum/magic_enum_all.hpp memory mfem/fem/doftrans.hpp 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/geom.hpp mfem/fem/gridfunc.hpp mfem/fem/intrules.hpp mfem/fem/pfespace.hpp mfem/linalg/densemat.hpp mfem/linalg/operator.hpp mfem/linalg/sparsemat.hpp mfem/linalg/vector.hpp mfem/mesh/mesh.hpp nameof.hpp numeric optional span stdexcept string string_view tuple type_traits unistd.h unordered_map utility variant vector version src/compile_time_checks/serif/model/state.cpp 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::model::computations::PressureBoundaryContribution serif serif::model serif::model::computations #pragmaonce #include<cmath> #include<limits> #include"serif/discretization/blocks/system_layout.hpp" #include"serif/backend/mfem/geometry/surface/boundary_inputs.hpp" #include"serif/backend/mfem/model/computations/interior.hpp" #include"serif/backend/mfem/geometry/surface/pressure_field_block.hpp" namespaceserif::model::computations { structPressureBoundaryContribution { }; template<geometry::surface::backend::mfem::PressureBoundaryFormulationFormulation,discretization::blocks::SystemLayoutTypeLayout> requiresutils::dependencies::ContainsKey<typename Formulation::CarrierTerm, typename Layout::Residuals> [[nodiscard]]autopressure_boundary_contribution(Layoutlayout) { usingnamespaceutils::dependencies; usingnamespaceutils::errors; usingnamespacediscretization::blocks; namespacesurface=geometry::surface; namespaceexecution=operators::backend::mfem; usingCarrier=typenameFormulation::CarrierTerm; serif_ensure<SERiFErrorCode::INVALID_DOF_COUNT>(layout.templatetotal_size<RepresentedValueKind::value>()<=static_cast<std::size_t>(std::numeric_limits<int>::max())&&layout.templatetotal_size<RepresentedValueKind::residual>()<=static_cast<std::size_t>(std::numeric_limits<int>::max()),"Pressure-boundarysystemcoordinatesmustfitMFEM'sintegerindexrange."); returnmake_computation<PressureBoundaryContribution,InputList<surface::keys::Linearization,surface::backend::mfem::keys::BoundaryDiscretization>>( [layout](conststd::vector<surface::SurfaceLinearization<Formulation>>&linearizations,constauto&boundary) { serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(boundary!=nullptr,"Pressure-boundaryassemblyrequiresitscompiledtraceowner."); boundary->validate(); surface::backend::mfem::detail::validate_boundary_samples(*boundary,linearizations); serif_ensure<SERiFErrorCode::MATRIX_SHAPE_ERROR>(layout.templatesize<RepresentedValueKind::residual,Carrier>()==static_cast<std::size_t>(boundary->selected_rows->binding().size()),"Thepressurecarrier'sresiduallayoutmustmatchitsfinite-elementrowbinding."); constintrow_offset=static_cast<int>(layout.templateoffset<RepresentedValueKind::residual,Carrier>()); backend::mfem::RowReplacementreplacement; replacement.selected_rows.reserve(linearizations.size()); replacement.residual.reserve(linearizations.size()); for(intindex=0;index<boundary->selected_rows->size();++index) { serif_ensure<SERiFErrorCode::NON_FINITE_VALUE_ERROR>(std::isfinite(linearizations[index].residual),"Apressurereplacementmustcontainfiniteacceptedresiduals."); replacement.selected_rows.push_back(row_offset+boundary->selected_rows->reduced_dofs()[index]); replacement.residual.push_back(linearizations[index].residual); } [&]<typename...Fields>(KeyList<Fields...>) { constautoadd_column=[&]<typenameField>() { ifconstexpr(ContainsKey<Field,typenameLayout::Values>) { automatrix=surface::backend::mfem::pressure_field_trace<Field>(*boundary,linearizations); serif_ensure<SERiFErrorCode::MATRIX_SHAPE_ERROR>(layout.templatesize<RepresentedValueKind::value,Field>()==static_cast<std::size_t>(matrix->Width()),"Asolvedpressureinput'scoordinatecountmustmatchitsfinite-elementtrace."); replacement.jacobian.push_back({.column_offset=static_cast<int>(layout.templateoffset<RepresentedValueKind::value,Field>()),.block=execution::linearized_block(std::move(matrix))}); } }; (add_column.templateoperator()<Fields>(),...); }(typenameFormulation::StateTerms{}); Contributionsresult; result.replacements.push_back(std::move(replacement)); returnresult; }); } }//namespaceserif::model::computations