On 6/5/2012 11:06 AM, Lukas Reichlin wrote:
On 05.06.2012, at 10:28, marco atzeri wrote:

On 6/5/2012 7:16 AM, Lukas Reichlin wrote:


Marco, Alexander

Could it be that you are using BLAS from ATLAS and LAPACK from Reference LAPACK 
3.4.1?
What I meant was using no ATLAS at all (and no Accelerate/Veclib). There is a BLAS included in 
Reference LAPACK 3.4.1. You have to specify "make blaslib" as it is not included in 
"all:" of the Top Level Makefile for LAPACK.

Regards,
Lukas



Lukas,
it is was a pure lapack + blas reference from 3.4.1 .

As cygwin package maintainer for Octave and Lapack (and something more)
I am pretty sure of the blas implementation.

Please find attached the diary for the two cases:

atlas-3.9.72 fails 3 tests
blas/lapack-3.4.1 fails 2 tests.


Regards
Marco
<diary_atlas_3.9.72><diary_lapack_3.4.1>

Marco, could you please run the tests from test_control separately? For example
        test @lti/feedback
        test hnamodred
Otherwise it's quite hard for me to figure out which test was failing. To make 
it easy for you, just copy-paste the block [1] below to the terminal.

Regards,
Lukas


[1] test_control.m
http://octave.svn.sourceforge.net/viewvc/octave/trunk/octave-forge/main/control/inst/test_control.m?revision=10558&view=markup

## test collection
test ltimodels

## LTI methods
test @lti/c2d
test @lti/d2c
test @lti/feedback
test @lti/horzcat
test @lti/inv
test @lti/minreal
test @lti/mtimes
test @lti/norm
test @lti/plus
test @lti/prescale
test @lti/sminreal
test @lti/zero

## robust control
test h2syn
test hinfsyn
test ncfsyn

## ARE solvers
test care
test dare
test kalman

## Lyapunov
test covar
test dlyap
## test dlyapchol  # TODO: add tests
test gram
test lyap
test lyapchol

## model order reduction
test bstmodred
test btamodred
test hnamodred
## test spamodred  # TODO: create test case

## controller order reduction
test btaconred
test cfconred
test fwcfconred
## test spaconred  # TODO: create test case

## identification
test fitfrd

## various oct-files
test ctrbf
test hsvd
test place

## various m-files
test ctrb
test filt
test initial
test issample
test margin
test obsv
test sigma



echo on should be enough...

attached

Regards
Marco
+ ## test collection
+ test ('ltimodels')
  ***** assert (ac, ac_e, 1e-4);
!!!!! test failed
assert (ac,ac_e,1e-4) expected
   0.00000   0.00000   0.00000   0.00000  -1.26270   0.43340   0.46660
   0.00000   2.00000   0.00000  -3.74170  -0.85200   0.29240  -0.43420
   0.00000   0.00000   1.78620   0.37800  -0.26510  -0.77230   0.00000
   0.00000   0.00000   0.00000   3.74170   0.85200  -0.29240   0.43420
   0.00000   0.00000   0.00000   0.00000  -1.55400   0.53340   0.57420
   0.00000   0.00000   0.00000   0.00000  -0.65330   0.22420   0.24140
   0.00000   0.00000   0.00000   0.00000  -0.58920   0.20220   0.21770
but got
  -0.00000   0.00000  -0.00000   0.00000  -1.33893   0.22482   0.39586
   0.00000   2.00000   0.00000  -3.74166  -0.86467   0.14519  -0.48091
   0.00000   0.00000   1.78619   0.37796  -0.13521  -0.80522  -0.00000
   0.00000   0.00000   0.00000   3.74166   0.86467  -0.14519   0.48091
   0.00000   0.00000   0.00000   0.00000  -1.78058   0.29898   0.52644
   0.00000   0.00000   0.00000   0.00000  -0.36617   0.06148   0.10826
   0.00000   0.00000   0.00000   0.00000  -0.53004   0.08900   0.15671
