Hi,

I am trying to use GotoBLAS2 on R 3.0 on Unix. I downloaded GotoBLAS2
source code from TACC web site, compiled it, and replaced libRblas.so with
libgoto2.so, following the instructions at the link
http://www.rochester.edu/college/gradstudents/jolmsted/files/computing/BLAS.pdf.
The simple matrix operations in R like "determinant" are 20 times faster
than before (I am using huge matrices), which is good. However, I cannot
use many cores in parallel now.

For example, below code runs forever. But if I use commented out "for"
instead of "foreach", it takes just a second. When I was using R's default
BLAS library, I could run below code (using many cores) (but it took more
time since BLAS was not optimized, of course)..

library("foreach")
library("doParallel")

registerDoParallel(cores=2)set.seed(100)
foreach (i = 1:2) %dopar% {# for (i in 1:2) {
a = replicate(1000, rnorm(1000))
d = determinant(a)

So, is it possible to use many cores at the same time with GotoBLAS2, do
you have any ideas?

Thanks a lot in advance.

-- 
-safiye

	[[alternative HTML version deleted]]

