?sbmv

向量与对称带状矩阵乘积，即 $\text{[math]}$ 。

alpha和beta是乘法系数，x和y是含有n个元素的向量，A是n*n的对称带状矩阵。

接口定义

C interface:

void cblas_ssbmv(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, const BLASINT K, const float alpha, const float *A, const BLASINT lda, const float *X, const BLASINT incX, const float beta, float *Y, const BLASINT incY);

void cblas_dsbmv(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, const BLASINT K, const double alpha, const double *A, const BLASINT lda, const double *X, const BLASINT incX, const double beta, double *Y, const BLASINT incY);

Fortran interface:

CALL SSBMV(UPLO, N, K, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)

CALL DSBMV(UPLO, N, K, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)

参数

参数名	类型	描述	输入/输出
order	枚举类型CBLAS_ORDER	表示矩阵是行主序或列主序。	输入
Uplo	枚举类型CBLAS_UPLO	表示使用矩阵A的上三角或下三角。如果Uplo=CblasUpper，则使用A的上三角部分。如果Uplo=CblasLower，则使用A的下三角部分。	输入
N	整型数	表示矩阵A的阶数，要求N大于或等于零。	输入
K	整型数	表示矩阵A的超对角阶数，要求K大于或等于零。	输入
alpha	在ssbmv中，alpha是单精度浮点类型。在dsbmv中，alpha是双精度浮点类型。	表示乘法系数。	输入
A	在ssbmv中，A是单精度浮点类型。在dsbmv中，A是双精度浮点类型。	对称带状矩阵，矩阵规模是lda*n。	输入
lda	整型数	矩阵A中每列存储的地址距离，要求lda大于或等于(k+1)。	输入
X	在ssbmv中，X是单精度浮点类型。在dsbmv中，X是双精度浮点类型。	向量X，向量规模至少是(1+(N-1)*abs(incX))。	输入
incX	整型数	表示X中向量增长步长，不能为零。	输入
beta	在ssbmv中，beta是单精度浮点类型。在dsbmv中，beta是双精度浮点类型。	乘法系数。	输入
Y	在ssbmv中，Y是单精度浮点类型。在dsbmv中，Y是双精度浮点类型。	向量Y，向量规模至少是(1+(N-1)*abs(incY))。	输入/输出
incY	整型数	表示Y中向量增长步长，不能为零。	输入

依赖

#include "kblas.h"

示例

C interface:

    int n = 4, k = 2; 
    float alpha = 2.0; 
    /** 
     *        |  .    .   2.0  2.0 | 
     *        |  .   1.0  2.0  3.0 | 
     *        | 1.0  2.0  2.0  4.0 | 
     *  A  =  | 1.0  3.0  2.0  5.0 | 
     *        | 1.0  4.0  2.0  6.0 | 
     *        | 1.0  5.0  2.0  2.0 | 
     *        |  .    .    .    .  | 
     *        |  .    .    .    .  | 
     */ 
    float a[32] = {0, 0, 1.0, 1.0, 1.0, 1.0, 0, 0, 
                   0, 1.0, 2.0, 3.0, 4.0, 5.0, 0, 0, 
                   2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 0, 0, 
                   2.0, 3.0, 4.0, 5.0, 6.0, 2.0, 0, 0}; 
    float x[4] = {2.0, 2.0, 3.0, 4.0}; 
    float y[4] = {1.0, 1.0, 13.0, 4.0}; 
 
    cblas_ssbmv(CblasColMajor,CblasUpper, n, k, alpha, a, 8, x, 1, 1.0, y, 1); 
    /** 
     * Output y = |21.0, 41.0, 65.0, 62.0| 
     */

Fortran interface:

      INTEGER :: N=4, K=2 
      REAL(4) :: ALPHA=2.0 
      REAL(4) :: BETA=1.0 
      REAL(4) :: A(32) 
      DATA A/0, 0, 1.0, 1.0, 1.0, 1.0, 0, 0, 
     $       0, 1.0, 2.0, 3.0, 4.0, 5.0, 0, 0, 
     $       2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 0, 0, 
     $       2.0, 3.0, 4.0, 5.0, 6.0, 2.0, 0, 0/ 
      REAL(4) :: X(4), Y(4) 
      DATA X/2.0, 2.0, 3.0, 4.0/ 
      DATA Y/1.0, 1.0, 13.0, 4.0/ 
      EXTERNAL SSBMV 
      CALL SSBMV('L', N, K, ALPHA, A, 8, X, 1, BETA, Y, 1) 
 
*     Output Y = |21.0, 41.0, 65.0, 62.0|

父主题： KML_BLAS Level 2函数