Changeset 20839 for branches/eam_branch_20081124/Ohana/src/opihi/cmd.astro/imfit-qrgauss.c

Timestamp:

Nov 25, 2008, 4:51:39 PM (18 years ago)

Author:

eugene

Message:

convert opihi vectors to double type, add option for int type as well

File:

• branches/eam_branch_20081124/Ohana/src/opihi/cmd.astro/imfit-qrgauss.c (modified) (4 diffs)

branches/eam_branch_20081124/Ohana/src/opihi/cmd.astro/imfit-qrgauss.c

@@ -1 +1 @@
 # include "imfit.h"
 
-float qrgaussTD (float, float, float *, int, float *);
+opihi_flt qrgaussTD (opihi_flt, opihi_flt, opihi_flt *, int, opihi_flt *);
 void  qrgaussCL ();
 
@@ -14 +14 @@
 
   /* allocate free and fixed parameters */
-  ALLOCATE (par, float, MAX (Npar, 1));
-  bzero (par, Npar*sizeof(float));
-  ALLOCATE (fpar, float, MAX (Nfpar, 1));
-  bzero (fpar, Nfpar*sizeof(float));
+  ALLOCATE (par, opihi_flt, MAX (Npar, 1));
+  bzero (par, Npar*sizeof(opihi_flt));
+  ALLOCATE (fpar, opihi_flt, MAX (Nfpar, 1));
+  bzero (fpar, Nfpar*sizeof(opihi_flt));
 
   par[0]  = get_variable_default ("Xg", 0);
@@ -44 +44 @@
 }
 
-float qrgaussTD (float x, float y, float *par, int Npar, float *dpar) {
+opihi_flt qrgaussTD (opihi_flt x, opihi_flt y, opihi_flt *par, int Npar, opihi_flt *dpar) {
 
-  float X, Y, px, py;
-  float z, r, q, f;
+  opihi_flt X, Y, px, py;
+  opihi_flt z, r, q, f;
 
   X = x - par[0];
@@ -77 +77 @@
 # if (0)
 /* one component, two slopes: (1 + z^M + z^N)^(-1) -- x, y, sx, sy, sxy, I, sky, sr */
-float qrgaussTD (float x, float y, float *par, int Npar, float *dpar) {
+opihi_flt qrgaussTD (opihi_flt x, opihi_flt y, opihi_flt *par, int Npar, opihi_flt *dpar) {
 
-  float X, Y, px, py;
-  float z, r, q, f;
+  opihi_flt X, Y, px, py;
+  opihi_flt z, r, q, f;
 
   /* if (par[7] < 1.0) par[7] = 1.0; */