core.h

Go to the documentation of this file.

#ifndef CORE_H_
#define CORE_H_
#include <math.h>
 
struct float_triplet {
  double x, y, z;
};
 
inline double sqr(double x) { 
  return x*x; 
}
inline double cub(double x) {
  return x*x*x;
}
inline double min(double x1, double x2, double x3) {
  return fmin(x1, fmin(x2, x3));
}
inline double max(double x1, double x2, double x3) {
  return fmax(x1, fmax(x2, x3));
}
inline double allmax(double *a, int n) {
  double v = a[0];
  for (int i = 1; i < n; ++i)
    if (v < a[i]) v = a[i];
  return v;
}
inline double allmin(double *a, int n) {
  double v = a[0];
  for (int i = 1; i < n; ++i)
    if (v > a[i]) v = a[i];
  return v;
}
inline int sign(double x) {
  return x < 0. ? -1 : 1;
}
 
int fcmp(double x1, double x2);
 
 
// unit conversion
 
inline double rad2deg(double phi) { return phi*180./M_PI; }
inline double deg2rad(double phi) { return phi*M_PI/180.; }
 
inline double km2m(double x) { return x*1.E3; }
inline double m2km(double x) { return x/1.E3; }
 
inline double day2sec(double x) { return x*86400.; }
inline double sec2day(double x) { return x/86400.; }
 
inline double au2m(double x) { return x*1.495978707E11; }
inline double m2au(double x) { return x/1.495978707E11; }
 
#endif