home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2001 May / SGI IRIX Base Documentation 2001 May.iso / usr / share / catman / p_man / cat3 / complib / ctptrs.z / ctptrs
Encoding:
Text File  |  1998-10-30  |  3.0 KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCTTTTPPPPTTTTRRRRSSSS((((3333FFFF))))                                                          CCCCTTTTPPPPTTTTRRRRSSSS((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CTPTRS - solve a triangular system of the form   A * X = B, A**T * X = B,
  10.      or A**H * X = B,
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE CTPTRS( UPLO, TRANS, DIAG, N, NRHS, AP, B, LDB, INFO )
  14.  
  15.          CHARACTER      DIAG, TRANS, UPLO
  16.  
  17.          INTEGER        INFO, LDB, N, NRHS
  18.  
  19.          COMPLEX        AP( * ), B( LDB, * )
  20.  
  21. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  22.      CTPTRS solves a triangular system of the form
  23.  
  24.      where A is a triangular matrix of order N stored in packed format, and B
  25.      is an N-by-NRHS matrix.  A check is made to verify that A is nonsingular.
  26.  
  27.  
  28. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  29.      UPLO    (input) CHARACTER*1
  30.              = 'U':  A is upper triangular;
  31.              = 'L':  A is lower triangular.
  32.  
  33.      TRANS   (input) CHARACTER*1
  34.              Specifies the form of the system of equations:
  35.              = 'N':  A * X = B     (No transpose)
  36.              = 'T':  A**T * X = B  (Transpose)
  37.              = 'C':  A**H * X = B  (Conjugate transpose)
  38.  
  39.      DIAG    (input) CHARACTER*1
  40.              = 'N':  A is non-unit triangular;
  41.              = 'U':  A is unit triangular.
  42.  
  43.      N       (input) INTEGER
  44.              The order of the matrix A.  N >= 0.
  45.  
  46.      NRHS    (input) INTEGER
  47.              The number of right hand sides, i.e., the number of columns of
  48.              the matrix B.  NRHS >= 0.
  49.  
  50.      AP      (input) COMPLEX array, dimension (N*(N+1)/2)
  51.              The upper or lower triangular matrix A, packed columnwise in a
  52.              linear array.  The j-th column of A is stored in the array AP as
  53.              follows:  if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
  54.              if UPLO = 'L', AP(i + (j-1)*(2*n-j)/2) = A(i,j) for j<=i<=n.
  55.  
  56.      B       (input/output) COMPLEX array, dimension (LDB,NRHS)
  57.              On entry, the right hand side matrix B.  On exit, if INFO = 0,
  58.              the solution matrix X.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCTTTTPPPPTTTTRRRRSSSS((((3333FFFF))))                                                          CCCCTTTTPPPPTTTTRRRRSSSS((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      LDB     (input) INTEGER
  75.              The leading dimension of the array B.  LDB >= max(1,N).
  76.  
  77.      INFO    (output) INTEGER
  78.              = 0:  successful exit
  79.              < 0:  if INFO = -i, the i-th argument had an illegal value
  80.              > 0:  if INFO = i, the i-th diagonal element of A is zero,
  81.              indicating that the matrix is singular and the solutions X have
  82.              not been computed.
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.