Bug #10178: norm returned error message: 'division by zero' for some sparse matrices.
[scilab.git] / scilab / modules / linear_algebra / tests / nonreg_tests / bug_10178.dia.ref
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2012 - Scilab Enterprises - Adeline CARNIS
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7 // <-- CLI SHELL MODE -->
8 // <-- Non-regression test for bug 10178 -->
9 //
10 // <-- Bugzilla URL -->
11 // http://bugzilla.scilab.org/show_bug.cgi?id=10178
12 //
13 // <-- Short Description -->
14 //    norm function falied for some sparse matrices.
15 // =============================================================================
16 A = [1 0 0 1 0;0 -1 -1 0 -1];
17 AS = sparse(A);
18 assert_checkalmostequal(norm(A), norm(AS));
19 assert_checkalmostequal(norm(A'), norm(AS'));
20 assert_checkequal(norm(A, 1), norm(AS, 1));
21 assert_checkequal(norm(A, %inf), norm(AS, %inf));
22 assert_checkequal(norm(A, 'fro'), norm(AS, 'fro'));
23 A = [1 1 1 1 1;-1 -1 -1 -1 -1];
24 AS = sparse(A);
25 assert_checkequal(norm(A), norm(AS));
26 assert_checkequal(norm(A)', norm(AS)');
27 assert_checkequal(norm(A, 1), norm(AS, 1));
28 assert_checkequal(norm(A, %inf), norm(AS, %inf));
29 assert_checkequal(norm(A, 'fro'), norm(AS, 'fro'));
30 A = [
31     0.    0.    0.    0.    0.  
32     0.    0.    0.    0.    0.  
33     0.    0.   -1.    0.    1.  
34     0.    0.    0.    0.    0.  
35     0.    0.    1.    0.   -1.  ];
36 AS = sparse(A);
37 assert_checkequal(norm(A), norm(AS));
38 assert_checkequal(norm(A, 1), norm(AS, 1));
39 assert_checkequal(norm(A, %inf), norm(AS, %inf));
40 assert_checkequal(norm(A, 'fro'), norm(AS, 'fro'));