Hello everyone, I'm using stats.zprob and need some more decimal points than standard python floats can give me, and I need the function to be fast. Can I get the standard scipy zprob or ndtr to return float96 in some way? Or do I need to modify something in special/cephes or something like that? cheers Martin
Hello everyone, I'm using stats.zprob and need some more decimal points than standard python floats can give me, and I need the function to be fast. Can I get the standard scipy zprob or ndtr to return float96 in some way? Or do I need to modify something in special/cephes or something like that? cheers Martin
Martin Blom wrote:
Hello everyone,
I'm using stats.zprob and need some more decimal points than standard python floats can give me, and I need the function to be fast. Can I get the standard scipy zprob or ndtr to return float96 in some way? Or do I need to modify something in special/cephes or something like that?
You would have to find or write a routine that used higher precision floats. We don't wrap any of them. Cephes does have quad-precision versions of many of their functions, but you would have to wrap them yourself. -- Robert Kern "I have come to believe that the whole world is an enigma, a harmless enigma that is made terrible by our own mad attempt to interpret it as though it had an underlying truth." -- Umberto Eco
Rrrrgh. That's what I feared! Well, thank you! Martin 2007/9/13, Robert Kern <robert.kern@gmail.com>:
Martin Blom wrote:
Hello everyone,
I'm using stats.zprob and need some more decimal points than standard python floats can give me, and I need the function to be fast. Can I get the standard scipy zprob or ndtr to return float96 in some way? Or do I need to modify something in special/cephes or something like that?
You would have to find or write a routine that used higher precision floats. We don't wrap any of them. Cephes does have quad-precision versions of many of their functions, but you would have to wrap them yourself.
-- Robert Kern
"I have come to believe that the whole world is an enigma, a harmless enigma that is made terrible by our own mad attempt to interpret it as though it had an underlying truth." -- Umberto Eco _______________________________________________ SciPy-user mailing list SciPy-user@scipy.org http://projects.scipy.org/mailman/listinfo/scipy-user
participants (2)
-
Martin Blom -
Robert Kern