I recently ran into this problem and came up with a solution using sparse matrices to compute the toepliz. I posted it at the torch forums. I hope it might help.