dtrsl
dtrsl - (obsolete) solve the linear system Ax = b for a triangular matrix A and vectors b
and x.
SUBROUTINE DTRSL( A, LDA, N, B, JOB, INFO)
INTEGER LDA, N, JOB, INFO
DOUBLE PRECISION A(LDA,*), B(*)
SUBROUTINE DTRSL_64( A, LDA, N, B, JOB, INFO)
INTEGER*8 LDA, N, JOB, INFO
DOUBLE PRECISION A(LDA,*), B(*)
#include <sunperf.h>
void dtrsl(double *a, int lda, int n, double *b, int job, int *info);
void dtrsl_64(double *a, long lda, long n, double *b, long job, long *info);
-
* A (input/output)
-
Matrix A.
-
* LDA (input)
-
Leading dimension of the array A as specified in a dimension or type statement. LDA >= max(1,N).
-
* N (input)
-
Order of the matrix A. N >= 0.
-
* B (input/output)
-
On entry, the right-hand side vector b.
On exit, the solution vector x.
-
* JOB (input)
-
Determines which operation the subroutine will perform:
- 00 Solve Ax = b, A lower triangular.
- 01 Solve Ax = b, A upper triangular.
- 10 Solve AHx = b, A lower triangular.
- 11 Solve AHx = b, A upper triangular.
- Note that ATx = AHx for real matrices.
-
* INFO (output)
-
On exit:
INFO = 0 Subroutine completed normally.
INFO > 0 Returns the index of the first zero diagonal element of A.