lapack_wrapper.hpp

Go to the documentation of this file.

#ifndef __LAPACK_WRAPPER_H__
#define __LAPACK_WRAPPER_H__
 
#include <stdint.h>
 
#include <complex>
#include <complex.h>
#include <algorithm>
#include "Type.hpp"
 
#ifdef UNI_MKL
  #include <mkl.h>
typedef MKL_INT blas_int;
 
#else
//#include <cblas.h>
//#ifdef BLAS_ILP64
//    typedef int64_t blas_int;
//#else
typedef int32_t blas_int;
  //#endif
  #include <lapacke.h>
  #include <cblas.h>
extern "C" {
 
// BLAS functions
void dgemm_(const char *transa, const char *transb, const blas_int *m, const blas_int *n,
            const blas_int *k, const double *alpha, const double *a, const blas_int *lda,
            const double *b, const blas_int *ldb, const double *beta, double *c,
            const blas_int *ldc);
void sgemm_(const char *transa, const char *transb, const blas_int *m, const blas_int *n,
            const blas_int *k, const float *alpha, const float *a, const blas_int *lda,
            const float *b, const blas_int *ldb, const float *beta, float *c, const blas_int *ldc);
void zgemm_(const char *transa, const char *transb, const blas_int *m, const blas_int *n,
            const blas_int *k, const std::complex<double> *alpha, const std::complex<double> *a,
            const blas_int *lda, const std::complex<double> *b, const blas_int *ldb,
            const std::complex<double> *beta, std::complex<double> *c, const blas_int *ldc);
void cgemm_(const char *transa, const char *transb, const blas_int *m, const blas_int *n,
            const blas_int *k, const std::complex<float> *alpha, const std::complex<float> *a,
            const blas_int *lda, const std::complex<float> *b, const blas_int *ldb,
            const std::complex<float> *beta, std::complex<float> *c, const blas_int *ldc);
 
// void   zgemm_(const char *transa, const char *transb, const blas_int *m, const blas_int *n, const
// blas_int *k,
//               const double _Complex *alpha, const double _Complex *a, const blas_int *lda, const
//               double _Complex *b, const blas_int *ldb, const double _Complex *beta, double
//               _Complex *c, const blas_int *ldc);
 
double dasum_(const blas_int *n, const double *x, const blas_int *incx);
 
void scopy_(const blas_int *n, const float *x, const blas_int *incx, float *y,
            const blas_int *incy);
void dcopy_(const blas_int *n, const double *x, const blas_int *incx, double *y,
            const blas_int *incy);
// void   zcopy_(const blas_int *n, const double *x, const blas_int *incx, std::complex<double> *y,
// const blas_int *incy);
 
void daxpy_(const blas_int *n, const double *alpha, const double *x, const blas_int *incx,
            double *y, const blas_int *incy);
void saxpy_(const blas_int *n, const float *alpha, const float *x, const blas_int *incx, float *y,
            const blas_int *incy);
void zaxpy_(const blas_int *n, const std::complex<double> *alpha, const std::complex<double> *x,
            const blas_int *incx, std::complex<double> *y, const blas_int *incy);
void caxpy_(const blas_int *n, const std::complex<float> *alpha, const std::complex<float> *x,
            const blas_int *incx, std::complex<float> *y, const blas_int *incy);
void dscal_(const blas_int *n, const double *a, double *x, const blas_int *incx);
void sscal_(const blas_int *n, const float *a, float *x, const blas_int *incx);
void zscal_(const blas_int *n, const std::complex<double> *a, std::complex<double> *x,
            const blas_int *incx);
void cscal_(const blas_int *n, const std::complex<float> *a, std::complex<float> *x,
            const blas_int *incx);
void zdscal_(const blas_int *n, const double *a, std::complex<double> *x, const blas_int *incx);
 
float snrm2_(const blas_int *n, const float *x, const blas_int *incx);
float scnrm2_(const blas_int *n, const std::complex<float> *x, const blas_int *incx);
double dnrm2_(const blas_int *n, const double *x, const blas_int *incx);
double dznrm2_(const blas_int *n, const std::complex<double> *x, const blas_int *incx);
 
void sgemv_(const char *trans, const blas_int *m, const blas_int *n, const float *alpha,
            const float *a, const blas_int *lda, const float *x, const blas_int *incx,
            const float *beta, const float *y, const blas_int *incy);
 
void dgemv_(const char *trans, const blas_int *m, const blas_int *n, const double *alpha,
            const double *a, const blas_int *lda, const double *x, const blas_int *incx,
            const double *beta, const double *y, const blas_int *incy);
 
void zgemv_(const char *trans, const blas_int *m, const blas_int *n,
            const std::complex<double> *alpha, const std::complex<double> *a, const blas_int *lda,
            const std::complex<double> *x, const blas_int *incx, const std::complex<double> *beta,
            const std::complex<double> *y, const blas_int *incy);
 
void cgemv_(const char *trans, const blas_int *m, const blas_int *n,
            const std::complex<float> *alpha, const std::complex<float> *a, const blas_int *lda,
            const std::complex<float> *x, const blas_int *incx, const std::complex<float> *beta,
            const std::complex<float> *y, const blas_int *incy);
 
double ddot_(const blas_int *n, const double *x, const blas_int *incx, const double *y,
             const blas_int *incy);
float sdot_(const blas_int *n, const float *x, const blas_int *incx, const float *y,
            const blas_int *incy);
 
  #ifndef FORTRAN_COMPLEX_FUNCTIONS_RETURN_VOID
std::complex<double> zdotc_(const blas_int *n, const std::complex<double> *x, const blas_int *incx,
                            const std::complex<double> *y, const blas_int *incy);
std::complex<double> zdotu_(const blas_int *n, const std::complex<double> *x, const blas_int *incx,
                            const std::complex<double> *y, const blas_int *incy);
std::complex<float> cdotc_(const blas_int *n, const std::complex<float> *x, const blas_int *incx,
                           const std::complex<float> *y, const blas_int *incy);
std::complex<float> cdotu_(const blas_int *n, const std::complex<float> *x, const blas_int *incx,
                           const std::complex<float> *y, const blas_int *incy);
  #else
void zdotc_(std::complex<double> *res, const blas_int *n, const std::complex<double> *x,
            const blas_int *incx, const std::complex<double> *y, const blas_int *incy);
void zdotu_(std::complex<double> *res, const blas_int *n, const std::complex<double> *x,
            const blas_int *incx, const std::complex<double> *y, const blas_int *incy);
void cdotc_(std::complex<float> *res, const blas_int *n, const std::complex<float> *x,
            const blas_int *incx, const std::complex<float> *y, const blas_int *incy);
void cdotu_(std::complex<float> *res, const blas_int *n, const std::complex<float> *x,
            const blas_int *incx, const std::complex<float> *y, const blas_int *incy);
  #endif
 
// D and Z LAPACK functions
/*
void dgeqp3_( const blas_int* m, const blas_int* n, double* a, const blas_int* lda,
              int* jpvt, double* tau, double* work, int* lwork, blas_int* info );
void zgeqp3_( const blas_int* m, const blas_int* n, std::complex<double>* a, const blas_int* lda,
              int* jpvt, std::complex<double>* tau, std::complex<double>* work,
              int* lwork, double* rwork, blas_int* info );
 
//void dgeqpf_( const blas_int* m, const blas_int* n, double* a, const blas_int* lda,
//              int* jpvt, double* tau, double* work, blas_int* info );
//void zgeqpf_( const blas_int* m, const blas_int* n, std::complex<double>* a, const blas_int* lda,
//              int* jpvt, std::complex<double>* tau, std::complex<double>* work,
//              double* rwork, blas_int* info );
//
 
void dgesvd_( const char* jobu, const char* jobvt, const blas_int* m,
              const blas_int* n, double* a, const blas_int* lda, double* s,
              double* u, const blas_int* ldu, double* vt, const blas_int* ldvt,
              double* work, const blas_int* lwork, blas_int* info );
 
void sgesvd_( const char* jobu, const char* jobvt, const blas_int* m,
              const blas_int* n, float* a, const blas_int* lda, float* s,
              float* u, const blas_int* ldu, float* vt, const blas_int* ldvt,
              float* work, const blas_int* lwork, blas_int* info );
 
void zgesvd_( const char* jobu, const char* jobvt, const blas_int* m,
              const blas_int* n, std::complex<double>* a, const blas_int* lda, double* s,
              std::complex<double>* u, const blas_int* ldu, std::complex<double>* vt, const
blas_int* ldvt, std::complex<double>* work, const blas_int* lwork, double* rwork, blas_int* info );
 
void cgesvd_( const char* jobu, const char* jobvt, const blas_int* m,
              const blas_int* n, std::complex<float>* a, const blas_int* lda, float* s,
              std::complex<float>* u, const blas_int* ldu, std::complex<float>* vt, const blas_int*
ldvt, std::complex<float>* work, const blas_int* lwork, float* rwork, blas_int* info );
 
void dgesdd_( const char* jobz, const blas_int* m, const blas_int* n, double* a,
              const blas_int* lda, double* s, double* u, const blas_int* ldu, double* vt, const
blas_int* ldvt, double* work, const blas_int* lwork, blas_int* iwork, blas_int* info ); void
zgesdd_( const char* jobz, const blas_int* m, const blas_int* n, std::complex<double>* a, const
blas_int* lda, double* s, std::complex<double>* u, const blas_int* ldu, std::complex<double>* vt,
const blas_int* ldvt, std::complex<double>* work, const blas_int* lwork, double* rwork, blas_int*
iwork, blas_int *info );
 
void ssyev_( const char* jobz, const char* uplo, const blas_int* n, float* a,
             const blas_int* lda, float* w, float* work, const blas_int* lwork,
             blas_int* info );
void dsyev_( const char* jobz, const char* uplo, const blas_int* n, double* a,
             const blas_int* lda, double* w, double* work, const blas_int* lwork,
             blas_int* info );
void cheev_( const char* jobz, const char* uplo, const blas_int* n, std::complex<float>* a,
             const blas_int* lda, float* w, std::complex<float>* work, const blas_int* lwork,
             const float* rwork, blas_int* info );
 
void zheev_( const char* jobz, const char* uplo, const blas_int* n, std::complex<double>* a,
             const blas_int* lda, double* w, std::complex<double>* work, const blas_int* lwork,
             const double* rwork, blas_int* info );
 
 
void zgeev_( const char* jobvl, const char* jobvr, const blas_int* n,  std::complex<double>* a,
             const blas_int* lda, std::complex<double>* w, std::complex<double> *vl, const blas_int
*ldvl, std::complex<double> *vr, const blas_int *ldvr, std::complex<double> *work, const blas_int*
lwork, double *rwork, blas_int* info ); void cgeev_( const char* jobvl, const char* jobvr, const
blas_int* n,  std::complex<float>* a, const blas_int* lda, std::complex<float>* w,
std::complex<float> *vl, const blas_int *ldvl, std::complex<float> *vr, const blas_int *ldvr,
std::complex<float> *work, const blas_int* lwork, float *rwork, blas_int* info ); void dgeev_( const
char* jobvl, const char* jobvr, const blas_int* n,  double* a, const blas_int* lda, double* wr,
double *wl, double *vl, const blas_int *ldvl, double *vr, const blas_int *ldvr, double *work, const
blas_int* lwork, blas_int* info ); void sgeev_( const char* jobvl, const char* jobvr, const
blas_int* n,  float* a, const blas_int* lda, float* wr, float *wl, float *vl, const blas_int *ldvl,
             float *vr, const blas_int *ldvr, float *work, const blas_int* lwork,
             blas_int* info );
 
void sstev_( const char* jobz, const blas_int* n, const float* d, const float* e, const float* z,
             const blas_int* ldaz, const float* work, blas_int* info );
void dstev_( const char* jobz, const blas_int* n, const double* d, const double* e, const double* z,
             const blas_int* ldaz, const double* work, blas_int* info );
 
void dgetrf_( const blas_int *m, const blas_int *n, const double *a,
              const blas_int *lda, const blas_int *ipiv, blas_int* info );
void sgetrf_( const blas_int *m, const blas_int *n, const float *a,
              const blas_int *lda, const blas_int *ipiv, blas_int* info );
 
void zgetrf_( const blas_int *m, const blas_int *n, const std::complex<double> *a,
              const blas_int *lda, const blas_int *ipiv, blas_int* info );
void cgetrf_( const blas_int *m, const blas_int *n, const std::complex<float> *a,
              const blas_int *lda, const blas_int *ipiv, blas_int* info );
 
void dgetri_( const blas_int *n, const double *a,  const blas_int *lda,
              const blas_int *ipiv, const double* work, const blas_int* lwork, blas_int* info );
void sgetri_( const blas_int *n, const float *a,  const blas_int *lda,
              const blas_int *ipiv, const float* work, const blas_int* lwork, blas_int* info );
void cgetri_( const blas_int *n, const std::complex<float> *a, const blas_int *lda,
              const blas_int *ipiv, const std::complex<float> *work, const blas_int *lwork, blas_int
*info ); void zgetri_( const blas_int *n, const std::complex<double> *a, const blas_int *lda, const
blas_int *ipiv, const std::complex<double> *work, const blas_int *lwork, blas_int *info ); void
sgelqf_( const blas_int* m, const blas_int* n, float* a, const blas_int* lda, float* tau, float*
work, const blas_int* lwork, blas_int* info );
 
void dgelqf_( const blas_int* m, const blas_int* n, double* a,
              const blas_int* lda, double* tau, double* work, const blas_int* lwork, blas_int* info
); void zgelqf_( const blas_int* m, const blas_int* n, std::complex<double>* a, const blas_int* lda,
std::complex<double>* tau, std::complex<double>* work, const blas_int* lwork, blas_int* info ); void
cgelqf_( const blas_int* m, const blas_int* n, std::complex<float>* a, const blas_int* lda,
std::complex<float>* tau, std::complex<float>* work, const blas_int* lwork, blas_int* info );
 
void dorglq_( const blas_int* m, const blas_int* n, const blas_int* k, double* a,
              const blas_int* lda, const double* tau, double* work, const blas_int* lwork, blas_int*
info ); void sorglq_( const blas_int* m, const blas_int* n, const blas_int* k, float* a, const
blas_int* lda, const float* tau, float* work, const blas_int* lwork, blas_int* info ); void zunglq_(
const blas_int* m, const blas_int* n, const blas_int* k, std::complex<double>* a, const blas_int*
lda, const std::complex<double>* tau, std::complex<double>* work, const blas_int* lwork, blas_int*
info ); void cunglq_( const blas_int* m, const blas_int* n, const blas_int* k, std::complex<float>*
a, const blas_int* lda, const std::complex<float>* tau, std::complex<float>* work, const blas_int*
lwork, blas_int* info ); void dgeqlf_( const blas_int* m, const blas_int* n, double* a, const
blas_int* lda, double* tau, double* work, const blas_int* lwork, blas_int* info ); void zgeqlf_(
const blas_int* m, const blas_int* n, std::complex<double>* a, const blas_int* lda,
std::complex<double>* tau, std::complex<double>* work, const blas_int* lwork, blas_int* info );
 
void dorgql_( const blas_int* m, const blas_int* n, const blas_int* k, double* a,
              const blas_int* lda, const double* tau, double* work, const blas_int* lwork, blas_int*
info ); void zungql_( const blas_int* m, const blas_int* n, const blas_int* k, std::complex<double>*
a, const blas_int* lda, const std::complex<double>* tau, std::complex<double>* work, const blas_int*
lwork, blas_int* info );
 
void dgeqrf_( const blas_int* m, const blas_int* n, double* a,
              const blas_int* lda, double* tau, double* work, const blas_int* lwork, blas_int* info
); void sgeqrf_( const blas_int* m, const blas_int* n, float* a, const blas_int* lda, float* tau,
float* work, const blas_int* lwork, blas_int* info ); void zgeqrf_( const blas_int* m, const
blas_int* n, std::complex<double>* a, const blas_int* lda, std::complex<double>* tau,
std::complex<double>* work, const blas_int* lwork, blas_int* info ); void cgeqrf_( const blas_int*
m, const blas_int* n, std::complex<float>* a, const blas_int* lda, std::complex<float>* tau,
std::complex<float>* work, const blas_int* lwork, blas_int* info ); void dorgqr_( const blas_int* m,
const blas_int* n, const blas_int* k, double* a, const blas_int* lda, const double* tau, double*
work, const blas_int* lwork, blas_int* info ); void sorgqr_( const blas_int* m, const blas_int* n,
const blas_int* k, float* a, const blas_int* lda, const float* tau, float* work, const blas_int*
lwork, blas_int* info );
 
void zungqr_( const blas_int* m, const blas_int* n, const blas_int* k, std::complex<double>* a,
              const blas_int* lda, const std::complex<double>* tau, std::complex<double>* work,
const blas_int* lwork, blas_int* info ); void cungqr_( const blas_int* m, const blas_int* n, const
blas_int* k, std::complex<float>* a, const blas_int* lda, const std::complex<float>* tau,
std::complex<float>* work, const blas_int* lwork, blas_int* info ); void dgerqf_( const blas_int* m,
const blas_int* n, double* a, const blas_int* lda, double* tau, double* work, const blas_int* lwork,
blas_int* info ); void zgerqf_( const blas_int* m, const blas_int* n, std::complex<double>* a, const
blas_int* lda, std::complex<double>* tau, std::complex<double>* work, const blas_int* lwork,
blas_int* info );
 
void dorgrq_( const blas_int* m, const blas_int* n, const blas_int* k, double* a,
              const blas_int* lda, const double* tau, double* work, const blas_int* lwork, blas_int*
info ); void sorgrq_( const blas_int* m, const blas_int* n, const blas_int* k, float* a, const
blas_int* lda, const float* tau, float* work, const blas_int* lwork, blas_int* info );
 
void zungrq_( const blas_int* m, const blas_int* n, const blas_int* k, std::complex<double>* a,
              const blas_int* lda, const std::complex<double>* tau, std::complex<double>* work,
const blas_int* lwork, blas_int* info );
 
void dstemr_( const char *jobz, const char *range, const blas_int *N, double *D, double *E,const
double *vl, const double *vu, const blas_int *il, const blas_int *iu, blas_int *M, double* W, double
*Z, const blas_int *ldz, const blas_int *nzc, blas_int *isuppz, bool *tryrac, double*work, blas_int*
lwork, blas_int*iwork, blas_int* liwork, blas_int* info);
*/
}  // extern C
 
