Add java options to picard tools

config.yaml

@@ -562,6 +562,14 @@ params:
     # Valid values: "processed", "merged"
     CollectMultipleMetrics-bams: "processed"
 
+    # For some specific error cases, it might be necessary to adjust java settings for the tools,
+    # such as `-Xmx10g` to increase the available memory within the Java virtual machine
+    # (because it limits itself otherwise, instead of using the available memory...).
+    # The last option, SortVcf-java-opts, is used by bcftools when using contig-group-size > 0.
+    MarkDuplicates-java-opts: ""
+    CollectMultipleMetrics-java-opts: ""
+    SortVcf-java-opts: ""
+
   # ----------------------------------------------------------------------
   #     dedup
   # ----------------------------------------------------------------------

rules/calling-bcftools.smk

@@ -65,16 +65,18 @@ if config["settings"].get("contig-group-size"):
                 " USE_JDK_DEFLATER=true USE_JDK_INFLATER=true"
                 if platform.system() == "Darwin"
                 else ""
-            )
+            ),
+            java_opts=config["params"]["picard"]["SortVcf-java-opts"]
         log:
-            "logs/vcflib/sort-genotyped.log"
+            "logs/picard/sort-genotyped.log"
         benchmark:
             "benchmarks/picard/sort-genotyped.bench.log"
         conda:
             "../envs/picard.yaml"
         shell:
             # Weird new picard syntax...
             "picard SortVcf "
+            "{params.java_opts} "
             "INPUT={input.vcf} "
             "OUTPUT={output.vcf} "
             "SEQUENCE_DICTIONARY={input.refdict} "

rules/duplicates-picard.smk

@@ -31,6 +31,7 @@ rule mark_duplicates:
         # and some system libraries, leading the JRE to exit with fatal error SIGSEGV caused by
         # libgkl_compression, see https://github.com/broadinstitute/picard/issues/1329.
         # Hence, on MacOS, we add the two settings recommended by the github issue.
+        config["params"]["picard"]["MarkDuplicates-java-opts"] + " " +
         config["params"]["picard"]["MarkDuplicates"] + (
             " USE_JDK_DEFLATER=true USE_JDK_INFLATER=true"
             if platform.system() == "Darwin"

rules/qc-bam.smk

@@ -191,6 +191,7 @@ rule picard_collectmultiplemetrics:
     log:
         "logs/picard/multiple_metrics/{sample}.log"
     params:
+        config["params"]["picard"]["CollectMultipleMetrics-java-opts"] + " " +
         config["params"]["picard"]["CollectMultipleMetrics-extra"] + (
             " USE_JDK_DEFLATER=true USE_JDK_INFLATER=true"
             if platform.system() == "Darwin"