D2a4. 連立一次方程式 (一般行列) (反復法ソルバー)

関数
void	bicg (int n, void(*matvec)(int, const double[], double[]), void(*matvectrans)(int, const double[], double[]), void(*psolve)(int, const double[], double[]), void(*psolvetrans)(int, const double[], double[]), void(*chkconv)(int, const double[], double, int, int *), const double b[], double x[], int maxiter, int *iter, double *res, int lwork, double work[], int *info)
	双共役勾配(BICG)法による連立一次方程式 Ax = b の解
void	bicg1 (int n, const double val[], const int rowptr[], const int colind[], const double b[], double x[], double tol, int maxiter, int *iter, double *res, int lwork, double work[], int *info)
	双共役勾配(BICG)法による連立一次方程式 Ax = b の解 (シンプルドライバ)
void	bicg_r (int n, const double b[], double x[], int maxiter, int *iter, double *res, int lwork, double work[], int *info, double xx[], double yy[], int *irev)
	双共役勾配(BICG)法による連立一次方程式 Ax = b の解 (リバースコミュニケーション版)
void	cgs (int n, void(*matvec)(int, const double[], double[]), void(*psolve)(int, const double[], double[]), void(*chkconv)(int, const double[], double, int, int *), const double b[], double x[], int maxiter, int *iter, double *res, int lwork, double work[], int *info)
	二乗共役勾配(CGS)法による連立一次方程式 Ax = b の解
void	cgs_r (int n, const double b[], double x[], int maxiter, int *iter, double *res, int lwork, double work[], int *info, double xx[], double yy[], int *irev)
	二乗共役勾配(CGS)法による連立一次方程式 Ax = b の解 (リバースコミュニケーション版)
void	diom (int n, void(*matvec)(int, const double[], double[]), void(*psolve)(int, const double[], double[]), void(*chkconv)(int, const double[], double, int, int *), const double b[], double x[], int m, int maxiter, int *iter, double *res, int lwork, double work[], int *info)
	不完全直交化法(DIOM)による連立一次方程式 Ax = b の解
void	diom_r (int n, const double b[], double x[], int m, int maxiter, int *iter, double *res, int lwork, double work[], int *info, double xx[], double yy[], int *irev)
	不完全直交化法(DIOM)による連立一次方程式 Ax = b の解 (リバースコミュニケーション版)
void	dqgmres (int n, void(*matvec)(int, const double[], double[]), void(*psolve)(int, const double[], double[]), void(*chkconv)(int, const double[], double, int, int *), const double b[], double x[], int m, int maxiter, int *iter, double *res, int lwork, double work[], int *info)
	疑似最小残差(DQGMRES)法による連立一次方程式 Ax = b の解
void	dqgmres_r (int n, const double b[], double x[], int m, int maxiter, int *iter, double *res, int lwork, double work[], int *info, double xx[], double yy[], int *irev)
	疑似最小残差(DQGMRES)法による連立一次方程式 Ax = b の解 (リバースコミュニケーション版)
void	fgmres (int n, void(*matvec)(int, const double[], double[]), void(*psolve)(int, const double[], double[]), void(*chkconv)(int, const double[], double, int, int *), const double b[], double x[], int m, int maxiter, int *iter, double *res, int lwork, double work[], int *info)
	最小残差(FGMRES)法による連立一次方程式 Ax = b の解
void	fgmres_r (int n, const double b[], double x[], int m, int maxiter, int *iter, double *res, int lwork, double work[], int *info, double xx[], double yy[], int *irev)
	最小残差(FGMRES)法による連立一次方程式 Ax = b の解 (リバースコミュニケーション版)
void	fom (int n, void(*matvec)(int, const double[], double[]), void(*psolve)(int, const double[], double[]), void(*chkconv)(int, const double[], double, int, int *), const double b[], double x[], int m, int maxiter, int *iter, double *res, int lwork, double work[], int *info)
	完全直交化法(FOM)による連立一次方程式 Ax = b の解
void	fom_r (int n, const double b[], double x[], int m, int maxiter, int *iter, double *res, int lwork, double work[], int *info, double xx[], double yy[], int *irev)
	完全直交化法(FOM)による連立一次方程式 Ax = b の解 (リバースコミュニケーション版)
