61 |
|
|
62 |
|
static int nabases = 3; /* current number of defined bases */ |
63 |
|
|
64 |
< |
#define FEQ(a,b) ((a)-(b) <= 1e-7 && (b)-(a) <= 1e-7) |
64 |
> |
#define FEQ(a,b) ((a)-(b) <= 1e-6 && (b)-(a) <= 1e-6) |
65 |
|
|
66 |
+ |
static int |
67 |
+ |
fequal(double a, double b) |
68 |
+ |
{ |
69 |
+ |
if (b != .0) |
70 |
+ |
a = a/b - 1.; |
71 |
+ |
return((a <= 1e-6) & (a >= -1e-6)); |
72 |
+ |
} |
73 |
+ |
|
74 |
|
// returns the name of the given tag |
75 |
|
#ifdef ezxml_name |
76 |
|
#undef ezxml_name |
235 |
|
if (!i) |
236 |
|
abase_list[nabases].lat[i].tmin = |
237 |
|
-abase_list[nabases].lat[i+1].tmin; |
238 |
< |
else if (!FEQ(atof(ezxml_txt(ezxml_child(ezxml_child(wbb, |
238 |
> |
else if (!fequal(atof(ezxml_txt(ezxml_child(ezxml_child(wbb, |
239 |
|
"ThetaBounds"), "LowerTheta"))), |
240 |
|
abase_list[nabases].lat[i].tmin)) |
241 |
|
error(WARNING, "theta values disagree in custom basis"); |