Parallel tools for solving incremental dense least squares problems. Application to space geodesy
Marc Baboulin, Luc Giraud, Serge Gratton, and Julien Langou
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  and PBLAS  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  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.
Published 2006-09-19 04:00:00 as ut-cs-06-582 (ID:140)