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.
AndersonE.BaiZ.BischofC.BlackfordS.DemmelJ.DongarraJ.Du CrozJ.GreenbaumA.HammarlingS.McKenneyA., and SorensenD., LAPACK user's guide, SIAM, 1999, Third edition.
2.
BaboulinM.GiraudL., and GrattonS., A parallel distributed solver for large dense symmetric systems: Applications to geodesy and electromagnetism problems, Int. J. of High Performance Computing Applications19 (2005), no. 4, 353–363.
3.
BaboulinM.GiraudL.GrattonS., and LangouJ., A distributed packed storage for large dense parallel in-core calculations, Technical Report TR/PA/05/30, CERFACS, Toulouse, France, 2005, To appear in Concurrency and Computation: Practice and Experience.
4.
BalminoG., The European GOCE Gravity Consortium (EGG-C), (April 2001), 7–12, Proceedings of the International GOCE User Workshop.
5.
BalminoG.BruinsmaS., and MartyJ-C., Numerical simulation of the gravity field recovery from GOCE mission data, (March 8–10, 2004), Proceedings of the Second International GOCE User Workshop “GOCE, The Geoid and Oceanography”, ESAESRIN, Frascati, Italy.
6.
BalminoG.CazenaveA.Comolet-TirmanA.HussonJ. C., and LefebvreM., Cours de gédésie dynamique et spatiale, ENSTA, 1982.
7.
BjörckA., Numerical methods for least squares problems, SIAM, 1996.
8.
BlackfordL.ChoiJ.ClearyA.D'AzevedoE.DemmelJ.DhillonI.DongarraJ.HammarlingS.HenryG.PetitetA.StanleyK.WalkerD., and WhaleyR., ScaLAPACK user's guide, SIAM, 1997.
9.
ChoiJ.DongarraJ.OstrouchovL.PetitetA.WalkerD., and WhaleyR., A proposal for a set of parallel basic linear algebra subprograms, Tech. report, 1995, LAPACK Working Note 100.
10.
ChoiJ.DongarraJ.OstrouchovL.PetitetA.WalkerD., and WhaleyR., The design and implementation of the ScaLAPACK LU, QR, and Cholesky factorization routines, Scientific Programming5 (1996), 173–184.
11.
DemmelJ. and DongarraJ., LAPACK 2005 prospectus: Reliable and scalable software for linear algebra computations on high end computers, (February 2005), LAPACK Working Note 164.
12.
GolubG. H.MannebackP., and TointPh. L., A comparison between some direct and iterative methods for certain large scale geodetic least squares problems, SIAM J. Scientific Computing7 (1986), no. 3, 799–816.
13.
GolubG. H. and van LoanC. F., Matrix computations, The Johns Hopkins University Press, 1996, Third edition.
14.
GunterB.ReileyW., and van de GeijnR., Parallel out-of-core Cholesky and QR factorizations with POOCLAPACK, IEEE Computer Society (2001), Proceedings of the 15th International Parallel and Distributed Processing Symposium (IPDPS).
15.
GunterB. and van de GeijnR., Parallel out-of-core computation and updating of the QR factorization, ACM Trans. Math. Softw.31 (2005), no. 1, 60–78.
16.
HighamN. J., Accuracy and stability of numerical algorithms, SIAM, 2002, Second edition.
17.
DennisJ. E.Jr. and SchnabelR. B., Numerical methods for unconstrained optimization and nonlinear equations, SIAM, 1996.
18.
KailathT.SayedA. H., and HassibiB., Linear estimation, Prentice-Hall, 2000.
19.
SünkelH., From Eötvös to milligal+, Final Report, ESA/ESTEC Contract No. 13392/98/NL/GD, Graz University of Technology, 2000.
20.
van de GeijnR., Using PLAPACK, The MIT Press, 1997.
