zgemv: Add a workaround for OpenMathLib#746

jeromerobert · Jan 21, 2016 · 72308e9 · 72308e9
1 parent 2f77aba
commit 72308e9
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/interface/zgemv.c b/interface/zgemv.c
@@ -244,7 +244,9 @@ void CNAME(enum CBLAS_ORDER order,
   // for alignment
   buffer_size = (buffer_size + 3) & ~3;
   STACK_ALLOC(buffer_size, FLOAT, buffer);
-
+  if(trans != CblasNoTrans)
+    // cgemv_t.S return NaN if there are NaN or Inf in the buffer (see bug #746)
+    memset(buffer, 0, sizeof(FLOAT) * buffer_size);
 #ifdef SMP
 
   nthreads_max = num_cpu_avail(2);