Changeset 35646 for branches/eam_branches/ipp-20130509/psModules/src/objects/pmPCMdata.c

Timestamp:

Jun 11, 2013, 9:23:27 AM (13 years ago)

Author:

eugene

Message:

new functions pmModelUseReff, pmModelAxesToParams, pmModelParamsToAxes; add useReff to pmModelSetShape; remove deprecated code (pmPSF_SXYtoModel, pmPSF_SXYtoModel); replace code which converts the PAR[SXX,SXY,SYY] values to shape and then axes with calls to pmModelParamsToAxes and pmModelAxesToParams

File:

• branches/eam_branches/ipp-20130509/psModules/src/objects/pmPCMdata.c (modified) (2 diffs)

branches/eam_branches/ipp-20130509/psModules/src/objects/pmPCMdata.c

@@ -291 +291 @@
     psAssert (modelPSF, "psf model must be defined");
     
-    psEllipseShape shape;
     psEllipseAxes axes;
-
-    shape.sx  = modelPSF->params->data.F32[PM_PAR_SXX];
-    shape.sy  = modelPSF->params->data.F32[PM_PAR_SYY];
-    shape.sxy = modelPSF->params->data.F32[PM_PAR_SXY];
-    axes = psEllipseShapeToAxes (shape, 20.0);
+    bool useReff = pmModelUseReff (modelPSF->type);
+    psF32 *PAR = modelPSF->params->data.F32;
+    pmModelParamsToAxes (&axes, PAR[PM_PAR_SXX], PAR[PM_PAR_SXY], PAR[PM_PAR_SYY], useReff);
     
-    float FWHM_MAJOR = 2*modelPSF->modelRadius (modelPSF->params, 0.5*modelPSF->params->data.F32[PM_PAR_I0]);
+    float FWHM_MAJOR = 2*modelPSF->modelRadius (modelPSF->params, 0.5*PAR[PM_PAR_I0]);
     float FWHM_MINOR = FWHM_MAJOR * (axes.minor / axes.major);
 
@@ -451 +448 @@
         psAssert (modelPSF, "psf model must be defined");
     
-        psEllipseShape shape;
         psEllipseAxes axes;
-
-        shape.sx  = modelPSF->params->data.F32[PM_PAR_SXX];
-        shape.sy  = modelPSF->params->data.F32[PM_PAR_SYY];
-        shape.sxy = modelPSF->params->data.F32[PM_PAR_SXY];
-        axes = psEllipseShapeToAxes (shape, 20.0);
+        bool useReff = pmModelUseReff (modelPSF->type);
+        psF32 *PAR = modelPSF->params->data.F32;
+        pmModelParamsToAxes (&axes, PAR[PM_PAR_SXX], PAR[PM_PAR_SXY], PAR[PM_PAR_SYY], useReff);
     
-        float FWHM_MAJOR = 2*modelPSF->modelRadius (modelPSF->params, 0.5*modelPSF->params->data.F32[PM_PAR_I0]);
+        float FWHM_MAJOR = 2*modelPSF->modelRadius (modelPSF->params, 0.5*PAR[PM_PAR_I0]);
         float FWHM_MINOR = FWHM_MAJOR * (axes.minor / axes.major);