r - 如何并行运行基本的 R 光栅操作？

Question

我有这些示例栅格：

s <- stack(replicate(10, {raster(matrix(runif(2160*4320), 2160),
xmn=-180, xmx=180, ymn=-90, ymx=90)}))

s2 <- stack(replicate(10, {raster(matrix(runif(2160*4320), 2160), 
xmn=-180, xmx=180, ymn=-90, ymx=90)}))

r <- raster(nrow=3000, ncol=7200, ymn=-60, ymx=90) 
values(r) <- 1:ncell(r)

我希望能够使用多核来运行 R 和 Raster 基本功能，例如：

s * s2
resample(s, r)