mathfunc.h

Go to the documentation of this file.

#ifndef mathfunc_h
#define mathfunc_h
 
#include "gadget.h"
 
const double pivalue = atan(1.0) * 4.0;
 
inline int max(int a, int b) {
  return ((a > b) ? a : b);
}
 
inline double max(double a, double b) {
  return ((a > b) ? a : b);
}
 
inline int min(int a, int b) {
  return ((a < b) ? a : b);
}
 
inline double min(double a, double b) {
  return ((a < b) ? a : b);
}
 
inline int isZero(double a) {
  return ((fabs(a) < verysmall) ? 1 : 0);
}
 
inline int isEqual(double a, double b) {
  return ((fabs(a - b) < verysmall) ? 1 : 0);
}
 
inline int isSmall(double a) {
  return ((fabs(a) < rathersmall) ? 1 : 0);
}
 
inline double logFactorial(double n) {
  return lgamma(n + 1.0);
}
 
inline double expRep(double n) {
  if (n > verysmall)
    return 1.0;
  else if (n < -25.0)
    return rathersmall;
  return exp(n);
}
 
inline double randomNumber() {
  int r = rand();
  double k = r % 32767;
  return (k / 32767.0);
}
 
inline double calcMortality(double pred, double pop, double t) {
  if (pred < verysmall)
    return 0.0;
  else if (pred > pop)
    return verybig;
  return (-log(1.0 - (pred / pop)) * t);
}
 
#endif