24 	con(con_), boxx(con_.boxx), boxy(con_.boxy), boxz(con_.boxz),  in voro_compute()
27 	id(con_.id), p(con_.p), co(con_.co), bxsq(boxx*boxx+boxy*boxy+boxz*boxz),  in voro_compute()
100 		mxs=boxx-fx;  in find_voronoi_cell()
347 	} else {m1=m2=0;gxs=boxx-fx;}  in compute_cell()
543 		xlo=(ei-i)*boxx-fx;xhi=xlo+boxx;  in compute_cell()
801 		xlo=di*boxx-fx;  in compute_min_max_radius()
809 				crs+=bxsq+2*(boxx*xlo+boxy*ylo+boxz*zlo);  in compute_min_max_radius()
813 				crs+=bxsq+2*(boxx*xlo+boxy*ylo-boxz*zlo);  in compute_min_max_radius()
816 				crs+=boxx*(2*xlo+boxx)+boxy*(2*ylo+boxy)+gzs;  in compute_min_max_radius()
824 				crs+=bxsq+2*(boxx*xlo-boxy*ylo+boxz*zlo);  in compute_min_max_radius()
828 				crs+=bxsq+2*(boxx*xlo-boxy*ylo-boxz*zlo);  in compute_min_max_radius()
831 				crs+=boxx*(2*xlo+boxx)+boxy*(-2*ylo+boxy)+gzs;  in compute_min_max_radius()
846 			crs+=gys+boxx*(2*xlo+boxx);  in compute_min_max_radius()
849 		xlo=(di+1)*boxx-fx;  in compute_min_max_radius()
857 				crs+=bxsq+2*(-boxx*xlo+boxy*ylo+boxz*zlo);  in compute_min_max_radius()
861 				crs+=bxsq+2*(-boxx*xlo+boxy*ylo-boxz*zlo);  in compute_min_max_radius()
864 				crs+=boxx*(-2*xlo+boxx)+boxy*(2*ylo+boxy)+gzs;  in compute_min_max_radius()
872 				crs+=bxsq+2*(-boxx*xlo-boxy*ylo+boxz*zlo);  in compute_min_max_radius()
876 				crs+=bxsq+2*(-boxx*xlo-boxy*ylo-boxz*zlo);  in compute_min_max_radius()
879 				crs+=boxx*(-2*xlo+boxx)+boxy*(-2*ylo+boxy)+gzs;  in compute_min_max_radius()
894 			crs+=gys+boxx*(-2*xlo+boxx);  in compute_min_max_radius()
951 	if(di>0) {t=di*boxx-fx;crs=t*t;}  in compute_min_radius()
952 	else if(di<0) {t=(di+1)*boxx-fx;crs=t*t;}  in compute_min_radius()