regenerated umfpack test ref 49/15549/1
Anais AUBERT [Mon, 17 Nov 2014 10:17:00 +0000 (11:17 +0100)]
Change-Id: I54e59b466e2b5bf27567c2c4a2c90a076ddd9bd0
(cherry picked from commit 11925fcb345e9f56b9c54bc88c3669a60b58063b)

scilab/modules/umfpack/tests/unit_tests/umfpack.dia.ref
scilab/modules/umfpack/tests/unit_tests/umfpack.tst

index 4583f19..59f4f17 100644 (file)
@@ -18,10 +18,10 @@ assert_checkfalse(execstr("umfpack(speye(3,3), ""\"", %s)"   ,"errcatch") == 0);
 refMsg = msprintf(_("%s: Wrong type for input argument #%d: A matrix expected.\n"), "umfpack", 3);
 assert_checkerror("umfpack(speye(3,3), ""\"", %s)", refMsg);
 A = sparse( [ 2  3  0  0  0;
-              3  0  4  0  6; 
-              0 -1 -3  2  0; 
-              0  0  1  0  0; 
-              0  4  2  0  1] );
+3  0  4  0  6;
+0 -1 -3  2  0;
+0  0  1  0  0;
+0  4  2  0  1] );
 b = [8 ; 45; -3; 3; 19];
 x = umfpack(A,"\",b);
 assert_checkalmostequal(x, [1 2 3 4 5]');
@@ -39,48 +39,30 @@ x = umfpack(b,"/",A);
 assert_checkfalse(norm(x*A - b) > %eps);
 // solve a complex system
 A = sparse( [ 2+%i  3+2*%i  0      0    0;
-              3-%i  0       4+%i   0    6-3*%i; 
-              0    -1+%i   -3+6*%i 2-%i 0; 
-              0     0       1-5*%i 0    0; 
-              0     4       2-%i   0    1] );
+3-%i  0       4+%i   0    6-3*%i;
+0    -1+%i   -3+6*%i 2-%i 0;
+0     0       1-5*%i 0    0;
+0     4       2-%i   0    1] );
 b = [ 3+13*%i ; 58+32*%i ; -19+13*%i ; 18-12*%i ; 22+16*%i ];
-x = umfpack(A,"\",b)  // x must be [1+i; 2+2i; 3+3i; 4 + 4i; 5+5i]
- x  =
-    1. + i    
-    2. + 2.i  
-    3. + 3.i  
-    4. + 4.i  
-    5. + 5.i  
-assert_checkfalse(norm(x - [1+%i; 2+2*%i; 3+3*%i; 4 + 4*%i; 5+5*%i]) < %eps)
- ans  =
-  T  
+x = umfpack(A,"\",b);  // x must be [1+i; 2+2i; 3+3i; 4 + 4i; 5+5i]
+assert_checkfalse(norm(x - [1+%i; 2+2*%i; 3+3*%i; 4 + 4*%i; 5+5*%i]) < %eps);
 A = sparse( [ 2  3  0  0  0;
-              3  0  4  0  6; 
-              0 -1 -3  2  0; 
-              0  0  1  0  0; 
-              0  4  2  0  1] );
+3  0  4  0  6;
+0 -1 -3  2  0;
+0  0  1  0  0;
+0  4  2  0  1] );
 Lup = umf_lufact(A);
-[OK, nrow, ncol, lnz, unz, udiag_nz, it] = umf_luinfo(Lup);  // OK must be %t, nrow=ncol = 5, 
+[OK, nrow, ncol, lnz, unz, udiag_nz, it] = umf_luinfo(Lup);  // OK must be %t, nrow=ncol = 5,
 [L,U,p,q,R] = umf_luget(Lup);
-nnz(L)  // must be equal to lnz
- ans  =
-    9.  
 assert_checkequal(nnz(L), lnz);
-nnz(U)  // must be equal to unz
- ans  =
-    9.  
 assert_checkequal(nnz(U), unz);
 umf_ludel(Lup); // clear memory
 // this is the test matrix from UMFPACK
 A = sparse( [ 2  3  0  0  0;
-              3  0  4  0  6; 
-              0 -1 -3  2  0; 
-              0  0  1  0  0; 
-              0  4  2  0  1] );
+3  0  4  0  6;
+0 -1 -3  2  0;
+0  0  1  0  0;
+0  4  2  0  1] );
 Lup = umf_lufact(A);
 [L,U,p,q,R] = umf_luget(Lup);
 B = A;