// Wrappers for BLAS and LAPACK functions
inline void dgemm(const char *transa, const char *transb, const blas_int *m, const blas_int *n,
                  const blas_int *k, const double *alpha, const double *a, const blas_int *lda,
                  const double *b, const blas_int *ldb, const double *beta, double *c,
                  const blas_int *ldc) {
  dgemm_(transa, transb, m, n, k, alpha, a, lda, b, ldb, beta, c, ldc);
}
inline void sgemm(const char *transa, const char *transb, const blas_int *m, const blas_int *n,
                  const blas_int *k, const float *alpha, const float *a, const blas_int *lda,
                  const float *b, const blas_int *ldb, const float *beta, float *c,
                  const blas_int *ldc) {
  sgemm_(transa, transb, m, n, k, alpha, a, lda, b, ldb, beta, c, ldc);
}
 
inline void zgemm(const char *transa, const char *transb, const blas_int *m, const blas_int *n,
                  const blas_int *k, const std::complex<double> *alpha,
                  const std::complex<double> *a, const blas_int *lda, const std::complex<double> *b,
                  const blas_int *ldb, const std::complex<double> *beta, std::complex<double> *c,
                  const blas_int *ldc) {
  zgemm_(transa, transb, m, n, k, (const std::complex<double> *)alpha,
         (const std::complex<double> *)a, lda, (const std::complex<double> *)b, ldb,
         (const std::complex<double> *)beta, (std::complex<double> *)c, ldc);
}
 
