scipy.sparse.linalg.

bicgstab#

scipy.sparse.linalg.bicgstab(A, b, x0=None, *, rtol=1e-05, atol=0.0, maxiter=None, M=None, callback=None)[Quelle]#

Löse Ax = b mit der BIConjugate Gradient STABilized Methode.

Parameter:

A{sparse array, ndarray, LinearOperator}: Die reelle oder komplexe N-by-N-Matrix des linearen Gleichungssystems. Alternativ kann A ein linearer Operator sein, der Ax und A^T x z. B. mit scipy.sparse.linalg.LinearOperator erzeugen kann.
bndarray: Rechte Seite des linearen Systems. Hat die Form (N,) oder (N,1).
x0ndarray: Startschätzung für die Lösung.
rtol, atolfloat, optional: Parameter für den Konvergenztest. Für Konvergenz sollte norm(b - A @ x) <= max(rtol*norm(b), atol) erfüllt sein. Der Standardwert ist atol=0. und rtol=1e-5.
maxiterinteger: Maximale Anzahl von Iterationen. Die Iteration wird nach maxiter-Schritten beendet, auch wenn die angegebene Toleranz nicht erreicht wurde.
M{sparse array, ndarray, LinearOperator}: Vorbedingung für A. Sie sollte eine Annäherung an die Inverse von A sein (siehe Hinweise). Eine effektive Vorbedingung verbessert die Konvergenzrate erheblich, was bedeutet, dass weniger Iterationen benötigt werden, um eine gegebene Fehlertoleranz zu erreichen.
callbackfunction: Vom Benutzer bereitgestellte Funktion, die nach jeder Iteration aufgerufen wird. Sie wird als callback(xk) aufgerufen, wobei xk der aktuelle Lösungsvektor ist.

Rückgabe:

xndarray

Die konvergierte Lösung.

infointeger

Liefert Konvergenzinformationen: 0 : erfolgreicher Abbruch >0 : Konvergenz zur Toleranz nicht erreicht, Anzahl der Iterationen <0 : Parameterzusammenbruch

Hinweise

Die Vorbedingung M sollte eine Matrix sein, so dass M @ A eine kleinere Konditionszahl hat als A, siehe [1].

Referenzen

[1]

„Vorkonditionierer“, Wikipedia, https://de.wikipedia.org/wiki/Vorkonditionierer

[2]

„Biconjugate gradient stabilized method“, Wikipedia, https://en.wikipedia.org/wiki/Biconjugate_gradient_stabilized_method

Beispiele

>>> import numpy as np
>>> from scipy.sparse import csc_array
>>> from scipy.sparse.linalg import bicgstab
>>> R = np.array([[4, 2, 0, 1],
...               [3, 0, 0, 2],
...               [0, 1, 1, 1],
...               [0, 2, 1, 0]])
>>> A = csc_array(R)
>>> b = np.array([-1, -0.5, -1, 2])
>>> x, exit_code = bicgstab(A, b, atol=1e-5)
>>> print(exit_code)  # 0 indicates successful convergence
0
>>> np.allclose(A.dot(x), b)
True