library(BiocParallel)

ll <- replicate(8, matrix(rnorm(1e6),1000), simplify=FALSE)
f <- function(x) mean(solve(x), trim=0.7)

(p <- SnowParam(4L))
bplapply(ll, f, BPPARAM = p)