Merge pull request #4694 from martin-frbg/issue3660

Add a minimum problem size for multithreading in GBMV

Author: martin-frbg

interface/gbmv.c

@@ -227,7 +227,10 @@ void CNAME(enum CBLAS_ORDER order,
   buffer = (FLOAT *)blas_memory_alloc(1);
 
 #ifdef SMP
-  nthreads = num_cpu_avail(2);
+  if (m * n < 250000 || kl+ku < 15 )
+    nthreads = 1;
+  else
+    nthreads = num_cpu_avail(2);
 
   if (nthreads == 1) {
 #endif

interface/zgbmv.c

@@ -251,7 +251,10 @@ void CNAME(enum CBLAS_ORDER order,
   buffer = (FLOAT *)blas_memory_alloc(1);
 
 #ifdef SMP
-  nthreads = num_cpu_avail(2);
+  if (m * n  < 125000 || ku + kl < 15) 
+    nthreads = 1;
+  else
+    nthreads = num_cpu_avail(2);
 
   if (nthreads == 1) {
 #endif