d08ed903aa4384e0b199021c9479dfbd73ff09e2
[scilab.git] / scilab / modules / scicos / tests / unit_tests / Solvers / DAE / DDaskrN_Platform.unix.dia.ref
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 //
10 // <-- XCOS TEST -->
11 //
12 ilib_verbose(0); //to remove ilib_* traces
13 // Import diagram
14 assert_checktrue(importXcosDiagram("SCI/modules/xcos/tests/unit_tests/Solvers/DAE/Platform.zcos"));
15 // Redefining messagebox() to avoid popup
16 prot = funcprot();
17 funcprot(0);
18 function messagebox(msg, msg_title)
19 endfunction
20 funcprot(prot);
21 Info = scicos_simulate(scs_m, list());
22 --------------------------------------------\
23  Main Modelica : TMPDIR/Platform_im.mo
24
25  Flat Modelica : TMPDIR/Platform_imf.mo
26  Simulation C code :TMPDIR/Platform_im.c
27 for i=1:3  // 'max step size' = 10^-i, precision
28     // Start by updating the clock block period (sampling)
29     Context.per = 5*10^-i;
30     // Modify solver + run DDaskr + save results
31     scs_m.props.tol(6) = 101;     // Solver
32     scicos_simulate(scs_m, Info, Context); // DDaskr
33     ddaskrval = res.values;       // Results
34     // Modify solver + run IDA + save results
35     scs_m.props.tol(6) = 100;     // Solver
36     scicos_simulate(scs_m, Info, Context); // IDA
37     idaval = res.values;          // Results
38     // Compare results
39     compa = abs(ddaskrval-idaval);
40     // Extract mean, standard deviation, maximum
41     mea = mean(compa);
42     [maxi, indexMaxi] = max(compa);
43     stdeviation = st_deviation(compa);
44     // Verifying closeness of the results
45     assert_checktrue(maxi <= 5*10^-(i+3));
46     assert_checktrue(mea <= 5*10^-(i+3));
47     assert_checktrue(stdeviation <= 5*10^-(i+3));
48 end
49 Warning : redefining function: dst                     . Use funcprot(0) to avoid this message
50
51 Warning : redefining function: dst                     . Use funcprot(0) to avoid this message
52
53 Warning : redefining function: dst                     . Use funcprot(0) to avoid this message
54
55 Warning : redefining function: dst                     . Use funcprot(0) to avoid this message
56
57 Warning : redefining function: dst                     . Use funcprot(0) to avoid this message
58
59 Warning : redefining function: dst                     . Use funcprot(0) to avoid this message