/* Note: THIS ROUTINE ASSUMES PERIODICITY!!! */
/* It also assumes that ALL points are included in the process. */
int process(                     /* Returns 0 if coalescence has */
	                         /* occured and 1 otherwise */
	    double time,         /* Value of $-M$ */
	    long length,         /* Number of data points */
	    double gamma,        /* Value of $\gamma$ */
	    double sigma,        /* Value of $\sigma$ */
	    double tau,          /* Value of $\tau$ */
	    double exponent,     /* The power we raise $J_{jk}$ to. */
	                         /* Determines tails of the prior */
	    location points[]    /* All info about all the points */
	    );
