IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Jan 29, 2009, 5:31:30 PM (17 years ago)
Author:
Paul Price
Message:

Changing 'weight' to 'variance'. Mostly propagating psModules changes.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/pap_branch_20090128/ppStack/src/ppStackThread.c

    r20711 r21236  
    3838        psFitsClose(stack->imageFits->data[i]);
    3939        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;
    4242    }
    4343    psFree(stack->imageFits);
    4444    psFree(stack->maskFits);
    45     psFree(stack->weightFits);
     45    psFree(stack->varianceFits);
    4646    return;
    4747}
    4848
    4949ppStackThreadData *ppStackThreadDataSetup(const psArray *cells, const psArray *imageNames,
    50                                           const psArray *maskNames, const psArray *weightNames,
     50                                          const psArray *maskNames, const psArray *varianceNames,
    5151                                          const pmConfig *config)
    5252{
     
    5454    PS_ASSERT_ARRAYS_SIZE_EQUAL(cells, imageNames, NULL);
    5555    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);
    5757
    5858    ppStackThreadData *stack = psAlloc(sizeof(ppStackThreadData)); // Thread data, to return
     
    6464    stack->imageFits  = psArrayAlloc(numInputs);
    6565    stack->maskFits   = psArrayAlloc(numInputs);
    66     stack->weightFits = psArrayAlloc(numInputs);
     66    stack->varianceFits = psArrayAlloc(numInputs);
    6767    for (int i = 0; i < numInputs; i++) {
    6868        if (!cells->data[i]) {
     
    7373        psString imageResolved = pmConfigConvertFilename(imageNames->data[i], config, false, false);
    7474        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);
    7676        stack->imageFits->data[i] = psFitsOpen(imageResolved, "r");
    7777        stack->maskFits->data[i] = psFitsOpen(maskResolved, "r");
    78         stack->weightFits->data[i] = psFitsOpen(weightResolved, "r");
     78        stack->varianceFits->data[i] = psFitsOpen(varianceResolved, "r");
    7979        psFree(imageResolved);
    8080        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]) {
    8383            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]);
    8585            return false;
    8686        }
     
    156156                // override the recorded last scan
    157157                ro->thisImageScan  = thread->firstScan;
    158                 ro->thisWeightScan = thread->firstScan;
     158                ro->thisVarianceScan = thread->firstScan;
    159159                ro->thisMaskScan   = thread->firstScan;
    160160                ro->lastImageScan  = thread->lastScan;
    161161                ro->lastMaskScan   = thread->lastScan;
    162                 ro->lastWeightScan = thread->lastScan;
     162                ro->lastVarianceScan = thread->lastScan;
    163163                ro->forceScan      = true;
    164164
    165165                psFits *imageFits  = stack->imageFits->data[i]; // FITS file for image
    166166                psFits *maskFits   = stack->maskFits->data[i]; // FITS file for mask
    167                 psFits *weightFits = stack->weightFits->data[i]; // FITS file for weight
     167                psFits *varianceFits = stack->varianceFits->data[i]; // FITS file for variance
    168168
    169169
     
    189189                }
    190190
    191                 if (pmReadoutMoreWeight(ro, weightFits, 0, rows, config)) {
     191                if (pmReadoutMoreVariance(ro, varianceFits, 0, rows, config)) {
    192192                    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",
    195196                                numChunk, i);
    196197                        *status = false;
Note: See TracChangeset for help on using the changeset viewer.