From f5fcc5baec1c5aea7dbd7a2a8fdd41ae8b422a6e Mon Sep 17 00:00:00 2001
From: Martin Kroeker <martin@ruby.chemie.uni-freiburg.de>
Date: Sat, 15 Aug 2020 13:30:29 +0200
Subject: [PATCH] Add trivial gemm test for multithread consistency

---
 cpp_thread_test/gemm64.cpp | 20 ++++++++++++++++++++
 1 file changed, 20 insertions(+)
 create mode 100644 cpp_thread_test/gemm64.cpp

diff --git a/cpp_thread_test/gemm64.cpp b/cpp_thread_test/gemm64.cpp
new file mode 100644
index 000000000..2c3442a2e
--- /dev/null
+++ b/cpp_thread_test/gemm64.cpp
@@ -0,0 +1,20 @@
+#include <iostream>
+#include <cblas.h>
+int main ( int argc, char* argv[] ) {
+    const long n = ((long)1 << 31) - 1;
+    std::cout << n <<std::endl;
+    float*  A = new float[n];
+    float*  B = new float[n];
+    float*  C = new float[1];
+    for(long i =0; i <n; i++){
+        A[i] = 1;
+        B[i] = 1;
+
+    }
+    cblas_sgemm(CblasRowMajor, CblasNoTrans, CblasNoTrans, 1, 1, n,  1.0, A, n, B, 1, 0.0, C, 1);
+    std::cout << *C <<std::endl;
+    delete[] A;
+    delete[] B;
+    delete[] C;
+    return 0;
+}