1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
|
[section:issues Known Issues, and Todo List]
This section lists those issues that are known about.
Predominantly this is a TODO list, or a list of possible
future enhancements. Items labled "High Priority" effect
the proper functioning of the component, and should be fixed
as soon as possible. Items labled "Medium Priority" are
desirable enhancements, often pertaining to the performance
of the component, but do not effect it's accuracy or functionality.
Items labled "Low Priority" should probably be investigated at
some point. Such classifications are obviously highly subjective.
If you don't see a component listed here, then we don't have any known
issues with it.
[h4 tgamma]
* Can the __lanczos be optimized any further? (low priority)
[h4 Incomplete Beta]
* Investigate Didonato and Morris' asymptotic expansion for large a and b
(medium priority).
[h4 Inverse Gamma]
* Investigate whether we can skip iteration altogether if the first approximation
is good enough (Medium Priority).
[h4 Polynomials]
* The Legendre and Laguerre Polynomials have surprisingly different error
rates on different platforms, considering they are evaluated with only
basic arithmetic operations. Maybe this is telling us something, or maybe not
(Low Priority).
[h4 Elliptic Integrals]
* Carlson's algorithms are essentially unchanged from Xiaogang Zhang's
Google Summer of Code student project, and are based on Carlson's
original papers. However, Carlson has revised his algorithms since then
(refer to the references in the elliptic integral docs for a list), to
improve performance and accuracy, we may be able to take advantage
of these improvements too (Low Priority).
* [para Carlson's algorithms (mainly R[sub J]) are somewhat prone to
internal overflow/underflow when the arguments are very large or small.
The homogeneity relations:]
[para R[sub F](ka, kb, kc) = k[super -1/2] R[sub F](a, b, c)]
[para and]
[para R[sub J](ka, kb, kc, kr) = k[super -3/2] R[sub J](a, b, c, r)]
[para could be used to sidestep trouble here: provided the problem domains
can be accurately identified. (Medium Priority).]
* Carlson's R[sub C] can be reduced to elementary funtions (asin and log),
would it be more efficient evaluated this way, rather than by Carlson's
algorithms? (Low Priority).
* Should we add an implementation of Carlson's R[sub G]? It's not
required for the Legendre form integrals, but some people may find it
useful (Low Priority).
* There are a several other integrals: D([phi], k), Z([beta], k),
[Lambda][sub 0]([beta], k) and Bulirsch's ['el] functions that could
be implemented using Carlson's integrals (Low Priority).
* The integrals K(k) and E(k) could be implemented using rational
approximations (both for efficiency and accuracy),
assuming we can find them. (Medium Priority).
* There is a sub-domain of __ellint_3 that is unimplemented (see the docs
for details), currently
it's not clear how to solve this issue, or if it's ever likely
to be an real problem in practice - especially as most other implementations
don't support this domain either (Medium Priority).
[h4 Inverse Hyperbolic Functions]
* These functions are inherited from previous Boost versions,
before __log1p became widely available. Would they be better expressed
in terms of this function? This is probably only an issue
for very high precision types (Low Priority).
[h4 Statistical distributions]
* Student's t Perhaps switch to normal distribution as a better approximation for very large degrees of freedom?
[endsect][/section:issues Known Issues, and Todo List]
[/
Copyright 2006 John Maddock and Paul A. Bristow.
Distributed under the Boost Software License, Version 1.0.
(See accompanying file LICENSE_1_0.txt or copy at
http://www.boost.org/LICENSE_1_0.txt).
]
|
