Changeset 10267

Timestamp:

Nov 28, 2006, 5:03:53 PM (20 years ago)

Author:

eugene

Message:

updates to psEllipse functions

Location:

trunk/psphot/src

Files:

• psphotChoosePSF.c (modified) (2 diffs)
• psphotModelTest.c (modified) (2 diffs)

trunk/psphot/src/psphotChoosePSF.c

@@ -288 +288 @@
     shape.sy  = modelPSF->params->data.F32[PM_PAR_SYY];
     shape.sxy = modelPSF->params->data.F32[PM_PAR_SXY];
-    axes = psEllipseShapeToAxes (shape);
+    axes = psEllipseShapeToAxes (shape, 20.0);
 
     psF64 FWHM_Y = FWHM_X * (axes.minor / axes.major);
@@ -330 +330 @@
         moments.xy = source->moments->Sxy;
 
-        axes = psEllipseMomentsToAxes (moments);
+        // limit axis ratio < 20.0
+        axes = psEllipseMomentsToAxes (moments, 20.0);
 
         FWHM_X += axes.major * 2.35;

trunk/psphot/src/psphotModelTest.c

@@ -107 +107 @@
     moments.y2 = source->moments->Sy;
     moments.xy = source->moments->Sxy;
-    psEllipseAxes axes = psEllipseMomentsToAxes (moments);
+    psEllipseAxes axes = psEllipseMomentsToAxes (moments, 20.0);
 
     fprintf (stderr, "axes: %f @ (%f, %f)\n", axes.theta*180/M_PI, axes.major, axes.minor);
@@ -147 +147 @@
     shape.sy  = 1.4 / model->params->data.F32[5];
     shape.sxy = model->params->data.F32[6];
-    axes = psEllipseShapeToAxes (shape);
+    axes = psEllipseShapeToAxes (shape, 20.0);
 
     fprintf (stderr, "guess: %f @ (%f, %f)\n", axes.theta*180/M_PI, axes.major, axes.minor);