Changeset 17556

Timestamp:

May 7, 2008, 10:54:46 AM (18 years ago)

Author:

eugene

Message:

add sum below peak

Location:

trunk/psastro/src

Files:

• psastro.h (modified) (1 diff)
• psastroLuminosityFunction.c (modified) (3 diffs)

trunk/psastro/src/psastro.h

@@ -23 +23 @@
     double mPeak;                       // mag of peak bin 
     int nPeak;                          // # of stars in peak bin
+    int sPeak;                          // sum of stars to peak bin
 } pmLumFunc;
 

trunk/psastro/src/psastroLuminosityFunction.c

@@ -22 +22 @@
   func->mPeak = 0.0;
   func->nPeak = 0;
+  func->sPeak = 0;
 
   return func;
@@ -62 +63 @@
   int iPeak = 0;
   int nPeak = 0;
+  int sPeak = 0;
   for (int i = 0; i < nMags->n; i++) {
     if (nMags->data.F32[i] < nPeak) continue;
     iPeak = i;
     nPeak = nMags->data.F32[i];
+    sPeak += nPeak;
   }
 
@@ -122 +125 @@
   lumFunc->mPeak = mMin + (iPeak + 0.5)*dMag;
   lumFunc->nPeak = nPeak;
+  lumFunc->sPeak = sPeak;
 
   psFree (lnMag);