Real-time implementation of an iterative solver for atmospheric tomography

by   Bernadett Stadler, et al.

The image quality of the new generation of earthbound Extremely Large Telescopes (ELTs) is heavily influenced by atmospheric turbulences. To compensate these optical distortions a technique called adaptive optics (AO) is used. Many AO systems require the reconstruction of the refractive index fluctuations in the atmosphere, called atmospheric tomography. The standard way of solving this problem is the Matrix Vector Multiplication, i.e., the direct application of a (regularized) generalized inverse of the system operator. However, over the last years the telescope sizes have increased significantly and the computational efficiency become an issue. Promising alternatives are iterative methods such as the Finite Element Wavelet Hybrid Algorithm (FEWHA), which is based on wavelets. Due to its efficient matrix-free representation of the underlying operators, the number of floating point operations and memory usage decreases significantly. In this paper, we focus on performance optimization techniques, such as parallel programming models, for the implementation of this iterative method on CPUs and GPUs. We evaluate the computational performance of our optimized, parallel version of FEWHA for ELT-sized test configurations.



There are no comments yet.


page 1

page 2

page 3

page 4


An augmented wavelet reconstructor for atmospheric tomography

Atmospheric tomography, i.e. the reconstruction of the turbulence profil...

Accelerating linear solvers for Stokes problems with C++ metaprogramming

The efficient solution of large sparse saddle point systems is very impo...

Algorithms and data structures for matrix-free finite element operators with MPI-parallel sparse multi-vectors

Traditional solution approaches for problems in quantum mechanics scale ...

Parallel two-scale finite element implementation of a system with varying microstructures

We propose a two-scale finite element method designed for heterogeneous ...

Learn an index operator by CNN for solving diffusive optical tomography: a deep direct sampling method

In this work, we investigate the diffusive optical tomography (DOT) prob...

Parallelization and scalability analysis of inverse factorization using the Chunks and Tasks programming model

We present three methods for distributed memory parallel inverse factori...

On the crucial impact of the coupling projector-backprojector in iterative tomographic reconstruction

The performance of an iterative reconstruction algorithm for X-ray tomog...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.