blob: ae471e34842295b52caec23a40701a982dfbf413 (
plain) (
blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
# FP_NONDEGENR -- If two floating point numbers are equivalent to within the
# machine epsilon, adjust their values until a nondegenerate range is obtained.
# The boolean function returns true if it has to MAKE the range nondegenerate,
# i.e., if it modifies their values.
bool procedure fp_nondegenr (x1, x2)
real x1, x2 # range to be adjusted
int n
bool fp_equalr()
begin
for (n=0; fp_equalr(x1,x2); n=n+1) {
x1 = x1 - max (abs(x1) * 0.01, 0.01)
x2 = x2 + max (abs(x2) * 0.01, 0.01)
}
return (n > 0)
end
|
