1 for (int c0
= 1; c0
<= floord(Ny
, 2) + 2; c0
+= 1)
2 for (int c1
= max(c0
- (c0
+ 1) / 2 + 1, c0
- 1); c1
<= min((Ny
+ 2 * c0
) / 4, c0
); 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
);
9 for (int c2
= 1; c2
<= 2; c2
+= 1)
10 S2((-Ny
+ 2 * c0
) / 4, (Ny
+ 2 * c0
) / 4, (-Ny
+ 2 * c0
) / 2, Ny
- 1, c2
);