centering.hpp array string_view serif/equations/contributions.hpp serif/equations/point.hpp src/compile_time_checks/serif/operators/computation_owners.cpp src/compile_time_checks/serif/operators/serial_backend.cpp src/include/serif/equations/equations.hpp serif/dimensions/dimensions.hpp serif/dimensions/quantities.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/physical_domains.hpp serif/discretization/domain/types.hpp serif/discretization/forms/operands.hpp serif/discretization/forms/value.hpp serif/discretization/functional_space/space.hpp src/include/serif/equations/centering.hpp serif/equations/contributions.hpp serif/equations/point.hpp serif/meta/enums.hpp serif/meta/names.hpp serif/utils/misc/concepts/numeric.hpp serif/utils/misc/std_helper/cleaning.hpp serif/utils/misc/std_helper/variadic_cleaning.hpp serif/utils/types/type_list.hpp array compare concepts cstddef format magic_enum/magic_enum_all.hpp nameof.hpp numeric span string_view tuple type_traits unordered_map utility src/compile_time_checks/serif/operators/computation_owners.cpp src/compile_time_checks/serif/operators/serial_backend.cpp src/include/serif/equations/centering.hpp src/include/serif/equations/equations.hpp serif::equations::CenteringConstraint serif serif::equations File containing the CenteringConstraint equation, which fixes the material point at the computational origin. Emily M. Boudreaux This is probably the simplest equation we have right now and therefore may serve as a useful example to learn from. #pragmaonce #include<array> #include<string_view> #include"serif/equations/contributions.hpp" #include"serif/equations/point.hpp" namespaceserif::equations{ structCenteringConstraintfinal:Contribution<discretization::forms::Operand<discretization::blocks::displacement::vector>>{ usingOperands=utils::types::TypeList<TestOperand>; staticconstexprstd::string_viewname="centeringconstraint"; staticconstexprboolposition_independent=true; staticconstexprstd::array<double,3>reference_position{}; [[nodiscard]]staticPointLinearization<Operands, discretization::forms::VectorValue>evaluate(constEquationPoint<Operands>&point){ PointLinearization<Operands, discretization::forms::VectorValue>result; result.residual=point.fields.get<Test>(); for(intaxis=0;axis<3;++axis){ result.field_derivatives.get<Test>()[axis][axis]=1.0; } returnresult; } }; }//namespaceserif::equations