Some code on circular convolution