sbeps=1.e-15;
localflag=false;
for (i=0;i<M;i++) {
n1=myran.int32p() % N;
n2=myran.int32p() % N;
double dis=tree3.disti(n1,n2);
// System.out.printf(dis << " %f\n", sqrt(SQR(kd3[n1].x[0]-kd3[n2].x[0])
// + SQR(kd3[n1].x[1]-kd3[n2].x[1])
// + SQR(kd3[n1].x[2]-kd3[n2].x[2])));
if (n1 == n2)
localflag = localflag || dis < 1.e99;