Bumping up clairs-to allocated memory and providing half of allocted threads

Author: skchronicles

config/cluster/slurm.json

@@ -76,7 +76,7 @@
         "threads": "24",
         "mem": "64G",
         "time": "1-18:00:00",
-        "gres": "lscratch:800"
+        "gres": "lscratch:1200"
     },
     "deepvariant": {
         "threads": "18",
@@ -297,7 +297,7 @@
         "threads": "24"
     },
     "clairs_tumor_only": {
-        "mem": "64G",
+        "mem": "120G",
         "time" : "1-00:00:00",
         "threads": "24"
     },

config/cluster/uge.json

@@ -287,7 +287,7 @@
     "threads": "8"
   },
   "clairs_tumor_only": {
-    "mem": "4G",
+    "mem": "5G",
     "partition": "",
     "threads": "8"
   },

workflow/rules/somatic.smk

@@ -50,6 +50,7 @@ def get_normal_pileup_table(wildcards):
         # Runs in tumor-only mode
         return []
 
+
 def get_somatic_tn_callers(wildcards):
     """Returns somatic variants found with tumor-normal variant
     callers. For tumor-normal samples, extra somatic callers 
@@ -706,7 +707,10 @@ rule clairs_tumor_only:
         genome  = config['references']['GENOME'],
         outdir = join(workpath, "clairs", "somatic", "{name}"),
     threads: 
-        int(allocated("threads", "clairs_tumor_only", cluster)),
+        # ClairS-TO over utilizes threads,
+        # testing has shown it over utilizes
+        # around 50% of the threads allocated
+        max(int(int(allocated("threads", "clairs_tumor_only", cluster))/2.0), 2),
     container: config['images']['clairs-to']
     envmodules: config['tools']['rlang']
     shell: """
@@ -789,7 +793,7 @@ rule deepsomatic:
     # mechanism for deletion on exit
     if [ ! -d "{params.tmpdir}" ]; then mkdir -p "{params.tmpdir}"; fi
     tmp=$(mktemp -d -p "{params.tmpdir}")
-    trap 'rm -rf "${{tmp}}"' EXIT
+    trap 'du -sh "${{tmp}}"; rm -rf "${{tmp}}"' EXIT
 
     run_deepsomatic \\
         --model_type={params.dv_model_type} \\