add ";" in pcg_list test
[scilab.git] / scilab / modules / sparse / tests / unit_tests / pcg_list.dia.ref
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2008 - INRIA - Michael Baudin
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7 // <-- CLI SHELL MODE -->
8 // Numerical tests
9 // Case where A is given as a function computing the right-hand side
10 mymatrix=[ 94  0   0   0    0   28  0   0   32  0
11 0   59  13  5    0   0   0   10  0   0
12 0   13  72  34   2   0   0   0   0   65
13 0   5   34  114  0   0   0   0   0   55
14 0   0   2   0    70  0   28  32  12  0
15 28  0   0   0    0   87  20  0   33  0
16 0   0   0   0    28  20  71  39  0   0
17 0   10  0   0    32  0   39  46  8   0
18 32  0   0   0    12  33  0   8   82  11
19 0   0   65  55   0   0   0   0   11  100];
20 b = [154.
21 87.
22 186.
23 208.
24 144.
25 168.
26 158.
27 135.
28 178.
29 231.];
30 function y=Atimesx(x,mymatrix)
31   y=mymatrix*x
32 endfunction
33 // With the default 10 iterations, the algorithm performs well
34 Alist = list(Atimesx,mymatrix);
35 [xcomputed, fail, err, iter, res]=pcg(Alist,b);
36 xexpected=ones(10,1);
37 if norm(xcomputed-xexpected)>10**3*%eps then bugmes();quit;end
38 if fail<>0 then bugmes();quit;end
39 if iter<>10 then bugmes();quit;end
40 if err > 10**3*%eps then bugmes();quit;end