[Octave-bug-tracker] [bug #61205] [octave forge] (control) Test failures with LAPACK 3.10.0

[Kai Torben Ohlhus]

Follow-up Comment #5, bug #61205 (project octave):

Regarding comment #4:

OS: openSUSE 15.3
Octave: 6.4.0 (Docker)
OpenBLAS: 0.3.19
control 3.4.0 release candidate


clear all

pkg load control

## extracted from bstmodred.m

A =  [ -0.04165  0.0000  4.9200  -4.9200  0.0000  0.0000  0.0000
       -5.2100  -12.500  0.0000   0.0000  0.0000  0.0000  0.0000
        0.0000   3.3300 -3.3300   0.0000  0.0000  0.0000  0.0000
        0.5450   0.0000  0.0000   0.0000 -0.5450  0.0000  0.0000
        0.0000   0.0000  0.0000   4.9200 -0.04165 0.0000  4.9200
        0.0000   0.0000  0.0000   0.0000 -5.2100 -12.500  0.0000
        0.0000   0.0000  0.0000   0.0000  0.0000  3.3300 -3.3300 ];

B =  [  0.0000   0.0000
        12.500   0.0000
        0.0000   0.0000
        0.0000   0.0000
        0.0000   0.0000
        0.0000   12.500
        0.0000   0.0000 ];

C =  [  1.0000   0.0000  0.0000   0.0000  0.0000  0.0000  0.0000
        0.0000   0.0000  0.0000   1.0000  0.0000  0.0000  0.0000
        0.0000   0.0000  0.0000   0.0000  1.0000  0.0000  0.0000 ];

D =  [  0.0000   0.0000
        0.0000   0.0000
        0.0000   0.0000 ];

G = ss (A, B, C, D, "scaled", true);

[Gr, Info] = bstmodred (G, "beta", 1.0, "tol1", 0.1, "tol2", 0.0);
[Ao, Bo, Co, Do] = ssdata (Gr);

Ae = [  1.2729   0.0000   6.5947   0.0000  -3.4229
        0.0000   0.8169   0.0000   2.4821   0.0000
       -2.9889   0.0000  -2.9028   0.0000  -0.3692
        0.0000  -3.3921   0.0000  -3.1126   0.0000
       -1.4767   0.0000  -2.0339   0.0000  -0.6107 ];

Be = [  0.1331  -0.1331
       -0.0862  -0.0862
       -2.6777   2.6777
       -3.5767  -3.5767
       -2.3033   2.3033 ];

Ce = [ -0.6907  -0.6882   0.0779   0.0958  -0.0038
        0.0676   0.0000   0.6532   0.0000  -0.7522
        0.6907  -0.6882  -0.0779   0.0958   0.0038 ];

De = [  0.0000   0.0000
        0.0000   0.0000
        0.0000   0.0000 ];

HSVe = [  0.8803   0.8506   0.8038   0.4494   0.3973   0.0214   0.0209 ].';

Mo = [Ao, Bo; Co, Do];
Me = [Ae, Be; Ce, De];

save ("-v7", "bstmodred_test_data.mat")



    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?61205>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/