NAG FL Interface
f07jrf computes the modified Cholesky factorization of a complex by Hermitian positive definite tridiagonal matrix .
|Integer, Intent (In)
|Integer, Intent (Out)
|Real (Kind=nag_wp), Intent (Inout)
|Complex (Kind=nag_wp), Intent (Inout)
|C Header Interface
f07jrf_ (const Integer *n, double d, Complex e, Integer *info)
The routine may be called by the names f07jrf, nagf_lapacklin_zpttrf or its LAPACK name zpttrf.
factorizes the matrix
is a unit lower bidiagonal matrix and
is a diagonal matrix with positive diagonal elements. The factorization may also be regarded as having the form
is a unit upper bidiagonal matrix.
On entry: , the order of the matrix .
– Real (Kind=nag_wp) array
the dimension of the array d
must be at least
On entry: must contain the diagonal elements of the matrix .
On exit: is overwritten by the diagonal elements of the diagonal matrix from the factorization of .
– Complex (Kind=nag_wp) array
the dimension of the array e
must be at least
On entry: must contain the subdiagonal elements of the matrix .
: is overwritten by the
subdiagonal elements of the lower bidiagonal matrix
can also be regarded as containing the
superdiagonal elements of the upper bidiagonal matrix
unless the routine detects an error (see Section 6
Error Indicators and Warnings
If , argument had an illegal value. An explanatory message is output, and execution of the program is terminated.
The leading minor of order is not positive definite,
the factorization could not be completed.
The leading minor of order is not positive definite, the factorization was
completed, but .
The computed factorization satisfies an equation of the form
is the machine precision
Following the use of this routine, f07jsf
can be used to solve systems of equations
, and f07juf
can be used to estimate the condition number of
Parallelism and Performance
f07jrf is not threaded in any implementation.
The total number of floating-point operations required to factorize the matrix is proportional to .
The real analogue of this routine is f07jdf
This example factorizes the Hermitian positive definite tridiagonal matrix