Floating point arithmetic.
From: Amit Bhatia (bhatia2_at_nospam.com)
Date: Sat, 10 Jul 2004 23:08:01 -0500
I am cross posting this on comp.lang.c as well: sorry for same.
The problem I am facing is as follows:
double a= 0.15;
const double VERYTINY =1.e-10;
I know b/a = 16 and hence the remainder is zero; but I am not
able to find any suitable thing to encode it into in c.
for example (fmod(b,a)>VERYTINY) returns true!
Now for this particular instance, (fmodf(b,a)>VERYTINY)
does return false.
But now if
then fmodf and fmod both don't help...
any suggestions on this?
I was pointed to a reference on floating point arithmetic, where they talk
of ulps etc, but is there a small function or fix to deal with this problem