octave-bug-tracker
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Octave-bug-tracker] [bug #65218] [octave-forge](control) test bstmodred


From: Dmitri A. Sergatskov
Subject: [Octave-bug-tracker] [bug #65218] [octave-forge](control) test bstmodred fails on aarch64
Date: Sat, 27 Jan 2024 00:43:23 -0500 (EST)

URL:
  <https://savannah.gnu.org/bugs/?65218>

                 Summary: [octave-forge](control) test bstmodred fails on
aarch64
                   Group: GNU Octave
               Submitter: dasergatskov
               Submitted: Sat 27 Jan 2024 05:43:23 AM UTC
                Category: Octave Forge Package
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: None
                  Status: None
             Assigned to: None
         Originator Name: 
        Originator Email: 
             Open/Closed: Open
                 Release: stable
         Discussion Lock: Any
        Operating System: GNU/Linux
           Fixed Release: None
         Planned Release: None


    _______________________________________________________

Follow-up Comments:


-------------------------------------------------------
Date: Sat 27 Jan 2024 05:43:23 AM UTC By: Dmitri A. Sergatskov <dasergatskov>

octave:3> test bstmodred verbose
>>>>>
/home/dima/.local/share/octave/api-v58/packages/control-4.0.0/bstmodred.m
***** shared Mo, Me, Info, HSVe
 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 ].';

 # Since bstmodred identifies the input/output behavior
 # only input/output behavior is tested. The state space
 # representation is not unique.
 [numo, deno] = tfdata (Gr, "vector");
 [nume, dene] = tfdata (ss (Ae,Be,Ce,De), "vector");

 Mo = [];
 Me = [];
 for iy = 1:size(Ce,1)
   for iu = 1:size(Be,2)
     d = max (abs(deno{iy,iu}));  # normalize to largest den coefficient
(numerical reasons)
     Mo = [ Mo ; numo{iy,iu}'/d ; deno{iy,iu}'/d ];
     d = max (abs(dene{iy,iu}));  # normalize to largest den coefficient
(numerical reasons)
     Me = [ Me ; nume{iy,iu}'/d ; dene{iy,iu}'/d ];
   endfor
 endfor

***** assert (Mo, Me, 1e-4);
!!!!! test failed
ASSERT errors for:  assert (Mo,Me,1e-4)

  Location  |  Observed  |  Expected  |  Reason
     .         O(66x1)      E(58x1)      Dimensions don't match
shared variables 
    Mo =

      -5.9695e-03
       1.2347e-01
       2.5897e-01
       1.1290e+00
       9.0385e-01
       1.0380e-02
       4.7088e-02
       2.3072e-01
       5.4155e-01
       1.0000e+00
       8.0998e-01
       9.0441e-05
       4.2763e-04
      -1.7740e-02
       3.2892e-01
       9.1610e-01
       1.0380e-02
       4.7088e-02
       2.3072e-01
       5.4155e-01
       1.0000e+00
       8.0998e-01
      -7.7034e-05
      -3.2912e-03
       6.6157e-02
       1.5103e-01
       4.3347e-01
       1.0380e-02
       4.7088e-02
       2.3072e-01
       5.4155e-01
       1.0000e+00
       8.0998e-01
       7.7034e-05
       3.2912e-03
      -6.6157e-02
      -1.5103e-01
      -4.3347e-01
       1.0380e-02
       4.7088e-02
       2.3072e-01
       5.4155e-01
       1.0000e+00
       8.0998e-01
       9.0441e-05
       4.2763e-04
      -1.7740e-02
       3.2892e-01
       9.1610e-01
       1.0380e-02
       4.7088e-02
       2.3072e-01
       5.4155e-01
       1.0000e+00
       8.0998e-01
      -5.9695e-03
       1.2347e-01
       2.5897e-01
       1.1290e+00
       9.0385e-01
       1.0380e-02
       4.7088e-02
       2.3072e-01
       5.4155e-01
       1.0000e+00
       8.0998e-01

    Me =

      -5.9698e-03
       1.2346e-01
       2.5895e-01
       1.1289e+00
       9.0382e-01
       1.0380e-02
       4.7089e-02
       2.3072e-01
       5.4155e-01
       1.0000e+00
       8.0996e-01
       8.7688e-05
       4.2734e-04
      -1.7746e-02
       3.2893e-01
       9.1605e-01
       1.0380e-02
       4.7089e-02
       2.3072e-01
       5.4155e-01
       1.0000e+00
       8.0996e-01
      -5.6743e-04
      -2.2611e-02
       5.3511e-01
       7.5318e-02
       1.6876e-01
       8.4399e-01
       1.0000e+00
       5.6743e-04
       2.2611e-02
      -5.3511e-01
       7.5318e-02
       1.6876e-01
       8.4399e-01
       1.0000e+00
       8.7688e-05
       4.2734e-04
      -1.7746e-02
       3.2893e-01
       9.1605e-01
       1.0380e-02
       4.7089e-02
       2.3072e-01
       5.4155e-01
       1.0000e+00
       8.0996e-01
      -5.9698e-03
       1.2346e-01
       2.5895e-01
       1.1289e+00
       9.0382e-01
       1.0380e-02
       4.7089e-02
       2.3072e-01
       5.4155e-01
       1.0000e+00
       8.0996e-01

    Info =

      scalar structure containing the fields:

        n = 7
        ns = 7
        hsv =

           0.880263
           0.850619
           0.803778
           0.449390
           0.397312
           0.021408
           0.020850

        nu = 0
        nr = 5

    HSVe =

       0.880300
       0.850600
       0.803800
       0.449400
       0.397300
       0.021400
       0.020900

octave:4> ver
----------------------------------------------------------------------
GNU Octave Version: 9.0.1 (hg id: 25784f22ead3)
GNU Octave License: GNU General Public License
Operating System: Linux 6.6.13-200.fc39.aarch64 #1 SMP PREEMPT_DYNAMIC Sat Jan
20 17:56:21 UTC 2024 aarch64
----------------------------------------------------------------------
Package Name  | Version | Installation directory
--------------+---------+-----------------------
     control *|   4.0.0 |
.../.local/share/octave/api-v58/packages/control-4.0.0
       optim  |   1.6.2 |
.../.local/share/octave/api-v58/packages/optim-1.6.2
      signal  |   1.4.5 |
.../.local/share/octave/api-v58/packages/signal-1.4.5
  statistics  |   1.6.1 | .../share/octave/api-v58/packages/statistics-1.6.1
      struct  |  1.0.18 |
.../.local/share/octave/api-v58/packages/struct-1.0.18


Plot of Mo and Me is attached.  








    _______________________________________________________
File Attachments:


-------------------------------------------------------
Name: MoMe.png  Size: 42KiB
<http://savannah.gnu.org/bugs/download.php?file_id=55621>

    AGPL NOTICE

These attachments are served by Savane. You can download the corresponding
source code of Savane at
https://git.savannah.nongnu.org/cgit/administration/savane.git/snapshot/savane-bd5dc6a791eeb858c14bff2c54b74c7dcaee57e4.tar.gz

    _______________________________________________________

Reply to this item at:

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

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




reply via email to

[Prev in Thread] Current Thread [Next in Thread]