?trsv

线性方程组求解，系数矩阵是三角矩阵。

$\text{[math]}$ 或 $\text{[math]}$ 或 $\text{[math]}$

其中，x和b为包含n个元素的向量，A是n*n的单元或非单元的三角矩阵。

接口定义

C interface:

void cblas_strsv(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const enum CBLAS_TRANSPOSE TransA, const enum CBLAS_DIAG Diag, const BLASINT N, const float *A, const BLASINT lda, float *X, const BLASINT incX);

void cblas_dtrsv(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const enum CBLAS_TRANSPOSE TransA, const enum CBLAS_DIAG Diag, const BLASINT N, const double *A, const BLASINT lda, double *X, const BLASINT incX);

void cblas_ctrsv(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const enum CBLAS_TRANSPOSE TransA, const enum CBLAS_DIAG Diag, const BLASINT N, const void *A, const BLASINT lda, void *X, const BLASINT incX);

void cblas_ztrsv(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const enum CBLAS_TRANSPOSE TransA, const enum CBLAS_DIAG Diag, const BLASINT N, const void *A, const BLASINT lda, void *X, const BLASINT incX);

Fortran interface:

CALL STRSV(UPLO, TRANS, DIAG, N, A, LDA, X, INCX)

CALL DTRSV(UPLO, TRANS, DIAG, N, A, LDA, X, INCX)

CALL CTRSV(UPLO, TRANS, DIAG, N, A, LDA, X, INCX)

CALL ZTRSV(UPLO, TRANS, DIAG, N, A, LDA, X, INCX)

参数

参数名	类型	描述	输入/输出
order	枚举类型CBLAS_ORDER	表示矩阵是行主序或列主序。	输入
Uplo	枚举类型CBLAS_UPLO	表示使用矩阵A的上三角或下三角。如果Uplo=CblasUpper，则使用A的上三角部分。如果Uplo=CblasLower，则使用A的下三角部分。	输入
TransA	枚举类型CBLAS_TRANSPOSE	矩阵A为常规矩阵，转置矩阵或共轭矩阵。如果trans=CblasNoTrans，则 $\text{[math]}$ 。如果trans=CblasTrans，则 $\text{[math]}$ 。如果trans=CblasConjTrans，则 $\text{[math]}$ 。	输入
Diag	枚举类型CBLAS_DIAG	表示A是否是单元三角阵。如果Diag=CblasUnit，则A是单元三角矩阵。如果Diag=CblasNonUnit，则A不是单元三角矩阵。	输入
N	整型数	表示矩阵A的阶数，要求N大于或等于零。	输入
A	在strsv中，A是单精度浮点类型。在dtrsv中，A是双精度浮点类型。在ctrsv中，A是单精度复数类型。在ztrsv中，A是双精度复数类型。	三角矩阵，矩阵规模至少是lda*N。	输入
lda	整型数	矩阵A中每列存储的地址距离，要求lda大于或等于max(1, N)。	输入
X	在strsv中，X是单精度浮点类型。在dtrsv中，X是双精度浮点类型。在ctrsv中，X是单精度复数类型。在ztrsv中，X是双精度复数类型。	向量X，向量规模至少是(1+(N-1)*abs(incX))。	输入/输出
incX	整型数	表示X中向量增长步长，不能为零。	输入

依赖

#include "kblas.h"

示例

C interface:

    int  n = 4, lda = 4,incx = 1; 
    float x[4] = {4.0, 3.0, 2.0, 1.0}; 
    /** 
     *        |  1.0 1.0  1.0  3.0 | 
     * A   =  |  .   3.0  3.0  2.0 | 
     *        |  .    .   2.0  3.0 | 
     *        |  .    .    .   1.0 | 
     */ 
    float ap[16] = {1.0, 0, 0, 0, 1.0, 3.0, 0, 0, 1.0, 3.0, 2.0, 0, 3.0, 2.0, 3.0, 1.0}; 
 
    cblas_strsv(CblasColMajor, CblasUpper, CblasTrans, CblasNonUnit, n, ap, lda, x, incx); 
    /** 
     * Output x = |4.000000, -0.333333, -0.500000, -8.833333| 
     */

Fortran interface:

      INTEGER :: N=4, INCX=1 
      INTEGER :: LDA=4 
      REAL(4) :: X(4), AP(16) 
      DATA X/4.0, 3.0, 2.0, 1.0/ 
      DATA AP/1.0, 0, 0, 0, 1.0, 3.0, 0, 0, 1.0, 3.0, 2.0, 0, 3.0, 
     $        2.0, 3.0, 1.0/ 
      EXTERNAL STRSV 
      CALL STRSV('U', 'T', 'N', N, A, LDA, X, INCX) 
 
*     Output X = |4.000000, -0.333333, -0.500000, -8.833333|

父主题： KML_BLAS Level 2函数