Matrices

The routines for matrix operations in this library are:

• PDA_DGEFA (SLATEC/CAMSUN)
  Factor a matrix using Gaussian elimination. This is needed before the determinant and inverse can be calculated by PDA_DGEDI.
• PDA_DGEDI (SLATEC/CAMSUN)
  Compute the determinant and inverse of a matrix using the factors computed by PDA_DGECO or PDA_DGEFA.
• PDA_DGEFS (SLATEC/CAMSUN)
  Solve a general system of linear equations. This solves the problem A * x = b. A is a square matrix, x and b are vectors. The problem A * X = B where all are matrices can be solved as several systems A * x = b. This routine supports such an undertaking since it is able to re-use a previous factorisation of A.
• PDA_DBOLS (SLATEC/CAMSUN)
  Solve the problem E * x = f (in the least squares sense) with bounds on selected x values. E is a matrix, x and f are vectors.
• PDA_LSQR
  Solves sparse unsymmetric, linear least squares and damped least squares problems