Version 1.0.0

Code Download

  • Single
  • Double

MC26 Sparse rectangular matrix: compute normal matrix

Given a sparse matrix A, this subroutine computes the sparse matrix ATA. Three forms of data storage are permitted for the input matrix: storage by columns, where the row indices and column pointers describe the matrix; storage by rows, where the column indices and row pointers describe the matrix; and the coordinate scheme, where both row and column indices describe the position of entries in the matrix.