Changeset 535d2509c2974c113e6c1bfa27ece5d5b8f99bd1

Timestamp:

03/06/11 18:30:15 (2 years ago)

Author:

Raymond Toy <toy.raymond@…>

Children:

c014c5f4e4b941aacf16a897a2cf689ed49b37ef

Parents:

258cd829289826458d6075caa279b253b8848749

git-committer:

Raymond Toy <toy.raymond@…> (03/06/11 18:30:15)

Message:

Add support for Jacobi elliptic functions

qd-elliptic.lisp:
o Implementation for Jacobi sn, cn, and dn functions.

qd-methods.lisp:
o Add EPSILON method to return the floating-point epsilon value for

the given float.

oct.asd:
o Build qd-elliptic.lisp.

Files:

• oct.asd (modified) (1 diff)
• qd-elliptic.lisp (added)
• qd-methods.lisp (modified) (1 diff)

oct.asd

@@ -59 +59 @@
    (:file "qd-complex"
           :depends-on ("qd-methods"))
+   (:file "qd-elliptic"
+          :depends-on ("qd-methods"))
    ))
 

qd-methods.lisp

@@ -1090 +1090 @@
 ;; %qd-real.
 (set-dispatch-macro-character #\# #\Q #'qd-class-reader)
+
+
+(defmethod epsilon ((m cl:float))
+  (etypecase m
+    (single-float single-float-epsilon)
+    (double-float double-float-epsilon)))
+
+(defmethod epsilon ((m cl:complex))
+  (epsilon (realpart m)))
+
+(defmethod epsilon ((m qd-real))
+  ;; What is the epsilon value for a quad-double?  This is complicated
+  ;; by the fact that things like (+ #q1 #q1q-100) is representable as
+  ;; a quad-double.  For most purposes we want epsilon to be close to
+  ;; the 212 bits of precision (4*53 bits) that we normally have with
+  ;; a quad-double.
+  (scale-float #q1 -212))
+
+(defmethod epsilon ((m qd-complex))
+  (epsilon (realpart m)))