maximum absolute error 0.287125 exceeds tolerance 0.0001
shared variables 
  scalar structure containing the fields:

    ac =

      -0.00000   0.00000  -0.00000   0.00000  -1.33893   0.22482   0.39586
       0.00000   2.00000   0.00000  -3.74166  -0.86467   0.14519  -0.48091
       0.00000   0.00000   1.78619   0.37796  -0.13521  -0.80522  -0.00000
       0.00000   0.00000   0.00000   3.74166   0.86467  -0.14519   0.48091
       0.00000   0.00000   0.00000   0.00000  -1.78058   0.29898   0.52644
       0.00000   0.00000   0.00000   0.00000  -0.36617   0.06148   0.10826
       0.00000   0.00000   0.00000   0.00000  -0.53004   0.08900   0.15671

    ec =

      -1.83254   1.00000   2.37525   0.00000  -0.94854   0.15927   1.75354
       0.48868  -0.00000   0.37702  -0.53452   0.09561   0.56938  -0.00000
      -0.17277   0.00000  -0.13330  -1.13389   0.06760   0.40261  -0.00000
       0.00000   0.00000   0.00000  -0.00000   0.86467  -0.14519   0.48091
       0.00000   0.00000   0.00000   0.00000  -1.00683  -0.08110   0.00000
       0.00000   0.00000   0.00000   0.00000   0.00000   1.21643  -0.00000
       0.00000   0.00000   0.00000   0.00000   0.00000   0.00000   1.00000

    bc =

       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 =

     Columns 1 through 5:

       1.3423e-16   1.0000e+00  -9.9893e-17   1.7834e-16  -1.3389e+00
       3.6651e-01  -9.6650e-17  -9.8026e-01  -1.6036e+00   9.5605e-02

     Columns 6 and 7:

       2.2482e-01   3.9586e-01
       5.6938e-01  -1.6989e-16

    q =

       0.00000   1.00000   0.00000   0.00000   0.00000   0.00000  -0.00000
       0.00000   0.00000   0.70711   0.00000   0.14244  -0.69261  -0.00000
       0.00000  -0.00000  -0.00000   0.00000   0.85880   0.17661  -0.48091
       0.00000   0.00000   0.00000  -1.00000   0.00000  -0.00000   0.00000
       0.00000   0.00000  -0.00000   0.00000   0.47105   0.09687   0.87677
       0.00000   0.00000   0.70711   0.00000  -0.14244   0.69261   0.00000
       1.00000   0.00000   0.00000   0.00000   0.00000   0.00000  -0.00000

    z =

       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.09561   0.56938  -0.00000
       0.00000  -0.00000   0.00000  -0.00000  -0.47427   0.07964   0.87677
       0.61085  -0.00000   0.47128   0.26726  -0.09561  -0.56938   0.00000
      -0.12217   0.00000  -0.09426  -0.80178  -0.09561  -0.56938   0.00000
      -0.00000   0.00000  -0.00000   0.00000  -0.86467   0.14519  -0.48091

    ncont =  3
    ac_e =

       0.00000   0.00000   0.00000   0.00000  -1.26270   0.43340   0.46660
       0.00000   2.00000   0.00000  -3.74170  -0.85200   0.29240  -0.43420
       0.00000   0.00000   1.78620   0.37800  -0.26510  -0.77230   0.00000
       0.00000   0.00000   0.00000   3.74170   0.85200  -0.29240   0.43420
       0.00000   0.00000   0.00000   0.00000  -1.55400   0.53340   0.57420
       0.00000   0.00000   0.00000   0.00000  -0.65330   0.22420   0.24140
       0.00000   0.00000   0.00000   0.00000  -0.58920   0.20220   0.21770

    ec_e =

      -1.83250   1.00000   2.37520   0.00000  -0.82140   0.28190   1.80160
       0.48870   0.00000   0.37700  -0.53450   0.18740   0.54610   0.00000
      -0.17280   0.00000  -0.13330  -1.13390   0.13250   0.38610   0.00000
       0.00000   0.00000   0.00000   0.00000   0.85200  -0.29240   0.43420
       0.00000   0.00000   0.00000   0.00000  -1.02600  -0.14960   0.00000
       0.00000   0.00000   0.00000   0.00000   0.00000   1.19370   0.00000
       0.00000   0.00000   0.00000   0.00000   0.00000   0.00000   1.00000

    bc_e =

       1   2   3
       2   1   0
       0   0   0
       0   0   0
       0   0   0
       0   0   0
       0   0   0

    cc_e =

       0.00000   1.00000   0.00000   0.00000  -1.26270   0.43340   0.46660
       0.36650   0.00000  -0.98030  -1.60360   0.18740   0.54610   0.00000

    q_e =

       0.00000   1.00000   0.00000   0.00000   0.00000   0.00000   0.00000
       0.00000   0.00000   0.70710   0.00000   0.27400  -0.65190   0.00000
       0.00000   0.00000   0.00000   0.00000   0.83040   0.34910  -0.43420
       0.00000   0.00000   0.00000  -1.00000   0.00000   0.00000   0.00000
       0.00000   0.00000   0.00000   0.00000   0.40030   0.16830   0.90080
       0.00000   0.00000   0.70710   0.00000  -0.27400   0.65190   0.00000
       1.00000   0.00000   0.00000   0.00000   0.00000   0.00000   0.00000

    z_e =

       0.00000   1.00000   0.00000   0.00000   0.00000   0.00000   0.00000
      -0.61080   0.00000   0.79170   0.00000   0.00000   0.00000   0.00000
       0.48870   0.00000   0.37700  -0.53450   0.18740   0.54610   0.00000
       0.00000   0.00000   0.00000   0.00000  -0.41070   0.14100   0.90080
       0.61080   0.00000   0.47130   0.26730  -0.18740  -0.54610   0.00000
      -0.12220   0.00000  -0.09430  -0.80180  -0.18740  -0.54610   0.00000
       0.00000   0.00000   0.00000   0.00000  -0.85200   0.29240  -0.43420

    ncont_e =  3
