bug 3025 fix
[scilab.git] / scilab / modules / sparse / tests / nonreg_tests / bug_3025.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 2277 -->
9 //
10 // <-- Bugzilla URL -->
11 // http://bugzilla.scilab.org/show_bug.cgi?id=2277
12 //
13 // <-- Short Description -->
14 //lusolve produces erroneous results with Scilab5 (it was ok on Scilab4)
15
16 A=sparse([2,2;2,4;3,3;3,7;4,2;4,4;6,6;6,8;7,3;7,7;8,6;8,8],..
17       [-1;1;-1;1;1;-1;-1;1;1;-1;1;-1],[10,10]) ;
18 [ptr,rkA]=lufact(A);
19 x=lusolve(ptr,zeros(10,1));ludel(ptr);
20 if norm(A*x)>1d-14 then pause,end
21
22