Index: trunk/archive/modules/src/pmAstrom.c
===================================================================
--- trunk/archive/modules/src/pmAstrom.c	(revision 5272)
+++ trunk/archive/modules/src/pmAstrom.c	(revision 5353)
@@ -223,2 +223,29 @@
     return (match);
 }
+
+pmAstromMatchedListFit (psArray *st1, psArray *st2, psArray *match, pmAstromOpt *opt) {
+
+  psVector *X = psVectorAlloc (match->n);
+  psVector *Y = psVectorAlloc (match->n);
+  psVector *x = psVectorAlloc (match->n);
+  psVector *y = psVectorAlloc (match->n);
+  
+  // take the matched stars, first fit 
+  for (i = 0; i < match->n; i++) {
+
+    pair = match->data[i];
+    ob1 = st1->data[pair->i1];
+    ob2 = st2->data[pair->i2];
+
+    X->data.F64[i] = ob1->P;
+    Y->data.F64[i] = ob1->Q;
+
+    x->data.F64[i] = ob2->P;
+    y->data.F64[i] = ob2->Q;
+
+    // use one or the other...
+    psPolynomial2D *xt = psVectorFitPolynomial2D (NULL, stats, NULL, 0, X, NULL, x, y)
+    psPolynomial2D *xt = psVectorClipFitPolynomial2D (NULL, stats, NULL, 0, X, NULL, x, y)
+
+    psPolynomial2D *yt = psVectorFitPolynomial2D (NULL, stats, NULL, 0, Y, NULL, x, y)
+    psPolynomial2D *yt = psVectorClipFitPolynomial2D (NULL, stats, NULL, 0, Y, NULL, x, y)