[Octave-bug-tracker] [bug #31479] Crash & bugs in eigs

[John W. Eaton]

Follow-up Comment #16, bug #31479 (project octave):

I checked again, and running


A = [1,0,0,-1;0,1,0,0;0,0,1,0;0,0,2,1]; 
for i = 1:1000
  i
  [vecs,vals,f] = eigs(A,1);
endfor


In my copy of Octave 3.3.90, which is built with ARPACK from the Debian
package, I see a crash after a few calls:


i =  1
i =  2
i =  3
i =  4
i =  5
i =  6
i =  7
*** glibc detected *** /usr/local/octave/3.3.90/bin/octave: free(): invalid
next size (fast): 0x00000000020b6760 ***
======= Backtrace: =========
/lib/libc.so.6(+0x71ad6)[0x7fba23d40ad6]
/lib/libc.so.6(cfree+0x6c)[0x7fba23d4584c]
/usr/local/octave/3.3.90/lib/octave-3.3.90/liboctinterp-3.3.90.so(_ZN5ArrayIdED2Ev+0x65)[0x7fba2a94e405]
/usr/local/octave/3.3.90/libexec/octave/3.3.90/oct/x86_64-unknown-linux-gnu/eigs.oct(_Z26EigsRealNonSymmetricMatrixI6MatrixEiRKT_SsiiRiR13ComplexMatrixR19ComplexColumnVectorS3_R12ColumnVectorSA_RSodbbii+0x156c)[0x7fba1ea9247c]
/usr/local/octave/3.3.90/libexec/octave/3.3.90/oct/x86_64-unknown-linux-gnu/eigs.oct(_Z5FeigsRK17octave_value_listi+0x4106)[0x7fba1ea89046]
/usr/local/octave/3.3.90/lib/octave-3.3.90/liboctinterp-3.3.90.so(_ZN14octave_builtin17do_multi_index_opEiRK17octave_value_listPKSt4listI13octave_lvalueSaIS4_EE+0x231)[0x7fba2ac799a1]
/usr/local/octave/3.3.90/lib/octave-3.3.90/liboctinterp-3.3.90.so(_ZN14octave_builtin7subsrefERKSsRKSt4listI17octave_value_listSaIS3_EEiPKS2_I13octave_lvalueSaIS8_EE+0x17c)[0x7fba2ac7896c]
/usr/local/octave/3.3.90/lib/octave-3.3.90/liboctinterp-3.3.90.so(_ZN12octave_value7subsrefERKSsRKSt4listI17octave_value_listSaIS3_EEiPKS2_I13octave_lvalueSaIS8_EE+0x16)[0x7fba2ad58916]
/usr/local/octave/3.3.90/lib/octave-3.3.90/liboctinterp-3.3.90.so(_ZN21tree_index_expression6rvalueEiPKSt4listI13octave_lvalueSaIS1_EE+0x2f9)[0x7fba2adf2209]
/usr/local/octave/3.3.90/lib/octave-3.3.90/liboctinterp-3.3.90.so(_ZN21tree_multi_assignment6rvalueEi+0x1d5)[0x7fba2add4545]
/usr/local/octave/3.3.90/lib/octave-3.3.90/liboctinterp-3.3.90.so(_ZN21tree_multi_assignment7rvalue1Ei+0x5e)[0x7fba2add40be]
/usr/local/octave/3.3.90/lib/octave-3.3.90/liboctinterp-3.3.90.so(_ZN14tree_evaluator15visit_statementER14tree_statement+0xfa)[0x7fba2ade389a]
/usr/local/octave/3.3.90/lib/octave-3.3.90/liboctinterp-3.3.90.so(_ZN14tree_evaluator20visit_statement_listER19tree_statement_list+0x58)[0x7fba2ade29e8]
/usr/local/octave/3.3.90/lib/octave-3.3.90/liboctinterp-3.3.90.so(_ZN14tree_evaluator24visit_simple_for_commandER23tree_simple_for_command+0x2ba)[0x7fba2ade3e9a]
/usr/local/octave/3.3.90/lib/octave-3.3.90/liboctinterp-3.3.90.so(_ZN14tree_evaluator15visit_statementER14tree_statement+0x58)[0x7fba2ade37f8]
/usr/local/octave/3.3.90/lib/octave-3.3.90/liboctinterp-3.3.90.so(_ZN14tree_evaluator20visit_statement_listER19tree_statement_list+0x58)[0x7fba2ade29e8]
/usr/local/octave/3.3.90/lib/octave-3.3.90/liboctinterp-3.3.90.so(_Z9main_loopv+0x14f)[0x7fba2ac1540f]
/usr/local/octave/3.3.90/lib/octave-3.3.90/liboctinterp-3.3.90.so(octave_main+0x12d5)[0x7fba2abc06d5]
/lib/libc.so.6(__libc_start_main+0xfd)[0x7fba23cedc4d]
/usr/local/octave/3.3.90/bin/octave[0x400769]
======= Memory map: ========
00400000-00401000 r-xp 00000000 fe:04 1766035                           
/usr/local/octave/3.3.90/bin/octave-3.3.90
00600000-00601000 rw-p 00000000 fe:04 1766035                           
/usr/local/octave/3.3.90/bin/octave-3.3.90
01b49000-02219000 rw-p 00000000 00:00 0                                 
[heap]
7fba18000000-7fba18021000 rw-p 00000000 00:00 0 
7fba18021000-7fba1c000000 ---p 00000000 00:00 0 
7fba1ea6c000-7fba1eabf000 r-xp 00000000 fe:04 1774291                   
/usr/local/octave/3.3.90/libexec/octave/3.3.90/oct/x86_64-unknown-linux-gnu/eigs.oct
7fba1eabf000-7fba1ecbe000 ---p 00053000 fe:04 1774291                   
/usr/local/octave/3.3.90/libexec/octave/3.3.90/oct/x86_64-unknown-linux-gnu/eigs.oct
7fba1ecbe000-7fba1ecc0000 rw-p 00052000 fe:04 1774291                   
/usr/local/octave/3.3.90/libexec/octave/3.3.90/oct/x86_64-unknown-linux-gnu/eigs.oct
7fba1ecc0000-7fba1ecc1000 rw-p 00000000 00:00 0 
7fba1ecc1000-7fba1ecdd000 r-xp 00000000 fe:04 1774327                   
/usr/local/octave/3.3.90/libexec/octave/3.3.90/oct/x86_64-unknown-linux-gnu/strfind.oct
7fba1ecdd000-7fba1eedd000 ---p 0001c000 fe:04 1774327                   
/usr/local/octave/3.3.90/libexec/octave/3.3.90/oct/x86_64-unknown-linux-gnu/strfind.oct
7fba1eedd000-7fba1eede000 rw-p 0001c000 fe:04 1774327                   
/usr/local/octave/3.3.90/libexec/octave/3.3.90/oct/x86_64-unknown-linux-gnu/strfind.oct
7fba1eede000-7fba20ee0000 rw-p 00000000 00:00 0 
7fba20ee0000-7fba20f36000 r-xp 00000000 fe:04 1774277                   
/usr/local/octave/3.3.90/libexec/octave/3.3.90/oct/x86_64-unknown-linux-gnu/cellfun.oct
7fba20f36000-7fba21135000 ---p 00056000 fe:04 1774277                   
/usr/local/octave/3.3.90/libexec/octave/3.3.90/oct/x86_64-unknown-linux-gnu/cellfun.oct
7fba21135000-7fba21137000 rw-p 00055000 fe:04 1774277                   
/usr/local/octave/3.3.90/libexec/octave/3.3.90/oct/x86_64-unknown-linux-gnu/cellfun.oct
7fba21137000-7fba21138000 rw-p 00000000 00:00 0 
7fba21138000-7fba21164000 r-xp 00000000 fe:04 1774297                   
/usr/local/octave/3.3.90/libexec/octave/3.3.90/oct/x86_64-unknown-linux-gnu/find.oct
7fba21164000-7fba21364000 ---p 0002c000 fe:04 1774297                   
/usr/local/octave/3.3.90/libexec/octave/3.3.90/oct/x86_64-unknown-linux-gnu/find.oct
7fba21364000-7fba21365000 rw-p 0002c000 fe:04 1774297                   
/usr/local/octave/3.3.90/libexec/octave/3.3.90/oct/x86_64-unknown-linux-gnu/find.oct
7fba21365000-7fba21366000 rw-p 00000000 00:00 0 
7fba21366000-7fba214db000 r--p 00000000 fe:04 721078                    
/usr/lib/locale/locale-archive
7fba214db000-7fba214e6000 r-xp 00000000 fe:00 350411                    
/lib/libnss_files-2.11.2.so
7fba214e6000-7fba216e5000 ---p 0000b000 fe:00 350411                    
/lib/libnss_files-2.11.2.so
7fba216e5000-7fba216e6000 r--p 0000a000 fe:00 350411                    
/lib/libnss_files-2.11.2.so
7fba216e6000-7fba216e7000 rw-p 0000b000 fe:00 350411                    
/lib/libnss_files-2.11.2.so
7fba216e7000-7fba216f1000 r-xp 00000000 fe:00 350408                    
/lib/libnss_nis-2.11.2.so
7fba216f1000-7fba218f0000 ---p 0000a000 fe:00 350408                    
/lib/libnss_nis-2.11.2.so
7fba218f0000-7fba218f1000 r--p 00009000 fe:00 350408                    
/lib/libnss_nis-2.11.2.so
7fba218f1000-7fba218f2000 rw-p 0000a000 fe:00 350408                    
/lib/libnss_nis-2.11.2.so
7fba218f2000-7fba21907000 r-xp 00000000 fe:00 350424                    
/lib/libnsl-2.11.2.so
7fba21907000-7fba21b06000 ---p 00015000 fe:00 350424                    
/lib/libnsl-2.11.2.so
7fba21b06000-7fba21b07000 r--p 00014000 fe:00 350424                    
/lib/libnsl-2.11.2.so
7fba21b07000-7fba21b08000 rw-p 00015000 fe:00 350424                    
/lib/libnsl-2.11.2.so
7fba21b08000-7fba21b0a000 rw-p 00000000 00:00 0 
7fba21b0a000-7fba21b11000 r-xp 00000000 fe:00 350421                    
/lib/libnss_compat-2.11.2.so
7fba21b11000-7fba21d10000 ---p 00007000 fe:00 350421                    
/lib/libnss_compat-2.11.2.so
7fba21d10000-7fba21d11000 r--p 00006000 fe:00 350421                    
/lib/libnss_compat-2.11.2.so
7fba21d11000-7fba21d12000 rw-p 00007000 fe:00 350421                    
/lib/libnss_compat-2.11.2.so
7fba21d12000-7fba21d17000 r-xp 00000000 fe:04 720379                    
/usr/lib/libXdmcp.so.6.0.0
7fba21d17000-7fba21f16000 ---p 00005000 fe:04 720379                    
/usr/lib/libXdmcp.so.6.0.0
7fba21f16000-7fba21f17000 rw-p 00004000 fe:04 720379                    
/usr/lib/libXdmcp.so.6.0.0
7fba21f17000-7fba21f19000 r-xp 00000000 fe:04 720377                    
/usr/lib/libXau.so.6.0.0
7fba21f19000-7fba22119000 ---p 00002000 fe:04 720377                    
/usr/lib/libXau.so.6.0.0
7fba22119000-7fba2211a000 rw-p 00002000 fe:04 720377                    
/usr/lib/libXau.so.6.0.0
7fba2211a000-7fba22121000 r-xp 00000000 fe:00 350416                    
/lib/librt-2.11.2.so
7fba22121000-7fba22320000 ---p 00007000 fe:00 350416                    
/lib/librt-2.11.2.sopanic: Aborted -- stopping myself...
attempting to save variables to `octave-core'...
save to `octave-core' complete
Aborted


With a copy of Octave built from the current sources (ARPACK sources moved to
Octave), I don't see the crash, but I do sometimes (every few thousand calls?)
see


error: eigs: error -8 in dneupd


Perhaps this is just because the algorithm is iterative and picks a starting
point using some random numbers?

I don't see either of the above errors in 10,000 calls if I set


A = rand (300);


so this seems somehow dependent on the small size of the test matrix in the
original bug report, which Rik has noted is not the type of matrix that eigs
was really designed for.

So it would be great if someone wants to try to debug the "error -8" problem,
but I don't think we should worry about this for the release.  The more
important thing for me is that the crash is gone and I'm no longer seeing
complaints from valgrind about uninitialized values.  Given that, I'd like to
close this report.  OK?

    _______________________________________________________

Reply to this item at:

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

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