Estimating the complex amplitude of a coherent optical wavefield is important in many metrology applications. One way to recover the phase information from measured intensity distributions are phase retrieval (PR) algorithms. Here we develop a numerical technique that accurately simulates the propagation of speckle fields and investigate the convergence properties of a PR algorithm.