A simple Parallel Program that uses ScaLAPACK package.

This program solves a linear system by calling ScaLAPACK routine PSGESV, for dense matrices. The matrix A is defined as A(I,J) = min(I,J). In this example, we solve linear system NxN ( N= 10) on processors grid NPROWxNPCOL (2x2), i.e. local array on each processor will be 5x5.

The program uses the following routines:

from BLACS library

  • BLACS_PINFO
  • BLACS_SETUP
  • BLACS_GET
  • BLACS_GRIDINIT
  • BLACS_GRIDINFO
  • BLACS_GRIDEXIT
  • BLACS_EXIT

    • from TOOLS library

  • DESCINIT
  • PSLAPRNT

    • from SCALAPACK library

  • PSGESV