octave-maintainers
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Small ARPACK problem


From: Robinson, Melvin D
Subject: Small ARPACK problem
Date: Sat, 7 Apr 2012 21:40:47 +0000

I am building Octave on another (older) system and am running into a small problem. 

make check yields:

  src/DLD-FUNCTIONS/eig.cc ............................... PASS   20/20
  src/DLD-FUNCTIONS/eigs.cc .............................. ** On entry to DLASCL parameter number  4 had an illegal value
make[1]: *** [check] Segmentation fault (core dumped)
make[1]: Leaving directory `/home/melvin/octave/test'
make: *** [check] Error 2

but from the command prompt gives:

octave:1> test eigs verbose
>>>>> /home/melvin/octave/src/DLD-FUNCTIONS/eigs.cc
  ***** shared n, k, A, d0, d2
 n = 20;
 k = 4;
 A = sparse ([3:n,1:n,1:(n-2)],[1:(n-2),1:n,3:n],[ones(1,n-2),4*ones(1,n),ones(1,n-2)]);
 d0 = eig (A);
 d2 = sort (d0);
 [~, idx] = sort (abs (d0));
 d0 = d0(idx);
 rand ("state", 42); # initialize generator to make eigs behavior reproducible
  ***** testif HAVE_ARPACK
 d1 = eigs (A, k);
 assert (d1, d0(end:-1:(end-k+1)), 1e-11);
!!!!! test failed
assert (d1,d0 (__end__:-1:(__end__ - k + 1)),1e-11) expected
   5.9190
   5.9190
   5.6825
   5.6825
but got
   5.9190
   5.9190
   5.6825
   5.6806
maximum absolute error 0.0019438 exceeds tolerance 1e-11

 

Now I know that I probably need to rebuild ARPACK again, but shouldn't Octave segfault from the command prompt too?

 

 


reply via email to

[Prev in Thread] Current Thread [Next in Thread]