Set svn:eol-style to native
[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 // <-- Non-regression test for bug 3263 -->
9 //
10 // <-- Bugzilla URL -->
11 // http://bugzilla.scilab.org/show_bug.cgi?id=3263
12 //
13 // <-- Short Description -->
14 // ludel(h) fails @ rev. 25935
15
16
17 a=[0.2,0.6,0.6,0.2,0.3;
18 0.8,0.8,0.7,0.2,0.9;
19 0,0.7,0.7,0.2,0.2;
20 0.3,0.9,0.2,0.9,0.3;
21 0.7,0.1,0.5,0.7,0.4];
22
23 b=[0.3;0.6;0.5;0.3;0.6];
24
25 A=sparse(a);
26
27 [h,rk]=lufact(A);
28
29 x=lusolve(h,b);
30
31 if norm(a*x-b)>1d-10 then pause,end
32
33 [P,L,U,Q]=luget(h);
34
35 if norm(P*L*U*Q-A)>1d-10 then pause,end
36
37 ierr = execstr('ludel(h);','errcatch');
38 if ierr <> 0 then pause,end;
39