Dear R HPC People:

I have the following function, which works:
> library(doMPI)
Loading required package: foreach
Loading required package: iterators
Loading required package: codetools
Loading required package: Rmpi
> cl <- startMPIcluster()
	58 slaves are spawned successfully. 0 failed.
> registerDoMPI(cl)
> xa <- seq(from=1,by=3,length=120)
> xb <- seq(from=3,by=3,length=120)
> zz <- foreach(i=1:10000,.combine=c) %dopar% {
+  x <- ts(arima.sim(list(order=c(1,0,0),ar=-0.9),n=360),start=1975,freq=12)
+  y <- aggregate(x,nfreq=4,sum)
+  arima(y,order=c(1,0,0))$model$phi
+  }
>      closeCluster(cl)
> 
>      
> 
[1] "Please use mpi.quit() to quit R"
> proc.time()
   user  system elapsed 
184.678   3.896 378.648 

But I thought that it might run faster.  Are my statements not parallelizable or have I done something wrong, please?
Finally, do I look at user time or elapsed time, please?

Thanks,
Erin


Erin M. Hodgess, PhD
Associate Professor
Department of Computer and Mathematical Sciences
University of Houston - Downtown
mailto: hodgesse@uhd.edu


	[[alternative HTML version deleted]]

