## File: Triangular-Systems.html

package info (click to toggle)
gsl-ref-html 2.3-1
• area: non-free
• in suites: bullseye, sid
• size: 6,876 kB
• ctags: 4,574
• sloc: makefile: 35
 file content (99 lines) | stat: -rw-r--r-- 5,655 bytes parent folder | download
 123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899  GNU Scientific Library – Reference Manual: Triangular Systems

Next: , Previous: Tridiagonal Systems, Up: Linear Algebra   [Index]

14.18 Triangular Systems

Function: int gsl_linalg_tri_upper_invert (gsl_matrix * T)
Function: int gsl_linalg_tri_lower_invert (gsl_matrix * T)
Function: int gsl_linalg_tri_upper_unit_invert (gsl_matrix * T)
Function: int gsl_linalg_tri_lower_unit_invert (gsl_matrix * T)

These functions calculate the in-place inverse of the triangular matrix T. When the upper prefix is specified, then the upper triangle of T is used, and when the lower prefix is specified, the lower triangle is used. If the unit prefix is specified, then the diagonal elements of the matrix T are taken as unity and are not referenced. Otherwise the diagonal elements are used in the inversion.

Function: int gsl_linalg_tri_upper_rcond (const gsl_matrix * T, double * rcond, gsl_vector * work)
Function: int gsl_linalg_tri_lower_rcond (const gsl_matrix * T, double * rcond, gsl_vector * work)

These functions estimate the reciprocal condition number, in the 1-norm, of the upper or lower N-by-N triangular matrix T. The reciprocal condition number is stored in rcond on output, and is defined by 1 / (||T||_1 \cdot ||T^{-1}||_1). Additional workspace of size 3 N is required in work.