[Octave-bug-tracker] [bug #29553] missing arpack: make check segfaults on eigs

[Petr Mikulik]

Follow-up Comment #11, bug #29553 (project octave):

1. I went through EIG::init(). Tracing and printing:

161       Array<double> wi (n, 1);
(gdb)
162       double *pwi = wi.fortran_vec ();
(gdb) p pwi
$24 = (double *) 0xb7fdf880
(gdb) p *pwi
$25 = 1.5251834411709584e-314
(gdb) p wi.rep
$26 = (class Array<double>::ArrayRep *) 0x83f6670
(gdb) next
...
174       F77_XFCN (dgeev, DGEEV, (F77_CONST_CHAR_ARG2 ("N", 1),
(gdb) p wi.rep
$27 = (class Array<double>::ArrayRep *) 0x83f6670
(gdb) 
...
185           double *pwork = work.fortran_vec ();
(gdb)
187           F77_XFCN (dgeev, DGEEV, (F77_CONST_CHAR_ARG2 ("N", 1),
(gdb) p wi.rep
$31 = (class Array<double>::ArrayRep *) 0x40580000
...
208           v.resize (nvr, nvr);
(gdb)
210           for (octave_idx_type j = 0; j < n; j++)
(gdb)
212               if (wi.elem (j) == 0.0)
(gdb)

Program received signal SIGSEGV, Segmentation fault.
0xb74de425 in Array<double>::make_unique (this=0xbfffc52c) at
../liboctave/Array.h:111
111           if (rep->count > 1)
(gdb) 


2. Blas and Lapack:
$ rpm -qa | grep blas
blas-3.2.0-2.4
libblas3-3.2.0-2.4
blas-man-3.2.0-2.4

$ rpm -qa | grep lapack
liblapack3-3.2.0-2.4
lapack-man-3.2.0-2.4
lapack-3.2.0-2.4



    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?29553>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/