--- trunk/SHAPES/visualizer.cpp 2004/06/24 15:56:05 1299 +++ trunk/SHAPES/visualizer.cpp 2004/06/25 18:33:17 1307 @@ -61,15 +61,22 @@ int main(int argc, char* argv[]){ // grid has a default value (default=51), so it is always given npts = args_info.grid_arg; - xmin = -5.0; - xmax = 5.0; + xmin = -10.0; + xmax = 10.0; - ymin = -5.0; - ymax = 5.0; + ymin = -10.0; + ymax = 10.0; - zmin = -5.0; - zmax = 5.0; + zmin = -10.0; + zmax = 10.0; + //sigmaProbe = 2.28; + //sProbe = 2.28; + //epsProbe = 0.020269601874; + sigmaProbe = 0.0; + sProbe = 0.0; + epsProbe = 1.0; + for (i = 0; i < npts; i++) { x = xmin + (xmax-xmin) * (double)i/(double)npts; @@ -80,8 +87,8 @@ int main(int argc, char* argv[]){ z = zmin + (zmax-zmin) * (double)k/(double)npts; r = sqrt(x*x + y*y + z*z); - theta = acos(z/r); - phi = atan(y/x); + costheta = z/r; + phi = atan2(y,x); sigmaShape = shape->getSigmaAt(costheta, phi); sShape = shape->getSAt(costheta, phi);