Changeset 9770 for trunk/psModules/src/objects/pmSourceIO_CMP.c

Timestamp:

Oct 28, 2006, 10:23:51 AM (20 years ago)

Author:

magnier

Message:

This is an important change to the PSF shape model.

• I have changed the definition of the elliptical contour parameters

(PM_PAR_SXX,SYY are now sigma_x,y/sqrt(2), rather than
sqrt(2)/sigma_x.

• I have changed the way in which the SXY term is modelled. rather

than fit SXY to a polynomial, I now fit SXY / (1/SXX^{2 + 1/SYY}2)^{2.
this representation is more likely to be well-described by a
polynomial. the other form is likely to vary as xy/r}N rather than xy
r^N.

these changes imply that the models all need to treat the SXX,SYY,SXY
terms in the same fashion (ie, fitting a contour of the form (x/SXX)<sup>2
+ (y/SYY)</sup>2 + SXY*x*y). If an arbitrary model uses some other meaning
for SXX, SYY, SXY, then it will have to define its own exceptions in
the functions pmPSFFromModels and pmModelFromPSF (the latter has an
implementation for each model already).

• I have cleaned up the pmModel_NAME.c functions to use the #define names

for the parameters (PM_PAR_XXX) everywhere, and to use #def names for
the functions. this latter change makes it easy to specify the model
function names in a single location.

File:

• trunk/psModules/src/objects/pmSourceIO_CMP.c (modified) (3 diffs)

trunk/psModules/src/objects/pmSourceIO_CMP.c

@@ -3 +3 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.17 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-10-19 21:16:49 $
+ *  @version $Revision: 1.18 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-10-28 20:23:51 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -35 +35 @@
 // XXX make sure the angle in correctly translated to/from degrees
 // XXX we lose all information from the 'type' field
+
+// XXX update this file is we convert to PAR[4] : SigmaX*sqrt(2) (not 1/SigmaX)
 
 // elixir-style pseudo FITS table (header + ascii list)
@@ -267 +269 @@
             shape = psEllipseAxesToShape (axes);
 
-            PAR[4] = shape.sx;
-            PAR[5] = shape.sy;
-            PAR[6] = shape.sxy;
+            PAR[PM_PAR_SXX] = shape.sx;
+            PAR[PM_PAR_SYY] = shape.sy;
+            PAR[PM_PAR_SXY] = shape.sxy;
 
             source->modelPSF = model;