KryPbcgs.c File Reference

Krylov subspace methods – Preconditioned BiCGstab. More...

#include <math.h>
#include <float.h>
#include "fasp.h"
#include "fasp_functs.h"
#include "KryUtil.inl"

Go to the source code of this file.

Functions
INT	fasp_solver_dcsr_pbcgs (dCSRmat *A, dvector *b, dvector *u, precond *pc, const REAL tol, const REAL abstol, const INT MaxIt, const SHORT StopType, const SHORT PrtLvl)
	Preconditioned BiCGstab method for solving Au=b for CSR matrix. More...
INT	fasp_solver_dbsr_pbcgs (dBSRmat *A, dvector *b, dvector *u, precond *pc, const REAL tol, const REAL abstol, const INT MaxIt, const SHORT StopType, const SHORT PrtLvl)
	Preconditioned BiCGstab method for solving Au=b for BSR matrix. More...
INT	fasp_solver_dblc_pbcgs (dBLCmat *A, dvector *b, dvector *u, precond *pc, const REAL tol, const REAL abstol, const INT MaxIt, const SHORT StopType, const SHORT PrtLvl)
	Preconditioned BiCGstab method for solving Au=b for BLC matrix. More...
INT	fasp_solver_dstr_pbcgs (dSTRmat *A, dvector *b, dvector *u, precond *pc, const REAL tol, const REAL abstol, const INT MaxIt, const SHORT StopType, const SHORT PrtLvl)
	Preconditioned BiCGstab method for solving Au=b for STR matrix. More...
INT	fasp_solver_pbcgs (mxv_matfree *mf, dvector *b, dvector *u, precond *pc, const REAL tol, const REAL abstol, const INT MaxIt, const SHORT StopType, const SHORT PrtLvl)
	Preconditioned BiCGstab method for solving Au=b. More...

Detailed Description

Krylov subspace methods – Preconditioned BiCGstab.

Note

This file contains Level-3 (Kry) functions. It requires: AuxArray.c, AuxMemory.c, AuxMessage.c, BlaArray.c, BlaSpmvBLC.c, BlaSpmvBSR.c, BlaSpmvCSR.c, and BlaSpmvSTR.c

This version is based on Matlab 2011a – Chunsheng Feng

See KrySPbcgs.c for a safer version

Reference: Y. Saad 2003 Iterative methods for sparse linear systems (2nd Edition), SIAM

---

Copyright (C) 2016–Present by the FASP team. All rights reserved.

Released under the terms of the GNU Lesser General Public License 3.0 or later.

TODO: abstol not used yet! –Chensong

Definition in file KryPbcgs.c.

Function Documentation

fasp_solver_dblc_pbcgs()

INT fasp_solver_dblc_pbcgs	(	dBLCmat *	A,
		dvector *	b,
		dvector *	u,
		precond *	pc,
		const REAL	tol,
		const REAL	abstol,
		const INT	MaxIt,
		const SHORT	StopType,
		const SHORT	PrtLvl
	)

Preconditioned BiCGstab method for solving Au=b for BLC matrix.

Parameters

A	Pointer to coefficient matrix
b	Pointer to dvector of right hand side
u	Pointer to dvector of DOFs
pc	Pointer to precond: structure of precondition
tol	Tolerance for relative residual
abstol	Tolerance for absolute residual
MaxIt	Maximal number of iterations
StopType	Stopping criteria type
PrtLvl	How much information to print out

Returns

Iteration number if converges; ERROR otherwise.

Author

Chunsheng Feng

Date

03/04/2016

Definition at line 705 of file KryPbcgs.c.

fasp_solver_dbsr_pbcgs()

INT fasp_solver_dbsr_pbcgs	(	dBSRmat *	A,
		dvector *	b,
		dvector *	u,
		precond *	pc,
		const REAL	tol,
		const REAL	abstol,
		const INT	MaxIt,
		const SHORT	StopType,
		const SHORT	PrtLvl
	)

Preconditioned BiCGstab method for solving Au=b for BSR matrix.

Parameters

A	Pointer to coefficient matrix
b	Pointer to dvector of right hand side
u	Pointer to dvector of DOFs
pc	Pointer to precond: structure of precondition
tol	Tolerance for relative residual
abstol	Tolerance for absolute residual
MaxIt	Maximal number of iterations
StopType	Stopping criteria type
PrtLvl	How much information to print out

Returns

Iteration number if converges; ERROR otherwise.

Author

Chunsheng Feng

Date

03/04/2016

Definition at line 383 of file KryPbcgs.c.

fasp_solver_dcsr_pbcgs()

INT fasp_solver_dcsr_pbcgs	(	dCSRmat *	A,
		dvector *	b,
		dvector *	u,
		precond *	pc,
		const REAL	tol,
		const REAL	abstol,
		const INT	MaxIt,
		const SHORT	StopType,
		const SHORT	PrtLvl
	)

Preconditioned BiCGstab method for solving Au=b for CSR matrix.

Parameters

A	Pointer to coefficient matrix
b	Pointer to dvector of right hand side
u	Pointer to dvector of DOFs
pc	Pointer to precond: structure of precondition
tol	Tolerance for relative residual
abstol	Tolerance for absolute residual
MaxIt	Maximal number of iterations
StopType	Stopping criteria type
PrtLvl	How much information to print out

Returns

Iteration number if converges; ERROR otherwise.

Author

Chunsheng Feng

Date

03/04/2016

Definition at line 62 of file KryPbcgs.c.

fasp_solver_dstr_pbcgs()

INT fasp_solver_dstr_pbcgs	(	dSTRmat *	A,
		dvector *	b,
		dvector *	u,
		precond *	pc,
		const REAL	tol,
		const REAL	abstol,
		const INT	MaxIt,
		const SHORT	StopType,
		const SHORT	PrtLvl
	)

Preconditioned BiCGstab method for solving Au=b for STR matrix.

Parameters

A	Pointer to coefficient matrix
b	Pointer to dvector of right hand side
u	Pointer to dvector of DOFs
pc	Pointer to precond: structure of precondition
tol	Tolerance for relative residual
abstol	Tolerance for absolute residual
MaxIt	Maximal number of iterations
StopType	Stopping criteria type
PrtLvl	How much information to print out

Returns

Iteration number if converges; ERROR otherwise.

Author

Chunsheng Feng

Date

03/04/2016

Definition at line 1027 of file KryPbcgs.c.

fasp_solver_pbcgs()

INT fasp_solver_pbcgs	(	mxv_matfree *	mf,
		dvector *	b,
		dvector *	u,
		precond *	pc,
		const REAL	tol,
		const REAL	abstol,
		const INT	MaxIt,
		const SHORT	StopType,
		const SHORT	PrtLvl
	)

Preconditioned BiCGstab method for solving Au=b.

Parameters

mf	Pointer to mxv_matfree: spmv operation
b	Pointer to dvector of right hand side
u	Pointer to dvector of DOFs
pc	Pointer to precond: structure of precondition
tol	Tolerance for relative residual
abstol	Tolerance for absolute residual
MaxIt	Maximal number of iterations
StopType	Stopping criteria type
PrtLvl	How much information to print out

Returns

Iteration number if converges; ERROR otherwise.

Author

Chunsheng Feng

Date

03/04/2016

Definition at line 1349 of file KryPbcgs.c.