Also implemented routine for calculating sum(a[i] * a[i-offset]>>shift)

Originally committed as revision 12974 to svn://svn.ffmpeg.org/ffmpeg/trunk