From: Joel VanderWerf <vjoel@...>
Date: 2001-12-16T19:50:55+09:00
Subject: [ruby-talk:28668] Re: Matrix#inverse broken? -- NOT

Joel VanderWerf wrote:
> 
> require 'matrix'
> m = Matrix[[1,2],[2,1]]
> m * (m.inv)
>   # ==> Matrix[[-1, -1], [-2, 1]]

m = Matrix[[1.0,2.0],[2.0,1.0]]
p m * (m.inv)
  # ==> Matrix[[1.0, 0.0], [1.110223025e-16, 1.0]]

It's me that's broken (and my FPU ain't perfect either).

--
Joel VanderWerf                          California PATH, UC Berkeley
mailto:vjoel@path.berkeley.edu                     Ph. (510) 231-9446
http://www.path.berkeley.edu                       FAX (510) 231-9512