Toggle navigation
Toggle navigation
This project
Loading...
Sign in
oct
/
oct
Go to a project
Toggle navigation
Toggle navigation pinning
Projects
Groups
Snippets
Help
Project
Activity
Repository
Pipelines
Graphs
Issues
0
Merge Requests
0
Wiki
Snippets
Network
Create a new issue
Builds
Commits
Issue Boards
Authored by
Raymond Toy
20120323 13:00:24 0700
Browse Files
Options
Browse Files
Download
Email Patches
Plain Diff
Commit
d80256d66d585a69813f277b115d039cc620e62d
d80256d6
1 parent
b1ae6953
Need to apply contagion for expintegrale.
Hide whitespace changes
Inline
Sidebyside
Showing
1 changed file
with
23 additions
and
20 deletions
qdgamma.lisp
qdgamma.lisp
View file @
d80256d
...
...
@@ 528,27 +528,30 @@
;; for arg(z) < pi.
;;
;;
(
cond
((
and
(
realp
v
)
(
minusp
v
))
;; E(v, z) = z^(v1)*incomplete_gamma_tail(v+1,z)
(
let
((
v
(

v
)))
(
let*
((
prec
(
floatcontagion
v
z
))
(
v
(
applycontagion
v
prec
))
(
z
(
applycontagion
z
prec
)))
(
cond
((
and
(
realp
v
)
(
minusp
v
))
;; E(v, z) = z^(v1)*incomplete_gamma_tail(v+1,z)
(
let
((
v
(

v
)))
(
*
(
expt
z
(

v
1
))
(
incompletegammatail
(
+
v
1
)
z
))))
((
<
(
abs
z
)
1
)
;; Use series for small z
(
sexpintegrale
v
z
))
((
>=
(
abs
(
phase
z
))
3.1
)
;; The continued fraction doesn't converge on the negative
;; real axis, and converges very slowly near the negative
;; real axis, so use the incompletegammatail function in
;; this region. "Closeness" to the negative real axis is
;; teken to mean that z is in a sector near the axis.
;;
;; E(v,z) = z^(v1)*incomplete_gamma_tail(1v,z)
(
*
(
expt
z
(

v
1
))
(
incompletegammatail
(
+
v
1
)
z
))))
((
<
(
abs
z
)
1
)
;; Use series for small z
(
sexpintegrale
v
z
))
((
>=
(
abs
(
phase
z
))
3.1
)
;; The continued fraction doesn't converge on the negative
;; real axis, and converges very slowly near the negative
;; real axis, so use the incompletegammatail function in
;; this region. "Closeness" to the negative real axis is
;; teken to mean that z is in a sector near the axis.
;;
;; E(v,z) = z^(v1)*incomplete_gamma_tail(1v,z)
(
*
(
expt
z
(

v
1
))
(
incompletegammatail
(

1
v
)
z
)))
(
t
;; Use continued fraction for everything else.
(
cfexpintegrale
v
z
))))
(
incompletegammatail
(

1
v
)
z
)))
(
t
;; Use continued fraction for everything else.
(
cfexpintegrale
v
z
)))))
;; Series for Fresnel S
;;
...
...
Please
register
or
login
to post a comment