boundary_computations.hpp serif/backend/mfem/geometry/surface/boundary_inputs.hpp serif/geometry/surface/computations.hpp serif/backend/mfem/shared/finite.hpp src/compile_time_checks/serif/geometry/surface_boundary.cpp src/include/serif/backend/mfem/geometry/surface/pressure_boundary.hpp src/include/serif/backend/mfem/model/graph.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/reconstruction.hpp serif/backend/mfem/discretization/term_binding.hpp serif/backend/mfem/discretization/term_coefficients.hpp src/include/serif/backend/mfem/geometry/surface/boundary_computations.hpp serif/backend/mfem/geometry/surface/boundary_data.hpp serif/backend/mfem/geometry/surface/boundary_inputs.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/evaluation.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/computations.hpp serif/geometry/surface/detail/batches.hpp serif/geometry/surface/detail/evaluation.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/utils/checks/finite.hpp serif/utils/dependencies/binding.hpp serif/utils/dependencies/computation.hpp serif/utils/dependencies/computation_graph.hpp serif/utils/dependencies/dependency_context.hpp serif/utils/dependencies/detail/graph.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 expected 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/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/geometry/surface_boundary.cpp src/compile_time_checks/serif/model/state.cpp src/include/serif/backend/mfem/geometry/surface/boundary_computations.hpp src/include/serif/backend/mfem/geometry/surface/pressure_boundary.hpp src/include/serif/backend/mfem/model/binding.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/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 serif serif::geometry serif::geometry::surface serif::geometry::surface::backend serif::geometry::surface::backend::mfem serif::geometry::surface::backend::mfem::detail #pragmaonce #include"serif/backend/mfem/geometry/surface/boundary_inputs.hpp" #include"serif/geometry/surface/computations.hpp" #include"serif/backend/mfem/shared/finite.hpp" namespaceserif::geometry::surface::backend::mfem { namespacedetail { template<typenameTerm,typenameBoundary,typenameInput> [[nodiscard]]std::expected<BoundaryValues<Term>,SurfaceFailure>sample_boundary_values(constBoundary&boundary,constInput&input) { usingnamespaceutils::errors; constauto&trace=boundary.traces.templateget<Term>(); serif_ensure<SERiFErrorCode::CACHE_FE_BINDING_ERROR, std::invalid_argument>(&trace.binding()==&input.binding,"Pressure-boundarycoefficientsmustusethebindingthatcompiledtheirtrace."); serif_ensure<SERiFErrorCode::INVALID_DOF_COUNT, std::invalid_argument>(input.values.Size()==trace.binding().size(),"Pressure-boundarycoefficientsmustmatchthesupportedfieldsize."); if(!utils::checks::is_finite(input.values)) { returnstd::unexpected(SurfaceFailure{.reason=SurfaceFailureReason::non_finite_state,.message="Apressure-boundaryfieldcontainsanon-finitecoefficient."}); } BoundaryValues<Term>result{.values=std::vector<double>(boundary.selected_rows->size())}; ::mfem::Vectorvalues(result.values.data(),static_cast<int>(result.values.size())); trace.interpolation().Mult(input.values,values); if(!utils::checks::is_finite(values)) { returnstd::unexpected(SurfaceFailure{.reason=SurfaceFailureReason::non_finite_state,.message="Afinite-elementboundarytraceproducedanon-finitevalue."}); } returnresult; } template<PressureBoundaryFormulationFormulation,typename...Terms> [[nodiscard]]autopressure_boundary_computations(utils::dependencies::KeyList<Terms...>) { usingnamespaceutils::dependencies; constautosample=[]<typenameTerm>() { returnmake_computation<keys::BoundarySamples<Term>,InputList<keys::BoundaryDiscretization,Term>>( [](constauto&boundary,constauto&input) { static_assert(std::same_as<typenamestd::remove_cvref_t<decltype(*boundary)>::Formulation,Formulation>,"Theboundarygraphmustusetheformulationthatcompileditstraces."); returnsample_boundary_values<Term>(*boundary,input); }); }; returnstd::tuple_cat(std::tuple{sample.templateoperator()<Terms>()...},make_pressure_surface_computations<Formulation, keys::BoundarySamples>()); } } template<PressureBoundaryFormulationFormulation> [[nodiscard]]automake_pressure_boundary_computations() { returndetail::pressure_boundary_computations<Formulation>(typenameFormulation::StateTerms{}); } template<PressureBoundaryFormulationFormulation> [[nodiscard]]automake_pressure_boundary_graph() { returnstd::apply([](auto...computations) { returnutils::dependencies::make_computation_graph(std::move(computations)...); },make_pressure_boundary_computations<Formulation>()); } }