# | Line 364 | Line 364 | int x, y; | |
---|---|---|
364 | xo = x-hl; | |
365 | else | |
366 | xo = 0; | |
367 | < | disc = 1. - (double)(xo*xo + y*y)/(sampdens*sampdens); |
368 | < | if (disc <= FTINY) |
367 | > | disc = 1. - (double)((long)xo*xo + (long)y*y)/((long)sampdens*sampdens); |
368 | > | if (disc <= FTINY*FTINY) |
369 | return(0.); | |
370 | return(1./(sampdens*sampdens*sqrt(disc))); | |
371 | } |
– | Removed lines |
+ | Added lines |
< | Changed lines |
> | Changed lines |