Changeset 235ac2
 Timestamp:
 04/12/12 21:37:27 (3 years ago)
 Branches:
 master
 Children:
 5daefc
 Parents:
 0823ac
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

qdbessel.lisp
r9fd2eb r235ac2 184 184 (let* ((iz (* #c(0 1) z)) 185 185 (i+ (exparci2 iz v))) 186 (cond (( = v (ftruncate v))186 (cond ((and (= v (ftruncate v)) (realp z)) 187 187 ;; We can simplify the result 188 (let ((c ( cis (* v (floatpi i+) 1/2))))188 (let ((c (exp (* v (floatpi i+) #c(0 1/2))))) 189 189 (/ (+ (* c i+) 190 190 (* (conjugate c) (conjugate i+))) … … 193 193 (t 194 194 (let ((i (exparci2 ( iz ) v))) 195 (/ (+ (* ( cis (* v (floatpi i+) 1/2))195 (/ (+ (* (exp (* v (floatpi i+) #c(0 1/2))) 196 196 i+) 197 (* ( cis (* v (floatpi i+) 1/2))197 (* (exp (* v (floatpi i+) #c(0 1/2))) 198 198 i)) 199 199 (floatpi i+) … … 395 395 (anclrhash) 396 396 (%bigaclrhash) 397 (cond (( = vv v)398 ;; v is an integer 397 (cond ((and (= vv v) (realp z)) 398 ;; v is an integer and z is real 399 399 (integerbesseljexparc v z)) 400 400 (t … … 403 403 (vpi (* v (floatpi (realpart z))))) 404 404 (+ (integerbesseljexparc v z) 405 (* z 406 (/ (sin vpi) vpi) 407 (+ (/ 1 z) 408 (sumab bign v z) 409 (sumbigia bign v z))))))))) 405 (if (= vv v) 406 0 407 (* z 408 (/ (sin vpi) vpi) 409 (+ (/ 1 z) 410 (sumab bign v z) 411 (sumbigia bign v z)))))))))) 410 412 411 413 ;; Bessel Y
Note: See TracChangeset
for help on using the changeset viewer.