[1] | 1 | head 1.3;
|
---|
| 2 | access;
|
---|
| 3 | symbols;
|
---|
| 4 | locks; strict;
|
---|
| 5 | comment @;;; @;
|
---|
| 6 |
|
---|
| 7 |
|
---|
| 8 | 1.3
|
---|
| 9 | date 2009.01.22.04.04.31; author marek; state Exp;
|
---|
| 10 | branches;
|
---|
| 11 | next 1.2;
|
---|
| 12 |
|
---|
| 13 | 1.2
|
---|
| 14 | date 2009.01.19.09.27.11; author marek; state Exp;
|
---|
| 15 | branches;
|
---|
| 16 | next 1.1;
|
---|
| 17 |
|
---|
| 18 | 1.1
|
---|
| 19 | date 2009.01.19.07.38.34; author marek; state Exp;
|
---|
| 20 | branches;
|
---|
| 21 | next ;
|
---|
| 22 |
|
---|
| 23 |
|
---|
| 24 | desc
|
---|
| 25 | @@
|
---|
| 26 |
|
---|
| 27 |
|
---|
| 28 | 1.3
|
---|
| 29 | log
|
---|
| 30 | @*** empty log message ***
|
---|
| 31 | @
|
---|
| 32 | text
|
---|
| 33 | @#|
|
---|
| 34 | $Id$
|
---|
| 35 | *--------------------------------------------------------------------------*
|
---|
| 36 | | Copyright (C) 1994, Marek Rychlik (e-mail: rychlik@@math.arizona.edu) |
|
---|
| 37 | | Department of Mathematics, University of Arizona, Tucson, AZ 85721 |
|
---|
| 38 | | |
|
---|
| 39 | | Everyone is permitted to copy, distribute and modify the code in this |
|
---|
| 40 | | directory, as long as this copyright note is preserved verbatim. |
|
---|
| 41 | *--------------------------------------------------------------------------*
|
---|
| 42 | |#
|
---|
| 43 | (defpackage "MONOM"
|
---|
| 44 | (:use "COMMON-LISP")
|
---|
| 45 | (:export monom-equal monom-divides-p monom-divisible-by-p monom-rel-prime
|
---|
| 46 | monom* nmonom* monom/
|
---|
| 47 | monom-lcm monom-gcd))
|
---|
| 48 |
|
---|
| 49 | (in-package "MONOM")
|
---|
| 50 |
|
---|
| 51 | #+debug(proclaim '(optimize (speed 0) (debug 3)))
|
---|
| 52 | #-debug(proclaim '(optimize (speed 3) (debug 0)))
|
---|
| 53 |
|
---|
| 54 | ;;----------------------------------------------------------------
|
---|
| 55 | ;; This package implements BASIC OPERATIONS ON MONOMIALS
|
---|
| 56 | ;;----------------------------------------------------------------
|
---|
| 57 | ;; DATA STRUCTURES: Monomials are represented as lists:
|
---|
| 58 | ;;
|
---|
| 59 | ;; monom: (n1 n2 ... nk) where ni are non-negative integers
|
---|
| 60 | ;;
|
---|
| 61 | ;;----------------------------------------------------------------
|
---|
| 62 | ;; EXAMPLES: Suppose that variables are x and y. Then
|
---|
| 63 | ;;
|
---|
| 64 | ;; Monom x*y^2 ---> (1 2)
|
---|
| 65 | ;;
|
---|
| 66 | ;;----------------------------------------------------------------
|
---|
| 67 |
|
---|
| 68 | (defun monom/ (m1 m2)
|
---|
| 69 | "Divide monomial M1 by monomial M2."
|
---|
| 70 | (mapcar #'- m1 m2))
|
---|
| 71 |
|
---|
| 72 | (defun monom* (m1 m2)
|
---|
| 73 | "Multiply monomial M1 by monomial M2."
|
---|
| 74 | (mapcar #'+ m1 m2))
|
---|
| 75 |
|
---|
| 76 | (defun nmonom* (m1 m2)
|
---|
| 77 | "Multiply monomials M1 and M2 - destructive version.
|
---|
| 78 | M1 is destructively modified, M2 is not modified."
|
---|
| 79 | (mapl #'(lambda (x y) (incf (car x) (car y))) m1 m2))
|
---|
| 80 |
|
---|
| 81 | (defun monom-divides-p (m1 m2)
|
---|
| 82 | "Returns T if monomial M1 divides monomial M2, NIL otherwise."
|
---|
| 83 | (every #'<= m1 m2))
|
---|
| 84 |
|
---|
| 85 | (defun monom-divisible-by-p (m1 m2)
|
---|
| 86 | "Returns T if monomial M1 is divisible by monomial M2, NIL otherwise."
|
---|
| 87 | (every #'>= m1 m2))
|
---|
| 88 |
|
---|
| 89 | (defun monom-rel-prime (m1 m2)
|
---|
| 90 | "Returns T if two monomials M1 and M2 are relatively prime (disjoint)."
|
---|
| 91 | (every #'(lambda (x y) (zerop (min x y))) m1 m2))
|
---|
| 92 |
|
---|
| 93 | (defun monom-equal (m1 m2)
|
---|
| 94 | "Returns T if two monomials M1 and M2 are equal."
|
---|
| 95 | (every #'= m1 m2))
|
---|
| 96 |
|
---|
| 97 | (defun monom-lcm (m1 m2)
|
---|
| 98 | "Returns least common multiple of monomials M1 and M2."
|
---|
| 99 | (mapcar #'max m1 m2))
|
---|
| 100 |
|
---|
| 101 | (defun monom-gcd (m1 m2)
|
---|
| 102 | "Returns greatest common divisor of monomials M1 and M2."
|
---|
| 103 | (mapcar #'min m1 m2))
|
---|
| 104 | @
|
---|
| 105 |
|
---|
| 106 |
|
---|
| 107 | 1.2
|
---|
| 108 | log
|
---|
| 109 | @*** empty log message ***
|
---|
| 110 | @
|
---|
| 111 | text
|
---|
| 112 | @d19 2
|
---|
| 113 | a20 2
|
---|
| 114 | ;;(proclaim '(optimize (speed 0) (debug 3)))
|
---|
| 115 | (proclaim '(optimize (speed 3) (debug 0)))
|
---|
| 116 | @
|
---|
| 117 |
|
---|
| 118 |
|
---|
| 119 | 1.1
|
---|
| 120 | log
|
---|
| 121 | @Initial revision
|
---|
| 122 | @
|
---|
| 123 | text
|
---|
| 124 | @d2 1
|
---|
| 125 | a2 1
|
---|
| 126 | $Id: monom.lisp,v 1.13 1997/12/02 23:41:29 marek Exp $
|
---|
| 127 | d19 2
|
---|
| 128 | a20 1
|
---|
| 129 | (proclaim '(optimize (speed 0) (debug 3)))
|
---|
| 130 | @
|
---|