D2a4. Solution of systems of linear equations (General matrices) (Iterative solvers)

Functions
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)
	Solution of linear system Ax = b using bi-conjugate gradient (BICG) method
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)
	Solution of linear system Ax = b using bi-conjugate gradient (BICG) method (Simple driver)
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)
	Solution of linear system Ax = b using bi-conjugate gradient (BICG) method (Reverse communication version)
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)
	Solution of linear system Ax = b using conjugate gradient squared (CGS) method
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)
	Solution of linear system Ax = b using conjugate gradient squared (CGS) method (Reverse communication version)
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)
	Solution of linear system Ax = b using direct incomplete orthogonalization method (DIOM)
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)
	Solution of linear system Ax = b using direct incomplete orthogonalization method (DIOM) (Reverse communication version)
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)
	Solution of linear system Ax = b using direct quasi generalized minimum residual (DQGMRES) method
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)
	Solution of linear system Ax = b using direct quasi generalized minimum residual (DQGMRES) method (Reverse communication version)
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)
	Solution of linear system Ax = b using generalized minimum residual (FGMRES) method
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)
	Solution of linear system Ax = b using generalized minimum residual (FGMRES) method (Reverse communication version)
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)
	Solution of linear system Ax = b using full orthogonalization method (FOM)
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)
	Solution of linear system Ax = b using full orthogonalization method (FOM) (Reverse communication version)
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)
	Solution of linear system Ax = b using generalized conjugate residual (GCR) method
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)
	Solution of linear system Ax = b using generalized conjugate residual (GCR) method (Reverse communication version)
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)
	Solution of linear system Ax = b using general product bi-conjugate gradient (GPBICG) method, bi-conjugate gradient stabilized (BICGSTAB) method or BICGSTAB2 method
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)
	Solution of linear system Ax = b using conjugate gradient squared (CGS) method (Reverse communication version)
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)
	Solution of linear system Ax = b using orthomin method
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)
	Solution of linear system Ax = b using orthomin method (Reverse communication version)
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)
	Solution of linear system Ax = b using quasi minimum residual (QMR) method
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)
	Solution of linear system Ax = b using quasi minimum residual (QMR) method (Reverse communication version)
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)
	Solution of linear system Ax = b using successive over-relaxation (SOR) method
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)
	Solution of linear system Ax = b using Successive over-relaxation (SOR) method (Reverse communication version)
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)
	Solution of linear system Ax = b using transpose free quasi minimum residual (TFQMR) method
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)
	Solution of linear system Ax = b using transpose free quasi minimum residual (TFQMR) method (Reverse communication version)

Detailed Description

This is the group of D2a4. Solution of systems of linear equations (General matrices) (Iterative solvers)