|
|
coherent
/*
* return 1 if args are within epsilon else 0
*/
double epsilon = 2.3e-16;
approx(a, b)
double a, b;
{
double e;
if (0.0 > (e = epsilon * b))
e = -e;
if (0.0 > (a -= b))
a = -a;
return (a <= e);
}
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.