1 double sqrt(double),fabs(double),sin(double);
9 for (l = 0; l < 9; l++)
12 int cx, cy, gx, gy, x, y;
13 gx = 2 > g / 3 ? 2 : g / 3;
14 gy = 2 > g / 3 ? 2 : g / 3;
15 for (y = 0 > cy - gy ? 0 : cy - gy; y <= (sys - 1 < cy + gy ? sys : cy + gy); y++)
17 int sx = 0 > cx - gx ? 0 : cx - gx;
18 short *ax = (short *) (y * sxs + sx);
20 for (x = sx; x <= (sxs - 1 < cx + gx ? sxs - 1 : cx + gx); x++)
22 double c=2.25, z=sqrt(fabs(1-c)), cz=(c>1?0.0:-10)>z?c>1?0:1:z;
26 for (l = sm; l <= sx; l++)
27 smap[l] = l > 0 ? 1 + foo(sin(.1 * l / sx)) : 1 - foo(sin(.1 * l / sm));