[Octave-bug-tracker] [bug #42097] 64bit Slackware, eig returning incorrect results

[Jordi Gutiérrez Hermoso]

Follow-up Comment #4, bug #42097 (project octave):

Can you use format long and come up with an example again? When I attempt your
example with the precision you provided, 


octave:10> A = [   1.53587   0.29450   0.15366
>    0.60792  -1.93114   0.67947
>   -0.31746   1.29087  -1.14817]
A =

   1.53587   0.29450   0.15366
   0.60792  -1.93114   0.67947
  -0.31746   1.29087  -1.14817

octave:11> A
A =

   1.53587   0.29450   0.15366
   0.60792  -1.93114   0.67947
  -0.31746   1.29087  -1.14817

octave:12> [V, D] = eig(A)
V =

   0.985834   0.027681  -0.122808
   0.163509  -0.738089   0.393531
  -0.037364   0.674136   0.911072

D =

Diagonal Matrix

   1.57889         0         0
         0  -2.57454         0
         0         0  -0.54780

octave:13> A*V - V*D
ans =

  -1.1102e-15   2.2204e-16   5.5511e-17
   1.1102e-16   1.1102e-15  -3.0531e-16
   2.0817e-17   0.0000e+00   5.5511e-16


which is as expected.

    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?42097>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/