Implementation of k-space PSTD method.