Initial commit

1 files changed, 34 insertions, 0 deletions

diff --git a/noao/onedspec/ecidentify/eccenter.x b/noao/onedspec/ecidentify/eccenter.x
new file mode 100644
index 00000000..730ad2a8
--- /dev/null
+++ b/noao/onedspec/ecidentify/eccenter.x
@@ -0,0 +1,34 @@
+include	"ecidentify.h"
+
+# EC_CENTER -- Locate the center of a feature.
+
+double procedure ec_center (ec, x, width, type)
+
+pointer	ec				# EC pointer
+double	x				# Initial guess
+real	width				# Feature width
+int	type				# Feature type
+
+double	dvalue
+real	value
+
+real	center1d()
+double	smw_c1trand()
+
+begin
+	if (IS_INDEFD(x))
+	    return (x)
+
+	dvalue = smw_c1trand (EC_PL(ec), x)
+	if (IS_INDEFD(dvalue))
+	    return (dvalue)
+
+	value = dvalue
+ 	value = center1d (value, IMDATA(ec,1), EC_NPTS(ec), width,
+	    abs (type), EC_CRADIUS(ec), EC_THRESHOLD(ec))
+
+	if (IS_INDEF(value))
+	    return (INDEFD)
+	else
+	    return (smw_c1trand (EC_LP(ec), double(value)))
+end