From 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 Mon Sep 17 00:00:00 2001 From: Joe Hunkeler Date: Tue, 11 Aug 2015 16:51:37 -0400 Subject: Repatch (from linux) of OSX IRAF --- sys/gio/fpndgr.x | 21 +++++++++++++++++++++ 1 file changed, 21 insertions(+) create mode 100644 sys/gio/fpndgr.x (limited to 'sys/gio/fpndgr.x') diff --git a/sys/gio/fpndgr.x b/sys/gio/fpndgr.x new file mode 100644 index 00000000..ae471e34 --- /dev/null +++ b/sys/gio/fpndgr.x @@ -0,0 +1,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 -- cgit