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]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
