Changeset 7903 for trunk/ppStats

Timestamp:

Jul 13, 2006, 7:20:26 PM (20 years ago)

Author:

Paul Price

Message:

Adding .cvsignore files. Code now compiles, but not yet tested.

Location:

trunk/ppStats

Files:

• . (modified) (1 prop)
• .cvsignore (added)
• src (modified) (1 prop)
• src/.cvsignore (added)
• src/ppStats.h (modified) (1 diff)
• src/ppStatsLoop.c (modified) (8 diffs)
• src/ppStatsSetup.c (modified) (2 diffs)

trunk/ppStats/src/ppStats.h

@@ -2 +2 @@
 #define PP_FRINGE_H
 
-#define RECIPE "PPSTATS"
+#define RECIPENAME "PPSTATS"
 #define TIMERNAME "PPSTATS"
 

trunk/ppStats/src/ppStatsLoop.c

@@ -14 +14 @@
     psString keyword;                    // Header from iteration
     while ((keyword = psListGetAndIncrement(keywordsIter))) {
-        psMetadataAddItem(target, psMetadataLookup(source, keyword));
+        psMetadataAddItem(target, psMetadataLookup(source, keyword), PS_LIST_TAIL, 0);
     }
     return;
@@ -31 +31 @@
     // Iterate through the FPA
     pmFPA *fpa = data->inFPA;           // The FPA of interest
-    if (psListLength(data->headers) > 0 && fpa->hdu && pmFPAHeader(fpa, data->inFile)) {
+    if (psListLength(data->headers) > 0 && fpa->hdu && pmFPAReadHeader(fpa, data->inFile)) {
         pmHDU *hdu = fpa->hdu;          // HDU for headers
         getHeaders(fpaResults, hdu->header, headersIter);
@@ -45 +45 @@
         psMetadata *chipResults = psMetadataAlloc(); // Metadata to hold the chip-level results
 
-        if (psListLength(data->headers) > 0 && chip->hdu && pmChipHeader(chip, data->inFile)) {
+        if (psListLength(data->headers) > 0 && chip->hdu && pmChipReadHeader(chip, data->inFile)) {
             pmHDU *hdu = chip->hdu;     // HDU for headers
             getHeaders(chipResults, hdu->header, headersIter);
@@ -59 +59 @@
             psMetadata *cellResults = psMetadataAlloc(); // Metadata to hold the cell-level results
 
-            if (psListLength(data->headers) > 0 && cell->hdu && pmCellHeader(cell, data->inFile)) {
+            if (psListLength(data->headers) > 0 && cell->hdu && pmCellReadHeader(cell, data->inFile)) {
                 pmHDU *hdu = cell->hdu;     // HDU for headers
                 getHeaders(cellResults, hdu->header, headersIter);
@@ -106 +106 @@
                 psImage *image = readout->image; // The image of interest
                 psImage *mask = readout->mask; // The mask image
-                int numSamples = sample * image->numCols * image->numRows; // Number of samples
-                int sampleSpace = 1.0 / sample; // Space between samples
+                int numSamples = data->sample * image->numCols * image->numRows; // Number of samples
+                int sampleSpace = 1.0 / data->sample; // Space between samples
                 psVector *sampleValues = psVectorAlloc(numSamples, PS_TYPE_F32); // Vector of samples
+                sampleValues->n = numSamples;
                 psVector *sampleMask = NULL;  // Corresponding mask
                 if (mask) {
                     sampleMask = psVectorAlloc(numSamples, PS_TYPE_U8);
+                    sampleMask->n = numSamples;
                 }
-                samples->n = numSamples;
                 for (int i = 0; i < numSamples; i++) {
                     int j = i * sampleSpace;
@@ -120 +121 @@
                     sampleValues->data.F32[i] = image->data.F32[y][x];
                     if (mask) {
-                        sampleMask->data.U8 = mask->data.U8[y][x];
+                        sampleMask->data.U8[i] = mask->data.U8[y][x];
                     }
                 }
@@ -136 +137 @@
 
             #define WRITE_STAT(SYMBOL, NAME, SOURCE) \
-            if (data->stats & SYMBOL) { \
+            if (data->stats->options & SYMBOL) { \
                 psMetadataAddF32(cellResults, PS_LIST_TAIL, NAME, 0, NULL, data->stats->SOURCE); \
             }
@@ -156 +157 @@
 
            // If we didn't add it before, add it now
-            if (!psMetadataLookup(&mdok, chipResults, cellName)) {
+            if (!psMetadataLookup(chipResults, cellName)) {
                 psMetadataAdd(chipResults, PS_LIST_TAIL, cellName, PS_DATA_METADATA,
                               "Results for cell", cellResults);

trunk/ppStats/src/ppStatsSetup.c

@@ -63 +63 @@
         psMetadataIterator *iterator = psMetadataIteratorAlloc(arguments, PS_LIST_HEAD, "^" FLAG "$"); \
         psMetadataItem *item; \
-        while ((item = psMetadataGetAndIncrement(chipIter))) { \
+        while ((item = psMetadataGetAndIncrement(iterator))) { \
             if (item->type != PS_DATA_STRING) { \
                 psLogMsg(__func__, PS_LOG_WARN, NAME " name is not of type STRING (%x) --- ignored.\n", \
@@ -86 +86 @@
 
     STATS_OPTION("-mean",     PS_STAT_SAMPLE_MEAN);
-    STATS_OPTION("-stdev",    PS_STAT_SAMPLE_STDDEV);
+    STATS_OPTION("-stdev",    PS_STAT_SAMPLE_STDEV);
     STATS_OPTION("-median",   PS_STAT_SAMPLE_MEDIAN);
     STATS_OPTION("-quartile", PS_STAT_SAMPLE_QUARTILE);