From: Philip N. <pr....@hc...> - 2012-09-03 10:30:22
|
PMFJI but, uhm.... Lukas Reichlin wrote: > On 03.09.2012, at 10:39, Rafael Laboissiere<ra...@la...> wrote: > >> * Rafael Laboissiere<ra...@la...> [2012-08-30 23:47]: >> >>> * Carnë Draug<car...@gm...> [2012-08-27 16:43]: >>> >>>> a new release of control package is out, version 2.3.53, by Lukas Reichlin. >>> >>> In exercising the tests in inst/@lti/minreal.m, I got the error below. >>> Is it normal? > Hi Rafael, > > You can check whether the observed and expected results are equivalent state-space models (i.e. state-transformation, see command prescale for formulae). This can be done, e.g., by inspection of the Hankel singular values (command hsvd), time response (step, impulse) or frequency response (sigma). > If they are the same, there should be nothing to worry about. If you want the same results, use Reference BLAS (and LAPACK) from www.netlib.org instead of ATLAS which you are probably using. The SLICOT authors recommend the use of the reference implementations. Correct results are more important than minor speed advantages of automatically tuned linear algebra software, aren't they? :-) ... shouldn't such an interpretation be reflected in the tests then, rather than a plain "assert" comparison that depends on what (allowed!) dependencies happen to be in place. Currently the test just seems to yield confusion. Philip |