From b407512054adf74695298dec3d53988b2404adb2 Mon Sep 17 00:00:00 2001
From: skchronicles <kuhnsa3@gmail.com>
Date: Thu, 31 Oct 2024 14:59:07 -0400
Subject: [PATCH] Exporting OpenBLAS thread pool variable prior to running
 deepsomatic

---
 workflow/rules/somatic.smk | 15 +++++++++++++++
 1 file changed, 15 insertions(+)

diff --git a/workflow/rules/somatic.smk b/workflow/rules/somatic.smk
index eda5f93..3c7398d 100644
--- a/workflow/rules/somatic.smk
+++ b/workflow/rules/somatic.smk
@@ -797,6 +797,21 @@ rule deepsomatic:
     tmp=$(mktemp -d -p "{params.tmpdir}")
     trap 'du -sh "${{tmp}}"; rm -rf "${{tmp}}"' EXIT
 
+    # Export OpenBLAS variable to
+    # control the number of threads
+    # in a thread pool. By setting
+    # this variable to 1, work is
+    # done in the thread that ran
+    # the operation, rather than
+    # disbatching the work to a
+    # thread pool. If this option
+    # is not provided, it can lead
+    # to nested parallelism.
+    # See this issue for more info:
+    # https://github.com/google/deepsomatic/issues/28
+    export OPENBLAS_NUM_THREADS=1
+
+    # Run deepsomatic
     run_deepsomatic \\
         --model_type={params.dv_model_type} \\
         --ref={params.genome} \\