NAG Library Routine Document
F07JDF (DPTTRF) computes the modified Cholesky factorization of a real by symmetric positive definite tridiagonal matrix .
The routine may be called by its
F07JDF (DPTTRF) 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.
- 1: – INTEGERInput
On entry: , the order of the matrix .
- 2: – REAL (KIND=nag_wp) arrayInput/Output
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 .
- 3: – REAL (KIND=nag_wp) arrayInput/Output
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
- 4: – INTEGEROutput
unless the routine detects an error (see Section 6
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, F07JEF (DPTTRS)
can be used to solve systems of equations
, and F07JGF (DPTCON)
can be used to estimate the condition number of
8 Parallelism and Performance
F07JDF (DPTTRF) is not threaded in any implementation.
The total number of floating-point operations required to factorize the matrix is proportional to .
The complex analogue of this routine is F07JRF (ZPTTRF)
This example factorizes the symmetric positive definite tridiagonal matrix
10.1 Program Text
Program Text (f07jdfe.f90)
10.2 Program Data
Program Data (f07jdfe.d)
10.3 Program Results
Program Results (f07jdfe.r)