for (int c1 = 0; c1 < M; c1 += 1)

  for (int c3 = 0; c3 < M; c3 += 1)

    S1(c1, c3);