D. Linear solve (Sparse matrices)

Functions
Sub	CsrDusmm (Trans As String, M As Long, N As Long, Nrhs As Long, Alpha As Double, Val() As Double, Rowptr() As Long, Colind() As Long, B() As Double, Beta As Double, C() As Double, Optional Info As Long, Optional Base As Long=-1)
	C <- αAB + βC or C <- αATB + βC (CSR)
Sub	CsrDusmv (Trans As String, M As Long, N As Long, Alpha As Double, Val() As Double, Rowptr() As Long, Colind() As Long, X() As Double, Beta As Double, Y() As Double, Optional Info As Long, Optional Base As Long=-1, Optional IncX As Long=1, Optional IncY As Long=1)
	y <- αAx + βy or y <- αATx + βy (CSR)
Sub	CsrDussm (Uplo As String, Trans As String, Diag As String, N As Long, Nrhs As Long, Val() As Double, Rowptr() As Long, Colind() As Long, X() As Double, Optional Info As Long, Optional Base As Long=-1)
	Solution of AX = B or ATX = B (Triangular matrices) (CSR)
Sub	CsrDussv (Uplo As String, Trans As String, Diag As String, N As Long, Val() As Double, Rowptr() As Long, Colind() As Long, X() As Double, Optional Info As Long, Optional Base As Long=-1, Optional IncX As Long=1, Optional Omega As Double=1)
	Solution of Ax = b or ATx = b (triangular matrices) (CSR)
Sub	SsrDusmv (Uplo As String, N As Long, Alpha As Double, Val() As Double, Rowptr() As Long, Colind() As Long, X() As Double, Beta As Double, Y() As Double, Optional Info As Long, Optional Base As Long=-1, Optional IncX As Long=1, Optional IncY As Long=1)
	y <- αAx + βy (CSR) (Symmetric matrix)
Sub	CooCsr (M As Long, N As Long, Nnz As Long, Val() As Double, Rowind() As Long, Colind() As Long, Val2() As Double, Rowptr2() As Long, Colind2() As Long, Optional Info As Long, Optional Base As Long=0, Optional Base2 As Long=0)
	COO -> CSR
Sub	CsrCoo (M As Long, N As Long, Val() As Double, Rowptr() As Long, Colind() As Long, Val2() As Double, Rowind2() As Long, Colind2() As Long, Optional Info As Long, Optional Nnz As Long, Optional Base As Long=-1, Optional Base2 As Long=0)
	CSR -> COO
Sub	CsrDense (M As Long, N As Long, Val() As Double, Rowptr() As Long, Colind() As Long, A() As Double, Optional Info As Long, Optional Base As Long=-1)
	CSR -> dense matrix
Sub	CsrSsr (Uplo As String, N As Long, Val() As Double, Rowptr() As Long, Colind() As Long, Val2() As Double, Rowptr2() As Long, Colind2() As Long, Optional Info As Long, Optional Base As Long=-1, Optional Base2 As Long=0)
	CSR (symmetric full matrix) -> SSR (CSR sparse matrix packed form)
Sub	DenseCsr (M As Long, N As Long, A() As Double, Val() As Double, Rowptr() As Long, Colind() As Long, Optional Info As Long, Optional Base As Long=0)
	Dense matrix -> CSR
Sub	SsrCsr (Uplo As String, N As Long, Val() As Double, Rowptr() As Long, Colind() As Long, Val2() As Double, Rowptr2() As Long, Colind2() As Long, Optional Info As Long, Optional Base As Long=-1, Optional Base2 As Long=0)
	SSR (CSR sparse matrix packed form) -> CSR (symmetric full matrix)
Sub	CsrTrans (M As Long, N As Long, Val() As Double, Rowptr() As Long, Colind() As Long, Val2() As Double, Colptr2() As Long, Rowind2() As Long, Optional Info As Long, Optional Base As Long=-1, Optional Base2 As Long=0)
	Transpose of sparse matrix (CSR)
Sub	Bicg1 (N As Long, Val() As Double, Rowptr() As Long, Colind() As Long, B() As Double, X() As Double, Optional Info As Long, Optional Iter As Long, Optional Res As Double, Optional MaxIter As Long=500, Optional Tol As Double=1.0E-10)
	Solution of linear system Ax = b using bi-conjugate gradient (BICG) method (Simple driver)
Sub	Cg1 (N As Long, Val() As Double, Rowptr() As Long, Colind() As Long, B() As Double, X() As Double, Optional Info As Long, Optional Iter As Long, Optional Res As Double, Optional MaxIter As Long=500, Optional Tol As Double=1.0E-10, Optional Uplo As String="F")
	Solution of linear system Ax = b using conjugate gradient (CG) method (Symmetric positive definite matrix) (Simple driver)

Detailed Description

This is the group of D. Linear solve programs. (Sparse matrices)