[Octave-bug-tracker] [bug #63052] [octave forge] (control) Inverse operation with transfer function variables resulting in wrong results

[Doug Stewart]

Follow-up Comment #6, bug #63052 (project octave):

I found that ,in __zero__.m
this line:
[zer, rank, infz, kronr, kronl] = __sl_ag08bd__ (sys.a, sys.e, sys.b, sys.c,
sys.d, sys.scaled);

is returning zer close to eps, when it should return 0;

I added:

if(abs(zer)<2*eps)
     zer=0;
   endif

This definitely reduced (removed) the spurious wrong numbers, that I was
seeing.

So the FORTRAN function __sl_ag08bd__  seems to have some round off problems.




    _______________________________________________________

Reply to this item at:

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

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