
double t_noedge(               /* Calculates T(h) without edge correction */
		coord *points,       /* An array of positions */
		int n,               /* Number of points */
		double a,            /* Length of box in x-direction */
		double b,            /* Length of box in y-direction */
		double h);           /* Parameter of t */

double t_tor(              /* Calculates T(h) on a torus */
	     coord *points,       /* An array of positions */
	     int n,               /* Number of points */
	     double a,            /* Length of box in x-direction */
	     double b,            /* Length of box in y-direction */
	     double h);           /* Parameter of t */