+ ## LTI methods
+ test ('@lti/c2d')
PASSES 5 out of 5 tests
+ test ('@lti/d2c')
PASSES 3 out of 3 tests
+ test ('@lti/feedback')
PASSES 5 out of 5 tests
+ test ('@lti/horzcat')
PASSES 1 out of 1 tests
+ test ('@lti/inv')
PASSES 1 out of 1 tests
+ test ('@lti/minreal')
  ***** assert (Ar, Ae, 1e-4);
!!!!! test failed
assert (Ar,Ae,1e-4) expected
   1.00000  -0.03930  -0.09800  -0.10660   0.07810  -0.23300   0.07770
   0.00000   1.03120   0.27170   0.26090  -0.15330   0.67580  -0.35530
   0.00000   0.00000   1.38870   0.66990  -0.42810   1.63890  -0.76150
   0.00000   0.00000   0.00000  -1.21470   0.24230  -0.97920   0.47880
   0.00000   0.00000   0.00000   0.00000  -1.05450   0.50350  -0.27880
   0.00000   0.00000   0.00000   0.00000   0.00000   1.63550  -0.43230
   0.00000   0.00000   0.00000   0.00000   0.00000   0.00000   1.00000
but got
   1.00000   0.03935   0.09802  -0.10656   0.07805  -0.23302   0.07773
   0.00000   1.03122   0.27169  -0.26093   0.15327  -0.67578   0.35529
   0.00000   0.00000   1.38869  -0.66991   0.42806  -1.63894   0.76153
   0.00000   0.00000   0.00000   1.21467  -0.24232   0.97919  -0.47884
   0.00000   0.00000   0.00000   0.00000   1.05454  -0.50348   0.27881
   0.00000   0.00000   0.00000   0.00000   0.00000   1.63547  -0.43234
   0.00000   0.00000   0.00000   0.00000   0.00000   0.00000   1.00000
maximum absolute error 3.27784 exceeds tolerance 0.0001
shared variables 
  scalar structure containing the fields:

    Ar =

       1.00000   0.03935   0.09802  -0.10656   0.07805  -0.23302   0.07773
       0.00000   1.03122   0.27169  -0.26093   0.15327  -0.67578   0.35529
       0.00000   0.00000   1.38869  -0.66991   0.42806  -1.63894   0.76153
       0.00000   0.00000   0.00000   1.21467  -0.24232   0.97919  -0.47884
       0.00000   0.00000   0.00000   0.00000   1.05454  -0.50348   0.27881
       0.00000   0.00000   0.00000   0.00000   0.00000   1.63547  -0.43234
       0.00000   0.00000   0.00000   0.00000   0.00000   0.00000   1.00000

    Br =

      -5.5968e-01   2.3631e-01
       4.8427e-01   4.9770e-02
       4.7273e-01   1.4915e-01
      -1.8019e-01  -1.1574e+00
      -5.9948e-01  -1.5560e-01
      -1.7294e-01  -3.9991e-01
      -3.3525e-17   2.5000e-01

    Cr =

       0.00000   0.00000   0.00000   0.00000   0.00000   0.00000   4.00000
       0.00000   0.00000   0.00000   0.00000   0.00000   3.15238  -1.75000

    Dr =

       0   0
       0   0

    Er =

       0.40999  -0.25895  -0.50798  -0.31085   0.07054   0.14290  -0.14769
       0.76292  -0.34645   0.09918   0.30065  -0.06189  -0.24831   0.01523
      -0.11198  -0.21245  -0.41840   0.12877  -0.05690   0.42134   0.61817
       0.00000   0.11215  -0.00387  -0.27713   0.07584  -0.09747  -0.39231
       0.00000   0.00000   0.37082   0.42899  -0.10064  -0.14015   0.26995
       0.00000   0.00000   0.00000   0.00000   0.94581  -0.22115   0.23779
       0.00000   0.00000   0.00000   0.57107   0.26478   0.59479  -0.50000

    Ae =

       1.00000  -0.03930  -0.09800  -0.10660   0.07810  -0.23300   0.07770
       0.00000   1.03120   0.27170   0.26090  -0.15330   0.67580  -0.35530
       0.00000   0.00000   1.38870   0.66990  -0.42810   1.63890  -0.76150
       0.00000   0.00000   0.00000  -1.21470   0.24230  -0.97920   0.47880
       0.00000   0.00000   0.00000   0.00000  -1.05450   0.50350  -0.27880
       0.00000   0.00000   0.00000   0.00000   0.00000   1.63550  -0.43230
       0.00000   0.00000   0.00000   0.00000   0.00000   0.00000   1.00000

    Be =

      -0.55970   0.23630
      -0.48430  -0.04980
      -0.47270  -0.14910
       0.18020   1.15740
       0.59950   0.15560
      -0.17290  -0.39990
       0.00000   0.25000

    Ce =

       0.00000   0.00000   0.00000   0.00000   0.00000   0.00000   4.00000
       0.00000   0.00000   0.00000   0.00000   0.00000   3.15240  -1.75000

    De =

       0   0
       0   0

    Ee =

       0.41000   0.25900   0.50800  -0.31090   0.07050   0.14290  -0.14770
      -0.76290  -0.34640   0.09920  -0.30070   0.06190   0.24830  -0.01520
       0.11200  -0.21240  -0.41840  -0.12880   0.05690  -0.42130  -0.61820
       0.00000   0.11220  -0.00390   0.27710  -0.07580   0.09750   0.39230
       0.00000   0.00000   0.37080  -0.42900   0.10060   0.14020  -0.26990
       0.00000   0.00000   0.00000   0.00000   0.94580  -0.22110   0.23780
       0.00000   0.00000   0.00000   0.57110   0.26480   0.59480  -0.50000

