On 8/23/2019 8:39 PM, Shaobo He wrote:
Hello everyone,
I'm in need of quad-precision floating-point type/operations in Racket
(the bigfloat module doesn't work for me). It appears there's none, am
I right?
If so, I think I will try to write one myself, which is going to be a
wrapper to libquadmath. Any suggestions about how I should approach
it? The first issue seems to be marshalling the C `__float128` type
into a Racket type.
Thanks,
Shaobo
I'm not sure exactly what you have in mind, and I am a total ignoramus
about Racket numerics, but it is possible to implement quad precision
floating point arithmetic using only double precision floating point,
while maintaining rigorous error bounds.
(... and so on, doubling precision precision at each step.)
These sources might be useful.
Design and Implementation of a High Precision Arithmetic with Rigorous
Error Bounds
Alexander WittigMSUHEP-081126December 2008
https://pdfs.semanticscholar.org/8450/3383403eca34dd7e92a12261a3ed7ea2340a.pdf
The Great Internet Mersenne Prime Search (GIMPS) project carries this to
extremes.
They use Fourier multiplication on integers with millions of bits - and
must have exact results.
https://www.mersenne.org/
Good luck with your project. I am interested in seeing the result!
--
Josh Rubin
jlru...@gmail.com
--
You received this message because you are subscribed to the Google Groups "Racket
Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to racket-users+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/racket-users/af8336d3-4845-1e49-4e9a-1e09ed97dd3a%40gmail.com.
