Changeset 18 for include/config-blas.h

Timestamp:

04/28/07 20:58:32 (2 years ago)

Author:

pernet

Message:

Introduction of the field Modular<float> and the corresponding specializations using the float BLAS.

Files:

• include/config-blas.h (modified) (8 diffs)

include/config-blas.h

@@ -1 +1 @@
 /* -*- mode: C++; tab-width: 8; indent-tabs-mode: t; c-basic-offset: 8 -*- */
-/* linbox/algorithms/lifting-container-base.h
+/* config-blas.h
  * Copyright (C) 2005  Pascal Giorgi
- *
+ *               2007  Clement Pernet
  * Written by Pascal Giorgi <pgiorgi@uwaterloo.ca>
  *
@@ -51 +51 @@
         // level 2 routines
         void dgemv_ (const char*, const int*, const int*, const double*, const double*, const int*, const double*, const int*, const double*, double*, const int*);
+        void sgemv_ (const char*, const int*, const int*, const float*, const float*, const int*, const float*, const int*, const float*, float*, const int*);
         void dger_  (const int*, const int*, const double*, const double*, const int*, const double*, const int*, double*, const int*);
 
         // level 3 routines
         void dtrsm_ (const char*, const char*, const char*, const char*, const int*, const int*, const double*, const double*, const int*, double*, const int*);
+        void strsm_ (const char*, const char*, const char*, const char*, const int*, const int*, const float*, const float*, const int*, float*, const int*);
         void dtrmm_ (const char*, const char*, const char*, const char*, const int*, const int*, const double*, const double*, const int*, double*, const int*);
+        void strmm_ (const char*, const char*, const char*, const char*, const int*, const int*, const float*, const float*, const int*, float*, const int*);
+        void sgemm_ (const char*, const char*, const int*, const int*, const int*, const float*, const float*, const int*, const float*, const int*, const float*, float*, const int*);
         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*);
 }
@@ -128 +132 @@
                         dgemv_ ( EXT_BLAS_TRANSPOSE(TransA), &M, &N, &alpha, A, &lda, X, &incX, &beta, Y, &incY);
         }
+        void cblas_sgemv(const enum CBLAS_ORDER Order, const enum CBLAS_TRANSPOSE TransA, const int M, const int N, const float alpha, 
+                         const float *A, const int lda, const float *X, const int incX, const float beta, float *Y, const int incY)
+        {
+                if (Order == CblasRowMajor)
+                        sgemv_ ( EXT_BLAS_TRANSPOSE_tr(TransA), &N, &M, &alpha, A, &lda, X, &incX, &beta, Y, &incY);
+                else
+                        sgemv_ ( EXT_BLAS_TRANSPOSE(TransA), &M, &N, &alpha, A, &lda, X, &incX, &beta, Y, &incY);
+        }
   
         void cblas_dger(const enum CBLAS_ORDER Order, const int M, const int N, const double alpha, const double *X, const int incX,
@@ -151 +163 @@
                         dtrsm_ ( EXT_BLAS_SIDE(Side), EXT_BLAS_UPLO(Uplo), EXT_BLAS_TRANSPOSE(TransA), EXT_BLAS_DIAG(Diag), &M, &N, &alpha, A, &lda, B, &ldb);
         }
+        void cblas_strsm(const enum CBLAS_ORDER Order, const enum CBLAS_SIDE Side, const enum CBLAS_UPLO Uplo, const enum CBLAS_TRANSPOSE TransA,
+                         const enum CBLAS_DIAG Diag, const int M, const int N, const float alpha, const float *A, const int lda,
+                         float *B, const int ldb)
+        {  
+                if (Order == CblasRowMajor) 
+                        strsm_ ( EXT_BLAS_SIDE_tr(Side), EXT_BLAS_UPLO_tr(Uplo), EXT_BLAS_TRANSPOSE(TransA), EXT_BLAS_DIAG(Diag), &N, &M, &alpha, A, &lda, B, &ldb); 
+                else
+                        strsm_ ( EXT_BLAS_SIDE(Side), EXT_BLAS_UPLO(Uplo), EXT_BLAS_TRANSPOSE(TransA), EXT_BLAS_DIAG(Diag), &M, &N, &alpha, A, &lda, B, &ldb);
+        }
   
         void cblas_dtrmm(const enum CBLAS_ORDER Order, const enum CBLAS_SIDE Side, const enum CBLAS_UPLO Uplo, const enum CBLAS_TRANSPOSE TransA,
@@ -160 +181 @@
                 else
                         dtrmm_ ( EXT_BLAS_SIDE(Side), EXT_BLAS_UPLO(Uplo), EXT_BLAS_TRANSPOSE(TransA), EXT_BLAS_DIAG(Diag), &M, &N, &alpha, A, &lda, B, &ldb);
+        }
+        void cblas_strmm(const enum CBLAS_ORDER Order, const enum CBLAS_SIDE Side, const enum CBLAS_UPLO Uplo, const enum CBLAS_TRANSPOSE TransA,
+                         const enum CBLAS_DIAG Diag, const int M, const int N, const float alpha, const float *A, const int lda,
+                         float *B, const int ldb)
+        {  
+                if (Order == CblasRowMajor)
+                        strmm_ ( EXT_BLAS_SIDE_tr(Side), EXT_BLAS_UPLO_tr(Uplo), EXT_BLAS_TRANSPOSE(TransA), EXT_BLAS_DIAG(Diag), &N, &M, &alpha, A, &lda, B, &ldb);
+                else
+                        strmm_ ( EXT_BLAS_SIDE(Side), EXT_BLAS_UPLO(Uplo), EXT_BLAS_TRANSPOSE(TransA), EXT_BLAS_DIAG(Diag), &M, &N, &alpha, A, &lda, B, &ldb);
         }
   
@@ -170 +200 @@
                 else
                         dgemm_ ( EXT_BLAS_TRANSPOSE(TransA), EXT_BLAS_TRANSPOSE(TransB), &M, &N, &K, &alpha, A, &lda, B, &ldb, &beta, C, &ldc);
+        }
+        void cblas_sgemm(const enum CBLAS_ORDER Order, const enum CBLAS_TRANSPOSE TransA, const enum CBLAS_TRANSPOSE TransB, const int M, const int N,
+                         const int K, const float alpha, const float *A, const int lda, const float *B, const int ldb,
+                         const float beta, float *C, const int ldc) 
+        {   
+                if (Order == CblasRowMajor)
+                        sgemm_ ( EXT_BLAS_TRANSPOSE(TransB), EXT_BLAS_TRANSPOSE(TransA), &N, &M, &K, &alpha, B, &ldb, A, &lda, &beta, C, &ldc);
+                else
+                        sgemm_ ( EXT_BLAS_TRANSPOSE(TransA), EXT_BLAS_TRANSPOSE(TransB), &M, &N, &K, &alpha, A, &lda, B, &ldb, &beta, C, &ldc);
         }
 
@@ -226 +265 @@
         void cblas_dgemv(const enum CBLAS_ORDER Order, const enum CBLAS_TRANSPOSE TransA, const int M, const int N, const double alpha, 
                          const double *A, const int lda, const double *X, const int incX, const double beta, double *Y, const int incY);
+
+        void cblas_sgemv(const enum CBLAS_ORDER Order, const enum CBLAS_TRANSPOSE TransA, const int M, const int N, const float alpha, 
+                         const float *A, const int lda, const float *X, const int incX, const float beta, float *Y, const int incY);
   
         void cblas_dger(const enum CBLAS_ORDER Order, const int M, const int N, const double alpha, const double *X, const int incX,
@@ -237 +279 @@
                          double *B, const int ldb);
   
+        void cblas_strsm(const enum CBLAS_ORDER Order, const enum CBLAS_SIDE Side, const enum CBLAS_UPLO Uplo, const enum CBLAS_TRANSPOSE TransA,
+                         const enum CBLAS_DIAG Diag, const int M, const int N, const float alpha, const float *A, const int lda,
+                         float *B, const int ldb);
+  
         void cblas_dtrmm(const enum CBLAS_ORDER Order, const enum CBLAS_SIDE Side, const enum CBLAS_UPLO Uplo, const enum CBLAS_TRANSPOSE TransA,
                          const enum CBLAS_DIAG Diag, const int M, const int N, const double alpha, const double *A, const int lda,
                          double *B, const int ldb);
   
+        void cblas_strmm(const enum CBLAS_ORDER Order, const enum CBLAS_SIDE Side, const enum CBLAS_UPLO Uplo, const enum CBLAS_TRANSPOSE TransA,
+                         const enum CBLAS_DIAG Diag, const int M, const int N, const float alpha, const float *A, const int lda,
+                         float *B, const int ldb);
+  
         void cblas_dgemm(const enum CBLAS_ORDER Order, const enum CBLAS_TRANSPOSE TransA, const enum CBLAS_TRANSPOSE TransB, const int M, const int N,
                          const int K, const double alpha, const double *A, const int lda, const double *B, const int ldb,
                          const double beta, double *C, const int ldc) ;
+        void cblas_sgemm(const enum CBLAS_ORDER Order, const enum CBLAS_TRANSPOSE TransA, const enum CBLAS_TRANSPOSE TransB, const int M, const int N,
+                         const int K, const float alpha, const float *A, const int lda, const float *B, const int ldb,
+                         const float beta, float *C, const int ldc) ;
 
         // LAPACK routines