+ test ('@lti/mtimes')
PASSES 1 out of 1 tests
+ test ('@lti/norm')
PASSES 4 out of 4 tests
+ test ('@lti/plus')
PASSES 1 out of 1 tests
+ test ('@lti/prescale')
PASSES 22 out of 22 tests
+ test ('@lti/sminreal')
PASSES 12 out of 12 tests
+ test ('@lti/zero')
PASSES 8 out of 8 tests
+ ## robust control
+ test ('h2syn')
PASSES 2 out of 2 tests
+ test ('hinfsyn')
PASSES 2 out of 2 tests
+ test ('ncfsyn')
PASSES 27 out of 27 tests
+ ## ARE solvers
+ test ('care')
PASSES 7 out of 7 tests
+ test ('dare')
PASSES 3 out of 3 tests
+ test ('kalman')
PASSES 3 out of 3 tests
+ ## Lyapunov
+ test ('covar')
PASSES 4 out of 4 tests
+ test ('dlyap')
PASSES 2 out of 2 tests
+ ## test dlyapchol  # TODO: add tests
+ test ('gram')
PASSES 5 out of 5 tests
+ test ('lyap')
PASSES 3 out of 3 tests
+ test ('lyapchol')
PASSES 3 out of 3 tests
+ ## model order reduction
+ test ('bstmodred')
  ***** assert (Mo, Me, 1e-4);
!!!!! test failed
assert (Mo,Me,1e-4) expected
   1.27290   0.00000   6.59470   0.00000  -3.42290   0.13310  -0.13310
   0.00000   0.81690   0.00000   2.48210   0.00000  -0.08620  -0.08620
  -2.98890   0.00000  -2.90280   0.00000  -0.36920  -2.67770   2.67770
   0.00000  -3.39210   0.00000  -3.11260   0.00000  -3.57670  -3.57670
  -1.47670   0.00000  -2.03390   0.00000  -0.61070  -2.30330   2.30330
  -0.69070  -0.68820   0.07790   0.09580  -0.00380   0.00000   0.00000
   0.06760   0.00000   0.65320   0.00000  -0.75220   0.00000   0.00000
   0.69070  -0.68820  -0.07790   0.09580   0.00380   0.00000   0.00000
but got
   1.27295  -0.00000  -6.59466   0.00000  -3.42287  -0.13307   0.13307
  -0.00000   0.81688   0.00000   2.48210   0.00000   0.08620   0.08620
   2.98890  -0.00000  -2.90283  -0.00000   0.36919  -2.67775   2.67775
  -0.00000  -3.39208   0.00000  -3.11263   0.00000   3.57669   3.57669
  -1.47666   0.00000   2.03393   0.00000  -0.61070   2.30328  -2.30328
   0.69073   0.68823   0.07791  -0.09576   0.00376   0.00000   0.00000
  -0.06755   0.00000   0.65316   0.00000   0.75223   0.00000   0.00000
  -0.69073   0.68823  -0.07791  -0.09576  -0.00376   0.00000   0.00000
