diff options
author | Joe Hunkeler <jhunkeler@gmail.com> | 2015-08-11 16:51:37 -0400 |
---|---|---|
committer | Joe Hunkeler <jhunkeler@gmail.com> | 2015-08-11 16:51:37 -0400 |
commit | 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 (patch) | |
tree | 4464880c571602d54f6ae114729bf62a89518057 /noao/onedspec/ecidentify/eccenter.x | |
download | iraf-osx-40e5a5811c6ffce9b0974e93cdd927cbcf60c157.tar.gz |
Repatch (from linux) of OSX IRAF
Diffstat (limited to 'noao/onedspec/ecidentify/eccenter.x')
-rw-r--r-- | noao/onedspec/ecidentify/eccenter.x | 34 |
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 |