---

mmath.h

Go to the documentation of this file.

/***************************************************************************
    copyright            : (C) 2000 by Marko Grönroos
    email                : magi@iki.fi
 ***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************
 *
 **/

#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

BEGIN_NAMESPACE (MagiC);


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);

#ifndef ROUND_DOUBLE_THRESHOLD
#define ROUND_DOUBLE_THRESHOLD 1e-15
#endif

inline bool isRoundZero (double x) {
    return fabs(x) < ROUND_DOUBLE_THRESHOLD;
}

END_NAMESPACE;

#endif

---