---

mmath.h

#ifndef __MATH_HPP__
#define __MATH_HPP__

#include <magic/mobject.h>
#include <magic/mpackarray.h>

// In SunOS/Solaris there is a fp-exception class in math.h.... SHINY BLOODY IDIOTS!
//#ifdef __P
//#undef __P
//#endif
//#ifdef SunOS
//#define exception exception_kludge
//#include <math.h>
//#undef exception
//#else
#include <math.h>
//#endif

int     rnd         (int range);

double  frnd        ();

double  gaussrnd    (double stdv);

inline double   sqr (double x) {return x*x;}

//double quadatan (double y, double x);

template <class T>
inline T sqr (T x) {
    return x*x;
}

inline double sigmoid (double x) {return 1/(1+exp(-x));}

// |x|
//inline double abs (double a) {
//  return a>=0? a:-a;
//}

inline void swap (double& x, double& y) {double tmp=x; x=y; y=tmp;}

#define UNDEFINED_FLOAT 1.2345E30

inline bool is_undef (double x) {
    return x>=1.2345E29;
}

typedef PackArray<double>   Vector;

double  sum         (const Vector& x);
double  min         (const Vector& x);
int     minIndex    (const Vector& x);
double  max         (const Vector& x);
int     maxIndex    (const Vector& x);
double  avg         (const Vector& x);
double  stddev      (const Vector& x);
double  stdmerr     (const Vector& x);

Vector  histogram   (const Vector& x, int n);

void    add (Vector& x, double m);

void    multiply (Vector& x, double m);

void    multiplyToUnity (Vector& x);

#endif

---