scipy.linalg.lapack.zhesvx#

scipy.linalg.lapack.zhesvx(a, b, af=None, ipiv=None, lwork=max(2 * n, 1), factored=0, lower=0, overwrite_a=0, overwrite_b=0) = <flapack function zhesvx>#

Solve a Hermitian indefinite system with condition estimation and error bounds (LAPACK zhesvx).

Parameters:

andarray: Hermitian matrix of shape (n, n).
bndarray: Right-hand side(s) of shape (n, nrhs).
afndarray, optional: Factorization to reuse when factored=1; otherwise it is computed.
ipivndarray, optional: Pivot indices to reuse when factored=1, 1-based; otherwise computed.
lworkint, optional: Size of the workspace. Default is max(2 * n, 1). Use zhesvx_lwork for the optimal value.
factoredint, optional: If nonzero, af and ipiv already hold the factorization and are reused. Default is 0.
lowerint, optional: If nonzero, the lower triangle of a is referenced; otherwise the upper. Default is 0.
overwrite_aint, optional: If nonzero, a may be overwritten in place. Default is 0.
overwrite_bint, optional: If nonzero, b may be overwritten in place. Default is 0.

Returns:

uduhndarray: The block diagonal D and the multipliers.
ipivndarray: Pivot indices describing the interchanges and the 1x1 / 2x2 block structure, 1-based. Pass them on to sytrs/sytri/sycon unchanged.
xndarray: Solution of the system.
rcondfloat: Estimate of the reciprocal condition number.
ferrndarray: Estimated forward error bound for each solution vector.
berrndarray: Componentwise relative backward error of each solution vector.
infoint: 0 on success; if negative, the -info-th argument had an illegal value; if 0 < info <= n, d[info-1, info-1] is exactly zero; if info = n+1, the matrix is singular to working precision.