Add the <-- JVM NOT MANDATORY --> TAG
[scilab.git] / scilab / modules / sparse / tests / nonreg_tests / bug_3263.tst
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2007-2008 - INRIA - Serge STEER <serge.steer@inria.fr>
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7
8 // <-- JVM NOT MANDATORY -->
9
10 // <-- Non-regression test for bug 3263 -->
11 //
12 // <-- Bugzilla URL -->
13 // http://bugzilla.scilab.org/show_bug.cgi?id=3263
14 //
15 // <-- Short Description -->
16 // ludel(h) fails @ rev. 25935
17
18
19 a=[0.2,0.6,0.6,0.2,0.3;
20 0.8,0.8,0.7,0.2,0.9;
21 0,0.7,0.7,0.2,0.2;
22 0.3,0.9,0.2,0.9,0.3;
23 0.7,0.1,0.5,0.7,0.4];
24
25 b=[0.3;0.6;0.5;0.3;0.6];
26
27 A=sparse(a);
28
29 [h,rk]=lufact(A);
30
31 x=lusolve(h,b);
32
33 if norm(a*x-b)>1d-10 then pause,end
34
35 [P,L,U,Q]=luget(h);
36
37 if norm(P*L*U*Q-A)>1d-10 then pause,end
38
39 ierr = execstr('ludel(h);','errcatch');
40 if ierr <> 0 then pause,end;
41