1*5971e316Smrg for (int c0 = 2; c0 < O; c0 += 1) 2*5971e316Smrg for (int c1 = 3; c1 < 2 * N - 2; c1 += 2) { 3*5971e316Smrg for (int c3 = 1; c3 <= M; c3 += 1) { 4*5971e316Smrg S1(c0, (c1 + 1) / 2, c3); 5*5971e316Smrg S2(c0, (c1 + 1) / 2, c3); 6*5971e316Smrg } 7*5971e316Smrg for (int c3 = 2; c3 < M; c3 += 1) 8*5971e316Smrg S3(c0, (c1 + 1) / 2, c3); 9*5971e316Smrg } 10