Collision catalyst handles stars with NAN coords/velocities

collide.c

@@ -54,8 +54,8 @@ point_t calculateCoM(RigidMass_t *stars, int nStars)
 	double totalY = 0;
 	for(int i = 0; i < nStars; i++)
 	{
-		totalX += stars[i].x;
-		totalY += stars[i].y;
+		if(!isnan(stars[i].x)) totalX += stars[i].x;
+		if(!isnan(stars[i].y)) totalY += stars[i].y;
 	}
 	point_t out;
 	out.x = totalX / nStars;
@@ -71,8 +71,8 @@ point_t calculateCoV(RigidMass_t *stars, int nStars)
 	double totalVy = 0;
 	for(int i = 0; i < nStars; i++)
 	{
-		totalVx += stars[i].vx;
-		totalVy += stars[i].vy;
+		if(!isnan(stars[i].vx)) totalVx += stars[i].vx;
+		if(!isnan(stars[i].vy)) totalVy += stars[i].vy;
 	}
 	point_t out;
 	out.x = totalVx / nStars;