regionMatrix() now has an mc.cores argument and uses mcmapply() instead of mapply(). Use with care with large fullCov objects -- have yet to test on a real data set