Changeset 36021 for branches/eam_branches/ipp-20130711/psModules/src/objects/pmModel_CentralPixel.c

Timestamp:

Aug 25, 2013, 12:08:56 PM (13 years ago)

Author:

eugene

Message:

convert SERSIC to use pmModelCP_SersicSubpix; clean up DEV and EXP versions; add pmPCMMakeModel, pmSersicKappa, pmSersicNorm functions

File:

• branches/eam_branches/ipp-20130711/psModules/src/objects/pmModel_CentralPixel.c (modified) (1 diff)

branches/eam_branches/ipp-20130711/psModules/src/objects/pmModel_CentralPixel.c

@@ -723 +723 @@
 }
 
+float pmSersicKappa (float Sindex) {
+    // this function is empirically derived from a fit to data for Sindex = 0.5 - 5.5
+    // constrain Sindex or kappa?
+    float kappa = -0.275552 + 1.972625*Sindex + 0.003487 * PS_SQR(Sindex);
+    return kappa;
+}
+
+float pmSersicNorm (float Sindex) {
+
+    float C0 = NAN;
+    float C1 = NAN;
+    float C2 = NAN;
+
+    // y = 0.201545 x^0 -0.950965 x^1 -0.315248 x^2 
+    // y = 0.402084 x^0 -1.357775 x^1 -0.105102 x^2 
+    // y = 0.619093 x^0 -1.591674 x^1 -0.041576 x^2 
+    // y = 0.770263 x^0 -1.696421 x^1 -0.023363 x^2 
+    // y = 0.885891 x^0 -1.755684 x^1 -0.015753 x^2 
+
+    if ((Sindex >= 0.0) && (Sindex < 1.0)) { 
+        C0 = 0.201545; C1 = -0.950965; C2 = -0.315248;
+        // y = 0.201545 x^0 -0.950965 x^1 -0.315248 x^2 
+    }
+    if ((Sindex >= 1.0) && (Sindex < 2.0)) { 
+        C0 = 0.402084; C1 = -1.357775; C2 = -0.105102;
+        // y = 0.402084 x^0 -1.357775 x^1 -0.105102 x^2 
+    }
+    if ((Sindex >= 2.0) && (Sindex < 3.0)) { 
+        C0 = 0.619093; C1 = -1.591674; C2 = -0.041576;
+        // y = 0.619093 x^0 -1.591674 x^1 -0.041576 x^2 
+    }
+    if ((Sindex >= 3.0) && (Sindex < 4.0)) { 
+        C0 = 0.770263; C1 = -1.696421; C2 = -0.023363;
+        // y = 0.770263 x^0 -1.696421 x^1 -0.023363 x^2 
+    }
+    if ((Sindex >= 4.0) && (Sindex < 5.5)) { 
+        C0 = 0.885891; C1 = -1.755684; C2 = -0.015753; 
+        // y = 0.885891 x^0 -1.755684 x^1 -0.015753 x^2 
+    }
+
+    if (isnan(C0)) return NAN;
+
+    float lnorm = C0 + C1*Sindex + C2*Sindex*Sindex;
+    float norm = exp(lnorm);
+    return norm;
+}
+
 # if (0)
 // create a vector containing only the unique entries in the input vector