Parallel Tools for Solving Incremental Dense Least Squares Problems: Application to Space Geodesy

We present a parallel distributed solver that enables us to solve incremental dense least squares arising in some parameter estimation problems. This solver is based on ScaLAPACK [8] and PBLAS [9] kernel routines. In the incremental process, the observations are collected periodically and the solver...

Full description

Bibliographic Details
Main Authors: Arc Baboulin, Luc Giraud, Serge Gratton, Julien Langou
Format: Article
Language:English
Published: SAGE Publishing 2009-03-01
Series:Journal of Algorithms & Computational Technology
Online Access:https://doi.org/10.1260/174830109787186541
Description
Summary:We present a parallel distributed solver that enables us to solve incremental dense least squares arising in some parameter estimation problems. This solver is based on ScaLAPACK [8] and PBLAS [9] kernel routines. In the incremental process, the observations are collected periodically and the solver updates the solution with new observations using a QR factorization algorithm. It uses a recently defined distributed packed format [3] that handles symmetric or triangular matrices in ScaLAPACK-based implementations. We provide performance analysis on IBM pSeries 690. We also present an example of application in the area of space geodesy for gravity field computations with some experimental results.
ISSN:1748-3018
1748-3026