Arnoldi tests should be available in nwni mode 39/15239/2
Bruno JOFRET [Thu, 18 Sep 2014 10:26:54 +0000 (12:26 +0200)]
Change-Id: I6c5f4a19f8a5c18a3e0bc0c579bbba1a22ba6322

14 files changed:
scilab/modules/arnoldi/tests/nonreg_tests/bug_dsaupd.dia.ref
scilab/modules/arnoldi/tests/nonreg_tests/bug_dsaupd.tst
scilab/modules/arnoldi/tests/unit_tests/dnaupd.dia.ref
scilab/modules/arnoldi/tests/unit_tests/dnaupd.tst
scilab/modules/arnoldi/tests/unit_tests/dneupd.dia.ref
scilab/modules/arnoldi/tests/unit_tests/dneupd.tst
scilab/modules/arnoldi/tests/unit_tests/dsaupd.dia.ref
scilab/modules/arnoldi/tests/unit_tests/dsaupd.tst
scilab/modules/arnoldi/tests/unit_tests/dseupd.dia.ref
scilab/modules/arnoldi/tests/unit_tests/dseupd.tst
scilab/modules/arnoldi/tests/unit_tests/znaupd.dia.ref
scilab/modules/arnoldi/tests/unit_tests/znaupd.tst
scilab/modules/arnoldi/tests/unit_tests/zneupd.dia.ref
scilab/modules/arnoldi/tests/unit_tests/zneupd.tst

index 5009b98..2d22910 100644 (file)
@@ -4,6 +4,7 @@
 //
 //  This file is distributed under the same license as the Scilab package.
 // =============================================================================
+// <-- CLI SHELL MODE -->
 // Tests for function dsaupd (added after a bug in the gateway: PutLhsVar missing)
 stacksize(300000);
 N = 1000;
index e62b021..7403fad 100644 (file)
@@ -5,6 +5,8 @@
 //  This file is distributed under the same license as the Scilab package.
 // =============================================================================
 
+// <-- CLI SHELL MODE -->
+
 // Tests for function dsaupd (added after a bug in the gateway: PutLhsVar missing)
 
 stacksize(300000);
index 5addde5..1ee3d1d 100644 (file)
@@ -4,6 +4,7 @@
 //
 //  This file is distributed under the same license as the Scilab package.
 // =============================================================================
+// <-- CLI SHELL MODE -->
 nx    = 10;
 nev   = 3;
 ncv   = 6;
@@ -46,7 +47,7 @@ while(ido <> 99)
         printf('Check the documentation of dnaupd\n\n');
     end
     if(ido == -1 | ido == 1)
-        // Perform matrix vector multiplication 
+        // Perform matrix vector multiplication
         workd(ipntr(2):ipntr(2) + nx -1) = A * workd(ipntr(1):ipntr(1) + nx - 1);
     end
 end
index b78a2ed..b964e28 100644 (file)
@@ -5,6 +5,8 @@
 //  This file is distributed under the same license as the Scilab package.
 // =============================================================================
 
+// <-- CLI SHELL MODE -->
+
 nx    = 10;
 
 nev   = 3;
@@ -61,7 +63,7 @@ while(ido <> 99)
     end
 
     if(ido == -1 | ido == 1)
-        // Perform matrix vector multiplication 
+        // Perform matrix vector multiplication
         workd(ipntr(2):ipntr(2) + nx -1) = A * workd(ipntr(1):ipntr(1) + nx - 1);
     end
 end
@@ -88,4 +90,3 @@ z(:,[c1, c2]) = [z(:,c1) + z(:,c2) * %i z(:,c1) - z(:,c2) * %i];
 z(:,$) = [];
 
 assert_checkalmostequal(A * z, z * d, sqrt(%eps), 1.e-10);
-
index bc265f2..8a58b92 100644 (file)
@@ -4,6 +4,7 @@
 //
 //  This file is distributed under the same license as the Scilab package.
 // =============================================================================
+// <-- CLI SHELL MODE -->
 // The following sets dimensions for this problem.
 nx    = 10;
 nev   = 3;
@@ -47,7 +48,7 @@ while(ido <> 99)
     printf('Check the documentation of dnaupd\n\n');
   end
   if(ido == -1 | ido == 1)
-    // Perform matrix vector multiplication 
+    // Perform matrix vector multiplication
     workd(ipntr(2):ipntr(2) + nx -1) = A * workd(ipntr(1):ipntr(1) + nx - 1);
   end
 end
@@ -68,4 +69,4 @@ if(modulo(nev + 1, 2) == 1)
 end
 z(:,[c1, c2]) = [z(:,c1) + z(:,c2) * %i z(:,c1) - z(:,c2) * %i];
 z(:,$) = [];
-assert_checkalmostequal(A * z, z * d, sqrt(%eps), 1.e-10); 
+assert_checkalmostequal(A * z, z * d, sqrt(%eps), 1.e-10);
index 7aee6a8..18abaaf 100644 (file)
@@ -5,6 +5,8 @@
 //  This file is distributed under the same license as the Scilab package.
 // =============================================================================
 
