[Rd] looks in liblapack.a not liblapack.so

José Matos jamatos at fc.up.pt
Wed Sep 28 16:12:42 CEST 2005


Peter Dalgaard wrote:

> 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 am sorry Peter, I am trying really hard to replicate this bug but I have
not been able to see the same result, no matter what I try.

  I have download the latest tar ball and then I run:

$ ./configure --with-lapack="-llapack -lblas"
...
R is now configured for x86_64-unknown-linux-gnu

  Source directory:          .
  Installation directory:    /usr/local

  C compiler:                gcc  -g -O2
  C++ compiler:              g++  -g -O2
  Fortran compiler:          gfortran  -g -O2

  Interfaces supported:      X11
  External libraries:        readline
  Additional capabilities:   PNG, JPEG, iconv, MBCS, NLS
  Options enabled:           R profiling

  Recommended packages:      yes

$ make -j8
$ make check

  It works.

  OTOH I am not sure that configure is accepting my options. Looking into
config.log I don't see that value being used, and I noticed that the lapack
module it is still being built.

  I read docs/manual/R-admin.html but without any difference. I have tried
different forms:

--with-lapack
--with-lapack="-llapack -lcblas"
--with-lapack="-llapack -lblas"

  What am I missing?

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

  Thanks,
-- 
José Abílio



More information about the R-devel mailing list