Index: trunk/ppTranslate/src/ppMopsDetections.c
===================================================================
--- trunk/ppTranslate/src/ppMopsDetections.c	(revision 28243)
+++ trunk/ppTranslate/src/ppMopsDetections.c	(revision 28623)
@@ -10,26 +10,46 @@
 static void mopsDetectionsFree(ppMopsDetections *det)
 {
-    psFree(det->component);
     psFree(det->raBoresight);
     psFree(det->decBoresight);
     psFree(det->filter);
-    psFree(det->header);
-    psFree(det->table);
     psFree(det->x);
     psFree(det->y);
     psFree(det->ra);
     psFree(det->dec);
-    psFree(det->deteffHeader);
-    psFree(det->deteffTable);
+    psFree(det->raErr);
+    psFree(det->decErr);
+    psFree(det->mag);
+    psFree(det->magErr);
+    psFree(det->chi2);
+    psFree(det->dof);
+    psFree(det->cr);
+    psFree(det->extended);
+    psFree(det->psfMajor);
+    psFree(det->psfMinor);
+    psFree(det->psfTheta);
+    psFree(det->quality);
+    psFree(det->numPix);
+    psFree(det->xxMoment);
+    psFree(det->xyMoment);
+    psFree(det->yyMoment);
+    psFree(det->flags);
+    psFree(det->diffSkyfileId);
+    psFree(det->naxis1);
+    psFree(det->naxis2);
+    psFree(det->mask);
+    psFree(det->nPos);
+    psFree(det->fPos);
+    psFree(det->nRatioBad);
+    psFree(det->nRatioMask);
+    psFree(det->nRatioAll);
 
     return;
 }
 
-ppMopsDetections *ppMopsDetectionsAlloc(void)
+ppMopsDetections *ppMopsDetectionsAlloc(long num)
 {
     ppMopsDetections *det = psAlloc(sizeof(ppMopsDetections)); // Detections, to return
     psMemSetDeallocator(det, (psFreeFunc)mopsDetectionsFree);
 
-    det->component = NULL;
     det->raBoresight = NULL;
     det->decBoresight = NULL;
@@ -43,13 +63,234 @@
     det->seeing = NAN;
     det->num = 0;
-    det->header = NULL;
-    det->table = NULL;
-    det->x = NULL;
-    det->y = NULL;
-    det->ra = NULL;
-    det->dec = NULL;
-    det->deteffHeader = NULL;
-    det->deteffTable = NULL;
+    det->x = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->y = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->ra = psVectorAllocEmpty(num, PS_TYPE_F64);
+    det->dec = psVectorAllocEmpty(num, PS_TYPE_F64);
+    det->raErr = psVectorAllocEmpty(num, PS_TYPE_F64);
+    det->decErr = psVectorAllocEmpty(num, PS_TYPE_F64);
+    det->mag = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->magErr = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->chi2 = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->dof = psVectorAllocEmpty(num, PS_TYPE_S32);
+    det->cr = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->extended = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->psfMajor = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->psfMinor = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->psfTheta = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->quality = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->numPix = psVectorAllocEmpty(num, PS_TYPE_S32);
+    det->xxMoment = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->xyMoment = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->yyMoment = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->flags = psVectorAllocEmpty(num, PS_TYPE_U32);
+    det->diffSkyfileId = psVectorAllocEmpty(num, PS_TYPE_S64);
+    det->naxis1 = psVectorAllocEmpty(num, PS_TYPE_S32);
+    det->naxis2 = psVectorAllocEmpty(num, PS_TYPE_S32);
+    det->mask = psVectorAllocEmpty(num, PS_TYPE_U8);
+    det->nPos = psVectorAllocEmpty(num, PS_TYPE_S32);
+    det->fPos = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->nRatioBad = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->nRatioMask = psVectorAllocEmpty(num, PS_TYPE_F32);
+    det->nRatioAll = psVectorAllocEmpty(num, PS_TYPE_F32);
 
     return det;
 }
