Version 1.1.0

Code Download

  • Single
  • Double

MI26 Unsymmetric system: BiCGStab (BiConjugate Gradient Stabilized) method

This routine uses the BiCGStab (BiConjugate Gradient Stabilized) method to solve the n × n unsymmetric linear system Ax = b, optionally using preconditioning. If PL, PR are the preconditioning matrices, the routine actually solves the preconditioned system

Āx̄ = b̄,

with Ā = PLAPR and b̄ = PLb and recovers the solution x = PRx̄. If PL = I, preconditioning is said to be from the right, if PR = I, it is said to be from the left, and otherwise it is from both sides. Reverse communication is used for preconditioning operations and matrix-vector products of the form Az.