surface_potential.hpp optional string_view serif/discretization/blocks/fields.hpp serif/equations/companion.hpp serif/equations/point.hpp serif/equations/rigid_rotation.hpp src/compile_time_checks/serif/model/barotropic_density.cpp src/compile_time_checks/serif/model/density_primal.cpp src/compile_time_checks/serif/model/execution_validation.cpp src/compile_time_checks/serif/model/normalization.cpp src/include/serif/backend/mfem/model/surface_potential.hpp src/include/serif/equations/equations.hpp src/include/serif/model/equilibrium_plan.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/companion.hpp serif/equations/point.hpp serif/equations/rigid_rotation.hpp src/include/serif/equations/surface_potential.hpp serif/meta/enums.hpp serif/meta/names.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/std_helper/cleaning.hpp serif/utils/misc/std_helper/variadic_cleaning.hpp serif/utils/misc/terminal/colors.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 nameof.hpp numeric optional span stdexcept string string_view tuple type_traits unistd.h unordered_map utility version src/compile_time_checks/serif/backend_independent.cpp src/compile_time_checks/serif/model/barotropic_density.cpp src/compile_time_checks/serif/model/density_primal.cpp src/compile_time_checks/serif/model/execution_validation.cpp src/compile_time_checks/serif/model/finite_elements.cpp src/compile_time_checks/serif/model/normalization.cpp src/compile_time_checks/serif/model/state.cpp src/include/serif/backend/mfem/model/binding.hpp src/include/serif/backend/mfem/model/computations/surface_potential.hpp src/include/serif/backend/mfem/model/equilibrium.hpp src/include/serif/backend/mfem/model/equilibrium_evaluator.hpp src/include/serif/backend/mfem/model/finite_elements.hpp src/include/serif/backend/mfem/model/graph.hpp src/include/serif/backend/mfem/model/normalization/computations.hpp src/include/serif/backend/mfem/model/surface_potential.hpp src/include/serif/backend/mfem/serialize/configuration.hpp src/include/serif/backend/mfem/serialize/discretization.hpp src/include/serif/backend/mfem/serialize/options.hpp src/include/serif/backend/mfem/solver/admissibility/geometry.hpp src/include/serif/backend/mfem/solver/equilibrium.hpp src/include/serif/backend/mfem/solver/nonlinear/stellar_problem.hpp src/include/serif/backend/mfem/solver/seed/state.hpp src/include/serif/equations/equations.hpp src/include/serif/equations/surface_potential.hpp src/include/serif/model/equilibrium_plan.hpp src/include/serif/model/execution_validation.hpp src/include/serif/model/finite_elements.hpp src/include/serif/model/model.hpp src/include/serif/model/normalization/plan.hpp src/include/serif/solver/continuation/equilibrium.hpp src/include/serif/solver/equilibrium.hpp src/include/serif/solver/seed/state.hpp src/include/serif/solver/solver.hpp src/include/serif/utils/serialize/equations.hpp src/include/serif/utils/serialize/model.hpp src/include/serif/utils/serialize/numerics.hpp src/include/serif/utils/serialize/serialize.hpp src/include/serif/utils/serialize/solver.hpp src/lib/serif/backend/mfem/model/finite_elements.cpp src/lib/serif/backend/mfem/model/surface_potential.cpp serif::equations::SurfacePotentialBalance serif serif::equations #pragmaonce #include<optional> #include<string_view> #include"serif/discretization/blocks/fields.hpp" #include"serif/equations/companion.hpp" #include"serif/equations/point.hpp" #include"serif/equations/rigid_rotation.hpp" namespaceserif::equations{ classSurfacePotentialBalancefinal{ public: usingPotential=discretization::blocks::gravity::potential; usingSurfaceLevel=discretization::blocks::barotropic_constant::mass_normalization; usingAngularVelocity=discretization::blocks::rotation::angular_velocity; usingOperands=utils::types::TypeList<discretization::forms::Operand<Potential>, discretization::forms::Operand<SurfaceLevel>,discretization::forms::Operand<AngularVelocity>>; usingResidualQuantity=dimensions::GravitationalPotential; staticconstexprstd::string_viewname="surfaceeffectivepotentialbalance"; explicitSurfacePotentialBalance(constRigidRotation::Vector&axis={0.0,0.0,1.0}, constRigidRotation::Vector&center={},std::optional<PointMassCompanion>companion=std::nullopt): m_rotation(RigidRotation::unit_rotation(axis,center)),m_companion(companion){} [[nodiscard]]PointLinearization<Operands>evaluate(constEquationPoint<Operands>&point)const{ constdoubleomega=point.fields.get<AngularVelocity>(); constdoublecentrifugal=m_rotation.centrifugal_potential(point.position); PointLinearization<Operands>result; result.residual=point.fields.get<Potential>()-omega*omega*centrifugal-point.fields.get<SurfaceLevel>(); result.field_derivatives.get<Potential>()=1.0; result.field_derivatives.get<SurfaceLevel>()=-1.0; result.field_derivatives.get<AngularVelocity>()=-2.0*omega*centrifugal; result.position_derivative=m_rotation.centrifugal_potential_gradient(point.position); for(double&component:result.position_derivative){ component*=-omega*omega; } if(m_companion){ result.residual+=m_companion->potential(point.position); constautogradient=m_companion->potential_gradient(point.position); for(intaxis=0;axis<3;++axis){result.position_derivative[axis]+=gradient[axis];} } returnresult; } [[nodiscard]]constRigidRotation&rotation()constnoexcept{returnm_rotation;} [[nodiscard]]conststd::optional<PointMassCompanion>&companion()constnoexcept{returnm_companion;} private: RigidRotationm_rotation; std::optional<PointMassCompanion>m_companion; }; }