|
|
1.1 root 1: /*
2: * return 1 if args are within epsilon else 0
3: */
4: double epsilon = 2.3e-16;
5:
6: approx(a, b)
7: double a, b;
8: {
9: double e;
10:
11: if (0.0 > (e = epsilon * b))
12: e = -e;
13:
14: if (0.0 > (a -= b))
15: a = -a;
16:
17: return (a <= e);
18: }
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.