nodal_evaluation.hpp cmath expected serif/operators/backend/mfem/nodal_data.hpp serif/operators/backend/mfem/rejection.hpp src/include/serif/operators/backend/mfem/nodal_computations.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/mapping/status.hpp serif/meta/enums.hpp serif/meta/names.hpp serif/operators/backend/mfem/nodal_data.hpp src/include/serif/operators/backend/mfem/nodal_evaluation.hpp serif/operators/backend/mfem/rejection.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 expected 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/compile_time_checks/serif/operators/computation_owners.cpp src/compile_time_checks/serif/operators/serial_backend.cpp 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/nodal_computations.hpp src/include/serif/operators/backend/mfem/nodal_evaluation.hpp serif serif::operators serif::operators::backend serif::operators::backend::mfem #pragmaonce #include<cmath> #include<expected> #include"serif/operators/backend/mfem/nodal_data.hpp" #include"serif/operators/backend/mfem/rejection.hpp" namespaceserif::operators::backend::mfem{ template<NodalValueEquationEquationT> [[nodiscard]]equations::EquationPoint<typenameEquationT::Operands>make_nodal_equation_point( const::mfem::Vector&values,constintfirst ){ usingValue=discretization::forms::OperandValue<equations::TestOperandOf<EquationT>>; equations::EquationPoint<typenameEquationT::Operands>point; for(intaxis=0;axis<discretization::forms::component_count<Value>;++axis){ discretization::forms::component(point.fields.templateget<typenameEquationT::Test>(),axis)=values[first+axis]; } returnpoint; } template<NodalValueEquationEquationT,discretization::domain::schema::IsSchemaSchemaT> [[nodiscard]]std::expected<void,PreparationRejection>store_nodal_linearization( constequations::PointLinearization<typenameEquationT::Operands, discretization::forms::OperandValue<equations::TestOperandOf<EquationT>>>&physics, constintfirst,NodalData<EquationT, SchemaT>&data ){ usingnamespaceutils::errors; usingValue=discretization::forms::OperandValue<equations::TestOperandOf<EquationT>>; constexprintcomponents=discretization::forms::component_count<Value>; for(intoutput=0;output<components;++output){ for(intaxis=0;axis<3;++axis){ serif_ensure<SERiFErrorCode::INVALID_ARGUMENT,std::logic_error>( (discretization::forms::derivative_component<Value, discretization::forms::VectorValue>(physics.position_derivative,output,axis)==0.0), "Anodalequationdeclaredposition-independentmusthavezeropositionderivative."); } data.residual[first+output]=discretization::forms::component(physics.residual,output); boolfinite=std::isfinite(data.residual[first+output]); for(intinput=0;input<components;++input){ constdoublepartial=discretization::forms::derivative_component<Value, Value>( physics.field_derivatives.templateget<typenameEquationT::Test>(),output,input); data.derivative(first+output,input)=partial; finite=finite&&std::isfinite(partial); } if(!finite){ returnstd::unexpected(PreparationRejection{ .reason=PreparationRejectionReason::NON_FINITE,.message="Non-finitenodalresidualorderivative."}); } } return{}; } template<NodalValueEquationEquationT,discretization::domain::schema::IsSchemaSchemaT> [[nodiscard]]std::expected<NodalData<EquationT,SchemaT>,PreparationRejection>evaluate_nodal_values( constEquationT&equation,constNodalDiscretization<EquationT, SchemaT>&discretization,const::mfem::Vector&values ){ usingValue=discretization::forms::OperandValue<equations::TestOperandOf<EquationT>>; constexprintcomponents=discretization::forms::component_count<Value>; NodalData<EquationT, SchemaT>data{ .discretization=discretization, .residual=::mfem::Vector(values.Size()), .derivative=::mfem::DenseMatrix(values.Size(),components) }; for(intfirst=0;first<values.Size();first+=components){ constautopoint=make_nodal_equation_point<EquationT>(values,first); equations::PointLinearization<typenameEquationT::Operands,Value>physics; try{ physics=equation.evaluate(point); }catch(conststd::domain_error&error){ returnstd::unexpected(PreparationRejection{.reason=PreparationRejectionReason::EQUATION,.message=error.what()}); } constautoaccepted=store_nodal_linearization<EquationT>(physics,first,data); if(!accepted){ returnstd::unexpected(accepted.error()); } } returndata; } }