[Octave-bug-tracker] [bug #55256] Wrong result raise a non-diagonalizable matrix by a fractional exponent

[A.R. Burgers]

Follow-up Comment #2, bug #55256 (project octave):

Here is what matlab does:


>> version
ans = '9.4.0.813654 (R2018a)'

>> A=[1 1 1; 0 1 1; 0 0 1];

>> A^0.5
Warning: Matrix is close to singular or badly scaled. Results may be
inaccurate. RCOND = 4.930381e-32. 
> In matlab.internal.math.mpower.viaEig (line 10)
 
ans =

     1     0     0
     0     1     0
     0     0     1



    _______________________________________________________

Reply to this item at:

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

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