D. 線形計算 (疎行列)

関数
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 または 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 または 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)
	AX = B または ATX = B の解 (三角行列) (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)
	Ax = b または ATx = b の解 (三角行列) (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) (対称行列)
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 -> 密行列
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 (対称なフル行列) -> SSR (CSR 対称行列形式)
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)
	密行列 -> 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 対称行列形式) -> CSR (対称なフル行列)
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)
	疎行列の転置 (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)
	双共役勾配(BICG)法による連立一次方程式 Ax = b の解 (シンプルドライバ)
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")
	共役勾配(CG)法による連立一次方程式 Ax = b の解 (正定値対称行列) (シンプルドライバ)

詳解

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