core.h

Go to the documentation of this file.

#ifndef SRC_CLIMATH_CORE_H_
#define SRC_CLIMATH_CORE_H_
 
#ifdef __cplusplus
extern "C" {
#endif
 
#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; }
 
#ifdef __cplusplus
} /* extern "C" */
#endif
 
#endif  // SRC_CLIMATH_CORE_H_