maximum absolute error 13.1894 exceeds tolerance 0.0001
shared variables 
  scalar structure containing the fields:

    Mo =

       1.27295  -0.00000  -6.59466   0.00000  -3.42287  -0.13307   0.13307
      -0.00000   0.81688   0.00000   2.48210   0.00000   0.08620   0.08620
       2.98890  -0.00000  -2.90283  -0.00000   0.36919  -2.67775   2.67775
      -0.00000  -3.39208   0.00000  -3.11263   0.00000   3.57669   3.57669
      -1.47666   0.00000   2.03393   0.00000  -0.61070   2.30328  -2.30328
       0.69073   0.68823   0.07791  -0.09576   0.00376   0.00000   0.00000
      -0.06755   0.00000   0.65316   0.00000   0.75223   0.00000   0.00000
      -0.69073   0.68823  -0.07791  -0.09576  -0.00376   0.00000   0.00000

    Me =

       1.27290   0.00000   6.59470   0.00000  -3.42290   0.13310  -0.13310
       0.00000   0.81690   0.00000   2.48210   0.00000  -0.08620  -0.08620
      -2.98890   0.00000  -2.90280   0.00000  -0.36920  -2.67770   2.67770
       0.00000  -3.39210   0.00000  -3.11260   0.00000  -3.57670  -3.57670
      -1.47670   0.00000  -2.03390   0.00000  -0.61070  -2.30330   2.30330
      -0.69070  -0.68820   0.07790   0.09580  -0.00380   0.00000   0.00000
       0.06760   0.00000   0.65320   0.00000  -0.75220   0.00000   0.00000
       0.69070  -0.68820  -0.07790   0.09580   0.00380   0.00000   0.00000

    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

+ test ('btamodred')
PASSES 2 out of 2 tests
+ test ('hnamodred')
PASSES 2 out of 2 tests
+ ## test spamodred  # TODO: create test case

+ ## controller order reduction
+ test ('btaconred')
PASSES 2 out of 2 tests
+ test ('cfconred')
PASSES 2 out of 2 tests
+ test ('fwcfconred')
PASSES 2 out of 2 tests
+ ## test spaconred  # TODO: create test case

+ ## identification
+ test ('fitfrd')
PASSES 1 out of 1 tests
+ ## various oct-files
+ test ('ctrbf')
PASSES 4 out of 4 tests
+ test ('hsvd')
PASSES 1 out of 1 tests
+ test ('place')
PASSES 4 out of 4 tests
+ ## various m-files
+ test ('ctrb')
PASSES 1 out of 1 tests
+ test ('filt')
PASSES 4 out of 4 tests
+ test ('initial')
PASSES 2 out of 2 tests
+ test ('issample')
PASSES 30 out of 30 tests
+ test ('margin')
PASSES 2 out of 2 tests
+ test ('obsv')
PASSES 1 out of 1 tests
+ test ('sigma')
PASSES 2 out of 2 tests
+ endscript;
octave:2>  exit

octave:3> source prova.m 
+ ## test collection
+ test ('ltimodels')
  ***** assert (ac, ac_e, 1e-4);
!!!!! test failed
assert (ac,ac_e,1e-4) expected
   0.00000   0.00000   0.00000   0.00000  -1.26270   0.43340   0.46660
   0.00000   2.00000   0.00000  -3.74170  -0.85200   0.29240  -0.43420
   0.00000   0.00000   1.78620   0.37800  -0.26510  -0.77230   0.00000
   0.00000   0.00000   0.00000   3.74170   0.85200  -0.29240   0.43420
   0.00000   0.00000   0.00000   0.00000  -1.55400   0.53340   0.57420
   0.00000   0.00000   0.00000   0.00000  -0.65330   0.22420   0.24140
   0.00000   0.00000   0.00000   0.00000  -0.58920   0.20220   0.21770
but got
   0.00000   0.00000  -0.00000   0.00000  -1.41132  -0.00509  -0.09025
   0.00000   2.00000   0.00000  -3.74166  -0.70600  -0.70821   0.00000
   0.00000   0.00000   1.78619   0.37796   0.05219  -0.05202  -0.81316
   0.00000   0.00000   0.00000   3.74166   0.70600   0.70821  -0.00000
   0.00000   0.00000   0.00000   0.00000   1.98980   0.00718   0.12724
   0.00000   0.00000   0.00000   0.00000   0.01122   0.00004   0.00072
   0.00000   0.00000   0.00000   0.00000   0.15568   0.00056   0.00996
