Added device ger(c), geru

Author: WestonJB

CMakeLists.txt

@@ -358,6 +358,7 @@ add_library(
     src/device_error.cc
     src/device_gemm.cc
     src/device_gemv.cc
+    src/device_ger.cc
     src/device_hemm.cc
     src/device_hemv.cc
     src/device_her2k.cc

include/blas/device_blas.hh

@@ -580,6 +580,80 @@ void trsv(
     std::complex<double>*       x, int64_t incx,
     blas::Queue& queue );
 
+//------------------------------------------------------------------------------
+void ger(
+    blas::Layout layout,
+    int64_t m, int64_t n,
+    float alpha,
+    float const* x, int64_t incx,
+    float const* y, int64_t incy,
+    float*       A, int64_t lda,
+    blas::Queue& queue );
+
+void ger(
+    blas::Layout layout,
+    int64_t m, int64_t n,
+    double alpha,
+    double const* x, int64_t incx,
+    double const* y, int64_t incy,
+    double*       A, int64_t lda,
+    blas::Queue& queue );
+
+void ger(
+    blas::Layout layout,
+    int64_t m, int64_t n,
+    std::complex<float> alpha,
+    std::complex<float> const* x, int64_t incx,
+    std::complex<float> const* y, int64_t incy,
+    std::complex<float>*       A, int64_t lda,
+    blas::Queue& queue );
+
+void ger(
+    blas::Layout layout,
+    int64_t m, int64_t n,
+    std::complex<double> alpha,
+    std::complex<double> const* x, int64_t incx,
+    std::complex<double> const* y, int64_t incy,
+    std::complex<double>*       A, int64_t lda,
+    blas::Queue& queue );
+
+//------------------------------------------------------------------------------
+void geru(
+    blas::Layout layout,
+    int64_t m, int64_t n,
+    float alpha,
+    float const* x, int64_t incx,
+    float const* y, int64_t incy,
+    float*       A, int64_t lda,
+    blas::Queue& queue );
+
+void geru(
+    blas::Layout layout,
+    int64_t m, int64_t n,
+    double alpha,
+    double const* x, int64_t incx,
+    double const* y, int64_t incy,
+    double*       A, int64_t lda,
+    blas::Queue& queue );
+
+void geru(
+    blas::Layout layout,
+    int64_t m, int64_t n,
+    std::complex<float> alpha,
+    std::complex<float> const* x, int64_t incx,
+    std::complex<float> const* y, int64_t incy,
+    std::complex<float>*       A, int64_t lda,
+    blas::Queue& queue );
+
+void geru(
+    blas::Layout layout,
+    int64_t m, int64_t n,
+    std::complex<double> alpha,
+    std::complex<double> const* x, int64_t incx,
+    std::complex<double> const* y, int64_t incy,
+    std::complex<double>*       A, int64_t lda,
+    blas::Queue& queue );
+
 //==============================================================================
 // Level 3 BLAS
 

src/cublas_wrappers.cc

@@ -1348,6 +1348,162 @@ void trsv(
             (cuDoubleComplex*) dx, incdx ) );
 }
 