inline void cgemm(const char *transa, const char *transb, const blas_int *m, const blas_int *n,
                  const blas_int *k, const std::complex<float> *alpha, const std::complex<float> *a,
                  const blas_int *lda, const std::complex<float> *b, const blas_int *ldb,
                  const std::complex<float> *beta, std::complex<float> *c, const blas_int *ldc) {
  cgemm_(transa, transb, m, n, k, (const std::complex<float> *)alpha,
         (const std::complex<float> *)a, lda, (const std::complex<float> *)b, ldb,
         (const std::complex<float> *)beta, (std::complex<float> *)c, ldc);
}
 
inline double dasum(const blas_int *n, const double *x, const blas_int *incx) {
  return dasum_(n, x, incx);
}
 
inline void dcopy(const blas_int &n, const double *x, const blas_int &incx, double *y,
                  const blas_int &incy) {
  dcopy_(&n, x, &incx, y, &incy);
}
inline void scopy(const blas_int &n, const float *x, const blas_int &incx, float *y,
                  const blas_int &incy) {
  scopy_(&n, x, &incx, y, &incy);
}
 
inline void daxpy(const blas_int *n, const double *alpha, const double *x, const blas_int *incx,
                  double *y, const blas_int *incy) {
  daxpy_(n, alpha, x, incx, y, incy);
}
inline void saxpy(const blas_int *n, const float *alpha, const float *x, const blas_int *incx,
                  float *y, const blas_int *incy) {
  saxpy_(n, alpha, x, incx, y, incy);
}
 