+// <-- CLI SHELL MODE -->
+
 // The following sets dimensions for this problem.
 
 nx    = 10;
@@ -56,14 +58,14 @@ while(ido <> 99)
   // either convergence is indicated or maxitr has been exceeded.
 
   [ido, resid, v, iparam, ipntr, workd, workl, info_dnaupd] = dnaupd(ido, bmat, nx, which, nev, tol, resid, ncv, v, iparam, ipntr, workd, workl, info_dnaupd);
-  
+
   if(info_dnaupd < 0)
     printf('\nError with dnaupd, info = %d\n',info_dnaupd);
     printf('Check the documentation of dnaupd\n\n');
   end
-  
+
   if(ido == -1 | ido == 1)
-    // Perform matrix vector multiplication 
+    // Perform matrix vector multiplication
     workd(ipntr(2):ipntr(2) + nx -1) = A * workd(ipntr(1):ipntr(1) + nx - 1);
   end
 end
@@ -89,5 +91,4 @@ end
 z(:,[c1, c2]) = [z(:,c1) + z(:,c2) * %i z(:,c1) - z(:,c2) * %i];
 z(:,$) = [];
 
-assert_checkalmostequal(A * z, z * d, sqrt(%eps), 1.e-10); 
-
+assert_checkalmostequal(A * z, z * d, sqrt(%eps), 1.e-10);
index 232574a..d9685ea 100644 (file)
@@ -4,6 +4,7 @@
 //
 //  This file is distributed under the same license as the Scilab package.
 // =============================================================================
+// <-- CLI SHELL MODE -->
 nx    = 10;
 nev   = 3;
 ncv   = 6;
@@ -15,9 +16,9 @@ ipntr   = zeros(14, 1);
 _select = zeros(ncv, 1);
 d       = zeros(nev, 1);
 z       = zeros(nx, nev);
-resid   = zeros(nx, 1); 
+resid   = zeros(nx, 1);
 v       = zeros(nx, ncv);
-workd   = zeros(3 * nx, 1); 
+workd   = zeros(3 * nx, 1);
 workl   = zeros(ncv * ncv + 8 * ncv, 1);
 // Build the symmetric test matrix
 A            = diag(10 * ones(nx,1));
@@ -43,7 +44,7 @@ while(ido <> 99)
     printf('Check the documentation of dsaupd\n\n');
   end
   if(ido == -1 | ido == 1)
-    // Perform matrix vector multiplication 
+    // Perform matrix vector multiplication
     workd(ipntr(2):ipntr(2) + nx - 1) = A * workd(ipntr(1):ipntr(1) + nx - 1);
   end
 end
@@ -57,4 +58,4 @@ if(info_dseupd < 0)
   printf('\nError with dseupd, info = %d\n', info_dseupd);
   printf('Check the documentation of dseupd.\n\n');
 end
-assert_checkalmostequal(A * z, z * diag(d), sqrt(%eps), 1.e-10); 
+assert_checkalmostequal(A * z, z * diag(d), sqrt(%eps), 1.e-10);
index e475269..73653b0 100644 (file)
@@ -5,6 +5,8 @@
 //  This file is distributed under the same license as the Scilab package.
 // =============================================================================
 
+// <-- CLI SHELL MODE -->
+
 nx    = 10;
 
 nev   = 3;
@@ -19,9 +21,9 @@ ipntr   = zeros(14, 1);
 _select = zeros(ncv, 1);
 d       = zeros(nev, 1);
 z       = zeros(nx, nev);
-resid   = zeros(nx, 1); 
+resid   = zeros(nx, 1);
 v       = zeros(nx, ncv);
-workd   = zeros(3 * nx, 1); 
+workd   = zeros(3 * nx, 1);
 workl   = zeros(ncv * ncv + 8 * ncv, 1);
 
 // Build the symmetric test matrix
@@ -51,14 +53,14 @@ while(ido <> 99)
   // either convergence is indicated or maxitr has been exceeded.
 
   [ido, resid, v, iparam, ipntr, workd, workl, info_dsaupd] = dsaupd(ido, bmat, nx, which, nev, tol, resid, ncv, v, iparam, ipntr, workd, workl, info_dsaupd);
-  
+
   if(info_dsaupd < 0)
     printf('\nError with dsaupd, info = %d\n',info_dsaupd);
     printf('Check the documentation of dsaupd\n\n');
   end
-  
+
   if(ido == -1 | ido == 1)
-    // Perform matrix vector multiplication 
+    // Perform matrix vector multiplication
     workd(ipntr(2):ipntr(2) + nx - 1) = A * workd(ipntr(1):ipntr(1) + nx - 1);
   end
 end
@@ -77,4 +79,4 @@ if(info_dseupd < 0)
 end
 
 
-assert_checkalmostequal(A * z, z * diag(d), sqrt(%eps), 1.e-10); 
+assert_checkalmostequal(A * z, z * diag(d), sqrt(%eps), 1.e-10);
index 2401b22..66890df 100644 (file)
@@ -4,6 +4,7 @@
 //
 //  This file is distributed under the same license as the Scilab package.
 // =============================================================================
