User talk:Hailholyghost: Difference between revisions
m
no edit summary
mNo edit summary |
mNo edit summary |
||
Line 13:
::You write the C code is trivial: really? How? You will need at least the incomplete beta function from [https://www.gnu.org/software/gsl/manual/html_node/Incomplete-Beta-Function.html GSL] (or IMSL, NAG or any library with it, even a Fortran-based library), and this will require more than one line (you may do more or less what is done here in Fortran or Python). Even getting a working GSL on Windows is not trivial at all. And what you call the non-trivial thing is really not about integration at all.
::[[User:Eoraptor|Eoraptor]] ([[User talk:Eoraptor|talk]]) 13:13, 8 December 2017 (UTC)
:::For instance, the R implementation uses a port of [https://dl.acm.org/citation.cfm?id=131776 TOMS 708] from the ACM. The original code is [http://calgo.acm.org/708.gz here]. Notice that, as usual with special functions, there are calls to different methods according the the values of the arguments. There are reasons not to reimplement this from scratch: it requires much research work to prove a given algorithm is correct, to find necessary coefficients with enough precision, and to get a correct and efficient implementation. Here you are not even trying to investigate the convergence of the integral, and it's an inefficient way to find an accurate answer. [[User:Eoraptor|Eoraptor]] ([[User talk:Eoraptor|talk]]) 13:30, 8 December 2017 (UTC)
| |||