core dump on 0.64.2 SPARC/Solaris 2.6 in eigen (PR#223)

Prof Brian D Ripley ripley@stats.ox.ac.uk
Mon, 12 Jul 1999 07:49:34 +0100 (BST)


On 12 Jul 1999, Peter Dalgaard BSA wrote:

> ripley@stats.ox.ac.uk writes:
> 
> > > Doug, could you 
> > > 
> > > - recompile eigen.o with -g 
> > > - run it with -d gdb
> > > - see if the crash resembles the one below
> > > 
> > > ...
> > > > eigen(cbind(-1,c(1:2,0),0:2)) # complex values
> > > 
> > > Program received signal SIGSEGV, Segmentation fault.
> > > hqr2_ (nm=0x1, n=0x70d99fc8, low=0xefffed94, igh=0xefffed90,
> > > h=0x70d99f80,
> > >     wr=0x70d99f68, wi=0x70d99f50, z=0x70d99f00, ierr=0x70d99ed0)
> > >     at eigen.f:2120
> > > 2120                   RA = RA + H(I,J) * H(J,NA)
> > > Current language:  auto; currently fortran
> > > (gdb) p j
> > > $1 = 0
> > > (gdb) p na
> > > $2 = 0
> > > (gdb) p en
> > > $3 = 3
> > > (gdb) p n
> > > $4 = (PTR TO -> ( integer )) 0x70d99fc8
> > > (gdb) p nm
> > > $5 = (PTR TO -> ( integer )) 0x1
> > 
> > Does for me (modulo addresses, of course).
> 
> Exactly the same?? (Doug's wasn't...)

Modulo addresses as I said, yes:

> eigen(cbind(-1,c(1:2,0),0:2)) # complex values

Program received signal SIGSEGV, Segmentation fault.
hqr2_ (nm=0x1, n=0xbe0a58, low=0xefffe414, igh=0xefffe410, h=0xbe0a10, 
    wr=0xbe09f8, wi=0xbe09e0, z=0xbe0990, ierr=0xbe0960) at eigen.f:2120
2120                   RA = RA + H(I,J) * H(J,NA)
Current language:  auto; currently fortran
(gdb) p j
$1 = 0
(gdb) p na
$2 = 0
(gdb) p en
$3 = 3
(gdb) p n
$4 = (PTR TO -> ( integer )) 0xbe0a58
(gdb) p nm
$5 = (PTR TO -> ( integer )) 0x1

which looks pretty much the same to me.


-- 
Brian D. Ripley,                  ripley@stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272860 (secr)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._