void	gcr (int n, void(*matvec)(int, const double[], double[]), void(*psolve)(int, const double[], double[]), void(*chkconv)(int, const double[], double, int, int *), const double b[], double x[], int m, int maxiter, int *iter, double *res, int lwork, double work[], int *info)
	一般化共役残差(GCR)法による連立一次方程式 Ax = b の解
void	gcr_r (int n, const double b[], double x[], int m, int maxiter, int *iter, double *res, int lwork, double work[], int *info, double xx[], double yy[], int *irev)
	一般化共役残差(GCR)法による連立一次方程式 Ax = b の解 (リバースコミュニケーション版)
void	gpbicg (int n, void(*matvec)(int, const double[], double[]), void(*psolve)(int, const double[], double[]), void(*chkconv)(int, const double[], double, int, int *), const double b[], double x[], int mode, int maxiter, int *iter, double *res, int lwork, double work[], int *info)
	積型双共役勾配(GPBICG)法, 安定化双共役勾配(BICGSTAB)法 または BICGSTAB2法による連立一次方程式 Ax = b の解
void	gpbicg_r (int n, const double b[], double x[], int mode, int maxiter, int *iter, double *res, int lwork, double work[], int *info, double xx[], double yy[], int *irev)
	積型双共役勾配(GPBICG)法, 安定化双共役勾配(BICGSTAB)法 または BICGSTAB2法による連立一次方程式 Ax = b の解 (リバースコミュニケーション版)
void	orthomin (int n, void(*matvec)(int, const double[], double[]), void(*psolve)(int, const double[], double[]), void(*chkconv)(int, const double[], double, int, int *), const double b[], double x[], int m, int maxiter, int *iter, double *res, int lwork, double work[], int *info)
	Orthomin法による連立一次方程式 Ax = b の解
void	orthomin_r (int n, const double b[], double x[], int m, int maxiter, int *iter, double *res, int lwork, double work[], int *info, double xx[], double yy[], int *irev)
	Orthomin法による連立一次方程式 Ax = b の解 (リバースコミュニケーション版)
void	qmr (int n, void(*matvec)(int, const double[], double[]), void(*matvectrans)(int, const double[], double[]), void(*psolve)(int, const double[], double[]), void(*psolvetrans)(int, const double[], double[]), void(*chkconv)(int, const double[], double, int, int *), const double b[], double x[], int maxiter, int *iter, double *res, int lwork, double work[], int *info)
	疑似最小残差(QMR)法による連立一次方程式 Ax = b の解
void	qmr_r (int n, const double b[], double x[], int maxiter, int *iter, double *res, int lwork, double work[], int *info, double xx[], double yy[], int *irev)
	疑似最小残差(QMR)法による連立一次方程式 Ax = b の解 (リバースコミュニケーション版)
void	sor (int n, void(*matvec)(int, const double[], double[]), void(*matsol)(int, const double[], double[]), void(*chkconv)(int, const double[], double, int, int *), const double b[], double x[], int maxiter, int *iter, double *res, int lwork, double work[], int *info)
	逐次的過剰緩和(SOR)法による連立一次方程式 Ax = b の解
void	sor_r (int n, const double b[], double x[], int maxiter, int *iter, double *res, int lwork, double work[], int *info, double xx[], double yy[], int *irev)
	逐次的過剰緩和(SOR)法による連立一次方程式 Ax = b の解 (リバースコミュニケーション版)
void	tfqmr (int n, void(*matvec)(int, const double[], double[]), void(*psolve)(int, const double[], double[]), void(*chkconv)(int, const double[], double, int, int *), const double b[], double x[], int maxiter, int *iter, double *res, int lwork, double work[], int *info)
	転置不要疑似最小残差(TFQMR)法による連立一次方程式 Ax = b の解
void	tfqmr_r (int n, const double b[], double x[], int maxiter, int *iter, double *res, int lwork, double work[], int *info, double xx[], double yy[], int *irev)
	転置不要疑似最小残差(TFQMR)法による連立一次方程式 Ax = b の解 (リバースコミュニケーション版)

詳解

D2a4. 連立一次方程式 (一般行列) (反復法ソルバー) プログラムを表示しています.