Blas.h

Go to the documentation of this file.

#ifndef BLAS_H
#define BLAS_H
#include <complex>
 
extern "C" {
 
/* Level 1 BLAS */
double dscal_(const int& n, const double&, double* x, const int&);
double zscal_(const int& n, const std::complex<double>&,
              std::complex<double>* x, const int&);
 
void zcopy_(const int& n, const std::complex<double>* x, const int& incx,
            std::complex<double>* y, const int& incy);
void dcopy_(const int& n, const double*, const int&, double*, const int&);
 
/* Level 2 BLAS */
 
void dgemv_(const char& trans, const int& nr, const int& nc,
            const double& alpha, const double* amat, const int& lda,
            const double* bv, const int& incx, const double& beta, double* cv,
            const int& incy);
 
void zgemv_(const char& trans, const int& nr, const int& nc,
            const std::complex<double>& alpha, const std::complex<double>* amat,
            const int& lda, const std::complex<double>* bv, const int& incx,
            const std::complex<double>& beta, std::complex<double>* cv,
            const int& incy);
 
void daxpy_(const int& n, const double& da, const double* dx, const int& incx,
            double* dy, const int& incy);
 
/* Level 3 BLAS */
void dgemm_(const char&, const char&, const int&, const int&, const int&,
            const double&, const double*, const int&, const double*, const int&,
            const double&, double*, const int&);
void zgemm_(const char&, const char&, const int&, const int&, const int&,
            const std::complex<double>&, const std::complex<double>*,
            const int&, const std::complex<double>*, const int&,
            const std::complex<double>&, std::complex<double>*, const int&);
 
// Lapack forward declarations
void dgeev_(char* JOBVL, char* JOBVR, int* N, double* A, int* LDA, double* WR,
            double* WI, double* VL, int* LDVL, double* VR, int* LDVR,
            double* WORK, int* LWORK, int* INFO);
void zgeev_(char* JOBVL, char* JOBVR, int* N, std::complex<double>* A, int* LDA,
            std::complex<double>* W, std::complex<double>* VL, int* LDVL,
            std::complex<double>* VR, int* LDVR, std::complex<double>* WORK,
            int* LWORK, double* RWORK, int* INFO);
void dgetrf_(int* M, int* N, double* A, int* LDA, int* IPIV, int* INFO);
void zgetrf_(int* M, int* N, std::complex<double>* A, int* LDA, int* IPIV,
             int* INFO);
double dlange_(char* NORM, int* M, int* N, double* A, int* LDA, double* work);
double zlange_(char* NORM, int* M, int* N, std::complex<double>* A, int* LDA,
               double* work);
 
void dsyev_(char* JOBZ, char* UPLO, int* N, double* A, int* LDA, double* W,
            double* WORK, int* LWORK, int* INFO);
void zheev_(char* JOBZ, char* UPLO, int* N, std::complex<double>* A, int* LDA,
            double* W, std::complex<double>* WORK, int* LWORK, double* RWORK,
            int* INFO);
}
 
#endif /* BLAS_H */