maximum absolute error 3.5438 exceeds tolerance 0.0001
shared variables 
  scalar structure containing the fields:

    ac =

       0.00000   0.00000  -0.00000   0.00000  -1.41132  -0.00509  -0.09025
       0.00000   2.00000   0.00000  -3.74166  -0.70600  -0.70821   0.00000
       0.00000   0.00000   1.78619   0.37796   0.05219  -0.05202  -0.81316
       0.00000   0.00000   0.00000   3.74166   0.70600   0.70821  -0.00000
       0.00000   0.00000   0.00000   0.00000   1.98980   0.00718   0.12724
       0.00000   0.00000   0.00000   0.00000   0.01122   0.00004   0.00072
       0.00000   0.00000   0.00000   0.00000   0.15568   0.00056   0.00996

    ec =

       1.83254   1.00000   2.37525   0.00000  -1.41064   1.40624  -0.18050
      -0.48868  -0.00000   0.37702  -0.53452  -0.03690   0.03679   0.57499
       0.17277   0.00000  -0.13330  -1.13389  -0.02609   0.02601   0.40658
       0.00000   0.00000   0.00000  -0.00000   0.70600   0.70821   0.00000
       0.00000   0.00000   0.00000   0.00000   1.00102  -0.00203  -0.03180
       0.00000   0.00000   0.00000   0.00000   0.00000   1.00101   0.03163
       0.00000   0.00000   0.00000   0.00000   0.00000   0.00000   1.22226

    bc =

       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 =

     Columns 1 through 5:

       3.5153e-18   1.0000e+00  -9.6723e-17   1.8014e-16  -1.4113e+00
      -3.6651e-01  -4.4660e-17  -9.8026e-01  -1.6036e+00  -3.6902e-02

     Columns 6 and 7:

      -5.0932e-03  -9.0251e-02
       3.6787e-02   5.7499e-01

    q =

       0.00000   1.00000  -0.00000   0.00000  -0.00000  -0.00000   0.00000
       0.00000  -0.00000   0.70711   0.00000   0.05530  -0.05501  -0.70279
       0.00000  -0.00000  -0.00000   0.00000  -0.70528  -0.70893   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.70460   0.70098  -0.11031
       0.00000  -0.00000   0.70711   0.00000  -0.05530   0.05501   0.70279
       1.00000   0.00000   0.00000   0.00000  -0.00000   0.00000   0.00000

    z =

      -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.03690   0.03679   0.57499
      -0.00000  -0.00000   0.00000  -0.00000  -0.70532   0.70312  -0.09025
      -0.61085  -0.00000   0.47128   0.26726   0.03690  -0.03679  -0.57499
       0.12217   0.00000  -0.09426  -0.80178   0.03690  -0.03679  -0.57499
       0.00000   0.00000  -0.00000   0.00000  -0.70600  -0.70821   0.00000

    ncont =  3
    ac_e =

       0.00000   0.00000   0.00000   0.00000  -1.26270   0.43340   0.46660
       0.00000   2.00000   0.00000  -3.74170  -0.85200   0.29240  -0.43420
       0.00000   0.00000   1.78620   0.37800  -0.26510  -0.77230   0.00000
       0.00000   0.00000   0.00000   3.74170   0.85200  -0.29240   0.43420
       0.00000   0.00000   0.00000   0.00000  -1.55400   0.53340   0.57420
       0.00000   0.00000   0.00000   0.00000  -0.65330   0.22420   0.24140
       0.00000   0.00000   0.00000   0.00000  -0.58920   0.20220   0.21770

    ec_e =

      -1.83250   1.00000   2.37520   0.00000  -0.82140   0.28190   1.80160
       0.48870   0.00000   0.37700  -0.53450   0.18740   0.54610   0.00000
      -0.17280   0.00000  -0.13330  -1.13390   0.13250   0.38610   0.00000
       0.00000   0.00000   0.00000   0.00000   0.85200  -0.29240   0.43420
       0.00000   0.00000   0.00000   0.00000  -1.02600  -0.14960   0.00000
       0.00000   0.00000   0.00000   0.00000   0.00000   1.19370   0.00000
       0.00000   0.00000   0.00000   0.00000   0.00000   0.00000   1.00000

    bc_e =

       1   2   3
       2   1   0
       0   0   0
       0   0   0
       0   0   0
       0   0   0
       0   0   0

    cc_e =

       0.00000   1.00000   0.00000   0.00000  -1.26270   0.43340   0.46660
       0.36650   0.00000  -0.98030  -1.60360   0.18740   0.54610   0.00000

    q_e =

       0.00000   1.00000   0.00000   0.00000   0.00000   0.00000   0.00000
       0.00000   0.00000   0.70710   0.00000   0.27400  -0.65190   0.00000
       0.00000   0.00000   0.00000   0.00000   0.83040   0.34910  -0.43420
       0.00000   0.00000   0.00000  -1.00000   0.00000   0.00000   0.00000
       0.00000   0.00000   0.00000   0.00000   0.40030   0.16830   0.90080
       0.00000   0.00000   0.70710   0.00000  -0.27400   0.65190   0.00000
       1.00000   0.00000   0.00000   0.00000   0.00000   0.00000   0.00000

    z_e =

       0.00000   1.00000   0.00000   0.00000   0.00000   0.00000   0.00000
      -0.61080   0.00000   0.79170   0.00000   0.00000   0.00000   0.00000
       0.48870   0.00000   0.37700  -0.53450   0.18740   0.54610   0.00000
       0.00000   0.00000   0.00000   0.00000  -0.41070   0.14100   0.90080
       0.61080   0.00000   0.47130   0.26730  -0.18740  -0.54610   0.00000
      -0.12220   0.00000  -0.09430  -0.80180  -0.18740  -0.54610   0.00000
       0.00000   0.00000   0.00000   0.00000  -0.85200   0.29240  -0.43420

    ncont_e =  3
