Add 3d correlation based on blas matrix multiplication 

padD = 0?