gravity_force.hpp concepts string_view serif/discretization/blocks/fields.hpp serif/equations/contributions.hpp serif/equations/point.hpp src/compile_time_checks/serif/operators/plan.cpp src/include/serif/equations/displacement_force.hpp 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 serif/equations/contributions.hpp src/include/serif/equations/gravity_force.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/model/execution_validation.cpp src/compile_time_checks/serif/model/normalization.cpp src/compile_time_checks/serif/operators/plan.cpp src/include/serif/equations/displacement_force.hpp src/include/serif/equations/equations.hpp src/include/serif/equations/gravity_force.hpp serif::equations::GravityForce serif serif::equations File containing the gravitational force contribution to the mechanical balance equation (displacement_force.hpp). Emily M. Boudreaux #pragmaonce #include<concepts> #include<string_view> #include"serif/discretization/blocks/fields.hpp" #include"serif/equations/contributions.hpp" #include"serif/equations/point.hpp" namespaceserif::equations{ structGravityForce:Contribution<discretization::forms::Operand<discretization::blocks::displacement::vector>,discretization::domain::StellarDomains>{ usingDensity=discretization::blocks::density::mass; usingGradient=discretization::blocks::gravity::gradient; usingDensityValue=discretization::forms::Operand<Density>; usingGradientValue=discretization::forms::Operand<Gradient>; usingOperands=utils::types::TypeList<DensityValue, GradientValue>; usingResult=discretization::forms::VectorValue; staticconstexprstd::string_viewname="gravitationalforce"; [[nodiscard]]staticPointLinearization<Operands, Result>evaluate( constEquationPoint<Operands>&point ){ constdoubledensity=point.fields.get<Density>(); constauto&gradient=point.fields.get<Gradient>(); PointLinearization<Operands, Result>result; for(intaxis=0;axis<3;++axis){ result.residual[axis]=density*gradient[axis]; result.field_derivatives.get<Density>()[axis]=gradient[axis]; result.field_derivatives.get<Gradient>()[axis][axis]=density; } returnresult; } template<typenameOperandT> [[nodiscard]]staticintadditional_order(OperandT/*operand_type*/,int/*position_order*/){ ifconstexpr(std::same_as<OperandT,DensityValue>){ //TheregisteredRTfamilyoforderphasvectorcomponentsofdegreeuptop+1. returnGradient::FunctionalSpace::order+1; }else{ returnDensity::FunctionalSpace::order; } } }; }//namespaceserif::equations