Re: [R] Segfaults of eigen

2011-06-07 Thread Juergen Rose
Dear Prof. Ripley, Am Montag, den 21.02.2011, 11:40 + schrieb Prof Brian Ripley: > So there is very likely a bug in your system software or compiler. > There is no such problem on other x86_64 systems (including > Fedora 14 and 12, FreeBSD, Solaris and Windows). > > You can debug segfaults f

Re: [R] Segfaults of eigen with blas-atlas at x86_64-pc-linux-gnu systems

2011-02-21 Thread Karl Ove Hufthammer
Juergen Rose wrote: > What are the differences between blas-reference, blas-atlas and > blas-goto. They are different. The reference BLAS is (relatively) slow; the other ones are fast. If you’re happy with the performance of the reference BLAS, or don’t use much matrix algebra-dependent code, t

Re: [R] Segfaults of eigen with blas-atlas at x86_64-pc-linux-gnu systems

2011-02-21 Thread Juergen Rose
Am Montag, den 21.02.2011, 11:25 +0100 schrieb Karl Ove Hufthammer: > Juergen Rose wrote: > > >> eigen(D) > > > > *** caught segfault *** > > address (nil), cause 'unknown' > > > > Traceback: > > 1: .Call("La_rs", x, only.values, PACKAGE = "base") > > 2: eigen(D) > > > > All systems are Gentoo

Re: [R] Segfaults of eigen

2011-02-21 Thread Prof Brian Ripley
So there is very likely a bug in your system software or compiler. There is no such problem on other x86_64 systems (including Fedora 14 and 12, FreeBSD, Solaris and Windows). You can debug segfaults for yourself (see 'Writing R Extensions') or ask your vendor for support. On Mon, 21 Feb 201

Re: [R] Segfaults of eigen

2011-02-21 Thread Karl Ove Hufthammer
Juergen Rose wrote: >> eigen(D) > > *** caught segfault *** > address (nil), cause 'unknown' > > Traceback: > 1: .Call("La_rs", x, only.values, PACKAGE = "base") > 2: eigen(D) > > All systems are Gentoo systems, i.e. R-2.12.1 and blas-atlas-3.9.23-r4 > is installed by compiling the sources. Rec

Re: [R] Segfaults of eigen

2011-02-21 Thread Juergen Rose
Am Montag, den 21.02.2011, 11:11 +0100 schrieb Juergen Rose: > Hi, > > with small matrices eigen works as expected: > > > eigen(cbind(c(1,4),c(4,7)), only.values = TRUE) > $values > [1] 9 -1 > > $vectors > NULL > > > eigen(cbind(c(1,4),c(4,7))) > $values > [1] 9 -1 > > $vectors > [

[R] Segfaults of eigen

2011-02-21 Thread Juergen Rose
Hi, with small matrices eigen works as expected: > eigen(cbind(c(1,4),c(4,7)), only.values = TRUE) $values [1] 9 -1 $vectors NULL > eigen(cbind(c(1,4),c(4,7))) $values [1] 9 -1 $vectors [,1] [,2] [1,] 0.4472136 -0.8944272 [2,] 0.8944272 0.4472136 > eigen(cbind(c(1,-1),c(1,-