Add lower threshold for multithreading

diff --git a/interface/lapack/getrf.c b/interface/lapack/getrf.c
index 02bb124..323370e 100644 (file)
--- a/interface/lapack/getrf.c
+++ b/interface/lapack/getrf.c
@@ -95,7 +95,14 @@ int NAME(blasint *M, blasint *N, FLOAT *a, blasint *ldA, blasint *ipiv, blasint
 
 #ifdef SMP
   args.common = NULL;
-  args.nthreads = num_cpu_avail(4);
+#ifndef DOUBLE
+  if (args.m*args.n < 40000)
+#else
+  if (args.m*args.n < 10000)
+#endif
+       args.nthreads=1;
+  else
+       args.nthreads = num_cpu_avail(4);
 
   if (args.nthreads == 1) {
 #endif

diff --git a/interface/lapack/zgetrf.c b/interface/lapack/zgetrf.c
index 7f8db94..d03541f 100644 (file)
--- a/interface/lapack/zgetrf.c
+++ b/interface/lapack/zgetrf.c
@@ -95,7 +95,10 @@ int NAME(blasint *M, blasint *N, FLOAT *a, blasint *ldA, blasint *ipiv, blasint
 
 #ifdef SMP
   args.common = NULL;
-  args.nthreads = num_cpu_avail(4);
+  if (args.m*args.n <10000)
+       args.nthreads = 1;
+  else
+       args.nthreads = num_cpu_avail(4);
 
   if (args.nthreads == 1) {
 #endif