point_constraints.hpp limits memory tuple type_traits utility vector serif/discretization/blocks/system_layout.hpp serif/model/computations/interior.hpp serif/operators/backend/mfem/linearized_block.hpp serif/operators/backend/mfem/nodal_computations.hpp serif/operators/backend/mfem/point_border.hpp src/include/serif/model/backend/mfem/binding.hpp src/include/serif/model/graph.hpp serif/dimensions/dimensions.hpp serif/dimensions/quantities.hpp serif/discretization/backend/mfem/basis_evaluation.hpp serif/discretization/backend/mfem/coefficient_tracking.hpp serif/discretization/backend/mfem/dof_map.hpp serif/discretization/backend/mfem/dof_selection.hpp serif/discretization/backend/mfem/point_dofs.hpp serif/discretization/backend/mfem/tensor_basis.hpp serif/discretization/backend/mfem/term_binding.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/backend/mfem/concepts.hpp serif/discretization/quadrature/backend/mfem/configure.hpp serif/discretization/quadrature/backend/mfem/mfem_resolver.hpp serif/discretization/quadrature/backend/mfem/mfem_static_order.hpp serif/discretization/quadrature/backend/order.hpp serif/discretization/quadrature/policy.hpp serif/discretization/quadrature/query.hpp serif/discretization/quadrature/rules.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/mapping/compactification/compactification.hpp serif/mapping/compactification/input_result.hpp serif/mapping/compactification/kelvin.hpp serif/mapping/compactification/variation.hpp serif/mapping/context.hpp serif/mapping/data.hpp serif/mapping/domain_mapper.hpp serif/mapping/status.hpp serif/mapping/transformations.hpp serif/mapping/variation.hpp serif/mapping/volume_cache.hpp serif/meta/enums.hpp serif/meta/names.hpp serif/model/backend/mfem/actions.hpp serif/model/backend/mfem/deformation.hpp serif/model/backend/mfem/system.hpp serif/model/backend/mfem/volume.hpp serif/model/backend/mfem/volume_actions.hpp serif/model/computations/interior.hpp src/include/serif/model/computations/point_constraints.hpp serif/model/computations/shape_derivatives.hpp serif/operators/backend/mfem/element_coordinates.hpp serif/operators/backend/mfem/element_data.hpp serif/operators/backend/mfem/element_evaluation.hpp serif/operators/backend/mfem/element_preparation.hpp serif/operators/backend/mfem/element_transpose.hpp serif/operators/backend/mfem/linearized_block.hpp serif/operators/backend/mfem/nodal_computations.hpp serif/operators/backend/mfem/nodal_data.hpp serif/operators/backend/mfem/nodal_evaluation.hpp serif/operators/backend/mfem/nodal_inputs.hpp serif/operators/backend/mfem/nodal_rows.hpp serif/operators/backend/mfem/operand_evaluation.hpp serif/operators/backend/mfem/point_border.hpp serif/operators/backend/mfem/preparation.hpp serif/operators/backend/mfem/quadrature.hpp serif/operators/backend/mfem/rejection.hpp serif/operators/backend/mfem/term_coefficients.hpp serif/operators/backend/mfem/volume_application.hpp serif/operators/backend/mfem/volume_computations.hpp serif/operators/backend/mfem/volume_data.hpp serif/operators/backend/mfem/volume_field_block.hpp serif/operators/backend/mfem/volume_inputs.hpp serif/operators/backend/mfem/volume_stages.hpp serif/operators/backend/mfem/workspace.hpp serif/operators/concepts.hpp serif/operators/contributions.hpp serif/operators/plan.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/computation_store.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/incremental_computation.hpp serif/utils/dependencies/key_list.hpp serif/utils/dependencies/key_list_operations.hpp serif/utils/dependencies/keys.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/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 cmath compare concepts cstddef cstdint cstdio cstdlib exception expected format functional initializer_list iostream limits magic_enum/magic_enum_all.hpp map 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/fe/fe_rt.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/include/serif/model/backend/mfem/binding.hpp src/include/serif/model/backend/mfem/equilibrium.hpp src/include/serif/model/computations/point_constraints.hpp src/include/serif/model/graph.hpp serif::model::computations::PointConstraintContribution serif serif::model serif::model::computations #pragmaonce #include<limits> #include<memory> #include<tuple> #include<type_traits> #include<utility> #include<vector> #include"serif/discretization/blocks/system_layout.hpp" #include"serif/model/computations/interior.hpp" #include"serif/operators/backend/mfem/linearized_block.hpp" #include"serif/operators/backend/mfem/nodal_computations.hpp" #include"serif/operators/backend/mfem/point_border.hpp" namespaceserif::model::computations { template<typenameScope> structPointConstraintContribution { }; template<operators::backend::mfem::AffinePointEquationEquation,typenameDiscretization> [[nodiscard]]autopoint_constraint_discretization(constEquation&,constDiscretization&discretization) { returnoperators::backend::mfem::make_point_border_discretization<Equation>(discretization.templatebinding<typenameEquation::Test>(),discretization.templatebinding<typenameEquation::Border>()); } template<typenameScope,typenameEquation,typenameSchema,discretization::blocks::SystemLayoutTypeLayout> [[nodiscard]]autopoint_constraint_contribution(Layoutlayout) { usingnamespaceutils::dependencies; usingnamespacediscretization::blocks; namespaceexecution=operators::backend::mfem; usingCarrier=typenameEquation::Test; usingBorder=typenameEquation::Border; usingInput=typenameexecution::nodal_keys::Scoped<Scope>::BorderLinearization; returnmake_computation<PointConstraintContribution<Scope>,InputList<Input>>([layout](constexecution::PointBorderData<Equation,Schema>&data) { usingnamespaceutils::errors; constintcarrier_size=data.nodal.discretization.residual_binding.size(); serif_ensure<SERiFErrorCode::INVALID_DOF_COUNT>(layout.templatesize<RepresentedValueKind::value,Carrier>()==static_cast<std::size_t>(carrier_size)&& layout.templatesize<RepresentedValueKind::residual,Carrier>()==static_cast<std::size_t>(carrier_size)&&layout.templatesize<RepresentedValueKind::value,Border>()==1&& layout.templatesize<RepresentedValueKind::residual,Border>()==1,"Themodellayoutmustmatchthepointconstraint'scarriersize{}andonescalarborderoneachside.",carrier_size); serif_ensure<SERiFErrorCode::INVALID_DOF_COUNT>(layout.templatetotal_size<RepresentedValueKind::value>()<=std::numeric_limits<int>::max()&& layout.templatetotal_size<RepresentedValueKind::residual>()<=std::numeric_limits<int>::max(),"Thepointconstraintlayoutexceedsthebackend'ssupportedintegercoefficientcount."); ::mfem::Vectorcarrier(carrier_size),phase(1); carrier=0.0; execution::add_point_border_residual(data,carrier,phase); constintpoint=data.nodal.discretization.selected_rows->reduced_dofs()[0]; constdoublederivative=data.nodal.derivative(0,0); autoforce=std::make_shared<::mfem::SparseMatrix>(carrier_size,1); autoconstraint=std::make_shared<::mfem::SparseMatrix>(1,carrier_size); force->Add(point,0,derivative); constraint->Add(0,point,derivative); force->Finalize(); constraint->Finalize(); constintcarrier_row=static_cast<int>(layout.templateoffset<RepresentedValueKind::residual,Carrier>()); constintborder_row=static_cast<int>(layout.templateoffset<RepresentedValueKind::residual,Border>()); constintcarrier_column=static_cast<int>(layout.templateoffset<RepresentedValueKind::value,Carrier>()); constintborder_column=static_cast<int>(layout.templateoffset<RepresentedValueKind::value,Border>()); returnContributions { .residual={{carrier_row,std::vector<double>(carrier.GetData(),carrier.GetData()+carrier_size)},{border_row,{phase[0]}}}, .jacobian={{carrier_row,border_column,execution::linearized_block(std::shared_ptr<const::mfem::SparseMatrix>{std::move(force)})}, {border_row,carrier_column,execution::linearized_block(std::shared_ptr<const::mfem::SparseMatrix>{std::move(constraint)})}}, .replacements={} }; }); } template<typenameScope=void,operators::backend::mfem::AffinePointEquationEquation,typenameSchema,discretization::blocks::SystemLayoutTypeLayout> requiresutils::dependencies::ContainsKey<typenameEquation::Test,typenameLayout::Values>&& utils::dependencies::ContainsKey<typenameEquation::Test,typenameLayout::Residuals>&& utils::dependencies::ContainsKey<typenameEquation::Border,typenameLayout::Values>&& utils::dependencies::ContainsKey<typenameEquation::Border,typenameLayout::Residuals> [[nodiscard]]automake_point_constraint_computations(Equationequation,constoperators::backend::mfem::PointBorderDiscretization<Equation,Schema>&discretization,Layoutlayout) { usingPointScope=std::conditional_t<std::same_as<Scope,void>,typenameEquation::Border,Scope>; returnstd::tuple_cat(operators::backend::mfem::make_point_border_computations<PointScope>(std::move(equation),discretization), std::tuple{point_constraint_contribution<PointScope, Equation, Schema>(std::move(layout))}); } }//namespaceserif::model::computations