diff options
Diffstat (limited to 'lib/matrix.rb')
-rw-r--r-- | lib/matrix.rb | 9 |
1 files changed, 6 insertions, 3 deletions
diff --git a/lib/matrix.rb b/lib/matrix.rb index 8c6b3d6e87..5ce130221f 100644 --- a/lib/matrix.rb +++ b/lib/matrix.rb @@ -956,8 +956,10 @@ class Matrix private :inverse_from # - # Matrix exponentiation. Currently implemented for integer powers only. + # Matrix exponentiation. # Equivalent to multiplying the matrix by itself N times. + # Non integer exponents will be handled by diagonalizing the matrix. + # # Matrix[[7,6], [3,9]] ** 2 # => 67 96 # 48 99 @@ -977,8 +979,9 @@ class Matrix return z if (other >>= 1).zero? x *= x end - when Float, Rational - Matrix.Raise ErrOperationNotImplemented, "**", self.class, other.class + when Numeric + v, d, v_inv = eigensystem + v * Matrix.diagonal(*d.each(:diagonal).map{|e| e ** other}) * v_inv else Matrix.Raise ErrOperationNotDefined, "**", self.class, other.class end |