inline void zaxpy(const blas_int *n, const std::complex<double> *alpha,
                  const std::complex<double> *x, const blas_int *incx, std::complex<double> *y,
                  const blas_int *incy) {
  zaxpy_(n, alpha, x, incx, y, incy);
}
inline void caxpy(const blas_int *n, const std::complex<float> *alpha, const std::complex<float> *x,
                  const blas_int *incx, std::complex<float> *y, const blas_int *incy) {
  caxpy_(n, alpha, x, incx, y, incy);
}
 
inline double dnrm2(const blas_int *n, const double *x, const blas_int *incx) {
  return dnrm2_(n, x, incx);
}
 
inline double dznrm2(const blas_int *n, const std::complex<double> *x, const blas_int *incx) {
  return dznrm2_(n, x, incx);
}
 
inline float snrm2(const blas_int *n, const float *x, const blas_int *incx) {
  return snrm2_(n, x, incx);
}
 
inline float scnrm2(const blas_int *n, const std::complex<float> *x, const blas_int *incx) {
  return scnrm2_(n, x, incx);
}
 
inline void dscal(const blas_int *n, const double *a, double *x, const blas_int *incx) {
  dscal_(n, a, x, incx);
}
inline void sscal(const blas_int *n, const float *a, float *x, const blas_int *incx) {
  sscal_(n, a, x, incx);
}
inline void zscal(const blas_int *n, const std::complex<double> *a, std::complex<double> *x,
                  const blas_int *incx) {
  zscal_(n, a, x, incx);
}
inline void cscal(const blas_int *n, const std::complex<float> *a, std::complex<float> *x,
                  const blas_int *incx) {
  cscal_(n, a, x, incx);
}
inline void zdscal(const blas_int *n, const double *a, std::complex<double> *x,
                   const blas_int *incx) {
  zdscal_(n, a, x, incx);
}
/*
inline void dsyev( const char* jobz, const char* uplo, const blas_int* n, double* a,
             const blas_int* lda, double* w, double* work, const blas_int* lwork,
             blas_int* info )
{ dsyev_(  jobz,  uplo,  n,  a, lda, w,  work,  lwork, info ); }
inline void ssyev( const char* jobz, const char* uplo, const blas_int* n, float* a,
             const blas_int* lda, float* w, float* work, const blas_int* lwork,
             blas_int* info )
{ ssyev_(  jobz,  uplo,  n,  a, lda, w,  work,  lwork, info ); }
 
inline void cheev( const char* jobz, const char* uplo, const blas_int* n, std::complex<float>* a,
             const blas_int* lda, float* w, std::complex<float>* work, const blas_int* lwork,
             const float* rwork, blas_int* info )
{ cheev_(  jobz,  uplo,  n,  a, lda, w,  work,  lwork, rwork, info ); }
 
inline void zheev( const char* jobz, const char* uplo, const blas_int* n, std::complex<double>* a,
             const blas_int* lda, double* w, std::complex<double>* work, const blas_int* lwork,
             const double* rwork, blas_int* info )
{ zheev_(  jobz,  uplo,  n,  a, lda, w,  work,  lwork, rwork, info ); }
 
inline void zgeev( const char* jobvl, const char* jobvr, const blas_int* n, std::complex<double>* a,
    const blas_int* lda, std::complex<double>* w, std::complex<double> *vl, const blas_int *ldvl,
    std::complex<double> *vr, const blas_int *ldvr, std::complex<double> *work, const blas_int*
lwork, double *rwork, blas_int* info )
{
  zgeev_(jobvl, jobvr, n, a, lda, w, vl, ldvl, vr, ldvr, work, lwork, rwork, info);
}
 
inline void cgeev( const char* jobvl, const char* jobvr, const blas_int* n, std::complex<float>* a,
    const blas_int* lda, std::complex<float>* w, std::complex<float> *vl, const blas_int *ldvl,
    std::complex<float> *vr, const blas_int *ldvr, std::complex<float> *work, const blas_int* lwork,
    float *rwork, blas_int* info )
{
  cgeev_(jobvl, jobvr, n, a, lda, w, vl, ldvl, vr, ldvr, work, lwork, rwork, info);
}
 
 
 
inline void dgesvd( const char* jobu, const char* jobvt, const blas_int* m,
              const blas_int* n, double* a, const blas_int* lda, double* s,
              double* u, const blas_int* ldu, double* vt, const blas_int* ldvt,
              double* work, const blas_int* lwork, blas_int* info )
{
  dgesvd_( jobu, jobvt, m, n, a, lda, s, u, ldu, vt, ldvt, work, lwork, info );
}
 
inline void sgesvd( const char* jobu, const char* jobvt, const blas_int* m,
              const blas_int* n, float* a, const blas_int* lda, float* s,
              float* u, const blas_int* ldu, float* vt, const blas_int* ldvt,
              float* work, const blas_int* lwork, blas_int* info )
{
  sgesvd_( jobu, jobvt, m, n, a, lda, s, u, ldu, vt, ldvt, work, lwork, info );
}
 
inline void zgesvd( const char* jobu, const char* jobvt, const blas_int* m,
              const blas_int* n, std::complex<double>* a, const blas_int* lda, double* s,
              std::complex<double>* u, const blas_int* ldu, std::complex<double>* vt, const
blas_int* ldvt, std::complex<double>* work, const blas_int* lwork, double* rwork, blas_int* info )
{
  zgesvd_( jobu, jobvt, m, n, a, lda, s, u, ldu, vt, ldvt, work, lwork, rwork, info );
}
 
inline void cgesvd( const char* jobu, const char* jobvt, const blas_int* m,
              const blas_int* n, std::complex<float>* a, const blas_int* lda, float* s,
              std::complex<float>* u, const blas_int* ldu, std::complex<float>* vt, const blas_int*
ldvt, std::complex<float>* work, const blas_int* lwork, float* rwork, blas_int* info )
{
  cgesvd_( jobu, jobvt, m, n, a, lda, s, u, ldu, vt, ldvt, work, lwork, rwork, info );
}
 
 
inline void dgesdd( const char* jobz, const blas_int* m, const blas_int* n, double* a,
              const blas_int* lda, double* s, double* u, const blas_int* ldu, double* vt, const
blas_int* ldvt, double* work, const blas_int* lwork, blas_int* iwork, blas_int* info )
{
  dgesdd_(jobz, m, n, a, lda, s, u, ldu, vt, ldvt, work, lwork, iwork, info);
}
 
inline void zgesdd( const char* jobz, const blas_int* m, const blas_int* n,
              std::complex<double>* a, const blas_int* lda, double* s,
              std::complex<double>* u, const blas_int* ldu,
              std::complex<double>* vt, const blas_int* ldvt,
              std::complex<double>* work, const blas_int* lwork,
              double* rwork, blas_int* iwork, blas_int *info )
{
  zgesdd_(jobz, m, n, a, lda, s, u, ldu, vt, ldvt, work, lwork, rwork, iwork, info);
}
*/
 
