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