Well, you can extend the trick above, if you have c channels, c divisible by k*l:
n,c,h,w = input.size()
res = input.reshape(n, c//k//l, k, l, l, h, w).sum(2).view(n, c//k, h, w)

This will sum k channels from an interleaving grid of spacing l.
Best regards
Thomas