blas.h

#ifndef MPQC_BLAS_H
#define MPQC_BLAS_H
 
#include <math/scmat/f77sym.h>
#include <mpqc_config.h>
#include <stdint.h>
 
#if defined(F77_INTEGER_WIDTH) && F77_INTEGER_WIDTH == 8
   typedef int64_t f77_integer_t;
#elif defined(F77_INTEGER_WIDTH) && F77_INTEGER_WIDTH == 4
   typedef int32_t f77_integer_t;
#else
# error "unknown F77_INTEGER_WIDTH"
#endif
 
typedef f77_integer_t blasint;
 
#ifdef __cplusplus
extern "C" {
#endif // __cplusplus
 
extern void F77_DGEMM(const char* transa, const char* transb, const blasint* nrows_opA,
                      const blasint* ncols_opB, const blasint* nrows_opB,
                      const double* alpha, const double* A, const blasint* lda,
                      const double* B, const blasint* ldb,
                      const double* beta,
                      double* C, const blasint* ldc);
 
extern void F77_DGEMV(const char* trans, const blasint* m, const blasint* n, const double* alpha,
                      const double* A, const blasint* lda, const double* X, const blasint* incX,
                      const double* beta, double* Y, const blasint* incY);
 
extern void F77_DAXPY(const blasint* n, const double* da, const double* dx,
                      const blasint* incx, double* dy, const blasint* incy);
 
extern double F77_DDOT(const blasint* n, const double* dx, const blasint* incx,
                       const double* dy, const blasint* incy);
 
extern void F77_DCOPY(const blasint *n, const double *dx, const blasint *incx, double *dy, const blasint *incy);
extern double F77_DNRM2(const blasint *n, const double *dx, const blasint *incx);
extern void F77_DSCAL(const blasint *n, const double *da, double *dx, const blasint *incx);
 
extern void F77_DSPMV(const char* uplo, const blasint* n, const double* alpha,
                      const double* A, const double* X, const blasint* incx,
                      const double* beta, double* Y,
                      const blasint* incy);
 
#ifdef __cplusplus
}
#endif // __cplusplus
 
#ifdef __cplusplus
namespace sc {
 
  //
  // taken from psi3/src/lib/libqt/blas_intfc.c
  //
 
  void C_DGEMM(char transa, char transb, int m, int n, int k, double alpha,
               const double *A, int nca, const double *B, int ncb, double beta,
               double *C, int ncc);
 
}
#endif // __cplusplus
 
#endif // MPQC_BLAS_H