A new, to my knowledge, procedure for retrieving the wave aberration from the point-spread function is presented. It uses the Levenberg-Marquardt optimization algorithm in a mutiresolution pyramidal scheme. The method, tested with simulated large aberrations without initial estimates, accelerates convergence and avoids stagnation in local minima.