Re: pkg test control errors
"Dmitri A. Sergatskov" <[email protected]>
| Newsgroups | gmane.comp.gnu.octave.maintainers |
|---|---|
| Message-ID | <CAO+XyQJtxtAggodfvLyaS0_4vOtd0o856shht4q4WzZKHAQxXw@mail.gmail.com> |
On Tue, Sep 1, 2020 at 12:09 AM Dmitri A. Sergatskov <[email protected]> wrote: > > > I got 12 FAILs all from ltimodels: > > Summary: > > PASS 355 > FAIL 12 > > The log file is attached. Dmitri.
ltimodels.txt
(text/plain, 88.5 KB)
>>>>> processing ltimodels
***** shared ltisys
ltisys = tf (12);
***** assert (ltisys.ts, -2);
***** assert (isct (ltisys));
***** assert (isdt (ltisys));
***** shared ltisys
ltisys = ss (17);
***** assert (ltisys.ts, -2);
***** assert (isct (ltisys));
***** assert (isdt (ltisys));
***** shared ltisys
ltisys = tf (1, [1 1]);
***** assert (ltisys.ts, 0);
***** assert (isct (ltisys));
***** assert (! isdt (ltisys));
***** shared ltisys, ts
ts = 0.1;
ltisys = ss (-1, 1, 1, 0, ts);
***** assert (ltisys.ts, ts);
***** assert (! isct (ltisys));
***** assert (isdt (ltisys));
***** shared Ac, Bc, Ace, Bce
A = [ 17.0 24.0 1.0 8.0 15.0
23.0 5.0 7.0 14.0 16.0
4.0 6.0 13.0 20.0 22.0
10.0 12.0 19.0 21.0 3.0
11.0 18.0 25.0 2.0 9.0 ];
B = [ -1.0 -4.0
4.0 9.0
-9.0 -16.0
16.0 25.0
-25.0 -36.0 ];
tol = 0;
A = A.'; # There's a little mistake in the example
# program of routine AB01OD in SLICOT 5.0
[Ac, Bc, U, ncont] = __sl_ab01od__ (A, B, tol);
Ace = [ 12.8848 3.2345 11.8211 3.3758 -0.8982
4.4741 -12.5544 5.3509 5.9403 1.4360
14.4576 7.6855 23.1452 26.3872 -29.9557
0.0000 1.4805 27.4668 22.6564 -0.0072
0.0000 0.0000 -30.4822 0.6745 18.8680 ];
Bce = [ 31.1199 47.6865
3.2480 0.0000
0.0000 0.0000
0.0000 0.0000
0.0000 0.0000 ];
***** assert <45008> (Ac, Ace, 1e-4);
***** assert (Bc, Bce, 1e-4);
***** shared ac, ec, bc, cc, q, z, ncont, ac_e, ec_e, bc_e, cc_e, q_e, z_e, ncont_e
a = [ 2 0 2 0 -1 3 1
0 1 0 0 1 0 0
0 0 0 1 0 0 1
0 0 2 0 -1 3 1
0 0 0 1 0 0 1
0 1 0 0 1 0 0
0 0 0 1 0 0 1 ];
e = [ 0 0 1 0 0 0 0
0 0 0 0 0 1 0
0 0 0 0 0 0 1
0 0 0 0 0 0 1
0 0 0 1 0 0 0
0 0 1 0 -1 0 0
1 3 0 2 0 0 0 ];
b = [ 2 1 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
1 2 3 ];
c = [ 1 0 0 1 0 0 1
0 -1 1 0 -1 1 0 ];
tol = 0;
[ac, ec, bc, cc, q, z, ncont] = __sl_tg01hd__ (a, e, b, c, tol);
ncont_e = 3;
ac_e = [ 0.0000 0.0000 0.0000 0.0000 -1.2627 0.4334 0.4666
0.0000 2.0000 0.0000 -3.7417 -0.8520 0.2924 -0.4342
0.0000 0.0000 1.7862 0.3780 -0.2651 -0.7723 0.0000
0.0000 0.0000 0.0000 3.7417 0.8520 -0.2924 0.4342
0.0000 0.0000 0.0000 0.0000 -1.5540 0.5334 0.5742
0.0000 0.0000 0.0000 0.0000 -0.6533 0.2242 0.2414
0.0000 0.0000 0.0000 0.0000 -0.5892 0.2022 0.2177 ];
ec_e = [ -1.8325 1.0000 2.3752 0.0000 -0.8214 0.2819 1.8016
0.4887 0.0000 0.3770 -0.5345 0.1874 0.5461 0.0000
-0.1728 0.0000 -0.1333 -1.1339 0.1325 0.3861 0.0000
0.0000 0.0000 0.0000 0.0000 0.8520 -0.2924 0.4342
0.0000 0.0000 0.0000 0.0000 -1.0260 -0.1496 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 1.1937 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 ];
bc_e = [ 1.0000 2.0000 3.0000
2.0000 1.0000 0.0000
0.0000 0.0000 0.0000
0.0000 0.0000 0.0000
0.0000 0.0000 0.0000
0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 ];
cc_e = [ 0.0000 1.0000 0.0000 0.0000 -1.2627 0.4334 0.4666
0.3665 0.0000 -0.9803 -1.6036 0.1874 0.5461 0.0000 ];
q_e = [ 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.7071 0.0000 0.2740 -0.6519 0.0000
0.0000 0.0000 0.0000 0.0000 0.8304 0.3491 -0.4342
0.0000 0.0000 0.0000 -1.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.4003 0.1683 0.9008
0.0000 0.0000 0.7071 0.0000 -0.2740 0.6519 0.0000
1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 ];
z_e = [ 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000
-0.6108 0.0000 0.7917 0.0000 0.0000 0.0000 0.0000
0.4887 0.0000 0.3770 -0.5345 0.1874 0.5461 0.0000
0.0000 0.0000 0.0000 0.0000 -0.4107 0.1410 0.9008
0.6108 0.0000 0.4713 0.2673 -0.1874 -0.5461 0.0000
-0.1222 0.0000 -0.0943 -0.8018 -0.1874 -0.5461 0.0000
0.0000 0.0000 0.0000 0.0000 -0.8520 0.2924 -0.4342 ];
***** assert <45008> (ac, ac_e, 1e-4);
***** assert <45008> (ec, ec_e, 1e-4);
***** assert <45008> (bc, bc_e, 1e-4);
***** assert <45008> (cc, cc_e, 1e-4);
***** assert <45008> (q, q_e, 1e-4);
***** assert <45008> (z, z_e, 1e-4);
***** assert (ncont, ncont_e);
***** shared ac, ec, bc, cc, q, z, ncont, ac_e2, ec_e2, bc_e2, cc_e2, q_e2, z_e2, ncont_e
a = [ 2 0 2 0 -1 3 1
0 1 0 0 1 0 0
0 0 0 1 0 0 1
0 0 2 0 -1 3 1
0 0 0 1 0 0 1
0 1 0 0 1 0 0
0 0 0 1 0 0 1 ];
e = [ 0 0 1 0 0 0 0
0 0 0 0 0 1 0
0 0 0 0 0 0 1
0 0 0 0 0 0 1
0 0 0 1 0 0 0
0 0 1 0 -1 0 0
1 3 0 2 0 0 0 ];
b = [ 2 1 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
1 2 3 ];
c = [ 1 0 0 1 0 0 1
0 -1 1 0 -1 1 0 ];
tol = 0;
[ac, ec, bc, cc, q, z, ncont] = __sl_tg01hd__ (a, e, b, c, tol);
ncont_e = 3;
ac_e2 = [ 0.00000 -0.00000 -0.00000 -0.00000 1.28956 0.49272 0.30702
0.00000 2.00000 -0.00000 -3.74166 -0.92273 0.38545 -0.00000
0.00000 0.00000 1.78619 -0.37796 0.09662 0.23131 -0.77706
0.00000 0.00000 0.00000 3.74166 0.92273 -0.38545 0.00000
0.00000 0.00000 0.00000 0.00000 1.65718 0.63317 0.39454
0.00000 0.00000 0.00000 0.00000 0.59595 0.22770 0.14188
0.00000 0.00000 0.00000 0.00000 0.47386 0.18105 0.11282 ];
ec_e2 = [ 1.83254 -1.00000 2.37525 -0.00000 0.73367 1.75633 0.61404
0.48868 -0.00000 -0.37702 -0.53452 0.06832 0.16356 -0.54947
0.17277 -0.00000 -0.13330 1.13389 -0.04831 -0.11566 0.38853
0.00000 0.00000 0.00000 -0.00000 0.92273 -0.38545 -0.00000
0.00000 0.00000 0.00000 0.00000 1.00350 0.01670 -0.05612
0.00000 0.00000 0.00000 0.00000 0.00000 1.01973 -0.13128
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.19687 ];
bc_e2 = [ -1.00000 -2.00000 -3.00000
2.00000 1.00000 0.00000
0.00000 0.00000 0.00000
0.00000 0.00000 0.00000
0.00000 0.00000 0.00000
0.00000 0.00000 0.00000
0.00000 0.00000 0.00000 ];
cc_e2 = [ -1.2948e-16 1.0000e+00 9.9893e-17 1.7834e-16 -1.2896e+00 -4.9272e-01 -3.0702e-01
3.6651e-01 -1.7636e-48 9.8026e-01 -1.6036e+00 6.8324e-02 1.6356e-01 -5.4947e-01 ];
q_e2 = [ 0.00000 1.00000 -0.00000 -0.00000 0.00000 -0.00000 0.00000
0.00000 0.00000 -0.70711 0.00000 -0.10213 -0.23892 0.65764
0.00000 -0.00000 0.00000 -0.00000 -0.91952 0.39305 -0.00000
0.00000 0.00000 -0.00000 -1.00000 0.00000 -0.00000 0.00000
0.00000 0.00000 0.00000 0.00000 -0.36556 -0.85519 -0.36746
0.00000 0.00000 -0.70711 0.00000 0.10213 0.23892 -0.65764
-1.00000 -0.00000 0.00000 0.00000 0.00000 -0.00000 0.00000 ];
z_e2 = [ 0.00000 1.00000 0.00000 0.00000 0.00000 -0.00000 0.00000
-0.61085 0.00000 -0.79175 0.00000 0.00000 -0.00000 0.00000
0.48868 -0.00000 -0.37702 -0.53452 0.06832 0.16356 -0.54947
0.00000 -0.00000 -0.00000 -0.00000 -0.36683 -0.87816 -0.30702
0.61085 -0.00000 -0.47128 0.26726 -0.06832 -0.16356 0.54947
-0.12217 0.00000 0.09426 -0.80178 -0.06832 -0.16356 0.54947
-0.00000 0.00000 0.00000 0.00000 -0.92273 0.38545 -0.00000 ];
***** assert (ac, ac_e2, 1e-4);
!!!!! test failed
ASSERT errors for: assert (ac,ac_e2,1e-4)
Location | Observed | Expected | Reason
(3,4) 0.37796 -0.37796 Abs err 0.75592 exceeds tol 0.0001 by 0.8
(1,5) -1.2627 1.2896 Abs err 2.5523 exceeds tol 0.0001 by 3
(2,5) -0.85202 -0.92273 Abs err 0.070714 exceeds tol 0.0001 by 0.07
(3,5) -0.26507 0.09662 Abs err 0.36169 exceeds tol 0.0001 by 0.4
(4,5) 0.85202 0.92273 Abs err 0.070714 exceeds tol 0.0001 by 0.07
(5,5) -1.554 1.6572 Abs err 3.2112 exceeds tol 0.0001 by 3
(6,5) -0.65327 0.59595 Abs err 1.2492 exceeds tol 0.0001 by 1
(7,5) -0.58916 0.47386 Abs err 1.063 exceeds tol 0.0001 by 1
(1,6) 0.4334 0.49272 Abs err 0.059315 exceeds tol 0.0001 by 0.06
(2,6) 0.29244 0.38545 Abs err 0.093011 exceeds tol 0.0001 by 0.09
(3,6) -0.77227 0.23131 Abs err 1.0036 exceeds tol 0.0001 by 1
(4,6) -0.29244 -0.38545 Abs err 0.093011 exceeds tol 0.0001 by 0.09
(5,6) 0.5334 0.63317 Abs err 0.099775 exceeds tol 0.0001 by 0.1
(6,6) 0.22422 0.2277 Abs err 0.0034762 exceeds tol 0.0001 by 0.003
(7,6) 0.20222 0.18105 Abs err 0.02117 exceeds tol 0.0001 by 0.02
(1,7) 0.46658 0.30702 Abs err 0.15956 exceeds tol 0.0001 by 0.2
(2,7) -0.43422 -0 Abs err 0.43422 exceeds tol 0.0001 by 0.4
(3,7) -2.0014e-16 -0.77706 Abs err 0.77706 exceeds tol 0.0001 by 0.8
(4,7) 0.43422 0 Abs err 0.43422 exceeds tol 0.0001 by 0.4
(5,7) 0.57423 0.39454 Abs err 0.17969 exceeds tol 0.0001 by 0.2
(6,7) 0.24139 0.14188 Abs err 0.099509 exceeds tol 0.0001 by 0.1
(7,7) 0.2177 0.11282 Abs err 0.10488 exceeds tol 0.0001 by 0.1
shared variables scalar structure containing the fields:
ac =
-0.0000 0.0000 -0.0000 0.0000 -1.2627 0.4334 0.4666
0 2.0000 0.0000 -3.7417 -0.8520 0.2924 -0.4342
0 0 1.7862 0.3780 -0.2651 -0.7723 -0.0000
0 0 0 3.7417 0.8520 -0.2924 0.4342
0 0 0 0 -1.5540 0.5334 0.5742
0 0 0 0 -0.6533 0.2242 0.2414
0 0 0 0 -0.5892 0.2022 0.2177
ec =
-1.8325 1.0000 2.3752 0.0000 -0.8214 0.2819 1.8016
0.4887 -0.0000 0.3770 -0.5345 0.1874 0.5461 0
-0.1728 0.0000 -0.1333 -1.1339 0.1325 0.3861 -0.0000
0 0 0 -0.0000 0.8520 -0.2924 0.4342
0 0 0 0 -1.0260 -0.1496 0.0000
0 0 0 0 0 1.1937 0.0000
0 0 0 0 0 0 1.0000
bc =
1 2 3
2 1 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
cc =
9.0388e-17 1.0000e+00 -8.1395e-17 1.8883e-16 -1.2627e+00 4.3340e-01 4.6658e-01
3.6651e-01 -7.1794e-17 -9.8026e-01 -1.6036e+00 1.8743e-01 5.4608e-01 -1.7044e-16
q =
0 1.0000 0 0 0 0 0
0 0 0.7071 0.0000 0.2740 -0.6519 0
0 -0.0000 -0.0000 -0.0000 0.8304 0.3491 -0.4342
0 -0.0000 0.0000 -1.0000 -0.0000 0.0000 -0.0000
0 0 0.0000 -0.0000 0.4003 0.1683 0.9008
0 -0.0000 0.7071 0.0000 -0.2740 0.6519 -0.0000
1.0000 0 0 0 0 0 0
z =
0.0000 1.0000 0 0 0 0 0
-0.6108 0.0000 0.7917 0.0000 0 0 0
0.4887 -0.0000 0.3770 -0.5345 0.1874 0.5461 -0.0000
-0.0000 0.0000 -0.0000 0.0000 -0.4107 0.1410 0.9008
0.6108 -0.0000 0.4713 0.2673 -0.1874 -0.5461 0.0000
-0.1222 0.0000 -0.0943 -0.8018 -0.1874 -0.5461 0
-0.0000 0.0000 -0.0000 0.0000 -0.8520 0.2924 -0.4342
ncont = 3
ac_e2 =
0 0 0 0 1.2896 0.4927 0.3070
0 2.0000 0 -3.7417 -0.9227 0.3855 0
0 0 1.7862 -0.3780 0.0966 0.2313 -0.7771
0 0 0 3.7417 0.9227 -0.3855 0
0 0 0 0 1.6572 0.6332 0.3945
0 0 0 0 0.5959 0.2277 0.1419
0 0 0 0 0.4739 0.1810 0.1128
ec_e2 =
1.8325 -1.0000 2.3752 0 0.7337 1.7563 0.6140
0.4887 0 -0.3770 -0.5345 0.0683 0.1636 -0.5495
0.1728 0 -0.1333 1.1339 -0.0483 -0.1157 0.3885
0 0 0 0 0.9227 -0.3855 0
0 0 0 0 1.0035 0.0167 -0.0561
0 0 0 0 0 1.0197 -0.1313
0 0 0 0 0 0 1.1969
bc_e2 =
-1 -2 -3
2 1 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
cc_e2 =
-1.2948e-16 1.0000e+00 9.9893e-17 1.7834e-16 -1.2896e+00 -4.9272e-01 -3.0702e-01
3.6651e-01 -1.7636e-48 9.8026e-01 -1.6036e+00 6.8324e-02 1.6356e-01 -5.4947e-01
q_e2 =
0 1.0000 0 0 0 0 0
0 0 -0.7071 0 -0.1021 -0.2389 0.6576
0 0 0 0 -0.9195 0.3931 0
0 0 0 -1.0000 0 0 0
0 0 0 0 -0.3656 -0.8552 -0.3675
0 0 -0.7071 0 0.1021 0.2389 -0.6576
-1.0000 0 0 0 0 0 0
z_e2 =
0 1.0000 0 0 0 0 0
-0.6109 0 -0.7917 0 0 0 0
0.4887 0 -0.3770 -0.5345 0.0683 0.1636 -0.5495
0 0 0 0 -0.3668 -0.8782 -0.3070
0.6109 0 -0.4713 0.2673 -0.0683 -0.1636 0.5495
-0.1222 0 0.0943 -0.8018 -0.0683 -0.1636 0.5495
0 0 0 0 -0.9227 0.3855 0
ncont_e = 3
***** assert (ec, ec_e2, 1e-4);
!!!!! test failed
ASSERT errors for: assert (ec,ec_e2,1e-4)
Location | Observed | Expected | Reason
(1,1) -1.8325 1.8325 Abs err 3.6651 exceeds tol 0.0001 by 4
(3,1) -0.17277 0.17277 Abs err 0.34554 exceeds tol 0.0001 by 0.3
(1,2) 1 -1 Abs err 2 exceeds tol 0.0001 by 2
(2,3) 0.37702 -0.37702 Abs err 0.75404 exceeds tol 0.0001 by 0.8
(3,4) -1.1339 1.1339 Abs err 2.2678 exceeds tol 0.0001 by 2
(1,5) -0.82141 0.73367 Abs err 1.5551 exceeds tol 0.0001 by 2
(2,5) 0.18743 0.06832 Abs err 0.11911 exceeds tol 0.0001 by 0.1
(3,5) 0.13253 -0.04831 Abs err 0.18084 exceeds tol 0.0001 by 0.2
(4,5) 0.85202 0.92273 Abs err 0.070714 exceeds tol 0.0001 by 0.07
(5,5) -1.026 1.0035 Abs err 2.0295 exceeds tol 0.0001 by 2
(1,6) 0.28193 1.7563 Abs err 1.4744 exceeds tol 0.0001 by 1
(2,6) 0.54608 0.16356 Abs err 0.38252 exceeds tol 0.0001 by 0.4
(3,6) 0.38614 -0.11566 Abs err 0.5018 exceeds tol 0.0001 by 0.5
(4,6) -0.29244 -0.38545 Abs err 0.093011 exceeds tol 0.0001 by 0.09
(5,6) -0.14964 0.0167 Abs err 0.16634 exceeds tol 0.0001 by 0.2
(6,6) 1.1937 1.0197 Abs err 0.17397 exceeds tol 0.0001 by 0.2
(1,7) 1.8016 0.61404 Abs err 1.1876 exceeds tol 0.0001 by 1
(2,7) 0 -0.54947 Abs err 0.54947 exceeds tol 0.0001 by 0.5
(3,7) -1.4463e-16 0.38853 Abs err 0.38853 exceeds tol 0.0001 by 0.4
(4,7) 0.43422 -0 Abs err 0.43422 exceeds tol 0.0001 by 0.4
(5,7) 3.9136e-16 -0.05612 Abs err 0.05612 exceeds tol 0.0001 by 0.06
(6,7) 7.1746e-17 -0.13128 Abs err 0.13128 exceeds tol 0.0001 by 0.1
(7,7) 1 1.1969 Abs err 0.19687 exceeds tol 0.0001 by 0.2
shared variables scalar structure containing the fields:
ac =
-0.0000 0.0000 -0.0000 0.0000 -1.2627 0.4334 0.4666
0 2.0000 0.0000 -3.7417 -0.8520 0.2924 -0.4342
0 0 1.7862 0.3780 -0.2651 -0.7723 -0.0000
0 0 0 3.7417 0.8520 -0.2924 0.4342
0 0 0 0 -1.5540 0.5334 0.5742
0 0 0 0 -0.6533 0.2242 0.2414
0 0 0 0 -0.5892 0.2022 0.2177
ec =
-1.8325 1.0000 2.3752 0.0000 -0.8214 0.2819 1.8016
0.4887 -0.0000 0.3770 -0.5345 0.1874 0.5461 0
-0.1728 0.0000 -0.1333 -1.1339 0.1325 0.3861 -0.0000
0 0 0 -0.0000 0.8520 -0.2924 0.4342
0 0 0 0 -1.0260 -0.1496 0.0000
0 0 0 0 0 1.1937 0.0000
0 0 0 0 0 0 1.0000
bc =
1 2 3
2 1 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
cc =
9.0388e-17 1.0000e+00 -8.1395e-17 1.8883e-16 -1.2627e+00 4.3340e-01 4.6658e-01
3.6651e-01 -7.1794e-17 -9.8026e-01 -1.6036e+00 1.8743e-01 5.4608e-01 -1.7044e-16
q =
0 1.0000 0 0 0 0 0
0 0 0.7071 0.0000 0.2740 -0.6519 0
0 -0.0000 -0.0000 -0.0000 0.8304 0.3491 -0.4342
0 -0.0000 0.0000 -1.0000 -0.0000 0.0000 -0.0000
0 0 0.0000 -0.0000 0.4003 0.1683 0.9008
0 -0.0000 0.7071 0.0000 -0.2740 0.6519 -0.0000
1.0000 0 0 0 0 0 0
z =
0.0000 1.0000 0 0 0 0 0
-0.6108 0.0000 0.7917 0.0000 0 0 0
0.4887 -0.0000 0.3770 -0.5345 0.1874 0.5461 -0.0000
-0.0000 0.0000 -0.0000 0.0000 -0.4107 0.1410 0.9008
0.6108 -0.0000 0.4713 0.2673 -0.1874 -0.5461 0.0000
-0.1222 0.0000 -0.0943 -0.8018 -0.1874 -0.5461 0
-0.0000 0.0000 -0.0000 0.0000 -0.8520 0.2924 -0.4342
ncont = 3
ac_e2 =
0 0 0 0 1.2896 0.4927 0.3070
0 2.0000 0 -3.7417 -0.9227 0.3855 0
0 0 1.7862 -0.3780 0.0966 0.2313 -0.7771
0 0 0 3.7417 0.9227 -0.3855 0
0 0 0 0 1.6572 0.6332 0.3945
0 0 0 0 0.5959 0.2277 0.1419
0 0 0 0 0.4739 0.1810 0.1128
ec_e2 =
1.8325 -1.0000 2.3752 0 0.7337 1.7563 0.6140
0.4887 0 -0.3770 -0.5345 0.0683 0.1636 -0.5495
0.1728 0 -0.1333 1.1339 -0.0483 -0.1157 0.3885
0 0 0 0 0.9227 -0.3855 0
0 0 0 0 1.0035 0.0167 -0.0561
0 0 0 0 0 1.0197 -0.1313
0 0 0 0 0 0 1.1969
bc_e2 =
-1 -2 -3
2 1 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
cc_e2 =
-1.2948e-16 1.0000e+00 9.9893e-17 1.7834e-16 -1.2896e+00 -4.9272e-01 -3.0702e-01
3.6651e-01 -1.7636e-48 9.8026e-01 -1.6036e+00 6.8324e-02 1.6356e-01 -5.4947e-01
q_e2 =
0 1.0000 0 0 0 0 0
0 0 -0.7071 0 -0.1021 -0.2389 0.6576
0 0 0 0 -0.9195 0.3931 0
0 0 0 -1.0000 0 0 0
0 0 0 0 -0.3656 -0.8552 -0.3675
0 0 -0.7071 0 0.1021 0.2389 -0.6576
-1.0000 0 0 0 0 0 0
z_e2 =
0 1.0000 0 0 0 0 0
-0.6109 0 -0.7917 0 0 0 0
0.4887 0 -0.3770 -0.5345 0.0683 0.1636 -0.5495
0 0 0 0 -0.3668 -0.8782 -0.3070
0.6109 0 -0.4713 0.2673 -0.0683 -0.1636 0.5495
-0.1222 0 0.0943 -0.8018 -0.0683 -0.1636 0.5495
0 0 0 0 -0.9227 0.3855 0
ncont_e = 3
***** assert (bc, bc_e2, 1e-4);
!!!!! test failed
ASSERT errors for: assert (bc,bc_e2,1e-4)
Location | Observed | Expected | Reason
(1,1) 1 -1 Abs err 2 exceeds tol 0.0001 by 2
(1,2) 2 -2 Abs err 4 exceeds tol 0.0001 by 4
(1,3) 3 -3 Abs err 6 exceeds tol 0.0001 by 6
shared variables scalar structure containing the fields:
ac =
-0.0000 0.0000 -0.0000 0.0000 -1.2627 0.4334 0.4666
0 2.0000 0.0000 -3.7417 -0.8520 0.2924 -0.4342
0 0 1.7862 0.3780 -0.2651 -0.7723 -0.0000
0 0 0 3.7417 0.8520 -0.2924 0.4342
0 0 0 0 -1.5540 0.5334 0.5742
0 0 0 0 -0.6533 0.2242 0.2414
0 0 0 0 -0.5892 0.2022 0.2177
ec =
-1.8325 1.0000 2.3752 0.0000 -0.8214 0.2819 1.8016
0.4887 -0.0000 0.3770 -0.5345 0.1874 0.5461 0
-0.1728 0.0000 -0.1333 -1.1339 0.1325 0.3861 -0.0000
0 0 0 -0.0000 0.8520 -0.2924 0.4342
0 0 0 0 -1.0260 -0.1496 0.0000
0 0 0 0 0 1.1937 0.0000
0 0 0 0 0 0 1.0000
bc =
1 2 3
2 1 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
cc =
9.0388e-17 1.0000e+00 -8.1395e-17 1.8883e-16 -1.2627e+00 4.3340e-01 4.6658e-01
3.6651e-01 -7.1794e-17 -9.8026e-01 -1.6036e+00 1.8743e-01 5.4608e-01 -1.7044e-16
q =
0 1.0000 0 0 0 0 0
0 0 0.7071 0.0000 0.2740 -0.6519 0
0 -0.0000 -0.0000 -0.0000 0.8304 0.3491 -0.4342
0 -0.0000 0.0000 -1.0000 -0.0000 0.0000 -0.0000
0 0 0.0000 -0.0000 0.4003 0.1683 0.9008
0 -0.0000 0.7071 0.0000 -0.2740 0.6519 -0.0000
1.0000 0 0 0 0 0 0
z =
0.0000 1.0000 0 0 0 0 0
-0.6108 0.0000 0.7917 0.0000 0 0 0
0.4887 -0.0000 0.3770 -0.5345 0.1874 0.5461 -0.0000
-0.0000 0.0000 -0.0000 0.0000 -0.4107 0.1410 0.9008
0.6108 -0.0000 0.4713 0.2673 -0.1874 -0.5461 0.0000
-0.1222 0.0000 -0.0943 -0.8018 -0.1874 -0.5461 0
-0.0000 0.0000 -0.0000 0.0000 -0.8520 0.2924 -0.4342
ncont = 3
ac_e2 =
0 0 0 0 1.2896 0.4927 0.3070
0 2.0000 0 -3.7417 -0.9227 0.3855 0
0 0 1.7862 -0.3780 0.0966 0.2313 -0.7771
0 0 0 3.7417 0.9227 -0.3855 0
0 0 0 0 1.6572 0.6332 0.3945
0 0 0 0 0.5959 0.2277 0.1419
0 0 0 0 0.4739 0.1810 0.1128
ec_e2 =
1.8325 -1.0000 2.3752 0 0.7337 1.7563 0.6140
0.4887 0 -0.3770 -0.5345 0.0683 0.1636 -0.5495
0.1728 0 -0.1333 1.1339 -0.0483 -0.1157 0.3885
0 0 0 0 0.9227 -0.3855 0
0 0 0 0 1.0035 0.0167 -0.0561
0 0 0 0 0 1.0197 -0.1313
0 0 0 0 0 0 1.1969
bc_e2 =
-1 -2 -3
2 1 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
cc_e2 =
-1.2948e-16 1.0000e+00 9.9893e-17 1.7834e-16 -1.2896e+00 -4.9272e-01 -3.0702e-01
3.6651e-01 -1.7636e-48 9.8026e-01 -1.6036e+00 6.8324e-02 1.6356e-01 -5.4947e-01
q_e2 =
0 1.0000 0 0 0 0 0
0 0 -0.7071 0 -0.1021 -0.2389 0.6576
0 0 0 0 -0.9195 0.3931 0
0 0 0 -1.0000 0 0 0
0 0 0 0 -0.3656 -0.8552 -0.3675
0 0 -0.7071 0 0.1021 0.2389 -0.6576
-1.0000 0 0 0 0 0 0
z_e2 =
0 1.0000 0 0 0 0 0
-0.6109 0 -0.7917 0 0 0 0
0.4887 0 -0.3770 -0.5345 0.0683 0.1636 -0.5495
0 0 0 0 -0.3668 -0.8782 -0.3070
0.6109 0 -0.4713 0.2673 -0.0683 -0.1636 0.5495
-0.1222 0 0.0943 -0.8018 -0.0683 -0.1636 0.5495
0 0 0 0 -0.9227 0.3855 0
ncont_e = 3
***** assert (cc, cc_e2, 1e-4);
!!!!! test failed
ASSERT errors for: assert (cc,cc_e2,1e-4)
Location | Observed | Expected | Reason
(2,3) -0.98026 0.98026 Abs err 1.9605 exceeds tol 0.0001 by 2
(1,5) -1.2627 -1.2896 Abs err 0.026881 exceeds tol 0.0001 by 0.03
(2,5) 0.18743 0.068324 Abs err 0.11911 exceeds tol 0.0001 by 0.1
(1,6) 0.4334 -0.49272 Abs err 0.92612 exceeds tol 0.0001 by 0.9
(2,6) 0.54608 0.16356 Abs err 0.38252 exceeds tol 0.0001 by 0.4
(1,7) 0.46658 -0.30702 Abs err 0.7736 exceeds tol 0.0001 by 0.8
(2,7) -1.7044e-16 -0.54947 Abs err 0.54947 exceeds tol 0.0001 by 0.5
shared variables scalar structure containing the fields:
ac =
-0.0000 0.0000 -0.0000 0.0000 -1.2627 0.4334 0.4666
0 2.0000 0.0000 -3.7417 -0.8520 0.2924 -0.4342
0 0 1.7862 0.3780 -0.2651 -0.7723 -0.0000
0 0 0 3.7417 0.8520 -0.2924 0.4342
0 0 0 0 -1.5540 0.5334 0.5742
0 0 0 0 -0.6533 0.2242 0.2414
0 0 0 0 -0.5892 0.2022 0.2177
ec =
-1.8325 1.0000 2.3752 0.0000 -0.8214 0.2819 1.8016
0.4887 -0.0000 0.3770 -0.5345 0.1874 0.5461 0
-0.1728 0.0000 -0.1333 -1.1339 0.1325 0.3861 -0.0000
0 0 0 -0.0000 0.8520 -0.2924 0.4342
0 0 0 0 -1.0260 -0.1496 0.0000
0 0 0 0 0 1.1937 0.0000
0 0 0 0 0 0 1.0000
bc =
1 2 3
2 1 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
cc =
9.0388e-17 1.0000e+00 -8.1395e-17 1.8883e-16 -1.2627e+00 4.3340e-01 4.6658e-01
3.6651e-01 -7.1794e-17 -9.8026e-01 -1.6036e+00 1.8743e-01 5.4608e-01 -1.7044e-16
q =
0 1.0000 0 0 0 0 0
0 0 0.7071 0.0000 0.2740 -0.6519 0
0 -0.0000 -0.0000 -0.0000 0.8304 0.3491 -0.4342
0 -0.0000 0.0000 -1.0000 -0.0000 0.0000 -0.0000
0 0 0.0000 -0.0000 0.4003 0.1683 0.9008
0 -0.0000 0.7071 0.0000 -0.2740 0.6519 -0.0000
1.0000 0 0 0 0 0 0
z =
0.0000 1.0000 0 0 0 0 0
-0.6108 0.0000 0.7917 0.0000 0 0 0
0.4887 -0.0000 0.3770 -0.5345 0.1874 0.5461 -0.0000
-0.0000 0.0000 -0.0000 0.0000 -0.4107 0.1410 0.9008
0.6108 -0.0000 0.4713 0.2673 -0.1874 -0.5461 0.0000
-0.1222 0.0000 -0.0943 -0.8018 -0.1874 -0.5461 0
-0.0000 0.0000 -0.0000 0.0000 -0.8520 0.2924 -0.4342
ncont = 3
ac_e2 =
0 0 0 0 1.2896 0.4927 0.3070
0 2.0000 0 -3.7417 -0.9227 0.3855 0
0 0 1.7862 -0.3780 0.0966 0.2313 -0.7771
0 0 0 3.7417 0.9227 -0.3855 0
0 0 0 0 1.6572 0.6332 0.3945
0 0 0 0 0.5959 0.2277 0.1419
0 0 0 0 0.4739 0.1810 0.1128
ec_e2 =
1.8325 -1.0000 2.3752 0 0.7337 1.7563 0.6140
0.4887 0 -0.3770 -0.5345 0.0683 0.1636 -0.5495
0.1728 0 -0.1333 1.1339 -0.0483 -0.1157 0.3885
0 0 0 0 0.9227 -0.3855 0
0 0 0 0 1.0035 0.0167 -0.0561
0 0 0 0 0 1.0197 -0.1313
0 0 0 0 0 0 1.1969
bc_e2 =
-1 -2 -3
2 1 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
cc_e2 =
-1.2948e-16 1.0000e+00 9.9893e-17 1.7834e-16 -1.2896e+00 -4.9272e-01 -3.0702e-01
3.6651e-01 -1.7636e-48 9.8026e-01 -1.6036e+00 6.8324e-02 1.6356e-01 -5.4947e-01
q_e2 =
0 1.0000 0 0 0 0 0
0 0 -0.7071 0 -0.1021 -0.2389 0.6576
0 0 0 0 -0.9195 0.3931 0
0 0 0 -1.0000 0 0 0
0 0 0 0 -0.3656 -0.8552 -0.3675
0 0 -0.7071 0 0.1021 0.2389 -0.6576
-1.0000 0 0 0 0 0 0
z_e2 =
0 1.0000 0 0 0 0 0
-0.6109 0 -0.7917 0 0 0 0
0.4887 0 -0.3770 -0.5345 0.0683 0.1636 -0.5495
0 0 0 0 -0.3668 -0.8782 -0.3070
0.6109 0 -0.4713 0.2673 -0.0683 -0.1636 0.5495
-0.1222 0 0.0943 -0.8018 -0.0683 -0.1636 0.5495
0 0 0 0 -0.9227 0.3855 0
ncont_e = 3
***** assert (q, q_e2, 1e-4);
!!!!! test failed
ASSERT errors for: assert (q,q_e2,1e-4)
Location | Observed | Expected | Reason
(7,1) 1 -1 Abs err 2 exceeds tol 0.0001 by 2
(2,3) 0.70711 -0.70711 Abs err 1.4142 exceeds tol 0.0001 by 1
(6,3) 0.70711 -0.70711 Abs err 1.4142 exceeds tol 0.0001 by 1
(2,5) 0.27402 -0.10213 Abs err 0.37615 exceeds tol 0.0001 by 0.4
(3,5) 0.83042 -0.91952 Abs err 1.7499 exceeds tol 0.0001 by 2
(5,5) 0.40029 -0.36556 Abs err 0.76585 exceeds tol 0.0001 by 0.8
(6,5) -0.27402 0.10213 Abs err 0.37615 exceeds tol 0.0001 by 0.4
(2,6) -0.65185 -0.23892 Abs err 0.41293 exceeds tol 0.0001 by 0.4
(3,6) 0.34908 0.39305 Abs err 0.043967 exceeds tol 0.0001 by 0.04
(5,6) 0.16827 -0.85519 Abs err 1.0235 exceeds tol 0.0001 by 1
(6,6) 0.65185 0.23892 Abs err 0.41293 exceeds tol 0.0001 by 0.4
(2,7) 0 0.65764 Abs err 0.65764 exceeds tol 0.0001 by 0.7
(3,7) -0.43422 -0 Abs err 0.43422 exceeds tol 0.0001 by 0.4
(5,7) 0.90081 -0.36746 Abs err 1.2683 exceeds tol 0.0001 by 1
(6,7) -2.1409e-32 -0.65764 Abs err 0.65764 exceeds tol 0.0001 by 0.7
shared variables scalar structure containing the fields:
ac =
-0.0000 0.0000 -0.0000 0.0000 -1.2627 0.4334 0.4666
0 2.0000 0.0000 -3.7417 -0.8520 0.2924 -0.4342
0 0 1.7862 0.3780 -0.2651 -0.7723 -0.0000
0 0 0 3.7417 0.8520 -0.2924 0.4342
0 0 0 0 -1.5540 0.5334 0.5742
0 0 0 0 -0.6533 0.2242 0.2414
0 0 0 0 -0.5892 0.2022 0.2177
ec =
-1.8325 1.0000 2.3752 0.0000 -0.8214 0.2819 1.8016
0.4887 -0.0000 0.3770 -0.5345 0.1874 0.5461 0
-0.1728 0.0000 -0.1333 -1.1339 0.1325 0.3861 -0.0000
0 0 0 -0.0000 0.8520 -0.2924 0.4342
0 0 0 0 -1.0260 -0.1496 0.0000
0 0 0 0 0 1.1937 0.0000
0 0 0 0 0 0 1.0000
bc =
1 2 3
2 1 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
cc =
9.0388e-17 1.0000e+00 -8.1395e-17 1.8883e-16 -1.2627e+00 4.3340e-01 4.6658e-01
3.6651e-01 -7.1794e-17 -9.8026e-01 -1.6036e+00 1.8743e-01 5.4608e-01 -1.7044e-16
q =
0 1.0000 0 0 0 0 0
0 0 0.7071 0.0000 0.2740 -0.6519 0
0 -0.0000 -0.0000 -0.0000 0.8304 0.3491 -0.4342
0 -0.0000 0.0000 -1.0000 -0.0000 0.0000 -0.0000
0 0 0.0000 -0.0000 0.4003 0.1683 0.9008
0 -0.0000 0.7071 0.0000 -0.2740 0.6519 -0.0000
1.0000 0 0 0 0 0 0
z =
0.0000 1.0000 0 0 0 0 0
-0.6108 0.0000 0.7917 0.0000 0 0 0
0.4887 -0.0000 0.3770 -0.5345 0.1874 0.5461 -0.0000
-0.0000 0.0000 -0.0000 0.0000 -0.4107 0.1410 0.9008
0.6108 -0.0000 0.4713 0.2673 -0.1874 -0.5461 0.0000
-0.1222 0.0000 -0.0943 -0.8018 -0.1874 -0.5461 0
-0.0000 0.0000 -0.0000 0.0000 -0.8520 0.2924 -0.4342
ncont = 3
ac_e2 =
0 0 0 0 1.2896 0.4927 0.3070
0 2.0000 0 -3.7417 -0.9227 0.3855 0
0 0 1.7862 -0.3780 0.0966 0.2313 -0.7771
0 0 0 3.7417 0.9227 -0.3855 0
0 0 0 0 1.6572 0.6332 0.3945
0 0 0 0 0.5959 0.2277 0.1419
0 0 0 0 0.4739 0.1810 0.1128
ec_e2 =
1.8325 -1.0000 2.3752 0 0.7337 1.7563 0.6140
0.4887 0 -0.3770 -0.5345 0.0683 0.1636 -0.5495
0.1728 0 -0.1333 1.1339 -0.0483 -0.1157 0.3885
0 0 0 0 0.9227 -0.3855 0
0 0 0 0 1.0035 0.0167 -0.0561
0 0 0 0 0 1.0197 -0.1313
0 0 0 0 0 0 1.1969
bc_e2 =
-1 -2 -3
2 1 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
cc_e2 =
-1.2948e-16 1.0000e+00 9.9893e-17 1.7834e-16 -1.2896e+00 -4.9272e-01 -3.0702e-01
3.6651e-01 -1.7636e-48 9.8026e-01 -1.6036e+00 6.8324e-02 1.6356e-01 -5.4947e-01
q_e2 =
0 1.0000 0 0 0 0 0
0 0 -0.7071 0 -0.1021 -0.2389 0.6576
0 0 0 0 -0.9195 0.3931 0
0 0 0 -1.0000 0 0 0
0 0 0 0 -0.3656 -0.8552 -0.3675
0 0 -0.7071 0 0.1021 0.2389 -0.6576
-1.0000 0 0 0 0 0 0
z_e2 =
0 1.0000 0 0 0 0 0
-0.6109 0 -0.7917 0 0 0 0
0.4887 0 -0.3770 -0.5345 0.0683 0.1636 -0.5495
0 0 0 0 -0.3668 -0.8782 -0.3070
0.6109 0 -0.4713 0.2673 -0.0683 -0.1636 0.5495
-0.1222 0 0.0943 -0.8018 -0.0683 -0.1636 0.5495
0 0 0 0 -0.9227 0.3855 0
ncont_e = 3
***** assert (z, z_e2, 1e-4);
!!!!! test failed
ASSERT errors for: assert (z,z_e2,1e-4)
Location | Observed | Expected | Reason
(2,3) 0.79175 -0.79175 Abs err 1.5835 exceeds tol 0.0001 by 2
(3,3) 0.37702 -0.37702 Abs err 0.75404 exceeds tol 0.0001 by 0.8
(5,3) 0.47128 -0.47128 Abs err 0.94256 exceeds tol 0.0001 by 0.9
(6,3) -0.094256 0.09426 Abs err 0.18852 exceeds tol 0.0001 by 0.2
(3,5) 0.18743 0.06832 Abs err 0.11911 exceeds tol 0.0001 by 0.1
(4,5) -0.4107 -0.36683 Abs err 0.043873 exceeds tol 0.0001 by 0.04
(5,5) -0.18743 -0.06832 Abs err 0.11911 exceeds tol 0.0001 by 0.1
(6,5) -0.18743 -0.06832 Abs err 0.11911 exceeds tol 0.0001 by 0.1
(7,5) -0.85202 -0.92273 Abs err 0.070714 exceeds tol 0.0001 by 0.07
(3,6) 0.54608 0.16356 Abs err 0.38252 exceeds tol 0.0001 by 0.4
(4,6) 0.14097 -0.87816 Abs err 1.0191 exceeds tol 0.0001 by 1
(5,6) -0.54608 -0.16356 Abs err 0.38252 exceeds tol 0.0001 by 0.4
(6,6) -0.54608 -0.16356 Abs err 0.38252 exceeds tol 0.0001 by 0.4
(7,6) 0.29244 0.38545 Abs err 0.093011 exceeds tol 0.0001 by 0.09
(3,7) -1.3635e-16 -0.54947 Abs err 0.54947 exceeds tol 0.0001 by 0.5
(4,7) 0.90081 -0.30702 Abs err 1.2078 exceeds tol 0.0001 by 1
(5,7) 3.4088e-17 0.54947 Abs err 0.54947 exceeds tol 0.0001 by 0.5
(6,7) 0 0.54947 Abs err 0.54947 exceeds tol 0.0001 by 0.5
(7,7) -0.43422 -0 Abs err 0.43422 exceeds tol 0.0001 by 0.4
shared variables scalar structure containing the fields:
ac =
-0.0000 0.0000 -0.0000 0.0000 -1.2627 0.4334 0.4666
0 2.0000 0.0000 -3.7417 -0.8520 0.2924 -0.4342
0 0 1.7862 0.3780 -0.2651 -0.7723 -0.0000
0 0 0 3.7417 0.8520 -0.2924 0.4342
0 0 0 0 -1.5540 0.5334 0.5742
0 0 0 0 -0.6533 0.2242 0.2414
0 0 0 0 -0.5892 0.2022 0.2177
ec =
-1.8325 1.0000 2.3752 0.0000 -0.8214 0.2819 1.8016
0.4887 -0.0000 0.3770 -0.5345 0.1874 0.5461 0
-0.1728 0.0000 -0.1333 -1.1339 0.1325 0.3861 -0.0000
0 0 0 -0.0000 0.8520 -0.2924 0.4342
0 0 0 0 -1.0260 -0.1496 0.0000
0 0 0 0 0 1.1937 0.0000
0 0 0 0 0 0 1.0000
bc =
1 2 3
2 1 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
cc =
9.0388e-17 1.0000e+00 -8.1395e-17 1.8883e-16 -1.2627e+00 4.3340e-01 4.6658e-01
3.6651e-01 -7.1794e-17 -9.8026e-01 -1.6036e+00 1.8743e-01 5.4608e-01 -1.7044e-16
q =
0 1.0000 0 0 0 0 0
0 0 0.7071 0.0000 0.2740 -0.6519 0
0 -0.0000 -0.0000 -0.0000 0.8304 0.3491 -0.4342
0 -0.0000 0.0000 -1.0000 -0.0000 0.0000 -0.0000
0 0 0.0000 -0.0000 0.4003 0.1683 0.9008
0 -0.0000 0.7071 0.0000 -0.2740 0.6519 -0.0000
1.0000 0 0 0 0 0 0
z =
0.0000 1.0000 0 0 0 0 0
-0.6108 0.0000 0.7917 0.0000 0 0 0
0.4887 -0.0000 0.3770 -0.5345 0.1874 0.5461 -0.0000
-0.0000 0.0000 -0.0000 0.0000 -0.4107 0.1410 0.9008
0.6108 -0.0000 0.4713 0.2673 -0.1874 -0.5461 0.0000
-0.1222 0.0000 -0.0943 -0.8018 -0.1874 -0.5461 0
-0.0000 0.0000 -0.0000 0.0000 -0.8520 0.2924 -0.4342
ncont = 3
ac_e2 =
0 0 0 0 1.2896 0.4927 0.3070
0 2.0000 0 -3.7417 -0.9227 0.3855 0
0 0 1.7862 -0.3780 0.0966 0.2313 -0.7771
0 0 0 3.7417 0.9227 -0.3855 0
0 0 0 0 1.6572 0.6332 0.3945
0 0 0 0 0.5959 0.2277 0.1419
0 0 0 0 0.4739 0.1810 0.1128
ec_e2 =
1.8325 -1.0000 2.3752 0 0.7337 1.7563 0.6140
0.4887 0 -0.3770 -0.5345 0.0683 0.1636 -0.5495
0.1728 0 -0.1333 1.1339 -0.0483 -0.1157 0.3885
0 0 0 0 0.9227 -0.3855 0
0 0 0 0 1.0035 0.0167 -0.0561
0 0 0 0 0 1.0197 -0.1313
0 0 0 0 0 0 1.1969
bc_e2 =
-1 -2 -3
2 1 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
cc_e2 =
-1.2948e-16 1.0000e+00 9.9893e-17 1.7834e-16 -1.2896e+00 -4.9272e-01 -3.0702e-01
3.6651e-01 -1.7636e-48 9.8026e-01 -1.6036e+00 6.8324e-02 1.6356e-01 -5.4947e-01
q_e2 =
0 1.0000 0 0 0 0 0
0 0 -0.7071 0 -0.1021 -0.2389 0.6576
0 0 0 0 -0.9195 0.3931 0
0 0 0 -1.0000 0 0 0
0 0 0 0 -0.3656 -0.8552 -0.3675
0 0 -0.7071 0 0.1021 0.2389 -0.6576
-1.0000 0 0 0 0 0 0
z_e2 =
0 1.0000 0 0 0 0 0
-0.6109 0 -0.7917 0 0 0 0
0.4887 0 -0.3770 -0.5345 0.0683 0.1636 -0.5495
0 0 0 0 -0.3668 -0.8782 -0.3070
0.6109 0 -0.4713 0.2673 -0.0683 -0.1636 0.5495
-0.1222 0 0.0943 -0.8018 -0.0683 -0.1636 0.5495
0 0 0 0 -0.9227 0.3855 0
ncont_e = 3
***** assert (ncont, ncont_e);
***** shared ao, eo, bo, co, q, z, nobsv, ao_e, eo_e, bo_e, co_e, q_e, z_e, nobsv_e
a = [ 2 0 0 0 0 0 0
0 1 0 0 0 1 0
2 0 0 2 0 0 0
0 0 1 0 1 0 1
-1 1 0 -1 0 1 0
3 0 0 3 0 0 0
1 0 1 1 1 0 1 ];
e = [ 0 0 0 0 0 0 1
0 0 0 0 0 0 3
1 0 0 0 0 1 0
0 0 0 0 1 0 2
0 0 0 0 0 -1 0
0 1 0 0 0 0 0
0 0 1 1 0 0 0 ];
b = [ 1 0
0 -1
0 1
1 0
0 -1
0 1
1 0 ];
c = [ 2 0 0 0 0 0 1
1 0 0 0 0 0 2
0 0 0 0 0 0 3 ];
tol = 0;
[ao, eo, bo, co, q, z, nobsv] = __sl_tg01id__ (a, e, b, c, tol);
nobsv_e = 3;
ao_e = [ 0.2177 0.2414 0.5742 0.4342 0.0000 -0.4342 0.4666
0.2022 0.2242 0.5334 -0.2924 -0.7723 0.2924 0.4334
-0.5892 -0.6533 -1.5540 0.8520 -0.2651 -0.8520 -1.2627
0.0000 0.0000 0.0000 3.7417 0.3780 -3.7417 0.0000
0.0000 0.0000 0.0000 0.0000 1.7862 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 2.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 ];
eo_e = [ 1.0000 0.0000 0.0000 0.4342 0.0000 0.0000 1.8016
0.0000 1.1937 -0.1496 -0.2924 0.3861 0.5461 0.2819
0.0000 0.0000 -1.0260 0.8520 0.1325 0.1874 -0.8214
0.0000 0.0000 0.0000 0.0000 -1.1339 -0.5345 0.0000
0.0000 0.0000 0.0000 0.0000 -0.1333 0.3770 2.3752
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000
0.0000 0.0000 0.0000 0.0000 -0.1728 0.4887 -1.8325 ];
bo_e = [ 0.4666 0.0000
0.4334 0.5461
-1.2627 0.1874
0.0000 -1.6036
0.0000 -0.9803
1.0000 0.0000
0.0000 0.3665 ];
co_e = [ 0.0000 0.0000 0.0000 0.0000 0.0000 2.0000 1.0000
0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 2.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 3.0000 ];
q_e = [ 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.7917 0.0000 -0.6108
0.0000 0.5461 0.1874 -0.5345 0.3770 0.0000 0.4887
0.9008 0.1410 -0.4107 0.0000 0.0000 0.0000 0.0000
0.0000 -0.5461 -0.1874 0.2673 0.4713 0.0000 0.6108
0.0000 -0.5461 -0.1874 -0.8018 -0.0943 0.0000 -0.1222
-0.4342 0.2924 -0.8520 0.0000 0.0000 0.0000 0.0000 ];
z_e = [ 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000
0.0000 -0.6519 0.2740 0.0000 0.7071 0.0000 0.0000
-0.4342 0.3491 0.8304 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 -1.0000 0.0000 0.0000 0.0000
0.9008 0.1683 0.4003 0.0000 0.0000 0.0000 0.0000
0.0000 0.6519 -0.2740 0.0000 0.7071 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 ];
***** assert <45008> (ao, ao_e, 1e-4);
***** assert <45008> (eo, eo_e, 1e-4);
***** assert <45008> (bo, bo_e, 1e-4);
***** assert <45008> (co, co_e, 1e-4);
***** assert <45008> (q, q_e, 1e-4);
***** assert <45008> (z, z_e, 1e-4);
***** assert (nobsv, nobsv_e);
***** shared ao, eo, bo, co, q, z, nobsv, ao_e, eo_e, bo_e, co_e, q_e, z_e, nobsv_e
a = [ 2 0 0 0 0 0 0
0 1 0 0 0 1 0
2 0 0 2 0 0 0
0 0 1 0 1 0 1
-1 1 0 -1 0 1 0
3 0 0 3 0 0 0
1 0 1 1 1 0 1 ];
e = [ 0 0 0 0 0 0 1
0 0 0 0 0 0 3
1 0 0 0 0 1 0
0 0 0 0 1 0 2
0 0 0 0 0 -1 0
0 1 0 0 0 0 0
0 0 1 1 0 0 0 ];
b = [ 1 0
0 -1
0 1
1 0
0 -1
0 1
1 0 ];
c = [ 2 0 0 0 0 0 1
1 0 0 0 0 0 2
0 0 0 0 0 0 3 ];
tol = 0;
[ao, eo, bo, co, q, z, nobsv] = __sl_tg01id__ (a, e, b, c, tol);
nobsv_e = 3;
ao_e = [ 0.11282 0.14188 0.39454 0.00000 -0.77706 -0.00000 0.30702
0.18105 0.22770 0.63317 -0.38545 0.23131 0.38545 0.49272
0.47386 0.59595 1.65718 0.92273 0.09662 -0.92273 1.28956
0.00000 0.00000 0.00000 3.74166 -0.37796 -3.74166 -0.00000
0.00000 0.00000 0.00000 0.00000 1.78619 -0.00000 -0.00000
0.00000 0.00000 0.00000 0.00000 0.00000 2.00000 -0.00000
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ];
eo_e = [ 1.19687 -0.13128 -0.05612 -0.00000 0.38853 -0.54947 0.61404
0.00000 1.01973 0.01670 -0.38545 -0.11566 0.16356 1.75633
0.00000 0.00000 1.00350 0.92273 -0.04831 0.06832 0.73367
0.00000 0.00000 0.00000 -0.00000 1.13389 -0.53452 -0.00000
0.00000 0.00000 0.00000 0.00000 -0.13330 -0.37702 2.37525
0.00000 0.00000 0.00000 0.00000 -0.00000 -0.00000 -1.00000
0.00000 0.00000 0.00000 0.00000 0.17277 0.48868 1.83254 ];
bo_e = [ -3.0702e-01 -5.4947e-01
-4.9272e-01 1.6356e-01
-1.2896e+00 6.8324e-02
1.7834e-16 -1.6036e+00
9.9893e-17 9.8026e-01
1.0000e+00 -1.7636e-48
-1.2948e-16 3.6651e-01 ];
co_e = [ 0.00000 0.00000 0.00000 0.00000 0.00000 2.00000 -1.00000
0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 -2.00000
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -3.00000 ];
q_e = [ 0.00000 -0.00000 0.00000 0.00000 0.00000 1.00000 0.00000
0.00000 -0.00000 0.00000 0.00000 -0.79175 0.00000 -0.61085
-0.54947 0.16356 0.06832 -0.53452 -0.37702 -0.00000 0.48868
-0.30702 -0.87816 -0.36683 -0.00000 -0.00000 -0.00000 0.00000
0.54947 -0.16356 -0.06832 0.26726 -0.47128 -0.00000 0.61085
0.54947 -0.16356 -0.06832 -0.80178 0.09426 0.00000 -0.12217
-0.00000 0.38545 -0.92273 0.00000 0.00000 0.00000 -0.00000 ];
z_e = [ 0.00000 -0.00000 0.00000 -0.00000 -0.00000 1.00000 0.00000
0.65764 -0.23892 -0.10213 0.00000 -0.70711 0.00000 0.00000
-0.00000 0.39305 -0.91952 -0.00000 0.00000 -0.00000 0.00000
0.00000 -0.00000 0.00000 -1.00000 -0.00000 0.00000 0.00000
-0.36746 -0.85519 -0.36556 0.00000 0.00000 0.00000 0.00000
-0.65764 0.23892 0.10213 0.00000 -0.70711 0.00000 0.00000
0.00000 -0.00000 0.00000 0.00000 0.00000 -0.00000 -1.00000 ];
***** assert (ao, ao_e, 1e-4);
!!!!! test failed
ASSERT errors for: assert (ao,ao_e,1e-4)
Location | Observed | Expected | Reason
(1,1) 0.2177 0.11282 Abs err 0.10488 exceeds tol 0.0001 by 0.1
(2,1) 0.20222 0.18105 Abs err 0.02117 exceeds tol 0.0001 by 0.02
(3,1) -0.58916 0.47386 Abs err 1.063 exceeds tol 0.0001 by 1
(1,2) 0.24139 0.14188 Abs err 0.099509 exceeds tol 0.0001 by 0.1
(2,2) 0.22422 0.2277 Abs err 0.0034762 exceeds tol 0.0001 by 0.003
(3,2) -0.65327 0.59595 Abs err 1.2492 exceeds tol 0.0001 by 1
(1,3) 0.57423 0.39454 Abs err 0.17969 exceeds tol 0.0001 by 0.2
(2,3) 0.5334 0.63317 Abs err 0.099775 exceeds tol 0.0001 by 0.1
(3,3) -1.554 1.6572 Abs err 3.2112 exceeds tol 0.0001 by 3
(1,4) 0.43422 0 Abs err 0.43422 exceeds tol 0.0001 by 0.4
(2,4) -0.29244 -0.38545 Abs err 0.093011 exceeds tol 0.0001 by 0.09
(3,4) 0.85202 0.92273 Abs err 0.070714 exceeds tol 0.0001 by 0.07
(1,5) -2.0014e-16 -0.77706 Abs err 0.77706 exceeds tol 0.0001 by 0.8
(2,5) -0.77227 0.23131 Abs err 1.0036 exceeds tol 0.0001 by 1
(3,5) -0.26507 0.09662 Abs err 0.36169 exceeds tol 0.0001 by 0.4
(4,5) 0.37796 -0.37796 Abs err 0.75592 exceeds tol 0.0001 by 0.8
(1,6) -0.43422 -0 Abs err 0.43422 exceeds tol 0.0001 by 0.4
(2,6) 0.29244 0.38545 Abs err 0.093011 exceeds tol 0.0001 by 0.09
(3,6) -0.85202 -0.92273 Abs err 0.070714 exceeds tol 0.0001 by 0.07
(1,7) 0.46658 0.30702 Abs err 0.15956 exceeds tol 0.0001 by 0.2
(2,7) 0.4334 0.49272 Abs err 0.059315 exceeds tol 0.0001 by 0.06
(3,7) -1.2627 1.2896 Abs err 2.5523 exceeds tol 0.0001 by 3
shared variables scalar structure containing the fields:
ao =
0.2177 0.2414 0.5742 0.4342 -0.0000 -0.4342 0.4666
0.2022 0.2242 0.5334 -0.2924 -0.7723 0.2924 0.4334
-0.5892 -0.6533 -1.5540 0.8520 -0.2651 -0.8520 -1.2627
0 0 0 3.7417 0.3780 -3.7417 0.0000
0 0 0 0 1.7862 0.0000 -0.0000
0 0 0 0 0 2.0000 0.0000
0 0 0 0 0 0 -0.0000
eo =
1.0000 0.0000 0.0000 0.4342 -0.0000 0 1.8016
0 1.1937 -0.1496 -0.2924 0.3861 0.5461 0.2819
0 0 -1.0260 0.8520 0.1325 0.1874 -0.8214
0 0 0 -0.0000 -1.1339 -0.5345 0.0000
0 0 0 0 -0.1333 0.3770 2.3752
0 0 0 0 0.0000 -0.0000 1.0000
0 0 0 0 -0.1728 0.4887 -1.8325
bo =
4.6658e-01 -1.7044e-16
4.3340e-01 5.4608e-01
-1.2627e+00 1.8743e-01
1.8883e-16 -1.6036e+00
-8.1395e-17 -9.8026e-01
1.0000e+00 -7.1794e-17
9.0388e-17 3.6651e-01
co =
0 0 0 0 0 2 1
0 0 0 0 0 1 2
0 0 0 0 0 0 3
q =
0 0 0 0 0 1.0000 0.0000
0 0 0 0.0000 0.7917 0.0000 -0.6108
-0.0000 0.5461 0.1874 -0.5345 0.3770 -0.0000 0.4887
0.9008 0.1410 -0.4107 0.0000 -0.0000 0.0000 -0.0000
0.0000 -0.5461 -0.1874 0.2673 0.4713 -0.0000 0.6108
0 -0.5461 -0.1874 -0.8018 -0.0943 0.0000 -0.1222
-0.4342 0.2924 -0.8520 0.0000 -0.0000 0.0000 -0.0000
z =
0 0 0 0 0 1.0000 0
0 -0.6519 0.2740 0.0000 0.7071 0 0
-0.4342 0.3491 0.8304 -0.0000 -0.0000 -0.0000 0
-0.0000 0.0000 -0.0000 -1.0000 0.0000 -0.0000 0
0.9008 0.1683 0.4003 -0.0000 0.0000 0 0
-0.0000 0.6519 -0.2740 0.0000 0.7071 -0.0000 0
0 0 0 0 0 0 1.0000
nobsv = 3
ao_e =
0.1128 0.1419 0.3945 0 -0.7771 0 0.3070
0.1810 0.2277 0.6332 -0.3855 0.2313 0.3855 0.4927
0.4739 0.5959 1.6572 0.9227 0.0966 -0.9227 1.2896
0 0 0 3.7417 -0.3780 -3.7417 0
0 0 0 0 1.7862 0 0
0 0 0 0 0 2.0000 0
0 0 0 0 0 0 0
eo_e =
1.1969 -0.1313 -0.0561 0 0.3885 -0.5495 0.6140
0 1.0197 0.0167 -0.3855 -0.1157 0.1636 1.7563
0 0 1.0035 0.9227 -0.0483 0.0683 0.7337
0 0 0 0 1.1339 -0.5345 0
0 0 0 0 -0.1333 -0.3770 2.3752
0 0 0 0 0 0 -1.0000
0 0 0 0 0.1728 0.4887 1.8325
bo_e =
-3.0702e-01 -5.4947e-01
-4.9272e-01 1.6356e-01
-1.2896e+00 6.8324e-02
1.7834e-16 -1.6036e+00
9.9893e-17 9.8026e-01
1.0000e+00 -1.7636e-48
-1.2948e-16 3.6651e-01
co_e =
0 0 0 0 0 2 -1
0 0 0 0 0 1 -2
0 0 0 0 0 0 -3
q_e =
0 0 0 0 0 1.0000 0
0 0 0 0 -0.7917 0 -0.6109
-0.5495 0.1636 0.0683 -0.5345 -0.3770 0 0.4887
-0.3070 -0.8782 -0.3668 0 0 0 0
0.5495 -0.1636 -0.0683 0.2673 -0.4713 0 0.6109
0.5495 -0.1636 -0.0683 -0.8018 0.0943 0 -0.1222
0 0.3855 -0.9227 0 0 0 0
z_e =
0 0 0 0 0 1.0000 0
0.6576 -0.2389 -0.1021 0 -0.7071 0 0
0 0.3931 -0.9195 0 0 0 0
0 0 0 -1.0000 0 0 0
-0.3675 -0.8552 -0.3656 0 0 0 0
-0.6576 0.2389 0.1021 0 -0.7071 0 0
0 0 0 0 0 0 -1.0000
nobsv_e = 3
***** assert (eo, eo_e, 1e-4);
!!!!! test failed
ASSERT errors for: assert (eo,eo_e,1e-4)
Location | Observed | Expected | Reason
(1,1) 1 1.1969 Abs err 0.19687 exceeds tol 0.0001 by 0.2
(1,2) 7.1746e-17 -0.13128 Abs err 0.13128 exceeds tol 0.0001 by 0.1
(2,2) 1.1937 1.0197 Abs err 0.17397 exceeds tol 0.0001 by 0.2
(1,3) 3.9136e-16 -0.05612 Abs err 0.05612 exceeds tol 0.0001 by 0.06
(2,3) -0.14964 0.0167 Abs err 0.16634 exceeds tol 0.0001 by 0.2
(3,3) -1.026 1.0035 Abs err 2.0295 exceeds tol 0.0001 by 2
(1,4) 0.43422 -0 Abs err 0.43422 exceeds tol 0.0001 by 0.4
(2,4) -0.29244 -0.38545 Abs err 0.093011 exceeds tol 0.0001 by 0.09
(3,4) 0.85202 0.92273 Abs err 0.070714 exceeds tol 0.0001 by 0.07
(1,5) -1.4463e-16 0.38853 Abs err 0.38853 exceeds tol 0.0001 by 0.4
(2,5) 0.38614 -0.11566 Abs err 0.5018 exceeds tol 0.0001 by 0.5
(3,5) 0.13253 -0.04831 Abs err 0.18084 exceeds tol 0.0001 by 0.2
(4,5) -1.1339 1.1339 Abs err 2.2678 exceeds tol 0.0001 by 2
(7,5) -0.17277 0.17277 Abs err 0.34554 exceeds tol 0.0001 by 0.3
(1,6) 0 -0.54947 Abs err 0.54947 exceeds tol 0.0001 by 0.5
(2,6) 0.54608 0.16356 Abs err 0.38252 exceeds tol 0.0001 by 0.4
(3,6) 0.18743 0.06832 Abs err 0.11911 exceeds tol 0.0001 by 0.1
(5,6) 0.37702 -0.37702 Abs err 0.75404 exceeds tol 0.0001 by 0.8
(1,7) 1.8016 0.61404 Abs err 1.1876 exceeds tol 0.0001 by 1
(2,7) 0.28193 1.7563 Abs err 1.4744 exceeds tol 0.0001 by 1
(3,7) -0.82141 0.73367 Abs err 1.5551 exceeds tol 0.0001 by 2
(6,7) 1 -1 Abs err 2 exceeds tol 0.0001 by 2
(7,7) -1.8325 1.8325 Abs err 3.6651 exceeds tol 0.0001 by 4
shared variables scalar structure containing the fields:
ao =
0.2177 0.2414 0.5742 0.4342 -0.0000 -0.4342 0.4666
0.2022 0.2242 0.5334 -0.2924 -0.7723 0.2924 0.4334
-0.5892 -0.6533 -1.5540 0.8520 -0.2651 -0.8520 -1.2627
0 0 0 3.7417 0.3780 -3.7417 0.0000
0 0 0 0 1.7862 0.0000 -0.0000
0 0 0 0 0 2.0000 0.0000
0 0 0 0 0 0 -0.0000
eo =
1.0000 0.0000 0.0000 0.4342 -0.0000 0 1.8016
0 1.1937 -0.1496 -0.2924 0.3861 0.5461 0.2819
0 0 -1.0260 0.8520 0.1325 0.1874 -0.8214
0 0 0 -0.0000 -1.1339 -0.5345 0.0000
0 0 0 0 -0.1333 0.3770 2.3752
0 0 0 0 0.0000 -0.0000 1.0000
0 0 0 0 -0.1728 0.4887 -1.8325
bo =
4.6658e-01 -1.7044e-16
4.3340e-01 5.4608e-01
-1.2627e+00 1.8743e-01
1.8883e-16 -1.6036e+00
-8.1395e-17 -9.8026e-01
1.0000e+00 -7.1794e-17
9.0388e-17 3.6651e-01
co =
0 0 0 0 0 2 1
0 0 0 0 0 1 2
0 0 0 0 0 0 3
q =
0 0 0 0 0 1.0000 0.0000
0 0 0 0.0000 0.7917 0.0000 -0.6108
-0.0000 0.5461 0.1874 -0.5345 0.3770 -0.0000 0.4887
0.9008 0.1410 -0.4107 0.0000 -0.0000 0.0000 -0.0000
0.0000 -0.5461 -0.1874 0.2673 0.4713 -0.0000 0.6108
0 -0.5461 -0.1874 -0.8018 -0.0943 0.0000 -0.1222
-0.4342 0.2924 -0.8520 0.0000 -0.0000 0.0000 -0.0000
z =
0 0 0 0 0 1.0000 0
0 -0.6519 0.2740 0.0000 0.7071 0 0
-0.4342 0.3491 0.8304 -0.0000 -0.0000 -0.0000 0
-0.0000 0.0000 -0.0000 -1.0000 0.0000 -0.0000 0
0.9008 0.1683 0.4003 -0.0000 0.0000 0 0
-0.0000 0.6519 -0.2740 0.0000 0.7071 -0.0000 0
0 0 0 0 0 0 1.0000
nobsv = 3
ao_e =
0.1128 0.1419 0.3945 0 -0.7771 0 0.3070
0.1810 0.2277 0.6332 -0.3855 0.2313 0.3855 0.4927
0.4739 0.5959 1.6572 0.9227 0.0966 -0.9227 1.2896
0 0 0 3.7417 -0.3780 -3.7417 0
0 0 0 0 1.7862 0 0
0 0 0 0 0 2.0000 0
0 0 0 0 0 0 0
eo_e =
1.1969 -0.1313 -0.0561 0 0.3885 -0.5495 0.6140
0 1.0197 0.0167 -0.3855 -0.1157 0.1636 1.7563
0 0 1.0035 0.9227 -0.0483 0.0683 0.7337
0 0 0 0 1.1339 -0.5345 0
0 0 0 0 -0.1333 -0.3770 2.3752
0 0 0 0 0 0 -1.0000
0 0 0 0 0.1728 0.4887 1.8325
bo_e =
-3.0702e-01 -5.4947e-01
-4.9272e-01 1.6356e-01
-1.2896e+00 6.8324e-02
1.7834e-16 -1.6036e+00
9.9893e-17 9.8026e-01
1.0000e+00 -1.7636e-48
-1.2948e-16 3.6651e-01
co_e =
0 0 0 0 0 2 -1
0 0 0 0 0 1 -2
0 0 0 0 0 0 -3
q_e =
0 0 0 0 0 1.0000 0
0 0 0 0 -0.7917 0 -0.6109
-0.5495 0.1636 0.0683 -0.5345 -0.3770 0 0.4887
-0.3070 -0.8782 -0.3668 0 0 0 0
0.5495 -0.1636 -0.0683 0.2673 -0.4713 0 0.6109
0.5495 -0.1636 -0.0683 -0.8018 0.0943 0 -0.1222
0 0.3855 -0.9227 0 0 0 0
z_e =
0 0 0 0 0 1.0000 0
0.6576 -0.2389 -0.1021 0 -0.7071 0 0
0 0.3931 -0.9195 0 0 0 0
0 0 0 -1.0000 0 0 0
-0.3675 -0.8552 -0.3656 0 0 0 0
-0.6576 0.2389 0.1021 0 -0.7071 0 0
0 0 0 0 0 0 -1.0000
nobsv_e = 3
***** assert (bo, bo_e, 1e-4);
!!!!! test failed
ASSERT errors for: assert (bo,bo_e,1e-4)
Location | Observed | Expected | Reason
(1,1) 0.46658 -0.30702 Abs err 0.7736 exceeds tol 0.0001 by 0.8
(2,1) 0.4334 -0.49272 Abs err 0.92612 exceeds tol 0.0001 by 0.9
(3,1) -1.2627 -1.2896 Abs err 0.026881 exceeds tol 0.0001 by 0.03
(1,2) -1.7044e-16 -0.54947 Abs err 0.54947 exceeds tol 0.0001 by 0.5
(2,2) 0.54608 0.16356 Abs err 0.38252 exceeds tol 0.0001 by 0.4
(3,2) 0.18743 0.068324 Abs err 0.11911 exceeds tol 0.0001 by 0.1
(5,2) -0.98026 0.98026 Abs err 1.9605 exceeds tol 0.0001 by 2
shared variables scalar structure containing the fields:
ao =
0.2177 0.2414 0.5742 0.4342 -0.0000 -0.4342 0.4666
0.2022 0.2242 0.5334 -0.2924 -0.7723 0.2924 0.4334
-0.5892 -0.6533 -1.5540 0.8520 -0.2651 -0.8520 -1.2627
0 0 0 3.7417 0.3780 -3.7417 0.0000
0 0 0 0 1.7862 0.0000 -0.0000
0 0 0 0 0 2.0000 0.0000
0 0 0 0 0 0 -0.0000
eo =
1.0000 0.0000 0.0000 0.4342 -0.0000 0 1.8016
0 1.1937 -0.1496 -0.2924 0.3861 0.5461 0.2819
0 0 -1.0260 0.8520 0.1325 0.1874 -0.8214
0 0 0 -0.0000 -1.1339 -0.5345 0.0000
0 0 0 0 -0.1333 0.3770 2.3752
0 0 0 0 0.0000 -0.0000 1.0000
0 0 0 0 -0.1728 0.4887 -1.8325
bo =
4.6658e-01 -1.7044e-16
4.3340e-01 5.4608e-01
-1.2627e+00 1.8743e-01
1.8883e-16 -1.6036e+00
-8.1395e-17 -9.8026e-01
1.0000e+00 -7.1794e-17
9.0388e-17 3.6651e-01
co =
0 0 0 0 0 2 1
0 0 0 0 0 1 2
0 0 0 0 0 0 3
q =
0 0 0 0 0 1.0000 0.0000
0 0 0 0.0000 0.7917 0.0000 -0.6108
-0.0000 0.5461 0.1874 -0.5345 0.3770 -0.0000 0.4887
0.9008 0.1410 -0.4107 0.0000 -0.0000 0.0000 -0.0000
0.0000 -0.5461 -0.1874 0.2673 0.4713 -0.0000 0.6108
0 -0.5461 -0.1874 -0.8018 -0.0943 0.0000 -0.1222
-0.4342 0.2924 -0.8520 0.0000 -0.0000 0.0000 -0.0000
z =
0 0 0 0 0 1.0000 0
0 -0.6519 0.2740 0.0000 0.7071 0 0
-0.4342 0.3491 0.8304 -0.0000 -0.0000 -0.0000 0
-0.0000 0.0000 -0.0000 -1.0000 0.0000 -0.0000 0
0.9008 0.1683 0.4003 -0.0000 0.0000 0 0
-0.0000 0.6519 -0.2740 0.0000 0.7071 -0.0000 0
0 0 0 0 0 0 1.0000
nobsv = 3
ao_e =
0.1128 0.1419 0.3945 0 -0.7771 0 0.3070
0.1810 0.2277 0.6332 -0.3855 0.2313 0.3855 0.4927
0.4739 0.5959 1.6572 0.9227 0.0966 -0.9227 1.2896
0 0 0 3.7417 -0.3780 -3.7417 0
0 0 0 0 1.7862 0 0
0 0 0 0 0 2.0000 0
0 0 0 0 0 0 0
eo_e =
1.1969 -0.1313 -0.0561 0 0.3885 -0.5495 0.6140
0 1.0197 0.0167 -0.3855 -0.1157 0.1636 1.7563
0 0 1.0035 0.9227 -0.0483 0.0683 0.7337
0 0 0 0 1.1339 -0.5345 0
0 0 0 0 -0.1333 -0.3770 2.3752
0 0 0 0 0 0 -1.0000
0 0 0 0 0.1728 0.4887 1.8325
bo_e =
-3.0702e-01 -5.4947e-01
-4.9272e-01 1.6356e-01
-1.2896e+00 6.8324e-02
1.7834e-16 -1.6036e+00
9.9893e-17 9.8026e-01
1.0000e+00 -1.7636e-48
-1.2948e-16 3.6651e-01
co_e =
0 0 0 0 0 2 -1
0 0 0 0 0 1 -2
0 0 0 0 0 0 -3
q_e =
0 0 0 0 0 1.0000 0
0 0 0 0 -0.7917 0 -0.6109
-0.5495 0.1636 0.0683 -0.5345 -0.3770 0 0.4887
-0.3070 -0.8782 -0.3668 0 0 0 0
0.5495 -0.1636 -0.0683 0.2673 -0.4713 0 0.6109
0.5495 -0.1636 -0.0683 -0.8018 0.0943 0 -0.1222
0 0.3855 -0.9227 0 0 0 0
z_e =
0 0 0 0 0 1.0000 0
0.6576 -0.2389 -0.1021 0 -0.7071 0 0
0 0.3931 -0.9195 0 0 0 0
0 0 0 -1.0000 0 0 0
-0.3675 -0.8552 -0.3656 0 0 0 0
-0.6576 0.2389 0.1021 0 -0.7071 0 0
0 0 0 0 0 0 -1.0000
nobsv_e = 3
***** assert (co, co_e, 1e-4);
!!!!! test failed
ASSERT errors for: assert (co,co_e,1e-4)
Location | Observed | Expected | Reason
(1,7) 1 -1 Abs err 2 exceeds tol 0.0001 by 2
(2,7) 2 -2 Abs err 4 exceeds tol 0.0001 by 4
(3,7) 3 -3 Abs err 6 exceeds tol 0.0001 by 6
shared variables scalar structure containing the fields:
ao =
0.2177 0.2414 0.5742 0.4342 -0.0000 -0.4342 0.4666
0.2022 0.2242 0.5334 -0.2924 -0.7723 0.2924 0.4334
-0.5892 -0.6533 -1.5540 0.8520 -0.2651 -0.8520 -1.2627
0 0 0 3.7417 0.3780 -3.7417 0.0000
0 0 0 0 1.7862 0.0000 -0.0000
0 0 0 0 0 2.0000 0.0000
0 0 0 0 0 0 -0.0000
eo =
1.0000 0.0000 0.0000 0.4342 -0.0000 0 1.8016
0 1.1937 -0.1496 -0.2924 0.3861 0.5461 0.2819
0 0 -1.0260 0.8520 0.1325 0.1874 -0.8214
0 0 0 -0.0000 -1.1339 -0.5345 0.0000
0 0 0 0 -0.1333 0.3770 2.3752
0 0 0 0 0.0000 -0.0000 1.0000
0 0 0 0 -0.1728 0.4887 -1.8325
bo =
4.6658e-01 -1.7044e-16
4.3340e-01 5.4608e-01
-1.2627e+00 1.8743e-01
1.8883e-16 -1.6036e+00
-8.1395e-17 -9.8026e-01
1.0000e+00 -7.1794e-17
9.0388e-17 3.6651e-01
co =
0 0 0 0 0 2 1
0 0 0 0 0 1 2
0 0 0 0 0 0 3
q =
0 0 0 0 0 1.0000 0.0000
0 0 0 0.0000 0.7917 0.0000 -0.6108
-0.0000 0.5461 0.1874 -0.5345 0.3770 -0.0000 0.4887
0.9008 0.1410 -0.4107 0.0000 -0.0000 0.0000 -0.0000
0.0000 -0.5461 -0.1874 0.2673 0.4713 -0.0000 0.6108
0 -0.5461 -0.1874 -0.8018 -0.0943 0.0000 -0.1222
-0.4342 0.2924 -0.8520 0.0000 -0.0000 0.0000 -0.0000
z =
0 0 0 0 0 1.0000 0
0 -0.6519 0.2740 0.0000 0.7071 0 0
-0.4342 0.3491 0.8304 -0.0000 -0.0000 -0.0000 0
-0.0000 0.0000 -0.0000 -1.0000 0.0000 -0.0000 0
0.9008 0.1683 0.4003 -0.0000 0.0000 0 0
-0.0000 0.6519 -0.2740 0.0000 0.7071 -0.0000 0
0 0 0 0 0 0 1.0000
nobsv = 3
ao_e =
0.1128 0.1419 0.3945 0 -0.7771 0 0.3070
0.1810 0.2277 0.6332 -0.3855 0.2313 0.3855 0.4927
0.4739 0.5959 1.6572 0.9227 0.0966 -0.9227 1.2896
0 0 0 3.7417 -0.3780 -3.7417 0
0 0 0 0 1.7862 0 0
0 0 0 0 0 2.0000 0
0 0 0 0 0 0 0
eo_e =
1.1969 -0.1313 -0.0561 0 0.3885 -0.5495 0.6140
0 1.0197 0.0167 -0.3855 -0.1157 0.1636 1.7563
0 0 1.0035 0.9227 -0.0483 0.0683 0.7337
0 0 0 0 1.1339 -0.5345 0
0 0 0 0 -0.1333 -0.3770 2.3752
0 0 0 0 0 0 -1.0000
0 0 0 0 0.1728 0.4887 1.8325
bo_e =
-3.0702e-01 -5.4947e-01
-4.9272e-01 1.6356e-01
-1.2896e+00 6.8324e-02
1.7834e-16 -1.6036e+00
9.9893e-17 9.8026e-01
1.0000e+00 -1.7636e-48
-1.2948e-16 3.6651e-01
co_e =
0 0 0 0 0 2 -1
0 0 0 0 0 1 -2
0 0 0 0 0 0 -3
q_e =
0 0 0 0 0 1.0000 0
0 0 0 0 -0.7917 0 -0.6109
-0.5495 0.1636 0.0683 -0.5345 -0.3770 0 0.4887
-0.3070 -0.8782 -0.3668 0 0 0 0
0.5495 -0.1636 -0.0683 0.2673 -0.4713 0 0.6109
0.5495 -0.1636 -0.0683 -0.8018 0.0943 0 -0.1222
0 0.3855 -0.9227 0 0 0 0
z_e =
0 0 0 0 0 1.0000 0
0.6576 -0.2389 -0.1021 0 -0.7071 0 0
0 0.3931 -0.9195 0 0 0 0
0 0 0 -1.0000 0 0 0
-0.3675 -0.8552 -0.3656 0 0 0 0
-0.6576 0.2389 0.1021 0 -0.7071 0 0
0 0 0 0 0 0 -1.0000
nobsv_e = 3
***** assert (q, q_e, 1e-4);
!!!!! test failed
ASSERT errors for: assert (q,q_e,1e-4)
Location | Observed | Expected | Reason
(3,1) -1.3635e-16 -0.54947 Abs err 0.54947 exceeds tol 0.0001 by 0.5
(4,1) 0.90081 -0.30702 Abs err 1.2078 exceeds tol 0.0001 by 1
(5,1) 3.4088e-17 0.54947 Abs err 0.54947 exceeds tol 0.0001 by 0.5
(6,1) 0 0.54947 Abs err 0.54947 exceeds tol 0.0001 by 0.5
(7,1) -0.43422 -0 Abs err 0.43422 exceeds tol 0.0001 by 0.4
(3,2) 0.54608 0.16356 Abs err 0.38252 exceeds tol 0.0001 by 0.4
(4,2) 0.14097 -0.87816 Abs err 1.0191 exceeds tol 0.0001 by 1
(5,2) -0.54608 -0.16356 Abs err 0.38252 exceeds tol 0.0001 by 0.4
(6,2) -0.54608 -0.16356 Abs err 0.38252 exceeds tol 0.0001 by 0.4
(7,2) 0.29244 0.38545 Abs err 0.093011 exceeds tol 0.0001 by 0.09
(3,3) 0.18743 0.06832 Abs err 0.11911 exceeds tol 0.0001 by 0.1
(4,3) -0.4107 -0.36683 Abs err 0.043873 exceeds tol 0.0001 by 0.04
(5,3) -0.18743 -0.06832 Abs err 0.11911 exceeds tol 0.0001 by 0.1
(6,3) -0.18743 -0.06832 Abs err 0.11911 exceeds tol 0.0001 by 0.1
(7,3) -0.85202 -0.92273 Abs err 0.070714 exceeds tol 0.0001 by 0.07
(2,5) 0.79175 -0.79175 Abs err 1.5835 exceeds tol 0.0001 by 2
(3,5) 0.37702 -0.37702 Abs err 0.75404 exceeds tol 0.0001 by 0.8
(5,5) 0.47128 -0.47128 Abs err 0.94256 exceeds tol 0.0001 by 0.9
(6,5) -0.094256 0.09426 Abs err 0.18852 exceeds tol 0.0001 by 0.2
shared variables scalar structure containing the fields:
ao =
0.2177 0.2414 0.5742 0.4342 -0.0000 -0.4342 0.4666
0.2022 0.2242 0.5334 -0.2924 -0.7723 0.2924 0.4334
-0.5892 -0.6533 -1.5540 0.8520 -0.2651 -0.8520 -1.2627
0 0 0 3.7417 0.3780 -3.7417 0.0000
0 0 0 0 1.7862 0.0000 -0.0000
0 0 0 0 0 2.0000 0.0000
0 0 0 0 0 0 -0.0000
eo =
1.0000 0.0000 0.0000 0.4342 -0.0000 0 1.8016
0 1.1937 -0.1496 -0.2924 0.3861 0.5461 0.2819
0 0 -1.0260 0.8520 0.1325 0.1874 -0.8214
0 0 0 -0.0000 -1.1339 -0.5345 0.0000
0 0 0 0 -0.1333 0.3770 2.3752
0 0 0 0 0.0000 -0.0000 1.0000
0 0 0 0 -0.1728 0.4887 -1.8325
bo =
4.6658e-01 -1.7044e-16
4.3340e-01 5.4608e-01
-1.2627e+00 1.8743e-01
1.8883e-16 -1.6036e+00
-8.1395e-17 -9.8026e-01
1.0000e+00 -7.1794e-17
9.0388e-17 3.6651e-01
co =
0 0 0 0 0 2 1
0 0 0 0 0 1 2
0 0 0 0 0 0 3
q =
0 0 0 0 0 1.0000 0.0000
0 0 0 0.0000 0.7917 0.0000 -0.6108
-0.0000 0.5461 0.1874 -0.5345 0.3770 -0.0000 0.4887
0.9008 0.1410 -0.4107 0.0000 -0.0000 0.0000 -0.0000
0.0000 -0.5461 -0.1874 0.2673 0.4713 -0.0000 0.6108
0 -0.5461 -0.1874 -0.8018 -0.0943 0.0000 -0.1222
-0.4342 0.2924 -0.8520 0.0000 -0.0000 0.0000 -0.0000
z =
0 0 0 0 0 1.0000 0
0 -0.6519 0.2740 0.0000 0.7071 0 0
-0.4342 0.3491 0.8304 -0.0000 -0.0000 -0.0000 0
-0.0000 0.0000 -0.0000 -1.0000 0.0000 -0.0000 0
0.9008 0.1683 0.4003 -0.0000 0.0000 0 0
-0.0000 0.6519 -0.2740 0.0000 0.7071 -0.0000 0
0 0 0 0 0 0 1.0000
nobsv = 3
ao_e =
0.1128 0.1419 0.3945 0 -0.7771 0 0.3070
0.1810 0.2277 0.6332 -0.3855 0.2313 0.3855 0.4927
0.4739 0.5959 1.6572 0.9227 0.0966 -0.9227 1.2896
0 0 0 3.7417 -0.3780 -3.7417 0
0 0 0 0 1.7862 0 0
0 0 0 0 0 2.0000 0
0 0 0 0 0 0 0
eo_e =
1.1969 -0.1313 -0.0561 0 0.3885 -0.5495 0.6140
0 1.0197 0.0167 -0.3855 -0.1157 0.1636 1.7563
0 0 1.0035 0.9227 -0.0483 0.0683 0.7337
0 0 0 0 1.1339 -0.5345 0
0 0 0 0 -0.1333 -0.3770 2.3752
0 0 0 0 0 0 -1.0000
0 0 0 0 0.1728 0.4887 1.8325
bo_e =
-3.0702e-01 -5.4947e-01
-4.9272e-01 1.6356e-01
-1.2896e+00 6.8324e-02
1.7834e-16 -1.6036e+00
9.9893e-17 9.8026e-01
1.0000e+00 -1.7636e-48
-1.2948e-16 3.6651e-01
co_e =
0 0 0 0 0 2 -1
0 0 0 0 0 1 -2
0 0 0 0 0 0 -3
q_e =
0 0 0 0 0 1.0000 0
0 0 0 0 -0.7917 0 -0.6109
-0.5495 0.1636 0.0683 -0.5345 -0.3770 0 0.4887
-0.3070 -0.8782 -0.3668 0 0 0 0
0.5495 -0.1636 -0.0683 0.2673 -0.4713 0 0.6109
0.5495 -0.1636 -0.0683 -0.8018 0.0943 0 -0.1222
0 0.3855 -0.9227 0 0 0 0
z_e =
0 0 0 0 0 1.0000 0
0.6576 -0.2389 -0.1021 0 -0.7071 0 0
0 0.3931 -0.9195 0 0 0 0
0 0 0 -1.0000 0 0 0
-0.3675 -0.8552 -0.3656 0 0 0 0
-0.6576 0.2389 0.1021 0 -0.7071 0 0
0 0 0 0 0 0 -1.0000
nobsv_e = 3
***** assert (z, z_e, 1e-4);
!!!!! test failed
ASSERT errors for: assert (z,z_e,1e-4)
Location | Observed | Expected | Reason
(2,1) 0 0.65764 Abs err 0.65764 exceeds tol 0.0001 by 0.7
(3,1) -0.43422 -0 Abs err 0.43422 exceeds tol 0.0001 by 0.4
(5,1) 0.90081 -0.36746 Abs err 1.2683 exceeds tol 0.0001 by 1
(6,1) -2.1409e-32 -0.65764 Abs err 0.65764 exceeds tol 0.0001 by 0.7
(2,2) -0.65185 -0.23892 Abs err 0.41293 exceeds tol 0.0001 by 0.4
(3,2) 0.34908 0.39305 Abs err 0.043967 exceeds tol 0.0001 by 0.04
(5,2) 0.16827 -0.85519 Abs err 1.0235 exceeds tol 0.0001 by 1
(6,2) 0.65185 0.23892 Abs err 0.41293 exceeds tol 0.0001 by 0.4
(2,3) 0.27402 -0.10213 Abs err 0.37615 exceeds tol 0.0001 by 0.4
(3,3) 0.83042 -0.91952 Abs err 1.7499 exceeds tol 0.0001 by 2
(5,3) 0.40029 -0.36556 Abs err 0.76585 exceeds tol 0.0001 by 0.8
(6,3) -0.27402 0.10213 Abs err 0.37615 exceeds tol 0.0001 by 0.4
(2,5) 0.70711 -0.70711 Abs err 1.4142 exceeds tol 0.0001 by 1
(6,5) 0.70711 -0.70711 Abs err 1.4142 exceeds tol 0.0001 by 1
(7,7) 1 -1 Abs err 2 exceeds tol 0.0001 by 2
shared variables scalar structure containing the fields:
ao =
0.2177 0.2414 0.5742 0.4342 -0.0000 -0.4342 0.4666
0.2022 0.2242 0.5334 -0.2924 -0.7723 0.2924 0.4334
-0.5892 -0.6533 -1.5540 0.8520 -0.2651 -0.8520 -1.2627
0 0 0 3.7417 0.3780 -3.7417 0.0000
0 0 0 0 1.7862 0.0000 -0.0000
0 0 0 0 0 2.0000 0.0000
0 0 0 0 0 0 -0.0000
eo =
1.0000 0.0000 0.0000 0.4342 -0.0000 0 1.8016
0 1.1937 -0.1496 -0.2924 0.3861 0.5461 0.2819
0 0 -1.0260 0.8520 0.1325 0.1874 -0.8214
0 0 0 -0.0000 -1.1339 -0.5345 0.0000
0 0 0 0 -0.1333 0.3770 2.3752
0 0 0 0 0.0000 -0.0000 1.0000
0 0 0 0 -0.1728 0.4887 -1.8325
bo =
4.6658e-01 -1.7044e-16
4.3340e-01 5.4608e-01
-1.2627e+00 1.8743e-01
1.8883e-16 -1.6036e+00
-8.1395e-17 -9.8026e-01
1.0000e+00 -7.1794e-17
9.0388e-17 3.6651e-01
co =
0 0 0 0 0 2 1
0 0 0 0 0 1 2
0 0 0 0 0 0 3
q =
0 0 0 0 0 1.0000 0.0000
0 0 0 0.0000 0.7917 0.0000 -0.6108
-0.0000 0.5461 0.1874 -0.5345 0.3770 -0.0000 0.4887
0.9008 0.1410 -0.4107 0.0000 -0.0000 0.0000 -0.0000
0.0000 -0.5461 -0.1874 0.2673 0.4713 -0.0000 0.6108
0 -0.5461 -0.1874 -0.8018 -0.0943 0.0000 -0.1222
-0.4342 0.2924 -0.8520 0.0000 -0.0000 0.0000 -0.0000
z =
0 0 0 0 0 1.0000 0
0 -0.6519 0.2740 0.0000 0.7071 0 0
-0.4342 0.3491 0.8304 -0.0000 -0.0000 -0.0000 0
-0.0000 0.0000 -0.0000 -1.0000 0.0000 -0.0000 0
0.9008 0.1683 0.4003 -0.0000 0.0000 0 0
-0.0000 0.6519 -0.2740 0.0000 0.7071 -0.0000 0
0 0 0 0 0 0 1.0000
nobsv = 3
ao_e =
0.1128 0.1419 0.3945 0 -0.7771 0 0.3070
0.1810 0.2277 0.6332 -0.3855 0.2313 0.3855 0.4927
0.4739 0.5959 1.6572 0.9227 0.0966 -0.9227 1.2896
0 0 0 3.7417 -0.3780 -3.7417 0
0 0 0 0 1.7862 0 0
0 0 0 0 0 2.0000 0
0 0 0 0 0 0 0
eo_e =
1.1969 -0.1313 -0.0561 0 0.3885 -0.5495 0.6140
0 1.0197 0.0167 -0.3855 -0.1157 0.1636 1.7563
0 0 1.0035 0.9227 -0.0483 0.0683 0.7337
0 0 0 0 1.1339 -0.5345 0
0 0 0 0 -0.1333 -0.3770 2.3752
0 0 0 0 0 0 -1.0000
0 0 0 0 0.1728 0.4887 1.8325
bo_e =
-3.0702e-01 -5.4947e-01
-4.9272e-01 1.6356e-01
-1.2896e+00 6.8324e-02
1.7834e-16 -1.6036e+00
9.9893e-17 9.8026e-01
1.0000e+00 -1.7636e-48
-1.2948e-16 3.6651e-01
co_e =
0 0 0 0 0 2 -1
0 0 0 0 0 1 -2
0 0 0 0 0 0 -3
q_e =
0 0 0 0 0 1.0000 0
0 0 0 0 -0.7917 0 -0.6109
-0.5495 0.1636 0.0683 -0.5345 -0.3770 0 0.4887
-0.3070 -0.8782 -0.3668 0 0 0 0
0.5495 -0.1636 -0.0683 0.2673 -0.4713 0 0.6109
0.5495 -0.1636 -0.0683 -0.8018 0.0943 0 -0.1222
0 0.3855 -0.9227 0 0 0 0
z_e =
0 0 0 0 0 1.0000 0
0.6576 -0.2389 -0.1021 0 -0.7071 0 0
0 0.3931 -0.9195 0 0 0 0
0 0 0 -1.0000 0 0 0
-0.3675 -0.8552 -0.3656 0 0 0 0
-0.6576 0.2389 0.1021 0 -0.7071 0 0
0 0 0 0 0 0 -1.0000
nobsv_e = 3
***** assert (nobsv, nobsv_e);
***** shared NUM, DEN, NUMe, DENe
A = [ -1.0 0.0 0.0
0.0 -2.0 0.0
0.0 0.0 -3.0 ];
B = [ 0.0 1.0 -1.0
1.0 1.0 0.0 ].';
C = [ 0.0 1.0 1.0
1.0 1.0 1.0 ];
D = [ 1.0 0.0
0.0 1.0 ];
[NUM, DEN] = tfdata (ss (A, B, C, D));
NUMe = {[1, 5, 7], [1]; [1], [1, 5, 5]};
DENe = {[1, 5, 6], [1, 2]; [1, 5, 6], [1, 3, 2]};
***** assert (NUM, NUMe, 1e-4);
***** assert (DEN, DENe, 1e-4);
***** shared NUM, DEN, NUMe, DENe
A = [ 0 ];
B = [ 1 ];
C = [ 1 ];
D = [ 0 ];
[NUM, DEN] = tfdata (ss (A, B, C, D));
NUMe = {[1]};
DENe = {[1, 0]};
***** assert (NUM, NUMe, 1e-4);
***** assert (DEN, DENe, 1e-4);
***** shared Mo, Me
INDEX = [ 3 3 ];
DCOEFF = [ 1.0 6.0 11.0 6.0
1.0 6.0 11.0 6.0 ];
UCOEFF = zeros (2, 2, 4);
u11 = [ 1.0 6.0 12.0 7.0 ];
u12 = [ 0.0 1.0 4.0 3.0 ];
u21 = [ 0.0 0.0 1.0 1.0 ];
u22 = [ 1.0 8.0 20.0 15.0 ];
UCOEFF(1,1,:) = u11;
UCOEFF(1,2,:) = u12;
UCOEFF(2,1,:) = u21;
UCOEFF(2,2,:) = u22;
[Ao, Bo, Co, Do] = __sl_td04ad__ (UCOEFF, DCOEFF, INDEX, 0);
Ae = [ 0.5000 -0.8028 0.9387
4.4047 -2.3380 2.5076
-5.5541 1.6872 -4.1620 ];
Be = [ -0.2000 -1.2500
0.0000 -0.6097
0.0000 2.2217 ];
Ce = [ 0.0000 -0.8679 0.2119
0.0000 0.0000 0.9002 ];
De = [ 1.0000 0.0000
0.0000 1.0000 ];
Mo = [Ao, Bo; Co, Do];
Me = [Ae, Be; Ce, De];
***** assert (Mo, Me, 1e-4);