Abstract

We present what we believe to be a new algorithm, FRactal Iterative Method (FRiM), aiming at the reconstruction of the optical wavefront from measurements provided by a wavefront sensor. As our application is adaptive optics on extremely large telescopes, our algorithm was designed with speed and best quality in mind. The latter is achieved thanks to a regularization that enforces prior statistics. To solve the regularized problem, we use the conjugate gradient method, which takes advantage of the sparsity of the wavefront sensor model matrix and avoids the storage and inversion of a huge matrix. The prior covariance matrix is, however, non-sparse, and we derive a fractal approximation to the Karhunen–Loève basis thanks to which the regularization by Kolmogorov statistics can be computed in $\mathcal{O}\left(N\right)$ operations, with N being the number of phase samples to estimate. Finally, we propose an effective preconditioning that also scales as $\mathcal{O}\left(N\right)$ and yields the solution in five to ten conjugate gradient iterations for any N. The resulting algorithm is therefore $\mathcal{O}\left(N\right)$. As an example, for a $128\times 128$ Shack–Hartmann wavefront sensor, the FRiM appears to be more than 100 times faster than the classical vector-matrix multiplication method.

© 2010 Optical Society of America

Performance comparison of wavefront reconstruction and control algorithms for Extremely Large Telescopes

I. Montilla, C. Béchet, M. Le Louarn, M. Reyes, and M. Tallon
J. Opt. Soc. Am. A 27(11) A9-A18 (2010)

Comparison of minimum-norm maximum likelihood and maximum a posteriori wavefront reconstructions for large adaptive optics systems

Clémentine Béchet, Michel Tallon, and Éric Thiébaut
J. Opt. Soc. Am. A 26(3) 497-508 (2009)

Cumulative Reconstructor: fast wavefront reconstruction algorithm for Extremely Large Telescopes

Matthias Rosensteiner
J. Opt. Soc. Am. A 28(10) 2132-2138 (2011)

Computational performance comparison of wavefront reconstruction algorithms for the European Extremely Large Telescope on multi-CPU architecture

Lu Feng, Enrico Fedrigo, Clémentine Béchet, Elisabeth Brunner, and Werther Pirani
Appl. Opt. 51(16) 3564-3583 (2012)

Finite element-wavelet hybrid algorithm for atmospheric tomography

Mykhaylo Yudytskiy, Tapio Helin, and Ronny Ramlau
J. Opt. Soc. Am. A 31(3) 550-560 (2014)