close Warning: Can't synchronize with repository "(default)" (The repository directory has changed, you should resynchronize the repository with: trac-admin $ENV repository resync '(default)'). Look in the Trac log for more information.

source: branches/f4grobner/ring.lisp@ 3006

Last change on this file since 3006 was 3006, checked in by Marek Rychlik, 9 years ago

* empty log message *

File size: 8.0 KB
Line 
1;;; -*- mode: lisp; package: maxima; syntax: common-lisp; base: 10 -*-
2;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
3;;;
4;;; copyright (c) 1999, 2002, 2009, 2015 marek rychlik <rychlik@u.arizona.edu>
5;;;
6;;; this program is free software; you can redistribute it and/or modify
7;;; it under the terms of the gnu general public license as published by
8;;; the free software foundation; either version 2 of the license, or
9;;; (at your option) any later version.
10;;;
11;;; this program is distributed in the hope that it will be useful,
12;;; but without any warranty; without even the implied warranty of
13;;; merchantability or fitness for a particular purpose. see the
14;;; gnu general public license for more details.
15;;;
16;;; you should have received a copy of the gnu general public license
17;;; along with this program; if not, write to the free software
18;;; foundation, inc., 59 temple place - suite 330, boston, ma 02111-1307, usa.
19;;;
20;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
21
22(defpackage "RING"
23 (:use :cl)
24 (:export "R-PARSE"
25 "UNIT-ELEMENT"
26 "R-ZEROP"
27 "R+"
28 "R-"
29 "R*"
30 "R+"
31 "R/"
32 "R-EXPT"
33 "R-LCM"
34 "R-EZGCD"
35 "R-GCD"
36 "R-TOTAL-DEGREE"
37 "R-DIMENSION"
38 "R-EXPONENTS"
39 "R-COEFF"
40 "R-SUGAR"
41 "R-DIVIDES-P"
42 "R-DIVIDES-LCM-P"
43 "R-LCM-DIVIDES-LCM-P"
44 "R-LCM-EQUAL-LCM-P"
45 "R-REL-PRIME-P"
46 "R-EQUALP"
47 "R-CLONE"
48 "R-ELT"
49 "R->LIST"
50 "R-DIVISIBLE-BY-P"
51 "R-REL-PRIME-P"
52 "R-DEPENDS-P"
53 "R-TENSOR-PRODUCT"
54 "R-CONTRACT"
55 "R-LENGTH"
56 "MULTIPLY-BY"
57 "DIVIDE-BY"
58 "ADD-TO"
59 "SUBTRACT-FROM"
60 "UNARY-MINUS"
61 "SCALAR"
62 "INSERT-ITEM"
63 "APPEND-ITEM"
64 "COPY-INSTANCE")
65 (:shadowing-import-from
66 #+openmcl-native-threads #:ccl
67 #+cmu #:pcl
68 #+sbcl #:sb-pcl
69 #+lispworks #:hcl
70 #+allegro #:mop
71 #+clisp #:clos
72 #:class-slots #:slot-definition-name)
73 (:documentation
74 "Implements ring operations. These are all operations that are
75performed on the coefficients by the package, and thus the coefficient
76ring can be changed by merely redefining these operations."))
77
78(in-package :ring)
79
80(defclass scalar ()
81 ((value :initarg :value :accessor value))
82 (:documentation "Wraps objects suitable as scalars/polynomial coefficients"))
83
84(defgeneric unit-element (class))
85
86(defgeneric r-zerop (object)
87 (:method ((self number)) (zerop self))
88 (:documentation "Tests whether a ring element is 0."))
89
90(defgeneric r+ (x y)
91 (:method ((x number) (y number)) (+ x y))
92 (:documentation "Adds ring elements."))
93
94(defgeneric r- (x y)
95 (:method ((x number) (y number)) (- x y))
96 (:documentation "Subtracts ring elements."))
97
98(defgeneric r* (x y)
99 (:method ((x number) (y number)) (* x y))
100 (:documentation "Multiplies ring elements."))
101
102(defgeneric r-tensor-product (x y)
103 (:documentation "Takes a tensor product of two objects."))
104
105(defgeneric r/ (x y)
106 (:method ((x number) (y number)) (/ x y))
107 (:documentation "Divides ring elements."))
108
109(defgeneric r-lcm (x y)
110 (:method ((x integer) (y integer)) (lcm x y))
111 (:documentation "Returns the least common multiple of ring elements."))
112
113(defgeneric r-expt (x y)
114 (:method ((x integer) (y integer)) (expt x y))
115 (:documentation "Raises X to power Y."))
116
117(defgeneric r-ezgcd (x y)
118 (:method ((x integer) (y integer)
119 &aux (c (gcd x y)))
120 (values c (/ x c) (/ y c)))
121 (:documentation "Solves the diophantine system: X=C*X1, Y=C*X2,
122C=GCD(X,Y). It returns C, X1 and Y1. The result may be obtained by
123the Euclidean algorithm."))
124
125(defgeneric r-gcd (x y)
126 (:method ((x integer) (y integer))
127 (gcd x y))
128 (:documentation "Returns GCD(X,Y)."))
129
130(defgeneric r-dimension (object))
131(defgeneric r-exponents (object))
132
133(defgeneric r-coeff (object))
134(defgeneric (setf r-coeff) (new-value object))
135
136(defgeneric r-total-degree (object &optional start end))
137
138(defgeneric r-divides-p (object1 object2)
139 (:method ((object1 integer) (object2 integer))
140 (zerop (rem object2 object1)))
141 (:documentation "Returns T if OBJECT1 divides OBJECT2"))
142
143(defgeneric r-divides-lcm-p (object1 object2 object3)
144 (:documentation "Returns T if OBJECT divides LCM(OBJECT2, OBJECT3), NIL otherwise."))
145
146(defgeneric r-lcm-divides-lcm-p (object1 object2 object3 object4)
147 (:documentation "Returns T if LCM(OBJECT1,OBJECT2) divides LCM(OBJECT3,OBJECT4), NIL otherwise."))
148
149(defgeneric r-lcm-equal-lcm-p (object1 object2 object3 object4)
150 (:documentation "Returns T if object LCM(OBJECT1,OBJECT2) equals LCM(OBJECT3,OBJECT4), NIL otherwise."))
151
152(defgeneric r-equalp (object1 object2)
153 (:method (object1 object2) (equalp object1 object2))
154 (:documentation "Equality using deep comparison of object slots."))
155
156(defgeneric r-elt (object index))
157
158(defgeneric (setf r-elt) (new-value object index))
159
160(defgeneric r-length (object))
161
162(defgeneric r->list (object))
163(defgeneric r-sugar (object))
164(defgeneric r-rel-prime-p (object1 object2))
165(defgeneric r-contract (object k))
166(defgeneric r-divisible-by-p (object1 object2))
167(defgeneric r-depends-p (object k))
168
169(defgeneric multiply-by (self other)
170 (:method (self other) (r* self other))
171 (:documentation "Multiply object SELF and OTHER and store the result
172into SELF. It returns SELF. For instances of a class, this operation
173may be destructive."))
174
175(defgeneric divide-by (self other)
176 (:method (self other) (r/ self other))
177 (:documentation "Divided object SELF by OTHER and store the result
178into SELF. It returns SELF. For instances of a class, this operation
179may be destructive."))
180
181(defgeneric add-to (self other)
182 (:documentation "Add to object SELF another object OTHER. For
183complex objects, it may destructively modify SELF and destructively
184modify/invalidate object OTHER. For standard classes implementing this
185method, the result should be an object which is EQ to SELF. For
186built-in classes, such as NUMBER, the returned object may not be EQ to
187the original, but it will be EQL to it.")
188 (:method (self other) (r+ self other)))
189
190(defgeneric subtract-from (self other)
191 (:documentation "Subtract from an object SELF another object OTHER.
192For complex objects, it may destructively modify SELF and
193destructively modify/invalidate object OTHER. For standard classes
194implementing this method, the result should be an object which is EQ
195to SELF. For built-in classes, such as NUMBER, the returned object may
196not be EQ to the original.")
197 (:method (self other) (r- self other)))
198
199(defgeneric unary-minus (self)
200 (:method ((x number)) (- x)))
201
202(defgeneric insert-item (self item))
203(defgeneric append-item (self item))
204
205;; Source: http://stackoverflow.com/questions/11067899/is-there-a-generic-method-for-cloning-clos-objects
206(defgeneric copy-instance (object &rest initargs &key &allow-other-keys)
207 (:documentation "Makes and returns a shallow copy of OBJECT.
208
209 An uninitialized object of the same class as OBJECT is allocated by
210 calling ALLOCATE-INSTANCE. For all slots returned by
211 CLASS-SLOTS, the returned object has the
212 same slot values and slot-unbound status as OBJECT.
213
214 REINITIALIZE-INSTANCE is called to update the copy with INITARGS.")
215 (:method ((object standard-object) &rest initargs &key &allow-other-keys)
216 (let* ((class (class-of object))
217 (copy (allocate-instance class)))
218 (dolist (slot-name (mapcar #'slot-definition-name (class-slots class)))
219 (when (slot-boundp object slot-name)
220 (setf (slot-value copy slot-name)
221 (slot-value object slot-name))))
222 (apply #'reinitialize-instance copy initargs))))
223
224#|
225;; Source: http://stackoverflow.com/questions/11067899/is-there-a-generic-method-for-cloning-clos-objects
226(defun shallow-copy-object (original)
227 (let* ((class (class-of original))
228 (copy (allocate-instance class)))
229 (dolist (slot (mapcar #'slot-definition-name (class-slots class)))
230 (when (slot-boundp original slot)
231 (setf (slot-value copy slot)
232 (slot-value original slot))))
233 copy))
234|#
Note: See TracBrowser for help on using the repository browser.