AnsweredAssumed Answered

error in BLAS - ACML

Question asked by lejeczek on May 21, 2015
Latest reply on Jun 13, 2015 by jtrudeau

it seems old (bug), using ACML 5.3.1

 

> ## PR#4582 %*% with NAs

> stopifnot(is.na(NA %*% 0), is.na(0 %*% NA))

> ## depended on the BLAS in use.

>

>

> ## found from fallback test in slam 0.1-15

> ## most likely indicates an inaedquate BLAS.

> x <- matrix(c(1, 0, NA, 1), 2, 2)

> y <- matrix(c(1, 0, 0, 2, 1, 0), 3, 2)

> (z <- tcrossprod(x, y))

     [,1] [,2] [,3]

[1,]   NA   NA    0

[2,]    2    1    0

> stopifnot(identical(z, x %*% t(y)))

Error: identical(z, x %*% t(y)) is not TRUE

 

here is RHEL bug I reported some (longish) time ago:

https://bugzilla.redhat.com/show_bug.cgi?id=1113657

 

Could AMD please comment on the issue? shed more light please?

regards

Outcomes