lapack and R.

Jose' Matos jamatos at fc.up.pt
Fri Sep 23 17:47:34 UTC 2005


Hi Tom,
  while reading the R-devel list I found a bug reported by Peter Dalgaard 
against the lapack package in FC4. The bug follows:

*************************************************************************
Hmm. Doesn't look like it is actually working, though. Install
lapack-devel, configure --with-lapack, and make check dies with

running code in 'base-Ex.R' ...make[4]: *** [base-Ex.Rout] Error 1
make[4]: Leaving directory `/home/pd/r-devel/BUILD/tests/Examples'
make[3]: *** [test-Examples-Base] Error 2
....
[pd at titmouse BUILD]$ tail tests/Examples/base-Ex.Rout.fail
> kappa(x2 <- cbind(x1,2:11))# high! [x2 is singular!]
[1] 8.351867e+16
>
> hilbert <- function(n) { i <- 1:n; 1 / outer(i - 1, i, "+") }
> sv9 <- svd(h9 <- hilbert(9))$ d
> kappa(h9)# pretty high!
[1] 728289254735
> kappa(h9, exact = TRUE) == max(sv9) / min(sv9)
Error in La.svd(x, nu, nv) : BLAS/LAPACK routine 'DGEBRD' gave error code -10
Execution halted

This happens on both x86_64 and x86 installs of FC4.

I have a strong sense of deja vu regarding this error.
*************************************************************************

It looks like the "dejá vu" feeling comes from a previous bad patch in Debian.

Should I add this to Bugzilla? And then under R or under lapack? :-)

All the best,
-- 
José Abílio




More information about the fedora-extras-list mailing list