?symv

向量与对称矩阵乘积

即 $\text{[math]}$ 。alpha和beta是乘法系数，x和y为包含n个元素向量，A为n*n的对称矩阵。

接口定义

C interface：

void cblas_ssymv(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, 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_dsymv(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, 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 SSYMV(UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)

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

参数

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

依赖

#include "kblas.h"

示例

C interface：

    int n = 3, lda = 3; 
    float alpha = 1.0, beta = 1.0; 
    int incx = 1, incy = 1; 
    /* 
     * A =  | 0.0   *    *  | 
     *      | 2.0  3.0   *  | 
     *      | 0.0  4.0  1.0 | 
     */ 
    float a[9] = {0, 2.0, 0, 0, 3.0, 4.0, 0, 0, 1.0}; 
    float x[3] = {5.0, 2.0, 1.0}; 
    float y[3] = {1.0, 1.0, 3.0}; 
 
    cblas_ssymv(CblasColMajor,CblasLower, n, alpha, a, lda, x, incx, beta, y, incy); 
    /* 
     * Output y: 5.000000, 21.000000, 12.000000 
     */

Fortran interface：

      INTEGER :: N=3, LDA=3 
      REAL(4) :: ALPHA=1.0, BETA=1.0 
      INTEGER :: INCX=1, INCY=1 
      REAL(4) :: A(9), X(3), Y(3) 
      DATA A/0, 2.0, 0, 0, 3.0, 4.0, 0, 0, 1.0/ 
      DATA X/5.0, 2.0, 1.0/ 
      DATA Y/1.0, 1.0, 3.0/ 
      EXTERNAL SSYMV 
      CALL SSYMV('L', N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY) 
 
*     Output Y: 5.000000, 21.000000, 12.000000

父主题： KML_BLAS Level 2函数