# 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)
R <- qr.R(qr(x))
xtx.inv <- chol2inv(R)
id <- xtx.inv
print(id)
finalize()
Run the code above in your browser using DataLab