added code for timing convolution; one small example comparing
abelian and nonabelian convolution seems to be working.