Add the <-- JVM NOT MANDATORY --> TAG
[scilab.git] / scilab / modules / sparse / tests / unit_tests / spisp.tst
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) ????-2008 - INRIA
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7
8 // <-- JVM NOT MANDATORY -->
9
10 m = 200;
11 n = 400;
12 mi = 30;
13 ni = 40;
14
15 A = sprand(n,m,0.1) + %i*sprand(n,m,0.1);
16 B = sprand(mi,ni,0.7) +%i*sprand(mi,ni,0.7);
17
18 ii = ["grand(mi,1,""uin"",1,m)";
19       "-sort(-grand(mi,1,""uin"",1,m))";
20       "1:mi";
21       "m:-1:(m-mi+1)";
22       "m-5:m+mi-6"];
23 jj = ["grand(ni,1,""uin"",1,n)";
24       "-sort(-grand(ni,1,""uin"",1,n))";
25       "1:ni";
26       "n:-1:(n-ni+1)";
27       "n-2:n+ni-3"];
28
29
30 switchA = ["real(A)" "A"];
31 switchB = ["real(B)" "B"];
32
33 for i=1:5
34    indi = evstr(ii(i));
35    for j=1:5
36       indj = evstr(jj(j));
37       for k=1:2
38          for l=1:2
39             AA = evstr(switchA(k));
40             Af = full(AA);
41             D = evstr(switchB(l));
42             Df = full(D);
43             AA(indi,indj) = D;
44             Af(indi,indj) = Df;      
45             if max(abs(Af - full(AA))) > 0 then
46                mprintf('\n\r test (%d,%d,%d,%d) fails',i,j,k,l)
47             end
48          end
49       end
50    end
51 end