dc28ac3a6944b6e06a0c6d2e744ee6e2ddf4d463
[scilab.git] / scilab / modules / cacsd / tests / unit_tests / sylv.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 // <-- ENGLISH IMPOSED -->
8 // Run with test_run('cacsd', 'sylv', ['no_check_error_output']);
9 A = rand(4, 4); C = rand(4, 3); B = rand(3, 3);
10 // Solve continuous-time Sylvester equation: A*X+X*B = C
11 X = sylv(A, B, C, "c");
12 expected_norm = 1.D-16;
13 if norm(A*X+X*B-C)-expected_norm > 1.D-15 then bugmes();quit;end
14 // Solve modified discrete-time Sylvester equation: A*X*B+X = C
15 X = sylv(A, B, C, "d");
16 if norm(A*X*B+X-C)-expected_norm > 1.D-15 then bugmes();quit;end
17 // Solve real discrete-time Sylvester equation: A*X*B+X = C
18 X = -sylv(-A, B, C, "d");
19 if norm(A*X*B-X-C)-expected_norm > 1.D-15 then bugmes();quit;end