Comparing execution times
Just a comment about your use of foreach:
On Mon, Apr 11, 2011 at 6:29 AM, Alaios <alaios at yahoo.com> wrote:
[snip]
C.Case. Foreach is considered to be easier to be applied to manycores.
foreach (i=1:dimz) %do% {
? ?print(sprintf('Creating the %d map',i));
? ?Shadowlist[,,i]<-f <- GaussRF(x=x, y=y, model=model, grid=TRUE,param=c(mean,variance,nugget,scale,Whit.alpha))
}
You are still running this sequentially.
To run in parallel, you need load the appropriate parallel backend,
and use %dopar%:
library(doMC)
foreach(i=1:dimz) %dopar% { ... }
Steve Lianoglou Graduate Student: Computational Systems Biology ?| Memorial Sloan-Kettering Cancer Center ?| Weill Medical College of Cornell University Contact Info: http://cbio.mskcc.org/~lianos/contact