From 496106642f77c7d319ee9966c067cdf0040462ab Mon Sep 17 00:00:00 2001
From: Martin Kroeker <martin@ruby.chemie.uni-freiburg.de>
Date: Wed, 27 Mar 2024 20:32:11 +0100
Subject: [PATCH] Cap the number of parallel threads

---
 interface/lapack/potrf.c | 16 ++++++++++------
 1 file changed, 10 insertions(+), 6 deletions(-)

diff --git a/interface/lapack/potrf.c b/interface/lapack/potrf.c
index 3abc80133..d4a56c94a 100644
--- a/interface/lapack/potrf.c
+++ b/interface/lapack/potrf.c
@@ -113,13 +113,17 @@ int NAME(char *UPLO, blasint *N, FLOAT *a, blasint *ldA, blasint *Info){
 #ifdef SMP
   args.common = NULL;
 #ifndef DOUBLE
-  if (args.n <128)
-#else
-  if (args.n <64)
-#endif
+  int nmax = 128;
+#else 
+  int nmax = 64;
+endif
+  if (args.n <nmax) {
     args.nthreads = 1;
-  else
-  args.nthreads = num_cpu_avail(4);
+  } else {
+    args.nthreads = num_cpu_avail(4);
+    if (args.n/args.nthreads <nmax)
+	    args.nthreads = args.n/nmax;
+  }
 
   if (args.nthreads == 1) {
 #endif