An augmented wavelet reconstructor for atmospheric tomography

Atmospheric tomography, i.e. the reconstruction of the turbulence profile in the atmosphere, is a challenging task for adaptive optics (AO) systems of the next generation of extremely large telescopes. Within the community of AO the first choice solver is the so called Matrix Vector Multiplication (MVM), which directly applies the (regularized) generalized inverse of the system operator to the data. For small telescopes this approach is feasible, however, for larger systems such as the European Extremely Large Telescope (ELT), the atmospheric tomography problem is considerably more complex and the computational efficiency becomes an issue. Iterative methods, such as the Finite Element Wavelet Hybrid Algorithm (FEWHA), are a promising alternative. FEWHA is a wavelet based reconstructor that uses the well-known iterative preconditioned conjugate gradient (PCG) method as a solver. The number of floating point operations and memory usage are decreased significantly by using a matrix-free representation of the forward operator. A crucial indicator for the real-time performance are the number of PCG iterations. In this paper, we propose an augmented version of FEWHA, where the number of iterations is decreased by $50\%$ using a Krylov subspace recycling technique. We demonstrate that a parallel implementation of augmented FEWHA allows the fulfilment of the real-time requirements of the ELT.

[1]  H. V. D. Vorst,et al.  The rate of convergence of Conjugate Gradients , 1986 .

[2]  Fabien Momey,et al.  Fractal iterative method for fast atmospheric tomography on extremely large telescopes , 2010, Astronomical Telescopes + Instrumentation.

[3]  B. Welsh,et al.  Imaging Through Turbulence , 1996 .

[4]  Ronny Ramlau,et al.  An efficient solution to the atmospheric turbulence tomography problem using Kaczmarz iteration , 2012 .

[5]  John B. Shoven,et al.  I , Edinburgh Medical and Surgical Journal.

[6]  L M Mugnier,et al.  Optimal wave-front reconstruction strategies for multiconjugate adaptive optics. , 2001, Journal of the Optical Society of America. A, Optics, image science, and vision.

[7]  R. Ramlau,et al.  A gradient-based method for atmospheric tomography , 2016 .

[8]  Roberto Biasi,et al.  Real-time implementation of an iterative solver for atmospheric tomography , 2020, ArXiv.

[10]  F. Natterer The Mathematics of Computerized Tomography , 1986 .

[11]  B. Ellerbroek,et al.  Split atmospheric tomography using laser and natural guide stars. , 2008, Journal of the Optical Society of America. A, Optics, image science, and vision.

[12]  T. Fusco,et al.  Coupling MOAO with integral field spectroscopy: specifications for the VLT and the E-ELT , 2008, 0808.1196.

[13]  Eric de Sturler,et al.  A survey of subspace recycling iterative methods , 2020, ArXiv.

[14]  Luc Gilles,et al.  Layer-oriented multigrid wavefront reconstruction algorithms for multiconjugate adaptive optics , 2003, SPIE Astronomical Telescopes + Instrumentation.

[15]  E. Diolaiti,et al.  MAORY: A Multi-conjugate Adaptive Optics RelaY for the E-ELT , 2010 .

[16]  Curtis R Vogel,et al.  Fourier domain preconditioned conjugate gradient algorithm for atmospheric tomography. , 2006, Applied optics.

[17]  Frédéric Guyomarc'h,et al.  An Augmented Conjugate Gradient Method for Solving Consecutive Symmetric Positive Definite Linear Systems , 2000, SIAM J. Matrix Anal. Appl..

[18]  Ruth C. Carter,et al.  Principles , 2003, Law’s Reality.

[19]  R. Ramlau,et al.  Efficient iterative tip/tilt reconstruction for atmospheric tomography , 2014 .

[20]  Francois Rigaut,et al.  Principles, limitations, and performance of multiconjugate adaptive optics , 2000, Astronomical Telescopes and Instrumentation.

[21]  R. Ramlau,et al.  A FREQUENCY DEPENDENT PRECONDITIONED WAVELET METHOD FOR ATMOSPHERIC TOMOGRAPHY , 2013 .

[22]  Y. Saad,et al.  On the Lánczos method for solving symmetric linear systems with several right-hand sides , 1987 .

[23]  M. E. Davison,et al.  The Ill-Conditioned Nature of the Limited Angle Tomography Problem , 1983 .

[24]  M. Kasper,et al.  Adaptive Optics for Astronomy , 2012, 1201.5741.

[25]  Tapio Helin,et al.  Wavelet methods in multi-conjugate adaptive optics , 2013, 1302.3734.

[26]  R. Ramlau,et al.  Optimal mirror deformation for multi conjugate adaptive optics systems , 2016 .

[27]  F. Courbin,et al.  The FALCON Concept: Multi-Object Spectroscopy Combined with MCAO in Near-IR , 2002 .

[28]  B. L. Ellerbroek,et al.  Inverse problems in astronomical adaptive optics , 2009 .

[29]  Ronny Ramlau,et al.  Finite element-wavelet hybrid algorithm for atmospheric tomography. , 2014, Journal of the Optical Society of America. A, Optics, image science, and vision.

[30]  Ronny Ramlau,et al.  A new temporal control approach for SCAO systems , 2019 .

[32]  Sylvain Oberti,et al.  Adaptive optics simulations for the European Extremely Large Telescope , 2006, Astronomical Telescopes + Instrumentation.

[33]  Michel Tallon,et al.  Fast minimum variance wavefront reconstruction for extremely large telescopes. , 2010, Journal of the Optical Society of America. A, Optics, image science, and vision.

[34]  Luc Gilles,et al.  A comparison of Multigrid V-cycle versus Fourier Domain Preconditioning for Laser Guide Star Atmospheric Tomography , 2007 .

[35]  Thierry Fusco,et al.  Tomographic wavefront error using multi-LGS constellation sensed with Shack-Hartmann wavefront sensors. , 2010, Journal of the Optical Society of America. A, Optics, image science, and vision.

[36]  Luc Gilles,et al.  Computationally efficient wavefront reconstructor for simulation of multiconjugate adaptive optics on giant telescopes , 2003, SPIE Astronomical Telescopes + Instrumentation.

[37]  Curtis R Vogel,et al.  Preconditioned conjugate gradient wave-front reconstructors for multiconjugate adaptive optics. , 2003, Applied optics.