Nicolas Vannieuwkerke
04/22/2026, 1:26 PMError in reducer$value.cache[[as.character(idx)]] <- values :
wrong args for environment subassignment
Calls: countSplitReads ... .bploop_impl -> .collect_result -> .reducer_add -> .reducer_add
In addition: Warning message:
In parallel::mccollect(wait = FALSE, timeout = 1) :
1 parallel job did not deliver a result
Execution halted
Anyone any ideas?Ata Jadid Ahari
04/22/2026, 1:27 PMNicolas Vannieuwkerke
04/22/2026, 1:27 PMNicolas Vannieuwkerke
04/22/2026, 2:20 PMAta Jadid Ahari
04/22/2026, 2:29 PMNicolas Vannieuwkerke
04/22/2026, 2:32 PMNicolas Vannieuwkerke
04/23/2026, 7:53 AMAta Jadid Ahari
04/23/2026, 9:54 AMAta Jadid Ahari
04/23/2026, 9:55 AMNicolas Vannieuwkerke
04/23/2026, 10:50 AMNicolas Vannieuwkerke
04/23/2026, 10:50 AMNicolas Vannieuwkerke
04/23/2026, 12:03 PMNicolas Vannieuwkerke
04/23/2026, 12:03 PMNicolas Vannieuwkerke
04/23/2026, 12:05 PMAta Jadid Ahari
04/23/2026, 12:07 PMNicolas Vannieuwkerke
04/23/2026, 12:12 PMNicolas Vannieuwkerke
04/23/2026, 12:13 PMAta Jadid Ahari
04/23/2026, 12:15 PMNicolas Vannieuwkerke
04/23/2026, 12:15 PMNicolas Vannieuwkerke
04/23/2026, 12:15 PMJash Trivedi
05/12/2026, 3:01 AMGENECOUNTS_COUNTREADS & SPLICECOUNTS_COUNTSPLITREADS tasks keep failing (exit code 140) on our HPC despite allocating what I thought was a generous amount of resources. For context, I am working with BAM files that have ~100M paired-end reads per sample and have total 370 samples. I’ve bumped the base memory to 128GB and set up dynamic retries in my hpc.config, but the jobs are still getting killed:
withName: '.*GENECOUNTS_COUNTREADS' {
memory = { 128.GB * task.attempt }
time = { 24.h * task.attempt }
errorStrategy = { task.exitStatus in [1, 104, 134, 137, 139, 140, 143, 255] ? 'retry' : 'finish' }
}
To troubleshoot, I also tried changing the number of cores/workers. When setting it to n-cpu=2 I ran into a Nextflow resource mismatch for the splicing step:
ERROR ~ Error executing process > 'DROP:ABERRANTSPLICING:SPLICECOUNTS_COUNTSPLITREADS (F216P003-100m_524015526)'
Caused by:
Process requirement exceeds available CPUs -- req: 3; avail: 2
When setting it to n-cpus=4, it throws this BiocParallel error (which I suspect is an underlying Out-Of-Memory kill being masked by R):
Error in reducer$value.cache[[as.character(idx)]] <- values :
wrong args for environment subassignment
Calls: countSplitReads ... .bploop_impl -> .collect_result -> .reducer_add -> .reducer_add
In addition: Warning message:
In parallel::mccollect(wait = FALSE, timeout = 1) :
What would be the best way forward here? Should I just keep throwing more memory/time at it, or is there a way to configure the underlying BiocParallel workers to chunk the BAM processing more efficiently?
On a related note regarding the handling of large datasets: has there been any discussion within the community about rebuilding the core counting logic or the FRASER/OUTRIDER packages in Rust instead of R? Given the memory overhead of `GenomicAlignments`/`BiocParallel`, I imagine a Rust-based implementation could drastically improve both memory efficiency and speed. With cohorts getting larger, the current BiocParallel + R memory model seems difficult to scale reliably on HPC environments. Would love to hear if this has been considered!Nicolas Vannieuwkerke
05/12/2026, 7:23 AMcpus directive for the process?Nicolas Vannieuwkerke
05/12/2026, 7:24 AMAta Jadid Ahari
05/13/2026, 9:05 AMAta Jadid Ahari
05/13/2026, 9:06 AMNicolas Vannieuwkerke
05/13/2026, 9:07 AMJash Trivedi
05/13/2026, 4:47 PMERROR ~ Error executing process > 'DROP:ABERRANTSPLICING:SPLICECOUNTS_COUNTSPLITREADS (F216P003-100m_524015526)'
Caused by:
Process requirement exceeds available CPUs -- req: 3; avail: 2Nicolas Vannieuwkerke
05/18/2026, 2:13 PMJash Trivedi
05/19/2026, 6:40 PMNicolas Vannieuwkerke
05/21/2026, 7:35 AM