@@ -90,10 +72,10 @@ assert_checkfalse(norm(B(p,q) - L*U) > %eps);
 umf_ludel(Lup);// clear memory
 // the same with a complex matrix
 A = sparse( [ 2+%i  3+2*%i  0      0    0;
-              3-%i  0       4+%i   0    6-3*%i; 
-              0    -1+%i   -3+6*%i 2-%i 0; 
-              0     0       1-5*%i 0    0; 
-              0     4       2-%i   0    1] );
+3-%i  0       4+%i   0    6-3*%i;
+0    -1+%i   -3+6*%i 2-%i 0;
+0     0       1-5*%i 0    0;
+0     4       2-%i   0    1] );
 Lup = umf_lufact(A);
 [L,U,p,q,R] = umf_luget(Lup);
 B = A;
index 351b70c..6f7a8ad 100644 (file)
@@ -24,10 +24,10 @@ refMsg = msprintf(_("%s: Wrong type for input argument #%d: A matrix expected.\n
 assert_checkerror("umfpack(speye(3,3), ""\"", %s)", refMsg);
 
 A = sparse( [ 2  3  0  0  0;
-              3  0  4  0  6; 
-              0 -1 -3  2  0; 
-              0  0  1  0  0; 
-              0  4  2  0  1] );
+3  0  4  0  6;
+0 -1 -3  2  0;
+0  0  1  0  0;
+0  4  2  0  1] );
 b = [8 ; 45; -3; 3; 19];
 x = umfpack(A,"\",b);
 
@@ -52,27 +52,24 @@ assert_checkfalse(norm(x*A - b) > %eps);
 
 // solve a complex system
 A = sparse( [ 2+%i  3+2*%i  0      0    0;
-              3-%i  0       4+%i   0    6-3*%i; 
-              0    -1+%i   -3+6*%i 2-%i 0; 
-              0     0       1-5*%i 0    0; 
-              0     4       2-%i   0    1] );
+3-%i  0       4+%i   0    6-3*%i;
+0    -1+%i   -3+6*%i 2-%i 0;
+0     0       1-5*%i 0    0;
+0     4       2-%i   0    1] );
 b = [ 3+13*%i ; 58+32*%i ; -19+13*%i ; 18-12*%i ; 22+16*%i ];
-x = umfpack(A,"\",b)  // x must be [1+i; 2+2i; 3+3i; 4 + 4i; 5+5i]
+x = umfpack(A,"\",b);  // x must be [1+i; 2+2i; 3+3i; 4 + 4i; 5+5i]
 
-assert_checkfalse(norm(x - [1+%i; 2+2*%i; 3+3*%i; 4 + 4*%i; 5+5*%i]) < %eps)
+assert_checkfalse(norm(x - [1+%i; 2+2*%i; 3+3*%i; 4 + 4*%i; 5+5*%i]) < %eps);
 
 A = sparse( [ 2  3  0  0  0;
-              3  0  4  0  6; 
-              0 -1 -3  2  0; 
-              0  0  1  0  0; 
-              0  4  2  0  1] );
+3  0  4  0  6;
+0 -1 -3  2  0;
+0  0  1  0  0;
+0  4  2  0  1] );
 Lup = umf_lufact(A);
-[OK, nrow, ncol, lnz, unz, udiag_nz, it] = umf_luinfo(Lup);  // OK must be %t, nrow=ncol = 5, 
+[OK, nrow, ncol, lnz, unz, udiag_nz, it] = umf_luinfo(Lup);  // OK must be %t, nrow=ncol = 5,
 [L,U,p,q,R] = umf_luget(Lup);
-nnz(L)  // must be equal to lnz
 assert_checkequal(nnz(L), lnz);
-
-nnz(U)  // must be equal to unz
 assert_checkequal(nnz(U), unz);
 
 umf_ludel(Lup); // clear memory
@@ -81,10 +78,10 @@ umf_ludel(Lup); // clear memory
 
 // this is the test matrix from UMFPACK
 A = sparse( [ 2  3  0  0  0;
-              3  0  4  0  6; 
-              0 -1 -3  2  0; 
-              0  0  1  0  0; 
-              0  4  2  0  1] );
+3  0  4  0  6;
+0 -1 -3  2  0;
+0  0  1  0  0;
+0  4  2  0  1] );
 Lup = umf_lufact(A);
 [L,U,p,q,R] = umf_luget(Lup);
 B = A;
@@ -96,10 +93,10 @@ umf_ludel(Lup);// clear memory
 
 // the same with a complex matrix
 A = sparse( [ 2+%i  3+2*%i  0      0    0;
-              3-%i  0       4+%i   0    6-3*%i; 
-              0    -1+%i   -3+6*%i 2-%i 0; 
-              0     0       1-5*%i 0    0; 
-              0     4       2-%i   0    1] );
+3-%i  0       4+%i   0    6-3*%i;
+0    -1+%i   -3+6*%i 2-%i 0;
+0     0       1-5*%i 0    0;
+0     4       2-%i   0    1] );
 Lup = umf_lufact(A);
 [L,U,p,q,R] = umf_luget(Lup);
 B = A;