Planets.hpp

Go to the documentation of this file.

// -------------------------------------------------------//
//
// SHAMROCK code for hydrodynamics
// Copyright (c) 2021-2026 Timothée David--Cléris <tim.shamrock@proton.me>
// SPDX-License-Identifier: CeCILL Free Software License Agreement v2.1
// Shamrock is licensed under the CeCILL 2.1 License, see LICENSE for more information
//
// -------------------------------------------------------//
 
#pragma once
 
#include "shambackends/sycl.hpp"
#include "shamunits/Constants.hpp"
#include "shamunits/UnitSystem.hpp"
 
namespace shamphys {
 
    template<class T>
    T hill_radius(T R, T m, T M) {
        return R * sycl::cbrt(m / (3 * M));
    }
 
    template<class T>
    T keplerian_speed(T G, T M, T R) {
        return sycl::sqrt(G * M / R);
    }
 
    template<class T, class Tu>
    T keplerian_speed(T M, T R, const shamunits::UnitSystem<Tu> usys = {}) {
        return keplerian_speed(shamunits::Constants{usys}.G(), M, R);
    }
 
} // namespace shamphys