+
+
+ppMopsDetections *ppMopsDetectionsRealloc(ppMopsDetections *det, long num)
+{
+    det->x = psVectorRealloc(det->x, num);
+    det->y = psVectorRealloc(det->y, num);
+    det->ra = psVectorRealloc(det->ra, num);
+    det->dec = psVectorRealloc(det->dec, num);
+    det->raErr = psVectorRealloc(det->raErr, num);
+    det->decErr = psVectorRealloc(det->decErr, num);
+    det->mag = psVectorRealloc(det->mag, num);
+    det->magErr = psVectorRealloc(det->magErr, num);
+    det->chi2 = psVectorRealloc(det->chi2, num);
+    det->dof = psVectorRealloc(det->dof, num);
+    det->cr = psVectorRealloc(det->cr, num);
+    det->extended = psVectorRealloc(det->extended, num);
+    det->psfMajor = psVectorRealloc(det->psfMajor, num);
+    det->psfMinor = psVectorRealloc(det->psfMinor, num);
+    det->psfTheta = psVectorRealloc(det->psfTheta, num);
+    det->quality = psVectorRealloc(det->quality, num);
+    det->numPix = psVectorRealloc(det->numPix, num);
+    det->xxMoment = psVectorRealloc(det->xxMoment, num);
+    det->xyMoment = psVectorRealloc(det->xyMoment, num);
+    det->yyMoment = psVectorRealloc(det->yyMoment, num);
+    det->flags = psVectorRealloc(det->flags, num);
+    det->diffSkyfileId = psVectorRealloc(det->diffSkyfileId, num);
+    det->naxis1 = psVectorRealloc(det->naxis1, num);
+    det->naxis2 = psVectorRealloc(det->naxis2, num);
+    det->mask = psVectorRealloc(det->mask, num);
+    det->nPos = psVectorRealloc(det->nPos, num);
+    det->fPos = psVectorRealloc(det->fPos, num);
+    det->nRatioBad = psVectorRealloc(det->nRatioBad, num);
+    det->nRatioMask = psVectorRealloc(det->nRatioMask, num);
+    det->nRatioAll = psVectorRealloc(det->nRatioAll, num);
+
+    return det;
+}
+
+
+bool ppMopsDetectionsAdd(ppMopsDetections *det, float x, float y, double ra, double dec,
+                         double raErr, double decErr, float mag, float magErr,
+                         float chi2, int dof, float cr, float extended, float psfMajor,
+                         float psfMinor, float psfTheta, float quality, int numPix,
+                         float xxMoment, float xyMoment, float yyMoment,
+                         psU32 flags, psS64 diffSkyfileId, int naxis1, int naxis2,
+                         int nPos, float fPos, float nRatioBad, float nRatioMask, float nRatioAll)
+{
+    psVectorAppend(det->x, x);
+    psVectorAppend(det->y, y);
+    psVectorAppend(det->ra, ra);
+    psVectorAppend(det->dec, dec);
+    psVectorAppend(det->raErr, raErr);
+    psVectorAppend(det->decErr, decErr);
+    psVectorAppend(det->mag, mag);
+    psVectorAppend(det->magErr, magErr);
+    psVectorAppend(det->chi2, chi2);
+    psVectorAppend(det->dof, dof);
+    psVectorAppend(det->cr, cr);
+    psVectorAppend(det->extended, extended);
+    psVectorAppend(det->psfMajor, psfMajor);
+    psVectorAppend(det->psfMinor, psfMinor);
+    psVectorAppend(det->psfTheta, psfTheta);
+    psVectorAppend(det->quality, quality);
+    psVectorAppend(det->numPix, numPix);
+    psVectorAppend(det->xxMoment, xxMoment);
+    psVectorAppend(det->xyMoment, xyMoment);
+    psVectorAppend(det->yyMoment, yyMoment);
+    psVectorAppend(det->flags, flags);
+    psVectorAppend(det->diffSkyfileId, diffSkyfileId);
+    psVectorAppend(det->naxis1, naxis1);
+    psVectorAppend(det->naxis2, naxis2);
+    psVectorAppend(det->mask, 0);
+    psVectorAppend(det->nPos, nPos);
+    psVectorAppend(det->fPos, fPos);
+    psVectorAppend(det->nRatioBad, nRatioBad);
+    psVectorAppend(det->nRatioMask, nRatioMask);
+    psVectorAppend(det->nRatioAll, nRatioAll);
+
+    return true;
+}
+
+
+bool ppMopsDetectionsCopySingle(ppMopsDetections *target, const ppMopsDetections *source, long index)
+{
+    psVectorAppend(target->x, source->x->data.F32[index]);
+    psVectorAppend(target->y, source->y->data.F32[index]);
+    psVectorAppend(target->ra, source->ra->data.F64[index]);
+    psVectorAppend(target->dec, source->dec->data.F64[index]);
+    psVectorAppend(target->raErr, source->raErr->data.F64[index]);
+    psVectorAppend(target->decErr, source->decErr->data.F64[index]);
+    psVectorAppend(target->mag, source->mag->data.F32[index]);
+    psVectorAppend(target->magErr, source->magErr->data.F32[index]);
+    psVectorAppend(target->chi2, source->chi2->data.F32[index]);
+    psVectorAppend(target->dof, source->dof->data.S32[index]);
+    psVectorAppend(target->cr, source->cr->data.F32[index]);
+    psVectorAppend(target->extended, source->extended->data.F32[index]);
+    psVectorAppend(target->psfMajor, source->psfMajor->data.F32[index]);
+    psVectorAppend(target->psfMinor, source->psfMinor->data.F32[index]);
+    psVectorAppend(target->psfTheta, source->psfTheta->data.F32[index]);
+    psVectorAppend(target->quality, source->quality->data.F32[index]);
+    psVectorAppend(target->numPix, source->numPix->data.S32[index]);
+    psVectorAppend(target->xxMoment, source->xxMoment->data.F32[index]);
+    psVectorAppend(target->xyMoment, source->xyMoment->data.F32[index]);
+    psVectorAppend(target->yyMoment, source->yyMoment->data.F32[index]);
+    psVectorAppend(target->flags, source->flags->data.U32[index]);
+    psVectorAppend(target->diffSkyfileId, source->diffSkyfileId->data.S64[index]);
+    psVectorAppend(target->naxis1, source->naxis1->data.S32[index]);
+    psVectorAppend(target->naxis2, source->naxis2->data.S32[index]);
+    psVectorAppend(target->mask, 0);
+    psVectorAppend(target->nPos, source->nPos->data.S32[index]);
+    psVectorAppend(target->fPos, source->fPos->data.F32[index]);
+    psVectorAppend(target->nRatioBad, source->nRatioBad->data.F32[index]);
+    psVectorAppend(target->nRatioMask, source->nRatioMask->data.F32[index]);
+    psVectorAppend(target->nRatioAll, source->nRatioAll->data.F32[index]);
+
+    target->num++;
+
+    return true;
+}
+
+
+bool ppMopsDetectionsPurge(ppMopsDetections *det)
+{
+    long num = 0;
+    for (long i = 0; i < det->num; i++) {
+        if (!det->mask->data.U8[i]) {
+            if (i == num) {
+                // No need to copy
+                num++;
+                continue;
+            }
+            det->x->data.F32[num] = det->x->data.F32[i];
+            det->y->data.F32[num] = det->y->data.F32[i];
+            det->ra->data.F64[num] = det->ra->data.F64[i];
+            det->dec->data.F64[num] = det->dec->data.F64[i];
+            det->raErr->data.F64[num] = det->raErr->data.F64[i];
+            det->decErr->data.F64[num] = det->decErr->data.F64[i];
+            det->mag->data.F32[num] = det->mag->data.F32[i];
+            det->magErr->data.F32[num] = det->magErr->data.F32[i];
+            det->chi2->data.F32[num] = det->chi2->data.F32[i];
+            det->dof->data.S32[num] = det->dof->data.S32[i];
+            det->cr->data.F32[num] = det->cr->data.F32[i];
+            det->extended->data.F32[num] = det->extended->data.F32[i];
+            det->psfMajor->data.F32[num] = det->psfMajor->data.F32[i];
+            det->psfMinor->data.F32[num] = det->psfMinor->data.F32[i];
+            det->psfTheta->data.F32[num] = det->psfTheta->data.F32[i];
+            det->quality->data.F32[num] = det->quality->data.F32[i];
+            det->numPix->data.S32[num] = det->numPix->data.S32[i];
+            det->xxMoment->data.F32[num] = det->xxMoment->data.F32[i];
+            det->xyMoment->data.F32[num] = det->xyMoment->data.F32[i];
+            det->yyMoment->data.F32[num] = det->yyMoment->data.F32[i];
+            det->flags->data.U32[num] = det->flags->data.U32[i];
+            det->diffSkyfileId->data.S64[num] = det->diffSkyfileId->data.S64[i];
+            det->naxis1->data.S32[num] = det->naxis1->data.S32[i];
+            det->naxis2->data.S32[num] = det->naxis2->data.S32[i];
+            det->mask->data.U8[num] = 0;
+            det->nPos->data.S32[num] = det->nPos->data.S32[i];
+            det->fPos->data.F32[num] = det->fPos->data.F32[i];
+            det->nRatioBad->data.F32[num] = det->nRatioBad->data.F32[i];
+            det->nRatioMask->data.F32[num] = det->nRatioMask->data.F32[i];
+            det->nRatioAll->data.F32[num] = det->nRatioAll->data.F32[i];
+            num++;
+        }
+    }
+    det->x->n = num;
+    det->y->n = num;
+    det->ra->n = num;
+    det->dec->n = num;
+    det->raErr->n = num;
+    det->decErr->n = num;
+    det->mag->n = num;
+    det->magErr->n = num;
+    det->chi2->n = num;
+    det->dof->n = num;
+    det->cr->n = num;
+    det->extended->n = num;
+    det->psfMajor->n = num;
+    det->psfMinor->n = num;
+    det->psfTheta->n = num;
+    det->quality->n = num;
+    det->numPix->n = num;
+    det->xxMoment->n = num;
+    det->xyMoment->n = num;
+    det->yyMoment->n = num;
+    det->flags->n = num;
+    det->diffSkyfileId->n = num;
+    det->naxis1->n = num;
+    det->naxis2->n = num;
+    det->mask->n = num;
+    det->num = num;
+    det->nPos->n = num;
+    det->fPos->n = num;
+    det->nRatioBad->n = num;
+    det->nRatioMask->n = num;
+    det->nRatioAll->n = num;
+
+    return true;
+}
+