Merge pull request #108 from SwayamInSync/fix-matmul-layout-89

adding F-contig layout support + other cleanups

Author: ngoldbaum

src/csrc/quadblas_interface.cpp

@@ -46,9 +46,22 @@ qblas_gemv(char layout, char trans, size_t m, size_t n,
            Sleef_quad *beta, Sleef_quad *y, size_t incy)
 {
     if (m == 0 || n == 0) {
+        int do_trans = (trans == 'T' || trans == 't' || trans == 'C' || trans == 'c');
+        size_t y_len = do_trans ? n : m;
+        if (y_len == 0) {
+            return 0;
+        }
+        if (!y || !beta) {
+            return -1;
+        }
+        Sleef_quad zero = Sleef_cast_from_doubleq1(0.0);
+        int beta_zero = Sleef_icmpeqq1(*beta, zero);
+        for (size_t i = 0; i < y_len; i++) {
+            y[i * incy] = beta_zero ? zero : Sleef_mulq1_u05(*beta, y[i * incy]);
+        }
         return 0;
     }
-    if (!alpha || !A || !x || !beta || !y) {
+    if (!alpha || !A || !x || !y) {
         return -1;
     }
     cblas_qgemv(to_layout(layout), to_trans(trans),
@@ -66,7 +79,7 @@ qblas_gemm(char layout, char transa, char transb,
            Sleef_quad *B, size_t ldb,
            Sleef_quad *beta, Sleef_quad *C, size_t ldc)
 {
-    if (m == 0 || n == 0 || k == 0) {
+    if (m == 0 || n == 0) {
         return 0;
     }
     if (!alpha || !A || !B || !beta || !C) {