+ ## LTI methods
+ test ('@lti/c2d')
PASSES 5 out of 5 tests
+ test ('@lti/d2c')
PASSES 3 out of 3 tests
+ test ('@lti/feedback')
PASSES 5 out of 5 tests
+ test ('@lti/horzcat')
PASSES 1 out of 1 tests
+ test ('@lti/inv')
PASSES 1 out of 1 tests
+ test ('@lti/minreal')
  ***** assert (Ar, Ae, 1e-4);
!!!!! test failed
assert (Ar,Ae,1e-4) expected
   1.00000  -0.03930  -0.09800  -0.10660   0.07810  -0.23300   0.07770
   0.00000   1.03120   0.27170   0.26090  -0.15330   0.67580  -0.35530
   0.00000   0.00000   1.38870   0.66990  -0.42810   1.63890  -0.76150
   0.00000   0.00000   0.00000  -1.21470   0.24230  -0.97920   0.47880
   0.00000   0.00000   0.00000   0.00000  -1.05450   0.50350  -0.27880
   0.00000   0.00000   0.00000   0.00000   0.00000   1.63550  -0.43230
   0.00000   0.00000   0.00000   0.00000   0.00000   0.00000   1.00000
but got
   1.00000   0.03935   0.09802  -0.10656   0.07805  -0.23302   0.07773
   0.00000   1.03122   0.27169  -0.26093   0.15327  -0.67578   0.35529
   0.00000   0.00000   1.38869  -0.66991   0.42806  -1.63894   0.76153
   0.00000   0.00000   0.00000   1.21467  -0.24232   0.97919  -0.47884
   0.00000   0.00000   0.00000   0.00000   1.05454  -0.50348   0.27881
   0.00000   0.00000   0.00000   0.00000   0.00000   1.63547  -0.43234
   0.00000   0.00000   0.00000   0.00000   0.00000   0.00000   1.00000
maximum absolute error 3.27784 exceeds tolerance 0.0001
shared variables 
  scalar structure containing the fields:

    Ar =

       1.00000   0.03935   0.09802  -0.10656   0.07805  -0.23302   0.07773
       0.00000   1.03122   0.27169  -0.26093   0.15327  -0.67578   0.35529
       0.00000   0.00000   1.38869  -0.66991   0.42806  -1.63894   0.76153
       0.00000   0.00000   0.00000   1.21467  -0.24232   0.97919  -0.47884
       0.00000   0.00000   0.00000   0.00000   1.05454  -0.50348   0.27881
       0.00000   0.00000   0.00000   0.00000   0.00000   1.63547  -0.43234
       0.00000   0.00000   0.00000   0.00000   0.00000   0.00000   1.00000

    Br =

      -5.5968e-01   2.3631e-01
       4.8427e-01   4.9770e-02
       4.7273e-01   1.4915e-01
      -1.8019e-01  -1.1574e+00
      -5.9948e-01  -1.5560e-01
      -1.7294e-01  -3.9991e-01
       1.1147e-17   2.5000e-01

    Cr =

       0.00000   0.00000   0.00000   0.00000   0.00000   0.00000   4.00000
       0.00000   0.00000   0.00000   0.00000   0.00000   3.15238  -1.75000

    Dr =

       0   0
       0   0

    Er =

       0.40999  -0.25895  -0.50798  -0.31085   0.07054   0.14290  -0.14769
       0.76292  -0.34645   0.09918   0.30065  -0.06189  -0.24831   0.01523
      -0.11198  -0.21245  -0.41840   0.12877  -0.05690   0.42134   0.61817
       0.00000   0.11215  -0.00387  -0.27713   0.07584  -0.09747  -0.39231
       0.00000   0.00000   0.37082   0.42899  -0.10064  -0.14015   0.26995
       0.00000   0.00000   0.00000   0.00000   0.94581  -0.22115   0.23779
       0.00000   0.00000   0.00000   0.57107   0.26478   0.59479  -0.50000

    Ae =

       1.00000  -0.03930  -0.09800  -0.10660   0.07810  -0.23300   0.07770
       0.00000   1.03120   0.27170   0.26090  -0.15330   0.67580  -0.35530
       0.00000   0.00000   1.38870   0.66990  -0.42810   1.63890  -0.76150
       0.00000   0.00000   0.00000  -1.21470   0.24230  -0.97920   0.47880
       0.00000   0.00000   0.00000   0.00000  -1.05450   0.50350  -0.27880
       0.00000   0.00000   0.00000   0.00000   0.00000   1.63550  -0.43230
       0.00000   0.00000   0.00000   0.00000   0.00000   0.00000   1.00000

    Be =

      -0.55970   0.23630
      -0.48430  -0.04980
      -0.47270  -0.14910
       0.18020   1.15740
       0.59950   0.15560
      -0.17290  -0.39990
       0.00000   0.25000

    Ce =

       0.00000   0.00000   0.00000   0.00000   0.00000   0.00000   4.00000
       0.00000   0.00000   0.00000   0.00000   0.00000   3.15240  -1.75000

    De =

       0   0
       0   0

    Ee =

       0.41000   0.25900   0.50800  -0.31090   0.07050   0.14290  -0.14770
      -0.76290  -0.34640   0.09920  -0.30070   0.06190   0.24830  -0.01520
       0.11200  -0.21240  -0.41840  -0.12880   0.05690  -0.42130  -0.61820
       0.00000   0.11220  -0.00390   0.27710  -0.07580   0.09750   0.39230
       0.00000   0.00000   0.37080  -0.42900   0.10060   0.14020  -0.26990
       0.00000   0.00000   0.00000   0.00000   0.94580  -0.22110   0.23780
       0.00000   0.00000   0.00000   0.57110   0.26480   0.59480  -0.50000

