Changeset 36026

Timestamp:

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

Author:

eugene

Message:

mv pmSersicNorm to psModules

File:

• branches/eam_branches/ipp-20130711/ppSim/src/ppSimMakeGalaxies.c (modified) (1 diff)

branches/eam_branches/ipp-20130711/ppSim/src/ppSimMakeGalaxies.c

@@ -202 +202 @@
     }
     return true;
-}
-
-float pmSersicNorm (float index) {
-
-    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 ((index >= 0.0) && (index < 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 ((index >= 1.0) && (index < 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 ((index >= 2.0) && (index < 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 ((index >= 3.0) && (index < 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 ((index >= 4.0) && (index < 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*index + C2*index*index;
-    float norm = exp(lnorm);
-    return norm;
 }