naginterfaces.library.matop.real_symm_posdef_fac¶
- naginterfaces.library.matop.real_symm_posdef_fac(k, a)[source]¶
real_symm_posdef_fac
performs a decomposition of a real symmetric positive definite band matrix.For full information please refer to the NAG Library document for f01bu
https://support.nag.com/numeric/nl/nagdoc_30.2/flhtml/f01/f01buf.html
- Parameters
- kint
, the change-over point in the decomposition.
- afloat, array-like, shape
The upper triangle of the symmetric band matrix , with the diagonal of the matrix stored in the th row of the array, and the superdiagonals within the band stored in the first rows of the array. Each column of the matrix is stored in the corresponding column of the array. For example, if and , the storage scheme is
Elements in the top left corner of the array are not used. The matrix elements within the band can be assigned to the correct elements of the array using the following code:
for j in range(A.shape[1]): for i in range(max(0, j+1-m1), j+1): a[i-j+m1-1, j] = A[i, j]
- Returns
- afloat, ndarray, shape
is overwritten by the corresponding elements of , and .
- Raises
- NagValueError
- (errno )
On entry, and .
Constraint: .
- (errno )
On entry, and .
Constraint: .
- (errno )
The matrix is not positive definite.
- (errno )
The matrix is not positive definite.
- Notes
No equivalent traditional C interface for this routine exists in the NAG Library.
The symmetric positive definite matrix , of order and bandwidth , is divided into the leading principal sub-matrix of order and its complement, where . A decomposition of the latter and an decomposition of the former are obtained by means of a sequence of elementary transformations, where is unit upper triangular, is unit lower triangular and is diagonal. Thus if , an decomposition of is obtained.
This function is specifically designed to precede
real_symm_posdef_geneig()
for the transformation of the symmetric-definite eigenproblem by the method of Crawford where and are of band form. In this context, is chosen to be close to and the decomposition is applied to the matrix .
- References
Wilkinson, J H, 1965, The Algebraic Eigenvalue Problem, Oxford University Press, Oxford
Wilkinson, J H and Reinsch, C, 1971, Handbook for Automatic Computation II, Linear Algebra, Springer–Verlag