companion.hpp cmath serif/equations/rigid_rotation.hpp serif/utils/error/errors.hpp src/include/serif/equations/barotropic_density.hpp src/include/serif/equations/equations.hpp src/include/serif/equations/surface_potential.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/companion.hpp serif/equations/rigid_rotation.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 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/compile_time_checks/serif/model/validation.cpp src/compile_time_checks/serif/operators/coordinate_identity_meanfield.cpp src/compile_time_checks/serif/operators/plan.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/serialize/top.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/barotropic_density.hpp src/include/serif/equations/companion.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 src/lib/serif/backend/mfem/serialize/top.cpp serif::equations::PointMassCompanion serif serif::equations #pragmaonce #include<cmath> #include"serif/equations/rigid_rotation.hpp" #include"serif/utils/error/errors.hpp" namespaceserif::equations{ classPointMassCompanionfinal{ public: usingVector=RigidRotation::Vector; PointMassCompanion( constdoublegravitational_mass, constVector&position, constbooltidal_only=false ): m_gravitational_mass(gravitational_mass), m_position(position), m_tidal_only(tidal_only){ usingnamespaceutils::errors; serif_ensure<SERiFErrorCode::INVALID_ARGUMENT,std::invalid_argument>(std::isfinite(gravitational_mass)&&gravitational_mass>0.0&& std::isfinite(position[0])&&std::isfinite(position[1])&&std::isfinite(position[2]), "Acompanionrequiresafinitepositivemassandafiniteposition."); } [[nodiscard]]doublepotential(constVector&x)const{ doubleresult=-m_gravitational_mass/std::hypot(x[0]-m_position[0],x[1]-m_position[1],x[2]-m_position[2]); if(m_tidal_only){ result+=uniform_scale()*(x[0]*m_position[0]+x[1]*m_position[1]+x[2]*m_position[2]); } returnresult; } [[nodiscard]]Vectorpotential_gradient(constVector&x)const{ constdoubledistance=std::hypot(x[0]-m_position[0],x[1]-m_position[1],x[2]-m_position[2]); constdoublescale=m_gravitational_mass/(distance*distance*distance); Vectorresult{scale*(x[0]-m_position[0]),scale*(x[1]-m_position[1]),scale*(x[2]-m_position[2])}; if(m_tidal_only){ for(intaxis=0;axis<3;++axis){result[axis]+=uniform_scale()*m_position[axis];} } returnresult; } [[nodiscard]]doublegravitational_mass()constnoexcept{returnm_gravitational_mass;} [[nodiscard]]constVector&position()constnoexcept{returnm_position;} [[nodiscard]]booltidal_only()constnoexcept{returnm_tidal_only;} private: [[nodiscard]]doubleuniform_scale()const{ constdoubledistance=std::hypot(m_position[0],m_position[1],m_position[2]); returnm_gravitational_mass/(distance*distance*distance); } doublem_gravitational_mass; Vectorm_position; boolm_tidal_only; }; }