23 inline T schwarzschild_radius(T M, T G, T c) {
24 return 2 * G * M / (c * c);
27 template<
class T,
class Tu>
28 T schwarzschild_radius(T M,
const shamunits::UnitSystem<Tu> usys = {}) {
29 return schwarzschild_radius(
30 M, shamunits::Constants{usys}.G(), shamunits::Constants{usys}.c());