From: Jay Belanger <jay.p.belanger@gmail.com>
Date: Wed, 9 Nov 2005 19:27:29 +0000 (+0000)
Subject: (math-pow-fancy):  Further expand product of square matrices.
X-Git-Tag: emacs-pretest-22.0.90~5992
X-Git-Url: http://git.eshelyaron.com/gitweb/?a=commitdiff_plain;h=2f884e83d44f32270309447d3439671898142182;p=emacs.git

(math-pow-fancy):  Further expand product of square matrices.
(math-check-known-square-matrixp): Do a more extensive check for
square matrices.
---

diff --git a/lisp/calc/calc-arith.el b/lisp/calc/calc-arith.el
index f8057c5f1b9..ba95ee7f713 100644
--- a/lisp/calc/calc-arith.el
+++ b/lisp/calc/calc-arith.el
@@ -374,6 +374,13 @@
          t)
         ((eq (car-safe a) '^)
          (math-check-known-square-matrixp (nth 1 a)))
+        ((or
+          (eq (car-safe a) '*)
+          (eq (car-safe a) '+)
+          (eq (car-safe a) '-))
+         (and
+          (math-check-known-square-matrixp (nth 1 a))
+          (math-check-known-square-matrixp (nth 2 a))))
         (t
          (let ((decl (if (eq (car a) 'var)
                          (or (assq (nth 2 a) math-decls-cache)
@@ -1960,7 +1967,8 @@
                   (if (and (= b -1)
                            (math-known-square-matrixp (nth 1 a))
                            (math-known-square-matrixp (nth 2 a)))
-                      (list '* (list '^ (nth 2 a) -1) (list '^ (nth 1 a) -1))
+                      (math-mul (math-pow-fancy (nth 2 a) -1) 
+                                (math-pow-fancy (nth 1 a) -1))
                     (list '^ a b)))
 		 ((and (eq (car-safe a) '*)
 		       (or (math-known-num-integerp b)