---

mj_math.h

Go to the documentation of this file.

/* -*- C++ -*- */

/* COPYRIGHT
 *
 * This file is part of Mustajuuri LGPL modules.
 *
 * Author: Tommi Ilmonen, 1999-2001.
 * Tommi.Ilmonen@hut.fi
 *
 * http://www.tml.hut.fi/~tilmonen/mustajuuri/
 *
 *
 * This file is licensed under the Lesser GNU Public License
 * (LGPL). Extra clauses can be found from the file MUSTAJUURI_LGPL and
 * further explanations from the file MUSTAJUURI_LICENSE. The license
 * text can be found in file "LGPL.txt".
 *
 * If you want a parallel license (for commercial reasons for example),
 * you should negotiate the matter with the author(s).
 *    
 */

#ifndef MJ_MATH_H
#define MJ_MATH_H

#include <math.h>



#define MJ_pi  3.1415926535897931

#define MJ_2pi 6.2831853071795862

#define MJ_piSquared 9.869604401089358

#define MJ_pi2 1.57079632679489660000

#define MJ_pi4 0.78539816339744828

#define MJ_ln2 0.693147180559945290

#define MJ_ln22 0.346573590279972640


/* Float wrappers. */

inline float MJ_cos(float v)
{
  return cosf(v);
}

inline float MJ_log(float v)
{
  return logf(v);
}

inline float MJ_exp(float v)
{
  return expf(v);
}

inline float MJ_pow(float v, float power)
{
  return powf(v, power);
}

inline float MJ_square(float v)
{
  return v * v;
}

inline float MJ_sin(float v)
{
  return sinf(v);
}

inline float MJ_sinh(float v)
{
  return sinhf(v);
}

inline float MJ_sqrt(float v)
{
  return sqrtf(v);
}

inline float MJ_tan(float v)
{
  return tanf(v);
}

inline float MJ_atan2(float v, float u)
{
  return atan2f(v,u);
}

/* Double wrappers. */

inline double MJ_cos(double v)
{
  return cos(v);
}

inline double MJ_log(double v)
{
  return log(v);
}

inline double MJ_exp(double v)
{
  return exp(v);
}

inline double MJ_pow(double v, double power)
{
  return pow(v, power);
}

inline double MJ_square(double v)
{
  return v * v;
}

inline double MJ_sin(double v)
{
  return sin(v);
}

inline double MJ_sinh(double v)
{
  return sinh(v);
}

inline double MJ_sqrt(double v)
{
  return sqrt(v);
}

inline double MJ_tan(double v)
{
  return tan(v);
}

inline double MJ_atan2(double v, double u)
{
  return atan2(v,u);
}

// Clamp (limit) to the given range

inline double MJ_clamp(double v, double min, double max)
{
  if(v < min) return min;
  if(v > max) return max;
  return v;
}

inline double MJ_clampD(double v, double min, double max)
{
  if(v < min) return min;
  if(v > max) return max;
  return v;
}

inline float MJ_clamp(float v, float min, float max)
{
  if(v < min) return min;
  if(v > max) return max;
  return v;
}

inline float MJ_clampF(float v, float min, float max)
{
  if(v < min) return min;
  if(v > max) return max;
  return v;
}

inline int MJ_clamp(int v, int min, int max)
{
  if(v < min) return min;
  if(v > max) return max;
  return v;
}


// Minmimum value

inline float MJ_min(float x1, float x2)
{
  if(x1 > x2) return x2;
  return x1;
}

inline double MJ_min(double x1, double x2)
{
  if(x1 > x2) return x2;
  return x1;
}

inline int MJ_min(int x1, int x2)
{
  if(x1 > x2) return x2;
  return x1;
}

inline long MJ_min(long x1, long x2)
{
  if(x1 > x2) return x2;
  return x1;
}

inline unsigned MJ_min(unsigned x1, unsigned x2)
{
  if(x1 > x2) return x2;
  return x1;
}

inline unsigned long MJ_min(unsigned long x1, unsigned long x2)
{
  if(x1 > x2) return x2;
  return x1;
}

// Maximum value

inline float MJ_max(float x1, float x2)
{
  if(x1 > x2) return x1;
  return x2;
}

inline double MJ_max(double x1, double x2)
{
  if(x1 > x2) return x1;
  return x2;
}

inline int MJ_max(int x1, int x2)
{
  if(x1 > x2) return x1;
  return x2;
}

inline unsigned MJ_max(unsigned x1, unsigned x2)
{
  if(x1 > x2) return x1;
  return x2;
}

inline unsigned long MJ_max(unsigned long x1, unsigned long x2)
{
  if(x1 > x2) return x1;
  return x2;
}


inline float MJ_abs(float x1)
{
  if(x1 > 0.0) return x1;
  else return -x1;
}

inline double MJ_abs(double x1)
{
  if(x1 > 0.0) return x1;
  else return -x1;
}

inline int MJ_abs(int x1)
{
  if(x1 > 0) return x1;
  else return -x1;
}

inline long MJ_abs(long x1)
{
  if(x1 > 0) return x1;
  else return -x1;
}

inline unsigned long MJ_abs(unsigned long x1)
{
  return x1;
}

//@name Decibel conversions

inline double MJ_absToDb(double d)
{ return (log10(d)) * 20.0; }
inline double MJ_squareToDb(double d)
{ return (log10(d)) * 10.0; }
inline double MJ_absToDb(double d, double clip)
{
  double tmp = (log10(d)) * 20.0;
  return tmp > clip ? tmp : clip;
}
inline double MJ_dbToAbs(double d)
{  return pow(10.0, d * 0.05); }
inline double MJ_dbToSquare(double d)
{ return pow(10.0, d * 0.1); }

inline float MJ_absToDb(float d)
{ return (log10f(d)) * 20.0; }
inline float MJ_squareToDb(float d)
{ return (log10f(d)) * 10.0; }
inline float MJ_dbToAbs(float d)
{ return powf(10.0, d * 0.05); }
inline float MJ_dbToSquare(float d)
{ return powf(10.0, d * 0.1); }

inline float MJ_centToMul(float cents)
{
  return powf(2.0, cents / 1200.0);
}
inline float MJ_mulToCent(float mul)
{ return logf(mul) / MJ_ln2 * 1200.0; }

inline double MJ_centToMul(double cents)
{ return pow(2.0, cents / 1200.0); }
inline double MJ_mulToCent(double mul)
{ return log(mul) / MJ_ln2 * 1200.0; }

#endif

---