+// <-- CLI SHELL MODE -->
 nx    = 10;
 nev   = 3;
 ncv   = 6;
@@ -15,9 +16,9 @@ ipntr   = zeros(14, 1);
 _select = zeros(ncv, 1);
 d       = zeros(nev, 1);
 z       = zeros(nx, nev);
-resid   = zeros(nx, 1); 
+resid   = zeros(nx, 1);
 v       = zeros(nx, ncv);
-workd   = zeros(3 * nx, 1); 
+workd   = zeros(3 * nx, 1);
 workl   = zeros(ncv * ncv + 8 * ncv, 1);
 // Build the symmetric test matrix
 A            = diag(10 * ones(nx,1));
@@ -43,7 +44,7 @@ while(ido <> 99)
     printf('Check the documentation of dsaupd\n\n');
   end
   if(ido == -1 | ido == 1)
-    // Perform matrix vector multiplication 
+    // Perform matrix vector multiplication
     workd(ipntr(2):ipntr(2) + nx - 1) = A * workd(ipntr(1):ipntr(1) + nx - 1);
   end
 end
@@ -53,4 +54,4 @@ howmany = 'A';
 info_dseupd = 0;
 [d, z, resid, v, iparam, ipntr, workd, workl, info_dseupd] = dseupd(rvec, howmany, _select, d, z, sigma, bmat, nx, which, nev, tol, resid, ncv, v, ...
                                                                     iparam, ipntr, workd, workl, info_dseupd);
-assert_checkalmostequal(A * z, z * diag(d), sqrt(%eps), 1.e-10); 
+assert_checkalmostequal(A * z, z * diag(d), sqrt(%eps), 1.e-10);
index e8fa73e..92369e1 100644 (file)
@@ -5,6 +5,8 @@
 //  This file is distributed under the same license as the Scilab package.
 // =============================================================================
 
+// <-- CLI SHELL MODE -->
+
 nx    = 10;
 
 nev   = 3;
@@ -19,9 +21,9 @@ ipntr   = zeros(14, 1);
 _select = zeros(ncv, 1);
 d       = zeros(nev, 1);
 z       = zeros(nx, nev);
-resid   = zeros(nx, 1); 
+resid   = zeros(nx, 1);
 v       = zeros(nx, ncv);
-workd   = zeros(3 * nx, 1); 
+workd   = zeros(3 * nx, 1);
 workl   = zeros(ncv * ncv + 8 * ncv, 1);
 
 // Build the symmetric test matrix
@@ -51,14 +53,14 @@ while(ido <> 99)
   // either convergence is indicated or maxitr has been exceeded.
 
   [ido, resid, v, iparam, ipntr, workd, workl, info_dsaupd] = dsaupd(ido, bmat, nx, which, nev, tol, resid, ncv, v, iparam, ipntr, workd, workl, info_dsaupd);
-  
+
   if(info_dsaupd < 0)
     printf('\nError with dsaupd, info = %d\n',info_dsaupd);
     printf('Check the documentation of dsaupd\n\n');
   end
-  
+
   if(ido == -1 | ido == 1)
-    // Perform matrix vector multiplication 
+    // Perform matrix vector multiplication
     workd(ipntr(2):ipntr(2) + nx - 1) = A * workd(ipntr(1):ipntr(1) + nx - 1);
   end
 end
@@ -71,4 +73,4 @@ info_dseupd = 0;
 [d, z, resid, v, iparam, ipntr, workd, workl, info_dseupd] = dseupd(rvec, howmany, _select, d, z, sigma, bmat, nx, which, nev, tol, resid, ncv, v, ...
                                                                     iparam, ipntr, workd, workl, info_dseupd);
 
-assert_checkalmostequal(A * z, z * diag(d), sqrt(%eps), 1.e-10); 
+assert_checkalmostequal(A * z, z * diag(d), sqrt(%eps), 1.e-10);
index 34a284f..5b69120 100644 (file)
@@ -4,6 +4,7 @@
 //
 //  This file is distributed under the same license as the Scilab package.
 // =============================================================================
+// <-- CLI SHELL MODE -->
 // <-- LINUX ONLY -->
 nx    = 10;
 nev   = 3;
index 579219f..4bc0677 100644 (file)
@@ -4,6 +4,8 @@
 //
 //  This file is distributed under the same license as the Scilab package.
 // =============================================================================
+
+// <-- CLI SHELL MODE -->
 // <-- LINUX ONLY -->
 nx    = 10;
 
index 34a284f..5b69120 100644 (file)
@@ -4,6 +4,7 @@
 //
 //  This file is distributed under the same license as the Scilab package.
 // =============================================================================
+// <-- CLI SHELL MODE -->
 // <-- LINUX ONLY -->
 nx    = 10;
 nev   = 3;
index 579219f..4bc0677 100644 (file)
@@ -4,6 +4,8 @@
 //
 //  This file is distributed under the same license as the Scilab package.
 // =============================================================================
+
+// <-- CLI SHELL MODE -->
 // <-- LINUX ONLY -->
 nx    = 10;