blas_lapack.h

Go to the documentation of this file.

#ifndef BLAS_LAPACK_H
#define BLAS_LAPACK_H


namespace gurls
{
    enum CBLAS_ORDER     {CblasRowMajor=9, CblasColMajor=10};

    enum CBLAS_TRANSPOSE {CblasNoTrans=0, CblasTrans=1, CblasConjTrans=2};

    enum CBLAS_UPLO      {CblasUpper=3, CblasLower=4};

    enum CBLAS_DIAG      {CblasNonUnit=5, CblasUnit=6};

    enum CBLAS_SIDE      {CblasLeft=7, CblasRight=8};

    class BlasUtils
    {
    public:

        static char charValue(int value)
        {
            static char chars[] = "NTCULNULR";

            return chars[value];
        }

    };

}

extern "C"
{


// ------ BLAS


float sdot_(int *n, float *x, int *incx, float *y, int *incy);

void sgemm_(char *transa, char *transb, int *m, int *n, int *k, float *alpha, float *a, int *lda, float *b, int *ldb, float *beta, float *c, int *ldc);

void sgemv_(char *trans, int *m, int *n, float *alpha, float *a, int *lda, float *x, int *incx, float *beta, float *y, int *incy);

void scopy_(int *n, float *x, int *incx, float *y, int *incy);

void saxpy_(int *n, float *alpha, float *x, int *incx, float *y, int *incy);

float snrm2_(int *n, float *x, int *incx);

void sscal_(int *n, float *a, float *x, int *incx);

void strsm_(char *side, char *uplo, char *transa, char *diag, int *m, int *n, float *alpha, float *a, int *lda, float *b, int *ldb);

void sswap_(int *n, float *sx, int *incx, float *sy, int *incy);

double ddot_(int *n, double *x, int *incx, double *y, int *incy);

void dgemm_(char *transa, char *transb, int *m, int *n, int *k, double *alpha, double *a, int *lda, double *b, int *ldb, double *beta, double *c, int *ldc);

void dgemv_(char *trans, int *m, int *n, double *alpha, double *a, int *lda, double *x, int *incx, double *beta, double *y, int *incy);

void dcopy_(int *n, double *x, int *incx, double *y, int *incy);

void daxpy_(int *n, double *alpha, double *x, int *incx, double *y, int *incy);

double dnrm2_(int *n, double *x, int *incx);

void dscal_(int *n, double *a, double *x, int *incx);

void dtrsm_(char *side, char *uplo, char *transa, char *diag, int *m, int *n, double *alpha, double *a, int *lda, double *b, int *ldb);

void dswap_(int *n, double *sx, int *incx, double *sy, int *incy);



// ------ LAPACK



int sgetrf_(int *m, int *n, float *a, int *lda, int *ipiv, int *info);

int sgesvd_(char *jobu, char *jobvt, int *m, int *n, float *a, int *lda, float *s, float *u, int *ldu, float *vt, int *ldvt, float *work, int *lwork, int *info);

int dgesvd_(char *jobu, char *jobvt, int *m, int *n, double *a, int *lda, double *s, double *u, int *ldu, double *vt, int *ldvt, double *work, int *lwork, int *info);

int sgeev_(char *jobvl, char *jobvr, int *n, float *a, int *lda, float *wr, float *wi, float *vl, int *ldvl, float *vr, int *ldvr, float *work, int *lwork, int *info);

int 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);

int sgetri_(int *n, float *a, int* lda, int *ipiv, float*work, int *lwork, int *info);

int spotrf_(char *UPLO, int *n, float *a, int *lda , int *info);

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

int sgelss_( int *m, int *n, int* nrhs, float *a, int *lda, float* b, int *ldb, float *s, float *rcond, int *rank, float *work, int *lwork, int *info);

int dgelss_( int *m, int *n, int* nrhs, double *a, int *lda, double* b, int *ldb, double *s, double *rcond, int *rank, double *work, int *lwork, int *info);

int ssyev_( char* jobz, char* uplo, int* n, float* a, int* lda, float* w, float* work, int* lwork, int* info );

int dsyev_( char* jobz, char* uplo, int* n, double* a, int* lda, double* w, double* work, int* lwork, int* info );

void sgeqp3_( int *m, int *n, float *A, int *lda, int *jpvt, float *tau, float *work, int *lwork, int *info);

void dgeqp3_( int *m, int *n, double *A, int *lda, int *jpvt, double *tau, double *work, int *lwork, int *info);

void sorgqr_(int *m, int *n, int *k, float *a, int *lda, float *tau, float *work, int *lwork, int *info);

void dorgqr_(int *m, int *n, int *k, double *a, int *lda, double *tau, double *work, int *lwork, int *info);

}

#include"gurls++/blas_lapack.hpp"

#endif //BLAS_LAPACK_H