* Bug #13058 fixed - Arpack-ng library has been updated and fixes this bug.
[scilab.git] / scilab / modules / arnoldi / tests / nonreg_tests / bug_13058.tst
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2014 - Scilab Enterprises - Adeline CARNIS
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7
8 // <-- CLI SHELL MODE -->
9
10 // <-- Non-regression test for bug 13058 -->
11 //
12 // <-- Bugzilla URL -->
13 // http://bugzilla.scilab.org/show_bug.cgi?id=13058
14 //
15 // <-- Short Description -->
16 //    The release 3.1.4 of arpack ng library fixes this bug.
17 // =============================================================================
18
19 A = [2.486880466472303  2.008746355685131  1.655976676384840  1.571428571428571  1.469387755102041  1.635568513119533  1.927113702623907
20 2.008746355685131  2.364431486880466  1.967930029154519  1.696793002915452  1.551020408163265  1.530612244897959  1.635568513119533
21 1.655976676384840  1.967930029154519  2.405247813411079  1.947521865889213  1.758017492711370  1.551020408163265  1.469387755102041
22 1.571428571428571  1.696793002915452  1.947521865889213  2.323615160349854  1.947521865889213  1.696793002915452  1.571428571428571
23 1.469387755102041  1.551020408163265  1.758017492711370  1.947521865889213  2.405247813411079  1.967930029154519  1.655976676384840
24 1.635568513119533  1.530612244897959  1.551020408163265  1.696793002915452  1.967930029154519  2.364431486880466  2.008746355685131
25 1.927113702623907  1.635568513119533  1.469387755102041  1.571428571428571  1.655976676384840  2.008746355685131  2.486880466472303];
26
27 d0 = spec(A);
28
29 B = eye(7, 7);
30 d1 = eigs(A, B, 4, "BE");
31 assert_checkalmostequal(d1, [d0(1:2); d0($-1:$)], 1.e-12);