The browser you are using is not supported by this website. All versions of Internet Explorer are no longer supported, either by us or Microsoft (read more here: https://www.microsoft.com/en-us/microsoft-365/windows/end-of-ie-support).

Please use a modern browser to fully experience our website, such as the newest versions of Edge, Chrome, Firefox or Safari etc.

Lennart Lindegren. Profile picture.

Lennart Lindegren

Professor

Lennart Lindegren. Profile picture.

Complexity of the Gaia astrometric least-squares problem and the (non-)feasibility of a direct solution method

Author

  • A. Bombrun
  • Lennart Lindegren
  • Berry Holl
  • S. Jordan

Summary, in English

The Gaia space astrometry mission (to be launched in 2012) will use a continuously spinning spacecraft to construct a global system of positions, proper motions and absolute parallaxes from relative position measurements made in an astrometric focal plane. This astrometric reduction can be cast as a classical least-squares problem, and the adopted baseline method for its solution uses a simple iteration algorithm. A potential weakness of this approach, as opposed to a direct solution, is that any finite number of iterations results in truncation errors that are difficult to quantify. Thus it is of interest to investigate alternative approaches, in particular the feasibility of a direct (non-iterative) solution. A simplified version of the astrometric reduction problem is studied in which the only unknowns are the astrometric parameters for a subset of the stars and the continuous three-axis attitude, thus neglecting further calibration issues. The specific design of the Gaia spacecraft and scanning law leads to an extremely large and sparse normal equations matrix. Elimination of the star parameters leads to a much smaller but less sparse system. We try different reordering schemes and perform symbolic Cholesky decomposition of this reduced normal matrix to study the fill-in for successively longer time span of simulated observations. Extrapolating to the full mission length, we conclude that a direct solution is not feasible with today's computational capabilities. Other schemes, e. g., eliminating the attitude parameters or orthogonalizing the observation equations, lead to similar or even worse problems. This negative result appears to be a consequence of the strong spatial and temporal connectivity among the unknowns achieved by two superposed fields of view and the scanning law, features that are in fact desirable and essential for minimizing large-scale systematic errors in the Gaia reference frame. We briefly consider also an approximate decomposition method a la Hipparcos, but conclude that it is either sub-optimal or effectively leads to an iterative solution.

Department/s

  • Lund Observatory - Undergoing reorganization

Publishing year

2010

Language

English

Publication/Series

Astronomy & Astrophysics

Volume

516

Document type

Journal article

Publisher

EDP Sciences

Topic

  • Astronomy, Astrophysics and Cosmology

Keywords

  • methods: numerical
  • methods: data analysis
  • space vehicles:
  • astrometry
  • instruments
  • reference systems

Status

Published

Research group

  • Observational and Theoretical Astrophysics

ISBN/ISSN/Other

  • ISSN: 0004-6361