;;; calc-arith.el --- arithmetic functions for Calc
;; Copyright (C) 1990, 1991, 1992, 1993, 2001, 2002, 2003, 2004,
-;; 2005 Free Software Foundation, Inc.
+;; 2005, 2006, 2007, 2008 Free Software Foundation, Inc.
;; Author: David Gillespie <daveg@synaptics.com>
-;; Maintainer: Jay Belanger <belanger@truman.edu>
+;; Maintainer: Jay Belanger <jay.p.belanger@gmail.com>
;; This file is part of GNU Emacs.
+;; GNU Emacs is free software: you can redistribute it and/or modify
+;; it under the terms of the GNU General Public License as published by
+;; the Free Software Foundation, either version 3 of the License, or
+;; (at your option) any later version.
+
;; GNU Emacs is distributed in the hope that it will be useful,
-;; but WITHOUT ANY WARRANTY. No author or distributor
-;; accepts responsibility to anyone for the consequences of using it
-;; or for whether it serves any particular purpose or works at all,
-;; unless he says so in writing. Refer to the GNU Emacs General Public
-;; License for full details.
-
-;; Everyone is granted permission to copy, modify and redistribute
-;; GNU Emacs, but only under the conditions described in the
-;; GNU Emacs General Public License. A copy of this license is
-;; supposed to have been given to you along with GNU Emacs so you
-;; can know your rights and responsibilities. It should be in a
-;; file named COPYING. Among other things, the copyright notice
-;; and this notice must be preserved on all copies.
+;; but WITHOUT ANY WARRANTY; without even the implied warranty of
+;; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+;; GNU General Public License for more details.
+
+;; You should have received a copy of the GNU General Public License
+;; along with GNU Emacs. If not, see <http://www.gnu.org/licenses/>.
;;; Commentary:
(real number)
(number)
(scalar)
+ (sqmatrix matrix vector)
(matrix vector)
(vector)
(const)))
(not (math-known-scalarp a t))))
(defun math-known-square-matrixp (a)
- (if (eq (car-safe a) '^)
- (math-known-square-matrixp (nth 1 a))
- (and (math-known-matrixp a)
- (or (math-square-matrixp a)
- (and (or
- (integerp calc-matrix-mode)
- (eq calc-matrix-mode 'square))
- (eq (car-safe a) 'var)
- (not (math-const-var a)))))))
-
+ (and (math-known-matrixp a)
+ (math-check-known-square-matrixp a)))
+
;;; Try to prove that A is a scalar (i.e., a non-vector).
(defun math-check-known-scalarp (a)
(cond ((Math-objectp a) t)
(let ((decl (if (eq (car a) 'var)
(or (assq (nth 2 a) math-decls-cache)
math-decls-all)
- (assq (car a) math-decls-cache))))
- (memq 'scalar (nth 1 decl))))))
+ (assq (car a) math-decls-cache)))
+ val)
+ (cond
+ ((memq 'scalar (nth 1 decl))
+ t)
+ ((and (eq (car a) 'var)
+ (symbolp (nth 2 a))
+ (boundp (nth 2 a))
+ (setq val (symbol-value (nth 2 a))))
+ (math-check-known-scalarp val))
+ (t
+ nil))))))
;;; Try to prove that A is *not* a scalar.
(defun math-check-known-matrixp (a)
(let ((decl (if (eq (car a) 'var)
(or (assq (nth 2 a) math-decls-cache)
math-decls-all)
- (assq (car a) math-decls-cache))))
- (memq 'vector (nth 1 decl))))))
-
+ (assq (car a) math-decls-cache)))
+ val)
+ (cond
+ ((memq 'matrix (nth 1 decl))
+ t)
+ ((and (eq (car a) 'var)
+ (symbolp (nth 2 a))
+ (boundp (nth 2 a))
+ (setq val (symbol-value (nth 2 a))))
+ (math-check-known-matrixp val))
+ (t
+ nil))))))
+
+;;; Given that A is a matrix, try to prove that it is a square matrix.
+(defun math-check-known-square-matrixp (a)
+ (cond ((math-square-matrixp a)
+ 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)
+ math-decls-all)
+ (assq (car a) math-decls-cache)))
+ val)
+ (cond
+ ((memq 'sqmatrix (nth 1 decl))
+ t)
+ ((and (eq (car a) 'var)
+ (boundp (nth 2 a))
+ (setq val (symbol-value (nth 2 a))))
+ (math-check-known-square-matrixp val))
+ ((and (or
+ (integerp calc-matrix-mode)
+ (eq calc-matrix-mode 'sqmatrix))
+ (eq (car-safe a) 'var))
+ t)
+ ((memq 'matrix (nth 1 decl))
+ nil)
+ (t
+ nil))))))
;;; Try to prove that A is a real (i.e., not complex).
(defun math-known-realp (a)
(and (math-known-scalarp b)
(math-add (nth 1 a) b))))
(and (eq (car-safe b) 'calcFunc-idn)
- (= (length a) 2)
+ (= (length b) 2)
(or (and (math-square-matrixp a)
(math-add a (math-mimic-ident (nth 1 b) a)))
(and (math-known-scalarp a)
(and (eq (car-safe b) '^)
(Math-looks-negp (nth 2 b))
(not (and (eq (car-safe a) '^) (Math-looks-negp (nth 2 a))))
+ (not (math-known-matrixp (nth 1 b)))
(math-div a (math-normalize
(list '^ (nth 1 b) (math-neg (nth 2 b))))))
(and (eq (car-safe a) '/)
(list 'calcFunc-idn (math-mul a (nth 1 b))))
(and (math-known-matrixp a)
(math-mul a (nth 1 b)))))
+ (and (math-identity-matrix-p a t)
+ (or (and (eq (car-safe b) 'calcFunc-idn)
+ (= (length b) 2)
+ (list 'calcFunc-idn (math-mul
+ (nth 1 (nth 1 a))
+ (nth 1 b))
+ (1- (length a))))
+ (and (math-known-scalarp b)
+ (list 'calcFunc-idn (math-mul
+ (nth 1 (nth 1 a)) b)
+ (1- (length a))))
+ (and (math-known-matrixp b)
+ (math-mul (nth 1 (nth 1 a)) b))))
+ (and (math-identity-matrix-p b t)
+ (or (and (eq (car-safe a) 'calcFunc-idn)
+ (= (length a) 2)
+ (list 'calcFunc-idn (math-mul (nth 1 a)
+ (nth 1 (nth 1 b)))
+ (1- (length b))))
+ (and (math-known-scalarp a)
+ (list 'calcFunc-idn (math-mul a (nth 1 (nth 1 b)))
+ (1- (length b))))
+ (and (math-known-matrixp a)
+ (math-mul a (nth 1 (nth 1 b))))))
(and (math-looks-negp b)
(math-mul (math-neg a) (math-neg b)))
(and (eq (car-safe b) '-)
(math-div-new-non-trig term))))
(defun math-div-symb-fancy (a b)
- (or (and math-simplify-only
+ (or (and (math-known-matrixp b)
+ (math-mul a (math-pow b -1)))
+ (and math-simplify-only
(not (equal a math-simplify-only))
(list '/ a b))
(and (Math-equal-int b 1) a)
(math-mul-zero b a))))
(list '/ a b)))
+;;; Division from the left.
+(defun calcFunc-ldiv (a b)
+ (if (math-known-scalarp a)
+ (math-div b a)
+ (math-mul (math-pow a -1) b)))
(defun calcFunc-mod (a b)
(math-normalize (list '% a b)))
(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)
(provide 'calc-arith)
-;;; arch-tag: 6c396b5b-14c6-40ed-bb2a-7cc2e8111465
+;; arch-tag: 6c396b5b-14c6-40ed-bb2a-7cc2e8111465
;;; calc-arith.el ends here