Utilities.hxx

Go to the documentation of this file.

/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *\
 * Copyright (c) 2025, Davide Stocco and Enrico Bertolazzi.                                      *
 *                                                                                               *
 * The Astro project is distributed under the BSD 2-Clause License.                              *
 *                                                                                               *
 * Davide Stocco                                                               Enrico Bertolazzi *
 * University of Trento                                                     University of Trento *
 * e-mail: davide.stocco@unitn.it                             e-mail: enrico.bertolazzi@unitn.it *
\* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
 
#pragma once
 
#ifndef ASTRO_UTILITIES_HXX
#define ASTRO_UTILITIES_HXX
 
namespace Astro {
 
  /*\
   |   ____
   |  |  _ \ _____      _____ _ __
   |  | |_) / _ \ \ /\ / / _ \ '__|
   |  |  __/ (_) \ V  V /  __/ |
   |  |_|   \___/ \_/\_/ \___|_|
   |
  \*/

  Real power2(Real x) {return x*x;}

  Real power3(Real x) {return x*x*x;}

  Real power4(Real x) {return x*x*x*x;}
 
  /*\
   |      _                _
   |     / \   _ __   __ _| | ___
   |    / _ \ | '_ \ / _` | |/ _ \
   |   / ___ \| | | | (_| | |  __/
   |  /_/   \_\_| |_|\__, |_|\___|
   |                 |___/
  \*/
 
  static Real const PI      = Real(3.141592653589793238462643383279502884197); 
  static Real const PIMUL2  = Real(6.283185307179586476925286766559005768394); 
  static Real const PIDIV2  = Real(1.570796326794896619231321691639751442098); 
  static Real const DEG2RAD = Real(0.017453292519943295769236907684886127134); 
  static Real const RAD2DEG = Real(57.29577951308232087679815481410517033240); 

  Real deg_to_rad(Real x) {return DEG2RAD*x;}

  Real rad_to_deg(Real x) {return RAD2DEG*x;}

  Real angle_in_range(Real x)
  {
    x = std::fmod(x, PIMUL2);
    while (x < Real(0.0)) {x += PIMUL2;}
    while (x > PIMUL2) {x -= PIMUL2;}
    return x;
  }

  Real angle_in_range_sym(Real x)
  {
    x = std::fmod(x, PIMUL2);
    while (x < -PI) {x += PIMUL2;}
    while (x > PI) {x -= PIMUL2;}
    return x;
  }
 
  /*\
   |   ____  _     _
   |  |  _ \‍(_)___| |_ __ _ _ __   ___ ___
   |  | | | | / __| __/ _` | '_ \ / __/ _ \
   |  | |_| | \__ \ || (_| | | | | (_|  __/
   |  |____/|_|___/\__\__,_|_| |_|\___\___|
   |
  \*/
 
  static Real const AU_TO_KM{1.49597870707e+08};   
  static Real const LY_TO_KM{9.4607304725808e+12}; 
  static Real const PC_TO_AU{6.48e+05/PI};         
  static Real const PC_TO_KM{PC_TO_AU*AU_TO_KM};   
  static Real const AU_TO_PC{1.0/PC_TO_AU};        
  static Real const LY_TO_PC{1.0/PC_TO_AU};        
  static Real const KM_TO_PC{1.0/PC_TO_KM};        
  static Real const PC_TO_LY{1.0/LY_TO_PC};        
  static Real const KM_TO_LY{1.0/LY_TO_KM};        
  static Real const LY_TO_AU{LY_TO_KM/AU_TO_KM};   
  static Real const AU_TO_LY{1.0/LY_TO_AU};        
  static Real const KM_TO_AU{1.0/AU_TO_KM};        
 
  static Real const KM_TO_M{1000.0};           
  static Real const M_TO_KM{1.0/KM_TO_M};      
  static Real const AU_TO_M{AU_TO_KM*KM_TO_M}; 
  static Real const LY_TO_M{LY_TO_KM*KM_TO_M}; 
  static Real const PC_TO_M{LY_TO_KM*KM_TO_M}; 
  static Real const M_TO_PC{LY_TO_KM*KM_TO_M}; 
  static Real const M_TO_LY{LY_TO_KM*KM_TO_M}; 
  static Real const M_TO_AU{LY_TO_KM*KM_TO_M}; 
 
 
  // Time units
  static Real const DAY_TO_SEC{86400.0};        
  static Real const SEC_TO_DAY{1.0/DAY_TO_SEC}; 
 
  static Real const KG_M_SEC2_TO_KG_AU_DAY2{M_TO_AU/(SEC_TO_DAY*SEC_TO_DAY)}; // Kg * m/s^2 => Kg * UA / day^2
 
  static Real const gravity_kg_m_s2{9.80665};
  static Real const gravity_kg_AU_DAY2{gravity_kg_m_s2*KG_M_SEC2_TO_KG_AU_DAY2};
 
 
  static Real const muSun_km3s2{1.32712440018E11}; // Km^3/s^2
  static Real const muSun_AU3DAY2{muSun_km3s2*(KM_TO_AU*KM_TO_AU*KM_TO_AU)/(SEC_TO_DAY*SEC_TO_DAY)}; // Km^3/s^2
//#define ASTRO_MU_SUN 1.32712440018e20
//#define ASTRO_MU_EARTH 398600441800000.0
//#define ASTRO_EARTH_VELOCITY 29784.6905
//#define ASTRO_EARTH_J2 1.08262668E-03
//#define ASTRO_EARTH_RADIUS 6378137
  // - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
 
  Real AU_to_KM(Real x) {return x * AU_TO_KM;}
 
  // - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
 
  Real AU_by_DAY_to_km_by_s(Real x) {return x * (AU_TO_KM/DAY_TO_SEC);}
 
  // - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
 
  Real AU_by_DAY2_to_km_by_s2(Real x) {return x * (AU_TO_KM/(DAY_TO_SEC*DAY_TO_SEC));}
 
  // - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
 
  Real AU_by_DAY3_to_km_by_s3(Real x) {return x * (AU_TO_KM/(DAY_TO_SEC*DAY_TO_SEC*DAY_TO_SEC));}
 
} // namespace Astro
 
#endif // ASTRO_UTILITIES_HXX