Changeset 1d9ec0 for qd-gamma.lisp


Ignore:
Timestamp:
03/22/12 05:24:48 (3 years ago)
Author:
Raymond Toy <toy.raymond@…>
Branches:
master
Children:
e10e40
Parents:
4b332e
Message:

Handle the case of negative v in exp-integral-e(v, z). Add tests.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • qd-gamma.lisp

    r4b332e r1d9ec0  
    528528  ;;
    529529  ;;
    530   (cond ((< (abs z) 1)
     530  (cond ((and (realp v) (minusp v))
     531         ;; E(-v, z) = z^(-v-1)*incomplete_gamma_tail(v+1,z)
     532         (let ((-v (- v)))
     533           (* (expt z (- v 1))
     534              (incomplete-gamma-tail (+ -v 1) z))))
     535        ((< (abs z) 1)
    531536         ;; Use series for small z
    532537         (s-exp-integral-e v z))
Note: See TracChangeset for help on using the changeset viewer.