+//------------------------------------------------------------------------------
+// ger
+//------------------------------------------------------------------------------
+void ger(
+    int64_t m, int64_t n,
+    float alpha,
+    float const* dx, int64_t incdx,
+    float const* dy, int64_t incdy,
+    float*       dA, int64_t ldda,
+    blas::Queue& queue )
+{
+    blas_dev_call(
+        cublasSger(
+            queue.handle(),
+            m, n,
+            &alpha,
+            dx, incdx,
+            dy, incdy,
+            dA, ldda ) );
+}
+
+//------------------------------------------------------------------------------
+void ger(
+    int64_t m, int64_t n,
+    double alpha,
+    double const* dx, int64_t incdx,
+    double const* dy, int64_t incdy,
+    double*       dA, int64_t ldda,
+    blas::Queue& queue )
+{
+    blas_dev_call(
+        cublasDger(
+            queue.handle(),
+            m, n,
+            &alpha,
+            dx, incdx,
+            dy, incdy,
+            dA, ldda ) );
+}
+
+//------------------------------------------------------------------------------
+void ger(
+    int64_t m, int64_t n,
+    std::complex<float> alpha,
+    std::complex<float> const* dx, int64_t incdx,
+    std::complex<float> const* dy, int64_t incdy,
+    std::complex<float>*       dA, int64_t ldda,
+    blas::Queue& queue )
+{
+    blas_dev_call(
+        cublasCgerc(
+            queue.handle(),
+            m, n,
+            (cuComplex*) &alpha,
+            (cuComplex*) dx, incdx,
+            (cuComplex*) dy, incdy,
+            (cuComplex*) dA, ldda ) );
+}
+
+//------------------------------------------------------------------------------
+void ger(
+    int64_t m, int64_t n,
+    std::complex<double> alpha,
+    std::complex<double> const* dx, int64_t incdx,
+    std::complex<double> const* dy, int64_t incdy,
+    std::complex<double>*       dA, int64_t ldda,
+    blas::Queue& queue )
+{
+    blas_dev_call(
+        cublasZgerc(
+            queue.handle(),
+            m, n,
+            (cuDoubleComplex*) &alpha,
+            (cuDoubleComplex*) dx, incdx,
+            (cuDoubleComplex*) dy, incdy,
+            (cuDoubleComplex*) dA, ldda ) );
+}
+
+//------------------------------------------------------------------------------
+// geru
+//------------------------------------------------------------------------------
+void geru(
+    int64_t m, int64_t n,
+    float alpha,
+    float const* dx, int64_t incdx,
+    float const* dy, int64_t incdy,
+    float*       dA, int64_t ldda,
+    blas::Queue& queue )
+{
+    blas_dev_call(
+        cublasSger(
+            queue.handle(),
+            m, n,
+            &alpha,
+            dx, incdx,
+            dy, incdy,
+            dA, ldda ) );
+}
+
+//------------------------------------------------------------------------------
+void geru(
+    int64_t m, int64_t n,
+    double alpha,
+    double const* dx, int64_t incdx,
+    double const* dy, int64_t incdy,
+    double*       dA, int64_t ldda,
+    blas::Queue& queue )
+{
+    blas_dev_call(
+        cublasDger(
+            queue.handle(),
+            m, n,
+            &alpha,
+            dx, incdx,
+            dy, incdy,
+            dA, ldda ) );
+}
+
+//------------------------------------------------------------------------------
+void geru(
+    int64_t m, int64_t n,
+    std::complex<float> alpha,
+    std::complex<float> const* dx, int64_t incdx,
+    std::complex<float> const* dy, int64_t incdy,
+    std::complex<float>*       dA, int64_t ldda,
+    blas::Queue& queue )
+{
+    blas_dev_call(
+        cublasCgeru(
+            queue.handle(),
+            m, n,
+            (cuComplex*) &alpha,
+            (cuComplex*) dx, incdx,
+            (cuComplex*) dy, incdy,
+            (cuComplex*) dA, ldda ) );
+}
+
+//------------------------------------------------------------------------------
+void geru(
+    int64_t m, int64_t n,
+    std::complex<double> alpha,
+    std::complex<double> const* dx, int64_t incdx,
+    std::complex<double> const* dy, int64_t incdy,
+    std::complex<double>*       dA, int64_t ldda,
+    blas::Queue& queue )
+{
+    blas_dev_call(
+        cublasZgeru(
+            queue.handle(),
+            m, n,
+            (cuDoubleComplex*) &alpha,
+            (cuDoubleComplex*) dx, incdx,
+            (cuDoubleComplex*) dy, incdy,
+            (cuDoubleComplex*) dA, ldda ) );
+}
+
 //==============================================================================
 // Level 3 BLAS - Device Interfaces