|
|
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.