On Thu, Oct 13, 2022 at 12:29 PM Henry Rich <henryhr...@gmail.com> wrote:
> I don't understand the question: what are you trying to compute?  Doesn't
> gmp support division?

Sure, and mpq_get_d() documented at
https://gmplib.org/manual/Rational-Conversions fails J's accuracy
tests, with off-by-least-significant-bit errors.

I've got most of those worked out, but I still fail on this test.

> Knuth has code for division IIRC.

It's not just division, it's the stability of the algorithm when
discarding precision.

Probably I just need to find a relevant reference on numerical analysis.

Thanks,

-- 
Raul
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm

Reply via email to