# NOT RUN {
# }
# NOT RUN {
# Save code in a file "demo.r" and run with 2 processors by
# > mpiexec -np 2 Rscript demo.r
library(pbdDMAT, quiet = TRUE)
init.grid()
comm.set.seed(diff=T)
x <- ddmatrix("rnorm", 3, 3, bldim=2)
R <- qr.R(qr(x))
xtx.inv <- chol2inv(R)
id <- as.matrix(xtx.inv %*% crossprod(x))
comm.print(id)
finalize()
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab