Index: trunk/psphot/src/pmSourceUtils.c
===================================================================
--- trunk/psphot/src/pmSourceUtils.c	(revision 4954)
+++ trunk/psphot/src/pmSourceUtils.c	(revision 4977)
@@ -1,44 +1,1 @@
 # include "psphot.h"
-
-bool pmSourceDefinePixels(pmSource *mySource, 
-			  const psImageData *imdata,
-			  psF32 x, 
-			  psF32 y,
-			  psF32 Radius)
-{
-    psRegion srcRegion;
-
-    // Grab a subimage of the original image of size (2 * outerRadius).
-    srcRegion = psRegionForSquare (x, y, Radius);
-    srcRegion = psRegionForImage (imdata->image, &srcRegion);
-
-    mySource->pixels = psImageSubset(imdata->image, srcRegion);
-    mySource->noise  = psImageSubset(imdata->noise, srcRegion);
-    mySource->mask   = psImageSubset(imdata->mask,  srcRegion);
-
-    return(mySource);
-}
-
-bool pmSourcePhotometry (float *fitMag, float *obsMag, pmModel *model, psImage *image, psImage *mask) {
-
-    float obsSum = 0;
-    float fitSum = 0;
-    float sky = model->params->data.F32[0];
-
-    pmModelFlux modelFluxFunc = pmModelFlux_GetFunction (model->type);
-    fitSum = modelFluxFunc (model->params);
-
-    for (int ix = 0; ix < image->numCols; ix++) {
-	for (int iy = 0; iy < image->numRows; iy++) {
-	    if (mask->data.U8[iy][ix]) continue;
-	    obsSum += image->data.F32[iy][ix] - sky;
-	}
-    }
-    if (obsSum <= 0) return false;
-    if (fitSum <= 0) return false;
-
-    *fitMag = -2.5*log10(fitSum);
-    *obsMag = -2.5*log10(obsSum);
-    return (true);
-} 
-