Changeset 16060 for trunk/Ohana/src/relastro

Timestamp:

Jan 14, 2008, 1:50:46 PM (19 years ago)

Author:

eugene

Message:

convert to new gauss-jordan code (non-Press); fixing Wall,Werror warnings

Location:

trunk/Ohana/src/relastro

Files:

• include/relastro.h (modified) (1 diff)
• src/ConfigInit.c (modified) (1 diff)
• src/FitPM.c (modified) (1 diff)
• src/FitPMandPar.c (modified) (1 diff)
• src/FitPar.c (modified) (2 diffs)
• src/ImageOps.c (modified) (1 diff)
• src/MosaicOps.c (modified) (1 diff)
• src/UpdateObjects.c (modified) (1 diff)
• src/fitpoly.c (modified) (1 diff)
• src/mkpolyterm.c (modified) (1 diff)
• src/plotstuff.c (modified) (1 diff)
• src/relastro.c (modified) (1 diff)

trunk/Ohana/src/relastro/include/relastro.h

@@ -207 +207 @@
 void unlock_image_db (FITS_DB *db);
 void create_image_db (FITS_DB *db);
+void save_catalogs (Catalog *catalog, int Ncatalog);
 
 int           main                PROTO((int argc, char **argv));

trunk/Ohana/src/relastro/src/ConfigInit.c

@@ -3 +3 @@
 void ConfigInit (int *argc, char **argv) {
 
-  double ZERO_POINT;
   char  *config, *file;
   char CatdirPhotcodeFile[256];

trunk/Ohana/src/relastro/src/FitPM.c

@@ -53 +53 @@
   B[3][0] = YT;
 
-  dgaussj (A, 4, B, 1);
+  dgaussjordan (A, B, 4, 1);
 
   fit[0].Ro = B[0][0];

trunk/Ohana/src/relastro/src/FitPMandPar.c

@@ -76 +76 @@
   B[4][0] = PRX + PDY;
 
-  dgaussj (A, 5, B, 1);
+  dgaussjordan (A, B, 5, 1);
 
   fit[0].Ro = B[0][0];

trunk/Ohana/src/relastro/src/FitPar.c

@@ -7 +7 @@
 
   double **A, **B;
-  double wx, wy, Wx, Wy, Tx, Ty, Tx2, Ty2, Xs, Ys, XT, YT;
-  double PR, PD, PRT, PDT, PRX, PDY, PR2, PD2;
+  double wx, wy, Wx, Wy, Xs, Ys;
+  double PR, PD, PRX, PDY, PR2, PD2;
 
   A = array_init (3, 3);
@@ -50 +50 @@
   B[2][0] = PRX + PDY;
 
-  dgaussj (A, 3, B, 1);
+  dgaussjordan (A, B, 3, 1);
 
   fit[0].Ro = B[0][0];

trunk/Ohana/src/relastro/src/ImageOps.c

@@ -289 +289 @@
       LM_to_XY (&ref[i].X, &ref[i].Y, ref[i].L, ref[i].M, &image[im].coords);
       break;
+      default:
+        fprintf (stderr, "invalid case");
+        abort();
     }
   }

trunk/Ohana/src/relastro/src/MosaicOps.c

@@ -22 +22 @@
 void initMosaics (Image *image, int Nimage) {
 
-  int i, j, status, found, NMOSAIC, *NIMLIST;
+  int i, j, found, NMOSAIC;
   unsigned int start, stop;
-  char *pname;
 
   Nmosaic = 0;

trunk/Ohana/src/relastro/src/UpdateObjects.c

@@ -179 +179 @@
 
         default:
-          fprintf (stderr, "programming error at %s, %s", __FILE__, __LINE__);
+          fprintf (stderr, "programming error at %s, %d", __FILE__, __LINE__);
           exit (2);
       }   

trunk/Ohana/src/relastro/src/fitpoly.c

@@ -137 +137 @@
   }     
 
-  dgaussj (matrix, fit[0].Nelems, vector, 2); 
+  dgaussjordan (matrix, vector, fit[0].Nelems, 2); 
 
   for (i = 0; i < fit[0].Nelems; i++) {

trunk/Ohana/src/relastro/src/mkpolyterm.c

@@ -61 +61 @@
       beta[1][0] = poly2d_eval (yfit, Nx, Ny, Xo, Yo);
 
-      dgaussj (alpha, 2, beta, 1);
+      dgaussjordan (alpha, beta, 2, 1);
 
       Xo -= beta[0][0];

trunk/Ohana/src/relastro/src/plotstuff.c

@@ -79 +79 @@
 
   float *values;
-  int i, Nbytes;
+  int i;
 
   if (Npts < 1) return;

trunk/Ohana/src/relastro/src/relastro.c

@@ -55 +55 @@
 
     default:
-      fprintf (stderr, "programming error at %s:%s", __FILE__, __LINE__);
+      fprintf (stderr, "programming error at %s:%d", __FILE__, __LINE__);
       exit (2);
   }