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@ 2513

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

* empty log message *

File size: 4.9 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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
23;;
24;; coefficient ring operations
25;;
26;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
27;;
28;; these are all operations that are performed on the coefficients by
29;; the package, and thus the coefficient ring can be changed by merely
30;; redefining these operations.
31;;
32;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
33
34(defpackage "RING"
35 (:use :cl)
36 (:export "R-PARSE"
37 "UNIT-ELEMENT"
38 "R-ZEROP"
39 "R+"
40 "R-"
41 "R*"
42 "R+"
43 "R/"
44 "R-EXPT"
45 "R-LCM"
46 "R-EZGCD"
47 "R-GCD"
48 "R-TOTAL-DEGREE"
49 "R-DIMENSION"
50 "R-EXPONENTS"
51 "R-COEFF"
52 "R-SUGAR"
53 "R-DIVIDES-P"
54 "R-DIVIDES-LCM-P"
55 "R-LCM-DIVIDES-LCM-P"
56 "R-LCM-EQUAL-LCM-P"
57 "R-REL-PRIME-P"
58 "R-EQUALP"
59 "R-ELT"
60 "R->LIST"
61 "R-DIVISIBLE-BY-P"
62 "R-REL-PRIME-P"
63 "R-DEPENDS-P"
64 "R-TENSOR-PRODUCT"
65 "R-CONTRACT"
66 "R-LENGTH"
67 "MULTIPLY-BY"
68 "ADD-TO"
69 "SUBTRACT-FROM"
70 "UNARY-UMINUS"
71 "SCALAR"
72 "INSERT-ITEM"
73 "APPEND-ITEM"
74))
75
76(in-package :ring)
77
78(defclass scalar ()
79 ((value :initarg :value :accessor value))
80 (:documentation "Wraps objects suitable as scalars/polynomial coefficients"))
81
82(defgeneric unit-element (class))
83
84(defgeneric r-zerop (object)
85 (:method ((self number)) (zerop self)))
86
87(defgeneric r+ (x y)
88 (:method ((x number) (y number)) (+ x y)))
89
90(defgeneric r- (x y)
91 (:method ((x number) (y number)) (- x y)))
92
93(defgeneric r* (x y)
94 (:method ((x number) (y number)) (* x y)))
95
96(defgeneric r-tensor-product (x y))
97
98(defgeneric r/ (x y)
99 (:method ((x number) (y number)) (/ x y)))
100
101(defgeneric r-lcm (x y)
102 (:method ((x integer) (y integer)) (lcm x y)))
103
104(defgeneric r-expt (x y)
105 (:method ((x integer) (y integer)) (expt x y)))
106
107(defgeneric r-ezgcd (x y)
108 (:method ((x integer) (y integer)
109 &aux (c (gcd x y)))
110 (values c (/ x c) (/ y c))))
111
112(defgeneric r-gcd (x y)
113 (:method ((x integer) (y integer))
114 (gcd x y)))
115
116(defgeneric r-dimension (object))
117(defgeneric r-exponents (object))
118(defgeneric r-coeff (object))
119
120(defgeneric r-total-degree (object &optional start end))
121
122(defgeneric r-divides-p (object1 object2)
123 (:method ((object1 integer) (object2 integer))
124 (zerop (rem object2 object1)))
125 (:documentation "Returns T if OBJECT1 divides OBJECT2"))
126
127(defgeneric r-divides-lcm-p (object1 object2 object3)
128 (:documentation "Returns T if OBJECT divides LCM(OBJECT2, OBJECT3), NIL otherwise."))
129
130(defgeneric r-lcm-divides-lcm-p (object1 object2 object3 object4)
131 (:documentation "Returns T if LCM(OBJECT1,OBJECT2) divides LCM(OBJECT3,OBJECT4), NIL otherwise."))
132
133(defgeneric r-lcm-equal-lcm-p (object1 object2 object3 object4)
134 (:documentation "Returns T if object LCM(OBJECT1,OBJECT2) equals LCM(OBJECT3,OBJECT4), NIL otherwise."))
135
136(defgeneric r-equalp (object1 object2)
137 (:method ((object1 t) (object2 t)) (equalp object1 object2))
138 (:documentation "Equality using deep comparison of object slots."))
139
140(defgeneric r-elt (object index))
141
142(defgeneric (setf r-elt) (new-value object index))
143
144(defgeneric r-length (object))
145
146(defgeneric r->list (object))
147(defgeneric r-sugar (object))
148(defgeneric r-rel-prime-p (object1 object2))
149(defgeneric r-contract (object k))
150(defgeneric r-divisible-by-p (object1 object2))
151(defgeneric r-depends-p (object k))
152
153(defgeneric multiply-by (self other)
154 (:method (self other) (r* self other)))
155
156(defgeneric add-to (self other)
157 (:method (self other) (r+ self other)))
158
159(defgeneric subtract-from (self other)
160 (:method (self other) (r- self other)))
161
162(defgeneric unary-uminus (self))
163
164(defgeneric insert-item (self item))
165(defgeneric append-item (self item))
166
Note: See TracBrowser for help on using the repository browser.