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

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

* empty log message *

File size: 4.6 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 "R-UNIT-FOR"
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))
72
73(in-package :ring)
74
75(defgeneric r-unit-for (object)
76 (:method ((self number)) 1))
77
78(defgeneric r-zerop (object)
79 (:method ((self number)) (zerop self)))
80
81(defgeneric r+ (x y)
82 (:method ((x number) (y number)) (+ x y)))
83
84(defgeneric r- (x y)
85 (:method ((x number) (y number)) (- x y)))
86
87(defgeneric r* (x y)
88 (:method ((x number) (y number)) (* x y)))
89
90(defgeneric r-tensor-product (x y))
91
92(defgeneric r/ (x y)
93 (:method ((x number) (y number)) (/ x y)))
94
95(defgeneric r-lcm (x y)
96 (:method ((x integer) (y integer)) (lcm x y)))
97
98(defgeneric r-expt (x y)
99 (:method ((x integer) (y integer)) (expt x y)))
100
101(defgeneric r-ezgcd (x y)
102 (:method ((x integer) (y integer)
103 &aux (c (gcd x y)))
104 (values c (/ x c) (/ y c))))
105
106(defgeneric r-gcd (x y)
107 (:method ((x integer) (y integer))
108 (gcd x y)))
109
110(defgeneric r-dimension (object))
111(defgeneric r-exponents (object))
112(defgeneric r-coeff (object))
113
114(defgeneric r-total-degree (object &optional start end))
115
116(defgeneric r-divides-p (object1 object2)
117 (:method ((object1 integer) (object2 integer))
118 (zerop (rem object2 object1)))
119 (:documentation "Returns T if OBJECT1 divides OBJECT2"))
120
121(defgeneric r-divides-lcm-p (object1 object2 object3)
122 (:documentation "Returns T if OBJECT divides LCM(OBJECT2, OBJECT3), NIL otherwise."))
123
124(defgeneric r-lcm-divides-lcm-p (object1 object2 object3 object4)
125 (:documentation "Returns T if LCM(OBJECT1,OBJECT2) divides LCM(OBJECT3,OBJECT4), NIL otherwise."))
126
127(defgeneric r-lcm-equal-lcm-p (object1 object2 object3 object4)
128 (:documentation "Returns T if object LCM(OBJECT1,OBJECT2) equals LCM(OBJECT3,OBJECT4), NIL otherwise."))
129
130(defgeneric r-equalp (object1 object2)
131 (:method ((object1 t) (object2 t)) (equalp object1 object2))
132 (:documentation "Equality using deep comparison of object slots."))
133
134(defgeneric r-elt (object index))
135
136(defgeneric (setf r-elt) (new-value object index))
137
138(defgeneric r-length (object))
139
140(defgeneric r->list (object))
141(defgeneric r-sugar (object))
142(defgeneric r-rel-prime-p (object1 object2))
143(defgeneric r-contract (object k))
144(defgeneric r-divisible-by-p (object1 object2))
145(defgeneric r-depends-p (object k))
146
147(defgeneric multiply-by (self other)
148 (:method (self other) (r* self other)))
149
150(defgeneric add-to (self other)
151 (:method (self other) (r+ self other)))
152
153(defgeneric subtract-from (self other)
154 (:method (self other) (r- self other)))
155
156(defmethod unary-uminus (self))
157
Note: See TracBrowser for help on using the repository browser.