14-blas-opt 363 Bytes
Newer Older
1
Further improvements in the reference BLAS routines.
2 3 4

DSYRK improved, which affects performance of R crossprod and tcrossprod.

Radford Neal's avatar
Radford Neal committed
5 6
DGEMM is improved for all cases except when the first but not second
operand is tranposed.
Radford Neal's avatar
Radford Neal committed
7 8 9 10 11 12 13

Also fixes the bug illustrate below:

  a <- matrix(c(NA,1,0,1),2,2,byrow=T)
  print(a%*%t(a))
  print(tcrossprod(a,a))
  print(tcrossprod(a))