include/math_utils.h

Go to the documentation of this file.

#if !defined (MATH_UTILS_H)
#define MATH_UTILS_H


/*Lapack solver for real symmetric matrices*/
extern "C"{

  int dsyev_(char* jobz, char* uplo, int *n, double *a,
             int *lda, double *w, double *work, int *lwork, 
             int *info);
  
  int dsymm_(char* side, char* uplo, int* m, int *n, double *alpha, 
             double* A, int *lda, double *B, int *ldb, double *beta, 
             double *C, int *ldc);

  int dpotrf_(char* uplo, int* n, double* a, int* lda, int* info);
 
}



namespace utils
{
  void get_eigensystem(double *eigenvalues, double* matrix, 
                       double* eigenvectors, int n);
  double* get_cholesky(double* matrix, int n) ;
  void mat_sqrt(double *A, int n);
  
  void mat_mult(double* A,double* B,double* C,int n);
  inline double max(double a,double b){return  (b<a)?a:b;}

}

#endif

---