1 Reply Latest reply on Jun 13, 2015 8:37 AM by jtrudeau

    error in BLAS - ACML

    lejeczek

      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