Scicos tests: fix SampleCLK blocks indices 73/18073/2
Paul Bignier [Tue, 12 Apr 2016 08:22:26 +0000 (10:22 +0200)]
Change-Id: I1a30f4a87b477b48e27e001a71fc1a91d4bb30d9

scilab/modules/scicos/tests/unit_tests/Solvers/ODE/DoPri.dia.ref
scilab/modules/scicos/tests/unit_tests/Solvers/ODE/DoPri.tst
scilab/modules/scicos/tests/unit_tests/Solvers/ODE/ExpRK.dia.ref
scilab/modules/scicos/tests/unit_tests/Solvers/ODE/ExpRK.tst
scilab/modules/scicos/tests/unit_tests/Solvers/ODE/ImpRK.dia.ref
scilab/modules/scicos/tests/unit_tests/Solvers/ODE/ImpRK.tst

index f83ca19..20e2b86 100644 (file)
@@ -14,8 +14,8 @@ assert_checktrue(importXcosDiagram("SCI/modules/xcos/tests/unit_tests/Solvers/OD
 Info = scicos_simulate(scs_m, list());
 for i=2:4  // 'max step size' = 5*10^-i, precision
  // Start by updating the clock block period (sampling)
- scs_m.objs(8).model.rpar(1) = 5*10^-i;
  scs_m.objs(9).model.rpar(1) = 5*10^-i;
+ scs_m.objs(10).model.rpar(1) = 5*10^-i;
  // Modify solver and 'max step size' + run DoPri + save results
  scs_m.props.tol(7) = 5*10^-i; scs_m.props.tol(6) = 5;           // 'max step size' + solver
  try scicos_simulate(scs_m, Info); catch disp(lasterror()); end;   // DoPri
index 46fe76c..fb91fb8 100644 (file)
@@ -16,31 +16,31 @@ Info = scicos_simulate(scs_m, list());
 
 for i=2:4  // 'max step size' = 5*10^-i, precision
 
- // Start by updating the clock block period (sampling)
- scs_m.objs(8).model.rpar(1) = 5*10^-i;
- scs_m.objs(9).model.rpar(1) = 5*10^-i;
-
- // Modify solver and 'max step size' + run DoPri + save results
- scs_m.props.tol(7) = 5*10^-i; scs_m.props.tol(6) = 5;           // 'max step size' + solver
- try scicos_simulate(scs_m, Info); catch disp(lasterror()); end;   // DoPri
- doprival = res.values;   // Results
-
- // Modify solver and 'max step size' + run CVode + save results
- scs_m.props.tol(7) = 0; scs_m.props.tol(6) = 1;
- try scicos_simulate(scs_m, Info); catch disp(lasterror()); end;
- cvval = res.values;
-
- // Compare results
- compa = abs(doprival-cvval);
-
- // Extract mean, standard deviation, maximum
- mea = mean(compa);
- [maxi, indexMaxi] = max(compa);
- stdeviation = stdev(compa);
-
- // Verifying closeness of the results
- assert_checktrue(maxi <= 10^-(i+2));
- assert_checktrue(mea <= 10^-(i+2));
- assert_checktrue(stdeviation <= 10^-(i+2));
+    // Start by updating the clock block period (sampling)
+    scs_m.objs(9).model.rpar(1) = 5*10^-i;
+    scs_m.objs(10).model.rpar(1) = 5*10^-i;
+
+    // Modify solver and 'max step size' + run DoPri + save results
+    scs_m.props.tol(7) = 5*10^-i; scs_m.props.tol(6) = 5;           // 'max step size' + solver
+    try scicos_simulate(scs_m, Info); catch disp(lasterror()); end;   // DoPri
+    doprival = res.values;   // Results
+
+    // Modify solver and 'max step size' + run CVode + save results
+    scs_m.props.tol(7) = 0; scs_m.props.tol(6) = 1;
+    try scicos_simulate(scs_m, Info); catch disp(lasterror()); end;
+    cvval = res.values;
+
+    // Compare results
+    compa = abs(doprival-cvval);
+
+    // Extract mean, standard deviation, maximum
+    mea = mean(compa);
+    [maxi, indexMaxi] = max(compa);
+    stdeviation = stdev(compa);
+
+    // Verifying closeness of the results
+    assert_checktrue(maxi <= 10^-(i+2));
+    assert_checktrue(mea <= 10^-(i+2));
+    assert_checktrue(stdeviation <= 10^-(i+2));
 
 end
index 1c92044..b0f3b20 100644 (file)
@@ -14,8 +14,8 @@ assert_checktrue(importXcosDiagram("SCI/modules/xcos/tests/unit_tests/Solvers/OD
 Info = scicos_simulate(scs_m, list());
 for i=2:4  // 'max step size' = 5*10^-i, precision
  // Start by updating the clock block period (sampling)
- scs_m.objs(7).model.rpar(1) = 5*10^-i;
- scs_m.objs(8).model.rpar(1) = 5*10^-i;
+ scs_m.objs(9).model.rpar(1) = 5*10^-i;
+ scs_m.objs(10).model.rpar(1) = 5*10^-i;
  // Modify solver and 'max step size' + run ExpRK + save results
  scs_m.props.tol(7) = 5*10^-i; scs_m.props.tol(6) = 6;     // 'max step size' + solver
  try scicos_simulate(scs_m, Info); catch disp(lasterror()); end;   // ExpRK
index 0ac6939..ec6f967 100644 (file)
@@ -16,31 +16,31 @@ Info = scicos_simulate(scs_m, list());
 
 for i=2:4  // 'max step size' = 5*10^-i, precision
 
- // Start by updating the clock block period (sampling)
- scs_m.objs(7).model.rpar(1) = 5*10^-i;
- scs_m.objs(8).model.rpar(1) = 5*10^-i;
-
- // Modify solver and 'max step size' + run ExpRK + save results
- scs_m.props.tol(7) = 5*10^-i; scs_m.props.tol(6) = 6;     // 'max step size' + solver
- try scicos_simulate(scs_m, Info); catch disp(lasterror()); end;   // ExpRK
- rkval = res.values;   // Results
-
- // Modify solver and 'max step size' + run CVode + save results
- scs_m.props.tol(7) = 0; scs_m.props.tol(6) = 1;
- try scicos_simulate(scs_m, Info); catch disp(lasterror()); end;
- cvval = res.values;
-
- // Compare results
- compa = abs(rkval-cvval);
-
- // Extract mean, standard deviation, maximum
- mea = mean(compa);
- [maxi, indexMaxi] = max(compa);
- stdeviation = stdev(compa);
-
- // Verifying closeness of the results
- assert_checktrue(maxi <= 10^-(i+1));
- assert_checktrue(mea <= 10^-(i+1));
- assert_checktrue(stdeviation <= 10^-(i+2));
+    // Start by updating the clock block period (sampling)
+    scs_m.objs(9).model.rpar(1) = 5*10^-i;
+    scs_m.objs(10).model.rpar(1) = 5*10^-i;
+
+    // Modify solver and 'max step size' + run ExpRK + save results
+    scs_m.props.tol(7) = 5*10^-i; scs_m.props.tol(6) = 6;     // 'max step size' + solver
+    try scicos_simulate(scs_m, Info); catch disp(lasterror()); end;   // ExpRK
+    rkval = res.values;   // Results
+
+    // Modify solver and 'max step size' + run CVode + save results
+    scs_m.props.tol(7) = 0; scs_m.props.tol(6) = 1;
+    try scicos_simulate(scs_m, Info); catch disp(lasterror()); end;
+    cvval = res.values;
+
+    // Compare results
+    compa = abs(rkval-cvval);
+
+    // Extract mean, standard deviation, maximum
+    mea = mean(compa);
+    [maxi, indexMaxi] = max(compa);
+    stdeviation = stdev(compa);
+
+    // Verifying closeness of the results
+    assert_checktrue(maxi <= 10^-(i+1));
+    assert_checktrue(mea <= 10^-(i+1));
+    assert_checktrue(stdeviation <= 10^-(i+2));
 
 end
index f62dd42..f17ce9a 100644 (file)
@@ -14,8 +14,8 @@ assert_checktrue(importXcosDiagram("SCI/modules/xcos/tests/unit_tests/Solvers/OD
 Info = scicos_simulate(scs_m, list());
 for i=2:4  // 'max step size' = 10^-i, precision
  // Start by updating the clock block period (sampling)
- scs_m.objs(7).model.rpar(1) = 5*10^-i;
- scs_m.objs(8).model.rpar(1) = 5*10^-i;
+ scs_m.objs(9).model.rpar(1) = 5*10^-i;
+ scs_m.objs(10).model.rpar(1) = 5*10^-i;
  // Modify solver and 'max step size' + run ImpRK + save results
  scs_m.props.tol(7) = 5*10^-i; scs_m.props.tol(6) = 7;        // 'max step size' + solver
  scs_m.props.tol(2) = 1d-12;                                   // reltol
index 7393f89..93ed683 100644 (file)
@@ -16,32 +16,32 @@ Info = scicos_simulate(scs_m, list());
 
 for i=2:4  // 'max step size' = 10^-i, precision
 
- // Start by updating the clock block period (sampling)
- scs_m.objs(7).model.rpar(1) = 5*10^-i;
- scs_m.objs(8).model.rpar(1) = 5*10^-i;
-
- // Modify solver and 'max step size' + run ImpRK + save results
- scs_m.props.tol(7) = 5*10^-i; scs_m.props.tol(6) = 7;        // 'max step size' + solver
- scs_m.props.tol(2) = 1d-12;                                   // reltol
- try scicos_simulate(scs_m, Info); catch disp(lasterror()); end  // ImpRK
- rkval = res.values;   // Results
-
- // Modify solver and reltol + run CVode + save results
- scs_m.props.tol(6) = 4; scs_m.props.tol(2) = 1d-15;
- try scicos_simulate(scs_m, Info); catch disp(lasterror()); end
- cvval = res.values;
-
- // Compare results
- compa = abs(rkval-cvval);
-
- // Extract mean, standard deviation, maximum
- mea = mean(compa);
- [maxi, indexMaxi] = max(compa);
- stdeviation = stdev(compa);
-
- // Verifying closeness of the results
- assert_checktrue(maxi <= 10^-i);
- assert_checktrue(mea <= 10^-i);
- assert_checktrue(stdeviation <= 10^-i);
+    // Start by updating the clock block period (sampling)
+    scs_m.objs(9).model.rpar(1) = 5*10^-i;
+    scs_m.objs(10).model.rpar(1) = 5*10^-i;
+
+    // Modify solver and 'max step size' + run ImpRK + save results
+    scs_m.props.tol(7) = 5*10^-i; scs_m.props.tol(6) = 7;        // 'max step size' + solver
+    scs_m.props.tol(2) = 1d-12;                                   // reltol
+    try scicos_simulate(scs_m, Info); catch disp(lasterror()); end  // ImpRK
+    rkval = res.values;   // Results
+
+    // Modify solver and reltol + run CVode + save results
+    scs_m.props.tol(6) = 4; scs_m.props.tol(2) = 1d-15;
+    try scicos_simulate(scs_m, Info); catch disp(lasterror()); end
+    cvval = res.values;
+
+    // Compare results
+    compa = abs(rkval-cvval);
+
+    // Extract mean, standard deviation, maximum
+    mea = mean(compa);
+    [maxi, indexMaxi] = max(compa);
+    stdeviation = stdev(compa);
+
+    // Verifying closeness of the results
+    assert_checktrue(maxi <= 10^-i);
+    assert_checktrue(mea <= 10^-i);
+    assert_checktrue(stdeviation <= 10^-i);
 
 end