D. 線形計算 (疎行列)

関数
def	csr_dusmm (trans, order, m, n, l, alpha, val, rowptr, colind, base, b, beta, c)
	C <- αAB + βC または C <- αATB + βC (CSR)
def	csr_dusmv (trans, m, n, alpha, val, rowptr, colind, base, x, incx, beta, y, incy)
	y <- αAx + βy または y <- αATx + βy (CSR)
def	csr_dussm (uplo, trans, diag, order, n, nrhs, val, rowptr, colind, base, x)
	AX = B または ATX = B の解 (三角行列) (CSR)
def	csr_dussv (uplo, trans, diag, n, val, rowptr, colind, base, x, incx)
	Ax = b または ATx = b の解 (三角行列) (CSR)
def	ssr_dusmv (uplo, n, alpha, val, rowptr, colind, base, x, incx, beta, y, incy)
	y <- αAx + βy (対称行列) (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 -> 密行列
def	csr_ssr (uplo, n, val, rowptr, colind, base, val2, rowptr2, colind2, base2)
	CSR (対称なフル行列) -> SSR (CSR 対称行列形式)
def	csr_trans (m, n, val, ptr, ind, base, val2, ptr2, ind2, base2)
	転置行列を求める (CSR)
def	dense_csr (m, n, a, maxnnz, val, rowptr, colind, base)
	密行列 -> CSR
def	ssr_csr (uplo, n, val, rowptr, colind, base, val2, rowptr2, colind2, base2)
	SSR (CSR 対称行列形式) -> CSR (対称なフル行列)
def	bicg1 (n, val, rowptr, colind, b, x, tol=1.0e-10, maxiter=500)
	双共役勾配(BICG)法による連立一次方程式 Ax = b の解 (シンプルドライバ)
def	cg1 (uplo, n, val, rowptr, colind, b, x, tol=1.0e-10, maxiter=500)
	共役勾配(CG)法による連立一次方程式 Ax = b の解 (正定値対称行列) (シンプルドライバ)

詳解

D. 線形計算 (疎行列) プログラムを表示しています.