inline void sgemv(const char *trans, const blas_int *m, const blas_int *n, const float *alpha,
                  const float *a, const blas_int *lda, const float *x, const blas_int *incx,
                  const float *beta, const float *y, const blas_int *incy) {
  sgemv_(trans, m, n, alpha, a, lda, x, incx, beta, y, incy);
}
inline void dgemv(const char *trans, const blas_int *m, const blas_int *n, const double *alpha,
                  const double *a, const blas_int *lda, const double *x, const blas_int *incx,
                  const double *beta, const double *y, const blas_int *incy) {
  dgemv_(trans, m, n, alpha, a, lda, x, incx, beta, y, incy);
}
inline void zgemv(const char *trans, const blas_int *m, const blas_int *n,
                  const std::complex<double> *alpha, const std::complex<double> *a,
                  const blas_int *lda, const std::complex<double> *x, const blas_int *incx,
                  const std::complex<double> *beta, const std::complex<double> *y,
                  const blas_int *incy) {
  zgemv_(trans, m, n, alpha, a, lda, x, incx, beta, y, incy);
}
inline void cgemv(const char *trans, const blas_int *m, const blas_int *n,
                  const std::complex<float> *alpha, const std::complex<float> *a,
                  const blas_int *lda, const std::complex<float> *x, const blas_int *incx,
                  const std::complex<float> *beta, const std::complex<float> *y,
                  const blas_int *incy) {
  cgemv_(trans, m, n, alpha, a, lda, x, incx, beta, y, incy);
}
 
