Changeset 2650 for trunk/psLib/src/dataManip/psMinimize.c

Timestamp:

Dec 7, 2004, 12:28:29 PM (22 years ago)

Author:

gusciora

Message:

...

File:

• trunk/psLib/src/dataManip/psMinimize.c (modified) (2 diffs)

trunk/psLib/src/dataManip/psMinimize.c

@@ -9 +9 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.92 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2004-12-01 19:56:06 $
+ *  @version $Revision: 1.93 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2004-12-07 22:28:29 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -455 +455 @@
 }
 
+/*
+XXX: from bug 230:
+ 
+We first perform a rotation:
+u = - (x-x0)*cos(theta) + (y-y0)*sin(theta)
+v = (x-x0)*cos(theta) + (y-y0)*sin(theta)
+ 
+Here u is the major axis, and v is the minor axis, x0,y0 is the centre, and
+theta is the position angle.
+ 
+Then the flux is
+ 
+flux = norm * exp(-( u*u/2.0/sigmau/sigmau + v*v/2.0/sigmav/sigmav)
+)/2.0/pi/sigmau/sigmav
+ 
+Here sigmau and sigmav are the widths of the major and minor axes.
+ 
+The "norm" parameter in the equation above corresponds to the normalisation.
+ 
+Suggest order:
+ 
+norm
+x0
+y0
+sigma_u
+sigma_v
+theta
+*/
 
 psVector *psMinimizeLMChi2Gauss2D(psImage *deriv,