/* ----------------------------------------------------------- *//*                                                             *//*                          ___                                *//*                       |_| | |_/   SPEECH                    *//*                       | | | | \   RECOGNITION               *//*                       =========   SOFTWARE                  */ /*                                                             *//*                                                             *//* ----------------------------------------------------------- *//*         Copyright: Microsoft Corporation                    *//*          1995-2000 Redmond, Washington USA                  *//*                    http://www.microsoft.com                 *//*                                                             *//*   Use of this software is governed by a License Agreement   *//*    ** See the file License for the Conditions of Use  **    *//*    **     This banner notice must not be removed      **    *//*                                                             *//* ----------------------------------------------------------- *//*         File: HMath.h:   Math Support                       *//* ----------------------------------------------------------- *//* !HVER!HMath:   3.2 [CUED 09/12/02] */#ifndef _HMATH_H_#define _HMATH_H_#ifdef __cplusplusextern "C" {#endif#ifdef PI#undef PI                /* PI is defined in Linux */#endif#define PI   3.14159265358979#define TPI  6.28318530717959     /* PI*2 */#define LZERO  (-1.0E10)   /* ~log(0) */#define LSMALL (-0.5E10)   /* log values < LSMALL are set to LZERO */#define MINEARG (-708.3)   /* lowest exp() arg  = log(MINLARG) */#define MINLARG 2.45E-308  /* lowest log() arg  = exp(MINEARG) *//* NOTE: On some machines it may be necessary to reduce the         values of MINEARG and MINLARG*/typedef float  LogFloat;   /* types just to signal log values */typedef double LogDouble;typedef enum {  /* Various forms of covariance matrix */   DIAGC,         /* diagonal covariance */   INVDIAGC,      /* inverse diagonal covariance */   FULLC,         /* inverse full rank covariance */   XFORMC,        /* arbitrary rectangular transform */   LLTC,          /* L' part of Choleski decomposition */   NULLC,         /* none - implies Euclidean in distance metrics */   NUMCKIND       /* DON'T TOUCH -- always leave as final element */} CovKind;typedef union {   SVector var;         /* if DIAGC or INVDIAGC */   STriMat inv;         /* if FULLC or LLTC */   SMatrix xform;       /* if XFORMC */} Covariance;/* ------------------------------------------------------------------- */void InitMath(void);/*   Initialise the module*//* ------------------ Vector Oriented Routines ----------------------- */void ZeroShortVec(ShortVec v);void ZeroIntVec(IntVec v);void ZeroVector(Vector v);void ZeroDVector(DVector v);/*   Zero the elements of v*/void CopyShortVec(ShortVec v1, ShortVec v2);void CopyIntVec(IntVec v1, IntVec v2);void CopyVector(Vector v1, Vector v2);void CopyDVector(DVector v1, DVector v2);/*   Copy v1 into v2; sizes must be the same*/Boolean ReadShortVec(Source *src, ShortVec v, Boolean binary);Boolean ReadIntVec(Source *src, IntVec v, Boolean binary);Boolean ReadVector(Source *src, Vector v, Boolean binary);/*   Read vector v from source in ascii or binary*/void WriteShortVec(FILE *f, ShortVec v, Boolean binary);void WriteIntVec(FILE *f, IntVec v, Boolean binary);void WriteVector(FILE *f, Vector v, Boolean binary);/*   Write vector v to stream f in ascii or binary*/void ShowShortVec(char * title, ShortVec v,int maxTerms);void ShowIntVec(char * title, IntVec v,int maxTerms);void ShowVector(char * title,Vector v,int maxTerms);void ShowDVector(char * title,DVector v,int maxTerms);/*   Print the title followed by upto maxTerms elements of v*//* Quadratic prod of a full square matrix C and an arbitry full matrix transform A */void LinTranQuaProd(Matrix Prod, Matrix A, Matrix C);/* ------------------ Matrix Oriented Routines ----------------------- */void ZeroMatrix(Matrix m);void ZeroDMatrix(DMatrix m);void ZeroTriMat(TriMat m);/*   Zero the elements of m*/void CopyMatrix (Matrix m1,  Matrix m2);void CopyDMatrix(DMatrix m1, DMatrix m2);void CopyTriMat (TriMat m1,  TriMat m2);/*   Copy matrix m1 to m2 which must have identical dimensions*/void Mat2DMat(Matrix m1,  DMatrix m2);void DMat2Mat(DMatrix m1, Matrix m2);void Mat2Tri (Matrix m1,  TriMat m2);void Tri2Mat (TriMat m1,  Matrix m2);/*   Convert matrix format from m1 to m2 which must have identical    dimensions*/Boolean ReadMatrix(Source *src, Matrix m, Boolean binary);Boolean ReadTriMat(Source *src, TriMat m, Boolean binary);/*   Read matrix from source into m using ascii or binary.   TriMat version expects m to be in upper triangular form   but converts to lower triangular form internally.*/   void WriteMatrix(FILE *f, Matrix m, Boolean binary);void WriteTriMat(FILE *f, TriMat m, Boolean binary);/*    Write matrix to stream in ascii or binary.  TriMat version     writes m in upper triangular form even though it is stored    in lower triangular form!*/void ShowMatrix (char * title,Matrix m, int maxCols,int maxRows);void ShowDMatrix(char * title,DMatrix m,int maxCols,int maxRows);void ShowTriMat (char * title,TriMat m, int maxCols,int maxRows);/*   Print the title followed by upto maxCols elements of upto   maxRows rows of m.*//* ------------------- Linear Algebra Routines ----------------------- */LogFloat CovInvert(TriMat c, Matrix invc);/*   Computes inverse of c in invc and returns the log of Det(c),   c must be positive definite.*/LogFloat CovDet(TriMat c);/*   Returns log of Det(c), c must be positive definite.*//* ------------- Singular Value Decomposition Routines --------------- */void SVD(DMatrix A, DMatrix U,  DMatrix V, DVector d);/*    Singular Value Decomposition (based on MESCHACH)   A is m x n ,  U is m x n,  W is diag N x 1, V is n x n*/void InvSVD(DMatrix A, DMatrix U, DVector W, DMatrix V, DMatrix Result);/*    Inverted Singular Value Decomposition (calls SVD)   A is m x n ,  U is m x n,  W is diag N x 1, V is n x n, Result is m x n *//* ------------------- Log Arithmetic Routines ----------------------- */LogDouble LAdd(LogDouble x, LogDouble y);/*   Return x+y where x and y are stored as logs,    sum < LSMALL is floored to LZERO */LogDouble LSub(LogDouble x, LogDouble y);/*   Return x-y where x and y are stored as logs,    diff < LSMALL is floored to LZERO */double L2F(LogDouble x);/*   Convert log(x) to real, result is floored to 0.0 if x < LSMALL *//* ------------------- Random Number Routines ------------------------ */void RandInit(int seed);/*    Initialise random number generators, if seed is -ve, then system    clock is used.  RandInit(-1) is called by InitMath.*/float RandomValue(void);/*   Return a random number in range 0.0->1.0 with uniform distribution*/float GaussDeviate(float mu, float sigma);/*   Return a random number with a N(mu,sigma) distribution*/#ifdef __cplusplus}#endif#endif  /* _HMATH_H_ *//* ------------------------- End of HMath.h -------------------------- */