From ad162055c41cbcb279705480c8405219688e34f2 Mon Sep 17 00:00:00 2001 From: marcandre Date: Sun, 13 Jan 2013 05:50:34 +0000 Subject: * lib/matrix/eigenvalue_decomposition: Fix eigensystem with complex eigenvectors. Patch by pypypy567. [Bug #7208] [ruby-dev:46251] [rubyspec:242f8e55bd] git-svn-id: svn+ssh://ci.ruby-lang.org/ruby/trunk@38799 b2dd03c8-39d4-4d8f-98ff-823fe69b080e --- ChangeLog | 6 ++++++ lib/matrix/eigenvalue_decomposition.rb | 4 ++-- 2 files changed, 8 insertions(+), 2 deletions(-) diff --git a/ChangeLog b/ChangeLog index f53ac70a41..f01befbe71 100644 --- a/ChangeLog +++ b/ChangeLog @@ -1,3 +1,9 @@ +Sun Jan 13 14:48:55 2013 Marc-Andre Lafortune + + * lib/matrix/eigenvalue_decomposition: Fix eigensystem with complex + eigenvectors. Patch by pypypy567. + [Bug #7208] [ruby-dev:46251] [rubyspec:242f8e55bd] + Sun Jan 13 14:06:00 2013 Zachary Scott * lib/irb.rb, lib/prime.rb: Typos in overview diff --git a/lib/matrix/eigenvalue_decomposition.rb b/lib/matrix/eigenvalue_decomposition.rb index ab5888ff37..1949406687 100644 --- a/lib/matrix/eigenvalue_decomposition.rb +++ b/lib/matrix/eigenvalue_decomposition.rb @@ -84,14 +84,14 @@ class Matrix private def build_eigenvectors # JAMA stores complex eigenvectors in a strange way - # See http://cio.nist.gov/esd/emaildir/lists/jama/msg01021.html + # See http://web.archive.org/web/20111016032731/http://cio.nist.gov/esd/emaildir/lists/jama/msg01021.html @e.each_with_index.map do |imag, i| if imag == 0 Array.new(@size){|j| @v[j][i]} elsif imag > 0 Array.new(@size){|j| Complex(@v[j][i], @v[j][i+1])} else - Array.new(@size){|j| Complex(@v[j][i], -@v[j][i-1])} + Array.new(@size){|j| Complex(@v[j][i-1], -@v[j][i])} end end end -- cgit v1.2.3