Changeset 34666 for branches/czw_branch/20120906/ppStack/src/ppStackCombineFinal.c

Timestamp:

Nov 15, 2012, 1:32:57 PM (14 years ago)

Author:

watersc1

Message:

Median code completes now, but is saving in the wrong output format.

File:

• branches/czw_branch/20120906/ppStack/src/ppStackCombineFinal.c (modified) (1 diff)

branches/czw_branch/20120906/ppStack/src/ppStackCombineFinal.c

@@ -3 +3 @@
 // This is the doomsday switch.
 // #define TESTING                         // Enable test output
-bool ppStackCombineBackground(ppStackThreadData *stack, ppStackOptions *options, pmConfig *config)
-{
-  psAssert(stack, "Require stack");
-  psAssert(options, "Require options");
-  psAssert(config, "Require configuration");
-
-  psTimerStart("PPSTACK_BKGMED");
-
-
-  pmReadout *bkgRO = options->bkgRO;
-  
-  psArray *inputs  = psArrayAlloc(options->num);
-  for (int i = 0; i < options->num; i++) {
-    ppStackFileActivationSingle(config, PPSTACK_FILES_BKG, true, i);
-    pmFPAfile *file = pmFPAfileSelectSingle(config->files, "PPSTACK.INPUT.BKGMODEL", i);
-    pmFPAview *view = ppStackFilesIterateDown(config);
-    pmReadout *ro = pmFPAviewThisReadout(view,file->fpa);
-    inputs->data[i] = ro;
-    pmFPAfileClose(file,view);
-  }
-  if (!pmStackSimpleMedianCombine(bkgRO,inputs)) {
-    psFree(inputs);
-    return(false);
-  }
-#if 0
-  if (!ppStackWriteImage("/tmp/test_forced.bkgmdl.fits",
-                         bkgRO->parent->parent->parent->hdu->header,
-                         bkgRO->image,
-                         config)) {
-    fprintf(stderr,"Failed to write image because fail.\n");
-  }
-#endif
-  for (int i = 0; i < options->num; i++) {
-    pmFPAfile *file = pmFPAfileSelectSingle(config->files, "PPSTACK.INPUT.BKGMODEL", i);
-    pmFPAview *view = ppStackFilesIterateDown(config);
-    bool success = pmFPAfileClose(file,view);
-    if (!success) {
-      psTrace("ppStack",5,"I failed at closing a file.\n");
-    }
-              
-    psFitsClose(file->fits);
-    file->fits = NULL;
-    file->header = NULL;
-    file->state = PM_FPA_STATE_CLOSED;
-    file->wrote_phu = false;
-    ppStackFileActivationSingle(config, PPSTACK_FILES_BKG, false, i);
-  }
-  psFree(inputs);
-  bkgRO->data_exists = true;
-  bkgRO->parent->data_exists = true;
-  bkgRO->parent->parent->data_exists = true;
-  
-  return(true);
-}
-  
-
-
 
 bool ppStackCombineFinal(ppStackThreadData *stack, psArray *covariances, ppStackOptions *options,