| 1 | for (int c0 = 1; c0 <= floord(Ny, 2) + 2; c0 += 1) |
|---|---|
| 2 | for (int c1 = max(c0 - 1, c0 / 2 + 1); c1 <= min(c0, (Ny + 2 * c0) / 4); c1 += 1) { |
| 3 | if (Ny + 2 * c0 >= 4 * c1 + 1) { |
| 4 | for (int c2 = 1; c2 <= 2; c2 += 1) { |
| 5 | S1(c0 - c1, c1, 2 * c0 - 2 * c1, -2 * c0 + 4 * c1, c2); |
| 6 | S2(c0 - c1, c1, 2 * c0 - 2 * c1, -2 * c0 + 4 * c1 - 1, c2); |
| 7 | } |
| 8 | } else { |
| 9 | for (int c2 = 1; c2 <= 2; c2 += 1) |
| 10 | S2((-Ny + 2 * c0) / 4, (Ny + 2 * c0) / 4, (-Ny / 2) + c0, Ny - 1, c2); |
| 11 | } |
| 12 | } |
| 13 |