- Timestamp:
- Jan 29, 2009, 5:31:30 PM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/pap_branch_20090128/ppStack/src/ppStackThread.c
r20711 r21236 38 38 psFitsClose(stack->imageFits->data[i]); 39 39 psFitsClose(stack->maskFits->data[i]); 40 psFitsClose(stack-> weightFits->data[i]);41 stack->imageFits->data[i] = stack->maskFits->data[i] = stack-> weightFits->data[i] = NULL;40 psFitsClose(stack->varianceFits->data[i]); 41 stack->imageFits->data[i] = stack->maskFits->data[i] = stack->varianceFits->data[i] = NULL; 42 42 } 43 43 psFree(stack->imageFits); 44 44 psFree(stack->maskFits); 45 psFree(stack-> weightFits);45 psFree(stack->varianceFits); 46 46 return; 47 47 } 48 48 49 49 ppStackThreadData *ppStackThreadDataSetup(const psArray *cells, const psArray *imageNames, 50 const psArray *maskNames, const psArray * weightNames,50 const psArray *maskNames, const psArray *varianceNames, 51 51 const pmConfig *config) 52 52 { … … 54 54 PS_ASSERT_ARRAYS_SIZE_EQUAL(cells, imageNames, NULL); 55 55 PS_ASSERT_ARRAYS_SIZE_EQUAL(cells, maskNames, NULL); 56 PS_ASSERT_ARRAYS_SIZE_EQUAL(cells, weightNames, NULL);56 PS_ASSERT_ARRAYS_SIZE_EQUAL(cells, varianceNames, NULL); 57 57 58 58 ppStackThreadData *stack = psAlloc(sizeof(ppStackThreadData)); // Thread data, to return … … 64 64 stack->imageFits = psArrayAlloc(numInputs); 65 65 stack->maskFits = psArrayAlloc(numInputs); 66 stack-> weightFits = psArrayAlloc(numInputs);66 stack->varianceFits = psArrayAlloc(numInputs); 67 67 for (int i = 0; i < numInputs; i++) { 68 68 if (!cells->data[i]) { … … 73 73 psString imageResolved = pmConfigConvertFilename(imageNames->data[i], config, false, false); 74 74 psString maskResolved = pmConfigConvertFilename(maskNames->data[i], config, false, false); 75 psString weightResolved = pmConfigConvertFilename(weightNames->data[i], config, false, false);75 psString varianceResolved = pmConfigConvertFilename(varianceNames->data[i], config, false, false); 76 76 stack->imageFits->data[i] = psFitsOpen(imageResolved, "r"); 77 77 stack->maskFits->data[i] = psFitsOpen(maskResolved, "r"); 78 stack-> weightFits->data[i] = psFitsOpen(weightResolved, "r");78 stack->varianceFits->data[i] = psFitsOpen(varianceResolved, "r"); 79 79 psFree(imageResolved); 80 80 psFree(maskResolved); 81 psFree( weightResolved);82 if (!stack->imageFits->data[i] || !stack->maskFits->data[i] || !stack-> weightFits->data[i]) {81 psFree(varianceResolved); 82 if (!stack->imageFits->data[i] || !stack->maskFits->data[i] || !stack->varianceFits->data[i]) { 83 83 psError(PS_ERR_UNKNOWN, false, "Unable to open convolved files %s, %s, %s", 84 (char*)imageNames->data[i], (char*)maskNames->data[i], (char*) weightNames->data[i]);84 (char*)imageNames->data[i], (char*)maskNames->data[i], (char*)varianceNames->data[i]); 85 85 return false; 86 86 } … … 156 156 // override the recorded last scan 157 157 ro->thisImageScan = thread->firstScan; 158 ro->this WeightScan = thread->firstScan;158 ro->thisVarianceScan = thread->firstScan; 159 159 ro->thisMaskScan = thread->firstScan; 160 160 ro->lastImageScan = thread->lastScan; 161 161 ro->lastMaskScan = thread->lastScan; 162 ro->last WeightScan = thread->lastScan;162 ro->lastVarianceScan = thread->lastScan; 163 163 ro->forceScan = true; 164 164 165 165 psFits *imageFits = stack->imageFits->data[i]; // FITS file for image 166 166 psFits *maskFits = stack->maskFits->data[i]; // FITS file for mask 167 psFits * weightFits = stack->weightFits->data[i]; // FITS file for weight167 psFits *varianceFits = stack->varianceFits->data[i]; // FITS file for variance 168 168 169 169 … … 189 189 } 190 190 191 if (pmReadoutMore Weight(ro, weightFits, 0, rows, config)) {191 if (pmReadoutMoreVariance(ro, varianceFits, 0, rows, config)) { 192 192 keepReading = true; 193 if (!pmReadoutReadChunkWeight(ro, weightFits, 0, rows, overlap, config)) { 194 psError(PS_ERR_IO, false, "Unable to read chunk %d for file PPSTACK.INPUT.WEIGHT %d", 193 if (!pmReadoutReadChunkVariance(ro, varianceFits, 0, rows, overlap, config)) { 194 psError(PS_ERR_IO, false, 195 "Unable to read chunk %d for file PPSTACK.INPUT.VARIANCE %d", 195 196 numChunk, i); 196 197 *status = false;
Note:
See TracChangeset
for help on using the changeset viewer.
