aboutsummaryrefslogtreecommitdiff
path: root/noao/onedspec/ecidentify/ecnearest.x
diff options
context:
space:
mode:
authorJoseph Hunkeler <jhunkeler@gmail.com>2015-07-08 20:46:52 -0400
committerJoseph Hunkeler <jhunkeler@gmail.com>2015-07-08 20:46:52 -0400
commitfa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 (patch)
treebdda434976bc09c864f2e4fa6f16ba1952b1e555 /noao/onedspec/ecidentify/ecnearest.x
downloadiraf-linux-fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4.tar.gz
Initial commit
Diffstat (limited to 'noao/onedspec/ecidentify/ecnearest.x')
-rw-r--r--noao/onedspec/ecidentify/ecnearest.x26
1 files changed, 26 insertions, 0 deletions
diff --git a/noao/onedspec/ecidentify/ecnearest.x b/noao/onedspec/ecidentify/ecnearest.x
new file mode 100644
index 00000000..7b061472
--- /dev/null
+++ b/noao/onedspec/ecidentify/ecnearest.x
@@ -0,0 +1,26 @@
+include <mach.h>
+include "ecidentify.h"
+
+# EC_NEAREST -- Find the nearest feature to a given coordinate.
+
+procedure ec_nearest (ec, fitnear)
+
+pointer ec # ID pointer
+double fitnear # Coordinate to find nearest feature
+
+int i, ec_next()
+double delta, delta1
+
+begin
+ EC_CURRENT(ec) = 0
+
+ i = 0
+ delta = MAX_REAL
+ while (ec_next (ec, i) != EOF) {
+ delta1 = abs (FIT(ec,i) - fitnear)
+ if (delta1 < delta) {
+ EC_CURRENT(ec) = i
+ delta = delta1
+ }
+ }
+end