Porting a sphere optimization program from LAPACK to ScaLAPACK

Abstract

The sphere optimization program sphopt was originally written as a sequential program using lapack, and was converted to use scalapack, primarily to overcome memory limitations. The conversion was relatively straightforward, using a small number of organizing principles which are widely applicable to the scalapack parallelization of serial code. The main innovation is the use of a compressed block cyclic storage scheme to store two symmetric matrices in little more than the storage needed for one matrix. The resulting psphopt program scales at least as well as the scalapack Cholesky decomposition routine which it uses. Show more