# | Line 61 | Line 61 | int main(int argc, char* argv[]){ | |
---|---|---|
61 | // grid has a default value (default=51), so it is always given | |
62 | npts = args_info.grid_arg; | |
63 | ||
64 | < | xmin = -5.0; |
65 | < | xmax = 5.0; |
64 | > | xmin = -10.0; |
65 | > | xmax = 10.0; |
66 | ||
67 | < | ymin = -5.0; |
68 | < | ymax = 5.0; |
67 | > | ymin = -10.0; |
68 | > | ymax = 10.0; |
69 | ||
70 | < | zmin = -5.0; |
71 | < | zmax = 5.0; |
70 | > | zmin = -10.0; |
71 | > | zmax = 10.0; |
72 | ||
73 | + | //sigmaProbe = 2.28; |
74 | + | //sProbe = 2.28; |
75 | + | //epsProbe = 0.020269601874; |
76 | + | sigmaProbe = 0.0; |
77 | + | sProbe = 0.0; |
78 | + | epsProbe = 1.0; |
79 | + | |
80 | for (i = 0; i < npts; i++) { | |
81 | x = xmin + (xmax-xmin) * (double)i/(double)npts; | |
82 | ||
# | Line 80 | Line 87 | int main(int argc, char* argv[]){ | |
87 | z = zmin + (zmax-zmin) * (double)k/(double)npts; | |
88 | ||
89 | r = sqrt(x*x + y*y + z*z); | |
90 | < | theta = acos(z/r); |
91 | < | phi = atan(y/x); |
90 | > | costheta = z/r; |
91 | > | phi = atan2(y,x); |
92 | ||
93 | sigmaShape = shape->getSigmaAt(costheta, phi); | |
94 | sShape = shape->getSAt(costheta, phi); |
– | Removed lines |
+ | Added lines |
< | Changed lines |
> | Changed lines |