Index: trunk/psastro/src/psastroLuminosityFunction.c
===================================================================
--- trunk/psastro/src/psastroLuminosityFunction.c	(revision 10880)
+++ trunk/psastro/src/psastroLuminosityFunction.c	(revision 11269)
@@ -49,5 +49,8 @@
   for (int i = 0; i < stars->n; i++) {
     star = stars->data[i];
+    if (!isfinite(star->Mag)) continue;
     int bin = (star->Mag - mMin) / dMag;
+    if (bin < 0) psAbort ("psastro", "bin cannot be negative!");
+    if (bin >= nBin) psAbort ("psastro", "bin cannot be > %d!", nBin);
     nMags->data.F32[bin] += 1.0;
   }