+ test ('@lti/mtimes')
PASSES 1 out of 1 tests
+ test ('@lti/norm')
PASSES 4 out of 4 tests
+ test ('@lti/plus')
PASSES 1 out of 1 tests
+ test ('@lti/prescale')
PASSES 22 out of 22 tests
+ test ('@lti/sminreal')
PASSES 12 out of 12 tests
+ test ('@lti/zero')
PASSES 8 out of 8 tests
+ ## robust control
+ test ('h2syn')
PASSES 2 out of 2 tests
+ test ('hinfsyn')
PASSES 2 out of 2 tests
+ test ('ncfsyn')
PASSES 27 out of 27 tests
+ ## ARE solvers
+ test ('care')
PASSES 7 out of 7 tests
+ test ('dare')
PASSES 3 out of 3 tests
+ test ('kalman')
PASSES 3 out of 3 tests
+ ## Lyapunov
+ test ('covar')
PASSES 4 out of 4 tests
+ test ('dlyap')
PASSES 2 out of 2 tests
+ ## test dlyapchol  # TODO: add tests
+ test ('gram')
PASSES 5 out of 5 tests
+ test ('lyap')
PASSES 3 out of 3 tests
+ test ('lyapchol')
PASSES 3 out of 3 tests
+ ## model order reduction
+ test ('bstmodred')
PASSES 2 out of 2 tests
+ test ('btamodred')
PASSES 2 out of 2 tests
+ test ('hnamodred')
PASSES 2 out of 2 tests
+ ## test spamodred  # TODO: create test case

+ ## controller order reduction
+ test ('btaconred')
PASSES 2 out of 2 tests
+ test ('cfconred')
PASSES 2 out of 2 tests
+ test ('fwcfconred')
PASSES 2 out of 2 tests
+ ## test spaconred  # TODO: create test case

+ ## identification
+ test ('fitfrd')
PASSES 1 out of 1 tests
+ ## various oct-files
+ test ('ctrbf')
PASSES 4 out of 4 tests
+ test ('hsvd')
PASSES 1 out of 1 tests
+ test ('place')
PASSES 4 out of 4 tests
+ ## various m-files
+ test ('ctrb')
PASSES 1 out of 1 tests
+ test ('filt')
PASSES 4 out of 4 tests
+ test ('initial')
PASSES 2 out of 2 tests
+ test ('issample')
PASSES 30 out of 30 tests
+ test ('margin')
PASSES 2 out of 2 tests
+ test ('obsv')
PASSES 1 out of 1 tests
+ test ('sigma')
PASSES 2 out of 2 tests
+ endscript;
octave:4> exit

------------------------------------------------------------------------------
Live Security Virtual Conference
Exclusive live event will cover all the ways today's security and 
threat landscape has changed and how IT managers can respond. Discussions 
will include endpoint security, mobile security and the latest in malware 
threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
_______________________________________________
Octave-dev mailing list
Octave-dev@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/octave-dev

Reply via email to