fix .combine in seqParallel()

NEWS

@@ -11,11 +11,18 @@ UTILITIES
     o new `seqGetData(, "$dosage_alt2")` and `seqGetData(, "$dosage_sp2")` for
       sex chromosomes, when the alleles are partially missing
 
+
+CHANGES IN VERSION 1.42.3
+-------------------------
+
 BUG FIXES
 
     o `seqGetData(, "$dosage_alt")` and `seqGetData(, "$dosage_sp")` work
       correctly when the ploidy is >2 and there are missing alleles
 
+    o fix a bug that `seqParallel()` does not call a user-defined '.combine'
+      when 'parallel=1'
+
 
 CHANGES IN VERSION 1.42.1
 -------------------------

R/Utilities.R

@@ -336,6 +336,14 @@ seqParallel <- function(cl=seqGetParallel(), gdsfile, FUN,
                 ans <- FUN(gdsfile, ...)
         }
 
+        if (is.function(.combine))
+        {
+            if (length(formals(args(.combine))) == 1L)
+                .combine(ans)
+            else
+                ans <- .combine(NULL, ans)
+        }
+
     } else if (inherits(cl, "cluster"))
     {
         ## multiple processes with a predefined cluster

README.md

@@ -18,7 +18,7 @@ The SeqArray package is built on top of [Genomic Data Structure (GDS)](https://g
 
 ## Bioconductor:
 
-Release Version: v1.42.1
+Release Version: v1.42.3
 
 [http://www.bioconductor.org/packages/release/bioc/html/SeqArray.html](http://www.bioconductor.org/packages/SeqArray)