point_border.hpp serif/operators/backend/mfem/nodal_rows.hpp src/include/serif/model/computations/point_constraints.hpp serif/dimensions/dimensions.hpp serif/dimensions/quantities.hpp serif/discretization/backend/mfem/dof_map.hpp serif/discretization/backend/mfem/dof_selection.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/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/equations/concepts.hpp serif/equations/contributions.hpp serif/equations/point.hpp serif/meta/enums.hpp serif/meta/names.hpp serif/operators/backend/mfem/nodal_data.hpp serif/operators/backend/mfem/nodal_rows.hpp src/include/serif/operators/backend/mfem/point_border.hpp serif/operators/concepts.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/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 iostream magic_enum/magic_enum_all.hpp memory mfem/fem/fespace.hpp mfem/fem/pfespace.hpp nameof.hpp numeric optional span stdexcept string string_view tuple type_traits unistd.h unordered_map utility 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 src/include/serif/operators/backend/mfem/point_border.hpp serif serif::operators serif::operators::backend serif::operators::backend::mfem #pragmaonce #include"serif/operators/backend/mfem/nodal_rows.hpp" namespaceserif::operators::backend::mfem{ template<AffinePointEquationEquationT,discretization::domain::schema::IsSchemaSchemaT> voidvalidate_point_border_outputs( constPointBorderData<EquationT, SchemaT>&data, const::mfem::Vector&carrier, const::mfem::Vector&phase ){ usingnamespaceutils::errors; validate_nodal_output(data.nodal,carrier); data.border_binding.validate(); serif_ensure<SERiFErrorCode::INVALID_DOF_COUNT,std::invalid_argument>( phase.Size()==1,"Pointphaseoutputmustcontainonescalar."); } template<AffinePointEquationEquationT,discretization::domain::schema::IsSchemaSchemaT> voidadd_point_border_residual( constPointBorderData<EquationT, SchemaT>&data, ::mfem::Vector&carrier, ::mfem::Vector&phase ){ validate_point_border_outputs(data,carrier,phase); constintpoint=data.nodal.discretization.selected_rows->reduced_dofs()[0]; carrier[point]+=data.border_coefficient*data.nodal.derivative(0,0); phase[0]=data.nodal.residual[0]; data.nodal.discretization.residual_binding.validate_values(carrier); data.border_binding.validate_values(phase); } template<AffinePointEquationEquationT,discretization::domain::schema::IsSchemaSchemaT,typenameCoefficientsT> voidapply_point_border_jacobian( constPointBorderData<EquationT, SchemaT>&data, constCoefficientsT&direction, ::mfem::Vector&carrier, ::mfem::Vector&phase ){ usingnamespaceutils::errors; validate_point_border_outputs(data,carrier,phase); constauto&input=direction.templateget<typenameEquationT::Test>(); constauto&border=direction.templateget<typenameEquationT::Border>(); serif_ensure<SERiFErrorCode::CACHE_FE_BINDING_ERROR,std::invalid_argument>( &input.binding==&data.nodal.discretization.residual_binding&&&border.binding==&data.border_binding, "Pointborderdirectionsmustusetheacceptedcoefficientbindings."); data.nodal.discretization.residual_binding.validate_values(input.values); data.border_binding.validate_values(border.values); constintpoint=data.nodal.discretization.selected_rows->reduced_dofs()[0]; constdoublecarrier_direction=input.values[point]; constdoubleborder_direction=border.values[0]; carrier[point]+=data.nodal.derivative(0,0)*border_direction; phase[0]=data.nodal.derivative(0,0)*carrier_direction; data.nodal.discretization.residual_binding.validate_values(carrier); data.border_binding.validate_values(phase); } template<AffinePointEquationEquationT,discretization::domain::schema::IsSchemaSchemaT> voidapply_point_border_jacobian_transpose( constPointBorderData<EquationT, SchemaT>&data, const::mfem::Vector&carrier_dual, const::mfem::Vector&phase_dual, ::mfem::Vector&carrier_result, ::mfem::Vector&border_result ){ validate_point_border_outputs(data,carrier_result,border_result); data.nodal.discretization.residual_binding.validate_values(carrier_dual); data.border_binding.validate_values(phase_dual); constintpoint=data.nodal.discretization.selected_rows->reduced_dofs()[0]; constdoublepoint_dual=carrier_dual[point]; constdoublescalar_dual=phase_dual[0]; carrier_result[point]+=data.nodal.derivative(0,0)*scalar_dual; border_result[0]+=data.nodal.derivative(0,0)*point_dual; data.nodal.discretization.residual_binding.validate_values(carrier_result); data.border_binding.validate_values(border_result); } }