Xcos solvers : ODEPACK's 'lsodar' implementation
[scilab.git] / scilab / modules / scicos / tests / unit_tests / LSodar.dia.ref
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2012 - Scilab Enterprises
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7 // <-- ENGLISH IMPOSED -->
8 // Execute with exec("SCI/modules/scicos/tests/unit_tests/LSodar.tst");
9 //  or test_run('scicos', 'LSodar', ['no_check_error_output']);
10 // Import diagram
11 loadScicos();
12 loadXcosLibs();
13 assert_checktrue(importXcosDiagram("SCI/modules/xcos/tests/unit_tests/LSodar_test.zcos"));
14 // Set solver to LSodar + run LSodar + save results
15 scs_m.props.tol(6) = 0;                                                                                 // Set solver to LSodar
16 try scicos_simulate(scs_m, 'nw'); catch disp(lasterror()); end; // Run LSodar
17 lsodarval = res.values;         // Results
18 time = res.time;                        // Time
19 // Set solver to CVode BDF/Newton + run + save results
20 scs_m.props.tol(6) = 1;
21 try scicos_simulate(scs_m, 'nw'); catch disp(lasterror()); end;
22 cvval = res.values;
23 // Compare results
24 compa = abs(lsodarval-cvval);
25 // Extract mean, standard deviation, maximum
26 mea = mean(compa);
27 [maxi, indexMaxi] = max(compa);
28 stdeviation = st_deviation(compa);
29 // Verifying closeness of the results
30 assert_checktrue(maxi <= 10^-8);
31 assert_checktrue(mea <= 10^-8);
32 assert_checktrue(stdeviation <= 10^-8);