D. Linear solve (Sparse matrices)

Functions
def	csr_dusmm (trans, order, m, n, l, alpha, val, rowptr, colind, base, b, beta, c)
	C <- αAB + βC or C <- αATB + βC (CSR)
def	csr_dusmv (trans, m, n, alpha, val, rowptr, colind, base, x, incx, beta, y, incy)
	y <- αAx + βy or y <- αATx + βy (CSR)
def	csr_dussm (uplo, trans, diag, order, n, nrhs, val, rowptr, colind, base, x)
	Solution of AX = B or ATX = B (Triangular matrices) (CSR)
def	csr_dussv (uplo, trans, diag, n, val, rowptr, colind, base, x, incx)
	Solution of Ax = b or ATx = b (triangular matrices) (CSR)
def	ssr_dusmv (uplo, n, alpha, val, rowptr, colind, base, x, incx, beta, y, incy)
	y <- αAx + βy (Symmetric matrices) (CSR)
def	coo_csr (m, n, nnz, val, rowind, colind, base, val2, rowptr2, colind2, base2)
	COO -> CSR
def	csr_coo (m, n, val, rowptr, colind, base, val2, rowind2, colind2, base2)
	CSR -> COO
def	csr_dense (m, n, val, rowptr, colind, base, a)
	CSR -> dense matrix
def	csr_ssr (uplo, n, val, rowptr, colind, base, val2, rowptr2, colind2, base2)
	CSR (symmetric full matrix) -> SSR (CSR sparse matrix packed form)
def	csr_trans (m, n, val, ptr, ind, base, val2, ptr2, ind2, base2)
	Transpose matrix (CSR)
def	dense_csr (m, n, a, maxnnz, val, rowptr, colind, base)
	Dense matrix -> CSR
def	ssr_csr (uplo, n, val, rowptr, colind, base, val2, rowptr2, colind2, base2)
	SSR (CSR sparse matrix packed form) -> CSR (symmetric full matrix)
def	bicg1 (n, val, rowptr, colind, b, x, tol=1.0e-10, maxiter=500)
	Solution of linear system Ax = b using bi-conjugate gradient (BICG) method (Simple driver)
def	cg1 (uplo, n, val, rowptr, colind, b, x, tol=1.0e-10, maxiter=500)
	Solution of linear system Ax = b using conjugate gradient (CG) method (Symmetric positive definite) (Simple driver)

Detailed Description

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