xi[0]=sin(phi)*cos(theta);
xi[1]=sin(phi)*sin(theta);
xi[2]=cos(phi);
line.p=p;
line.xi=xi;
f=line.linmin();
d=sqrt(SQR(line.p[0])+SQR(line.p[1])+SQR(line.p[2]));
// System.out.printf(abs(d-d0) << " " << abs(f-f0));
localflag = localflag || abs(d-d0) > sbeps1;
globalflag = globalflag || localflag;