pixlist = pt.query_disc(nside, vc, radius, nest, inclusive);
int nlist = pixlist.size();
for (int i = 0; i < nlist; i++) {
ipix = ((Long) pixlist.get(i)).longValue();
Vector3d v = pt.pix2vect_ring(nside,ipix);
double dist = pt.AngDist(v,vc);
assertTrue(dist<=2.*radius);
}
cpix = pt.ang2pix_nest(nside,theta,phi);
Vector3d vc1 = pt.pix2vect_nest(nside, cpix);