Changeset 20711 for trunk/ppStack/src/ppStackMatch.c

Timestamp:

Nov 12, 2008, 4:59:10 PM (17 years ago)

Author:

Paul Price

Message:

Can't figure out why the weightings which are being placed on readout->analysis aren't getting to the combination stage, so I'm putting them in a vector and passing them explicitly.

File:

• trunk/ppStack/src/ppStackMatch.c (modified) (2 diffs)

trunk/ppStack/src/ppStackMatch.c

@@ -47 +47 @@
 #endif
 
-bool ppStackMatch(pmReadout *readout, psArray **regions, psArray **kernels, float *chi2,
+bool ppStackMatch(pmReadout *readout, psArray **regions, psArray **kernels, float *chi2, float *weighting,
                   const psArray *sources, const pmPSF *psf, psRandom *rng, const pmConfig *config)
 {
@@ -384 +384 @@
         (void)psBinaryOp(readout->image, readout->image, "+",
                          psScalarAlloc(- psStatsGetValue(bg, PS_STAT_ROBUST_MEDIAN), PS_TYPE_F32));
+        *weighting = 1.0 / PS_SQR(psStatsGetValue(bg, PS_STAT_ROBUST_STDEV));
         psMetadataAddF32(readout->analysis, PS_LIST_TAIL, "PPSTACK.WEIGHTING", 0,
                          "Weighting by 1/noise^2 for stack",