Changeset 235ac2e4756b7c3f03815b267d292a26e8e1df37

Timestamp:

04/12/12 14:37:27 (14 months ago)

Author:

Raymond Toy <toy.raymond@…>

Children:

5daefcc43f8b88741634f7a982bdb8b976739a1c

Parents:

0823acf49a91cbdb9a3e7f1faae2efe128580038

git-committer:

Raymond Toy <toy.raymond@…> (04/12/12 14:37:27)

Message:

Fix for complex args.

Files:

• qd-bessel.lisp (modified) (4 diffs)

qd-bessel.lisp

@@ -184 +184 @@
   (let* ((iz (* #c(0 1) z))
          (i+ (exp-arc-i-2 iz v)))
-    (cond ((= v (ftruncate v))
+    (cond ((and (= v (ftruncate v)) (realp z))
            ;; We can simplify the result
-           (let ((c (cis (* v (float-pi i+) -1/2))))
+           (let ((c (exp (* v (float-pi i+) #c(0 -1/2)))))
              (/ (+ (* c i+)
                    (* (conjugate c) (conjugate i+)))
@@ -193 +193 @@
           (t
            (let ((i- (exp-arc-i-2 (- iz ) v)))
-             (/ (+ (* (cis (* v (float-pi i+) -1/2))
+             (/ (+ (* (exp (* v (float-pi i+) #c(0 -1/2)))
                       i+)
-                   (* (cis (* v (float-pi i+) 1/2))
+                   (* (exp (* v (float-pi i+) #c(0 1/2)))
                       i-))
                 (float-pi i+)
@@ -395 +395 @@
     (an-clrhash)
     (%big-a-clrhash)
-    (cond ((= vv v)
-           ;; v is an integer
+    (cond ((and (= vv v) (realp z))
+           ;; v is an integer and z is real
            (integer-bessel-j-exp-arc v z))
           (t
@@ -403 +403 @@
                  (vpi (* v (float-pi (realpart z)))))
              (+ (integer-bessel-j-exp-arc v z)
-                (* z
-                   (/ (sin vpi) vpi)
-                   (+ (/ -1 z)
-                      (sum-ab big-n v z)
-                      (sum-big-ia big-n v z)))))))))
+                (if (= vv v)
+                    0
+                    (* z
+                       (/ (sin vpi) vpi)
+                       (+ (/ -1 z)
+                          (sum-ab big-n v z)
+                          (sum-big-ia big-n v z))))))))))
 
 ;; Bessel Y