inline double ddot(const blas_int *n, const double *x, const blas_int *incx, const double *y,
                   const blas_int *incy) {
  return ddot_(n, x, incx, y, incy);
}
inline float sdot(const blas_int *n, const float *x, const blas_int *incx, const float *y,
                  const blas_int *incy) {
  return sdot_(n, x, incx, y, incy);
}
 
inline void zdotc(std::complex<double> *res, const blas_int *n, const std::complex<double> *x,
                  const blas_int *incx, const std::complex<double> *y, const blas_int *incy) {
  #ifndef FORTRAN_COMPLEX_FUNCTIONS_RETURN_VOID
  *res = zdotc_(n, x, incx, y, incy);
  #else
  zdotc_(res, n, x, incx, y, incy);
  #endif
}
inline void zdotu(std::complex<double> *res, const blas_int *n, const std::complex<double> *x,
                  const blas_int *incx, const std::complex<double> *y, const blas_int *incy) {
  #ifndef FORTRAN_COMPLEX_FUNCTIONS_RETURN_VOID
  *res = zdotu_(n, x, incx, y, incy);
  #else
  zdotu_(res, n, x, incx, y, incy);
  #endif
}
inline void cdotc(std::complex<float> *res, const blas_int *n, const std::complex<float> *x,
                  const blas_int *incx, const std::complex<float> *y, const blas_int *incy) {
  #ifndef FORTRAN_COMPLEX_FUNCTIONS_RETURN_VOID
  *res = cdotc_(n, x, incx, y, incy);
  #else
  cdotc_(res, n, x, incx, y, incy);
  #endif
}
inline void cdotu(std::complex<float> *res, const blas_int *n, const std::complex<float> *x,
                  const blas_int *incx, const std::complex<float> *y, const blas_int *incy) {
  #ifndef FORTRAN_COMPLEX_FUNCTIONS_RETURN_VOID
  *res = cdotu_(n, x, incx, y, incy);
  #else
  cdotu_(res, n, x, incx, y, incy);
  #endif
}
 
  /*
  inline void dstev( const char* jobz, const blas_int* n, const double* d, const double* e, const
  double* z, const blas_int* ldaz, const double* work, blas_int* info )
  {
    dstev_( jobz, n, d, e, z, ldaz, work, info );
  }
 
  inline void sstev( const char* jobz, const blas_int* n, const float* d, const float* e, const
  float* z, const blas_int* ldaz, const float* work, blas_int* info )
  {
    sstev_( jobz, n, d, e, z, ldaz, work, info );
  }
 
  inline void dgetrf( const blas_int *m, const blas_int *n, const double *a,  const blas_int *lda,
  const blas_int *ipiv, blas_int* info )
  {
    dgetrf_( m, n, a, lda, ipiv, info );
  }
  inline void sgetrf( const blas_int *m, const blas_int *n, const float *a,  const blas_int *lda,
  const blas_int *ipiv, blas_int* info )
  {
    sgetrf_( m, n, a, lda, ipiv, info );
  }
 
  inline void zgetrf( const blas_int *m, const blas_int *n, const std::complex<double> *a,  const
  blas_int *lda, const blas_int *ipiv, blas_int* info )
  {
    zgetrf_(m, n, a, lda, ipiv, info);
  }
  inline void cgetrf( const blas_int *m, const blas_int *n, const std::complex<float> *a,  const
  blas_int *lda, const blas_int *ipiv, blas_int* info )
  {
    cgetrf_(m, n, a, lda, ipiv, info);
  }
 
  inline void dgetri( const blas_int *n, const double *a,  const blas_int *lda, const blas_int
  *ipiv, const double* work, const blas_int* lwork, blas_int* info )
  {
    dgetri_(n, a, lda, ipiv, work, lwork, info);
  }
  inline void sgetri( const blas_int *n, const float *a,  const blas_int *lda, const blas_int *ipiv,
  const float* work, const blas_int* lwork, blas_int* info )
  {
    sgetri_(n, a, lda, ipiv, work, lwork, info);
  }
  inline void zgetri( const blas_int *n, const std::complex<double> *a, const blas_int *lda, const
  blas_int *ipiv, const std::complex<double> *work, const blas_int *lwork, blas_int *info )
  {
    zgetri_(n, a, lda, ipiv, work, lwork, info);
  }
  inline void cgetri( const blas_int *n, const std::complex<float> *a, const blas_int *lda, const
  blas_int *ipiv, const std::complex<float> *work, const blas_int *lwork, blas_int *info )
  {
    cgetri_(n, a, lda, ipiv, work, lwork, info);
  }
 
  //inline void dgeqpf( const blas_int* m, const blas_int* n, double* a, const blas_int* lda, int*
  jpvt, double* tau, double* work, blas_int* info ){
  //
  //dgeqpf_(m, n, a, lda, jpvt, tau, work, info);
  //
  //}
 
 
  inline void dgeqp3( const blas_int* m, const blas_int* n, double* a, const blas_int* lda, int*
  jpvt, double* tau, double* work, int* lwork, blas_int* info ){
 
    dgeqp3_(m, n, a, lda, jpvt, tau, work, lwork,info);
 
  }
 
  inline void dgelqf( const blas_int* m, const blas_int* n, double* a,
                      const blas_int* lda, double* tau, double* work, const blas_int* lwork,
  blas_int* info )
  {
    dgelqf_(m, n, a, lda, tau, work, lwork, info );
  }
  inline void sgelqf( const blas_int* m, const blas_int* n, float* a,
                      const blas_int* lda, float* tau, float* work, const blas_int* lwork, blas_int*
  info )
  {
    sgelqf_(m, n, a, lda, tau, work, lwork, info );
  }
  inline void dorglq( const blas_int* m, const blas_int* n, const blas_int* k, double* a,
                      const blas_int* lda, const double* tau, double* work, const blas_int* lwork,
  blas_int* info )
  {
    dorglq_(m, n, k, a, lda, tau, work, lwork, info );
  }
  inline void sorglq( const blas_int* m, const blas_int* n, const blas_int* k, float* a,
                      const blas_int* lda, const float* tau, float* work, const blas_int* lwork,
  blas_int* info )
  {
    sorglq_(m, n, k, a, lda, tau, work, lwork, info );
  }
 
 
 
  inline void dgeqlf( const blas_int* m, const blas_int* n, double* a,
                      const blas_int* lda, double* tau, double* work, const blas_int* lwork,
  blas_int* info )
  {
    dgeqlf_(m, n, a, lda, tau, work, lwork, info );
  }
 
  inline void dorgql( const blas_int* m, const blas_int* n, const blas_int* k, double* a,
                      const blas_int* lda, const double* tau, double* work, const blas_int* lwork,
  blas_int* info )
  {
    dorgql_(m, n, k, a, lda, tau, work, lwork, info );
  }
 
  inline void dgeqrf( const blas_int* m, const blas_int* n, double* a,
                      const blas_int* lda, double* tau, double* work, const blas_int* lwork,
  blas_int* info )
  {
    dgeqrf_(m, n, a, lda, tau, work, lwork, info );
  }
  inline void sgeqrf( const blas_int* m, const blas_int* n, float* a,
                      const blas_int* lda, float* tau, float* work, const blas_int* lwork, blas_int*
  info )
  {
    sgeqrf_(m, n, a, lda, tau, work, lwork, info );
  }
  inline void dorgqr( const blas_int* m, const blas_int* n, const blas_int* k, double* a,
                      const blas_int* lda, const double* tau, double* work, const blas_int* lwork,
  blas_int* info )
  {
    dorgqr_(m, n, k, a, lda, tau, work, lwork, info );
  }
  inline void sorgqr( const blas_int* m, const blas_int* n, const blas_int* k, float* a,
                      const blas_int* lda, const float* tau, float* work, const blas_int* lwork,
  blas_int* info )
  {
    sorgqr_(m, n, k, a, lda, tau, work, lwork, info );
  }
 
  inline void dgerqf( const blas_int* m, const blas_int* n, double* a,
                      const blas_int* lda, double* tau, double* work, const blas_int* lwork,
  blas_int* info )
  {
    dgerqf_(m, n, a, lda, tau, work, lwork, info );
  }
 
  inline void dorgrq(const blas_int* m, const blas_int* n, const blas_int* k, double* a,
                     const blas_int* lda, const double* tau, double* work, const blas_int* lwork,
  blas_int* info )
  {
    dorgrq_(m, n, k, a, lda, tau, work, lwork, info );
  }
  inline void sorgrq(const blas_int* m, const blas_int* n, const blas_int* k, float* a,
                     const blas_int* lda, const float* tau, float* work, const blas_int* lwork,
  blas_int* info )
  {
    sorgrq_(m, n, k, a, lda, tau, work, lwork, info );
  }
  inline void zgeqrf( const blas_int* m, const blas_int* n, std::complex<double>* a,
                      const blas_int* lda, std::complex<double>* tau, std::complex<double>* work,
  const blas_int* lwork, blas_int* info )
  {
    zgeqrf_(m, n, a, lda, tau, work, lwork, info );
  }
  inline void cgeqrf( const blas_int* m, const blas_int* n, std::complex<float>* a,
                      const blas_int* lda, std::complex<float>* tau, std::complex<float>* work,
  const blas_int* lwork, blas_int* info )
  {
    cgeqrf_(m, n, a, lda, tau, work, lwork, info );
  }
  inline void zungqr( const blas_int* m, const blas_int* n, const blas_int* k, std::complex<double>*
  a, const blas_int* lda, const std::complex<double>* tau, std::complex<double>* work, const
  blas_int* lwork, blas_int* info )
  {
    zungqr_(m, n, k, a, lda, tau, work, lwork, info );
  }
  inline void cungqr( const blas_int* m, const blas_int* n, const blas_int* k, std::complex<float>*
  a, const blas_int* lda, const std::complex<float>* tau, std::complex<float>* work, const blas_int*
  lwork, blas_int* info )
  {
    cungqr_(m, n, k, a, lda, tau, work, lwork, info );
  }
 
  inline void zgerqf( const blas_int* m, const blas_int* n, std::complex<double>* a,
                      const blas_int* lda, std::complex<double>* tau, std::complex<double>* work,
  const blas_int* lwork, blas_int* info )
  {
    zgerqf_(m, n, a, lda, tau, work, lwork, info );
  }
 
  inline void zungrq( const blas_int* m, const blas_int* n, const blas_int* k, std::complex<double>*
  a, const blas_int* lda, const std::complex<double>* tau, std::complex<double>* work, const
  blas_int* lwork, blas_int* info )
  {
    zungrq_(m, n, k, a, lda, tau, work, lwork, info );
  }
 
  inline void zgeqp3( const blas_int* m, const blas_int* n, std::complex<double>* a, const blas_int*
  lda, int* jpvt, std::complex<double>* tau, std::complex<double>* work, int* lwork, double* rwork,
  blas_int* info ){ zgeqp3_(m, n, a, lda, jpvt, tau, work, lwork,rwork, info);
  }
 
  //inline void zgeqpf( const blas_int* m, const blas_int* n, std::complex<double>* a, const
  blas_int* lda,
  //              int* jpvt, std::complex<double>* tau, std::complex<double>* work,
  //              double* rwork, blas_int* info ){
  //  zgeqpf_(m, n, a, lda, jpvt, tau, work, rwork, info);
  //}
 
  inline void zgelqf( const blas_int* m, const blas_int* n, std::complex<double>* a,
                      const blas_int* lda, std::complex<double>* tau, std::complex<double>* work,
  const blas_int* lwork, blas_int* info )
  {
    zgelqf_(m, n, a, lda, tau, work, lwork, info );
  }
  inline void cgelqf( const blas_int* m, const blas_int* n, std::complex<float>* a,
                      const blas_int* lda, std::complex<float>* tau, std::complex<float>* work,
  const blas_int* lwork, blas_int* info )
  {
    cgelqf_(m, n, a, lda, tau, work, lwork, info );
  }
 
  inline void zunglq( const blas_int* m, const blas_int* n, const blas_int* k, std::complex<double>*
  a, const blas_int* lda, const std::complex<double>* tau, std::complex<double>* work, const
  blas_int* lwork, blas_int* info )
  {
    zunglq_(m, n, k, a, lda, tau, work, lwork, info );
  }
  inline void cunglq( const blas_int* m, const blas_int* n, const blas_int* k, std::complex<float>*
  a, const blas_int* lda, const std::complex<float>* tau, std::complex<float>* work, const blas_int*
  lwork, blas_int* info )
  {
    cunglq_(m, n, k, a, lda, tau, work, lwork, info );
  }
  inline void zgeqlf( const blas_int* m, const blas_int* n, std::complex<double>* a,
                      const blas_int* lda, std::complex<double>* tau, std::complex<double>* work,
  const blas_int* lwork, blas_int* info )
  {
    zgeqlf_(m, n, a, lda, tau, work, lwork, info );
  }
 
  inline void zungql( const blas_int* m, const blas_int* n, const blas_int* k, std::complex<double>*
  a, const blas_int* lda, const std::complex<double>* tau, std::complex<double>* work, const
  blas_int* lwork, blas_int* info )
  {
    zungql_(m, n, k, a, lda, tau, work, lwork, info );
  }
  */
#endif  // UNI_MKL
 
#endif