fecd6c23ce2b49aa61e28741f00a9e8d3482c0ff
[scilab.git] / scilab / modules / scicos / tests / unit_tests / DDaskrN_Platform.tst
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2013 - Scilab Enterprises - Paul Bignier
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7
8 // <-- ENGLISH IMPOSED -->
9 // <-- XCOS TEST -->
10
11 // Import diagram
12 loadScicos();
13 loadXcosLibs();
14 assert_checktrue(importXcosDiagram("SCI/modules/xcos/tests/unit_tests/DDaskr_Platform_test.zcos"));
15
16 // Redefining messagebox() to avoid popup
17 prot = funcprot();
18 funcprot(0);
19 function messagebox(msg, msg_title)
20  disp(msg);
21 endfunction
22 funcprot(prot);
23
24 // looking for the CLOCK_f/EVTDLY_f to update period
25 for path_1=1:length(scs_m.objs)
26     if typeof(scs_m.objs(path_1))=="Block" & scs_m.objs(path_1).gui=="CLOCK_c" then
27         CLOCK_f_scs_m = scs_m.objs(path_1).model.rpar;
28         break;
29     end
30 end
31 for path_2=1:length(CLOCK_f_scs_m)
32     if typeof(CLOCK_f_scs_m.objs(path_2))=="Block" & CLOCK_f_scs_m.objs(path_2).gui=="EVTDLY_c" then
33         EVTDLY_f_blk = CLOCK_f_scs_m.objs(path_2);
34         break;
35     end
36 end
37
38 for i=1:3  // 'max step size' = 10^-i, precision
39
40     // Start by updating the clock block period (sampling)
41     scs_m.objs(path_1).model.rpar.objs(path_2).graphics.exprs = [string(5*(10^-i));"0"];
42
43     // Modify solver + run DDaskr + save results
44     scs_m.props.tol(6) = 101;     // Solver
45     scicos_simulate(scs_m, 'nw'); // DDaskr
46     ddaskrval = res.values;       // Results
47     time = res.time;              // Time
48
49     // Modify solver + run IDA + save results
50     scs_m.props.tol(6) = 100;     // Solver
51     scicos_simulate(scs_m, 'nw'); // IDA
52     idaval = res.values;          // Results
53
54     // Compare results
55     compa = abs(ddaskrval-idaval);
56
57     // Extract mean, standard deviation, maximum
58     mea = mean(compa);
59     [maxi, indexMaxi] = max(compa);
60     stdeviation = st_deviation(compa);
61
62     // Verifying closeness of the results
63     assert_checktrue(maxi <= 10^-(i+3));
64     assert_checktrue(mea <= 10^-(i+3));
65     assert_checktrue(stdeviation <= 10^-(i+3));
66
67 end