Pipelining Computational Stages of the Tomographic Reconstructor for Multi-Object Adaptive Optics on a Multi-GPU System

The European Extremely Large Telescope project (E-ELT) is one of Europe's highest priorities in ground-based astronomy. ELTs are built on top of a variety of highly sensitive and critical astronomical instruments. In particular, a new instrument called MOSAIC has been proposed to perform multi-object spectroscopy using the Multi-Object Adaptive Optics (MOAO) technique. The core implementation of the simulation lies in the intensive computation of a tomographic reconstruct or (TR), which is used to drive the deformable mirror in real time from the measurements. A new numerical algorithm is proposed (1) to capture the actual experimental noise and (2) to substantially speed up previous implementations by exposing more concurrency, while reducing the number of floating-point operations. Based on the Matrices Over Runtime System at Exascale numerical library (MORSE), a dynamic scheduler drives all computational stages of the tomographic reconstruct or simulation and allows to pipeline and to run tasks out-of order across different stages on heterogeneous systems, while ensuring data coherency and dependencies. The proposed TR simulation outperforms asymptotically previous state-of-the-art implementations up to 13-fold speedup. At more than 50000 unknowns, this appears to be the largest-scale AO problem submitted to computation, to date, and opens new research directions for extreme scale AO simulations.

[1]  Emmanuel Agullo,et al.  Towards an Efficient Tile Matrix Inversion of Symmetric Positive Definite Matrices on Multicore Architectures , 2010, VECPAR.

[2]  J. Cuby,et al.  ELT-MOS White Paper: Science Overview & Requirements , 2013, 1303.0029.

[3]  David E. Keyes,et al.  Toward Accelerating the Matrix Inversion Computation of Symmetric Positive-Definite Matrices on Heterogeneous GPU-Based Systems , 2011 .

[4]  Julien Langou,et al.  A Class of Parallel Tiled Linear Algebra Algorithms for Multicore Architectures , 2007, Parallel Comput..

[5]  Philippe Dierickx,et al.  E-ELT update of project and effect of change to 39m design , 2012, Other Conferences.

[6]  Robert A. van de Geijn,et al.  Supermatrix out-of-order scheduling of matrix operations for SMP and multi-core architectures , 2007, SPAA '07.

[7]  Alejandro Duran,et al.  Ompss: a Proposal for Programming Heterogeneous Multi-Core Architectures , 2011, Parallel Process. Lett..

[8]  John Shalf,et al.  The International Exascale Software Project roadmap , 2011, Int. J. High Perform. Comput. Appl..

[9]  G. Rousset,et al.  Tomography approach for multi-object adaptive optics. , 2010, Journal of the Optical Society of America. A, Optics, image science, and vision.

[10]  Mathieu Puech,et al.  Project overview of OPTIMOS-EVE: the fibre-fed multi-object spectrograph for the E-ELT , 2010, Astronomical Telescopes + Instrumentation.

[11]  A. Sevin,et al.  GPUs for adaptive optics: simulations and real-time control , 2012, Other Conferences.

[12]  A. Sevin,et al.  MOAO first on-sky demonstration with CANARY , 2011 .

[13]  J. Conan,et al.  Tomographic reconstruction for wide-field adaptive optics systems: Fourier domain analysis and fundamental limitations. , 2008, Journal of the Optical Society of America. A, Optics, image science, and vision.

[14]  Robert A. van de Geijn,et al.  The libflame Library for Dense Matrix Computations , 2009, Computing in Science & Engineering.

[15]  W. R. Sutherland,et al.  The on-line graphical specification of computer procedures , 1966 .

[16]  Francois Rigaut,et al.  Laser guide star in adaptive optics : the tilt determination problem , 1992 .

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

[18]  A. G. Basden,et al.  Monte-Carlo simulation of ELT scale multi-object adaptive optics deformable mirror requirements and tolerances , 2013 .

[19]  Jack Dongarra,et al.  QUARK Users' Guide: QUeueing And Runtime for Kernels , 2011 .

[20]  B. Ellerbroek,et al.  Sky coverage modeling for the whole sky for laser guide star multiconjugate adaptive optics. , 2012, Applied optics.

[21]  Miska Le Louarn Progress and prospects in AO simulation capabilities , 2010 .

[22]  Jack Dongarra,et al.  ScaLAPACK Users' Guide , 1987 .

[23]  T. Fusco,et al.  PSD-based simulation algorithm for Wide FoV AO design: application to ELT studies , 2008, Astronomical Telescopes + Instrumentation.

[24]  Eric Gendron,et al.  Astronomical adaptive optics. I. Modal control optimization. , 1994 .

[25]  Cédric Augonnet,et al.  StarPU: a unified platform for task scheduling on heterogeneous multicore architectures , 2011, Concurr. Comput. Pract. Exp..

[26]  David E. Keyes,et al.  High Performance Pseudo-analytical Simulation of Multi-Object Adaptive Optics over Multi-GPU Systems , 2014, Euro-Par.

[27]  Fabrice Madec,et al.  EAGLE: a MOAO fed multi-IFU NIR workhorse for E-ELT , 2010, Astronomical Telescopes + Instrumentation.

[28]  T. Fusco,et al.  EAGLE MOAO system conceptual design and related technologies , 2010, Astronomical Telescopes + Instrumentation.

[29]  Thomas Hérault,et al.  Flexible Development of Dense Linear Algebra Algorithms on Massively Parallel Architectures with DPLASMA , 2011, 2011 IEEE International Symposium on Parallel and Distributed Processing Workshops and Phd Forum.

[30]  Emmanuel Agullo,et al.  Comparative study of one-sided factorizations with multiple software packages on multi-core hardware , 2009, Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis.