diff options
| -rw-r--r-- | ChangeLog | 6 | ||||
| -rw-r--r-- | lib/matrix/eigenvalue_decomposition.rb | 4 |
2 files changed, 8 insertions, 2 deletions
@@ -1,3 +1,9 @@ +Sun Jan 13 14:48:55 2013 Marc-Andre Lafortune <ruby-core@marc-andre.ca> + + * 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 <zachary@zacharyscott.net> * 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 |