?symv
向量与对称矩阵乘积
即
。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的上三角或下三角。
|
输入 |
N |
整型数 |
表示矩阵A的阶数，要求N大于或等于零。 |
输入 |
alpha |
|
表示系数。 |
输入 |
A |
|
对称矩阵A(lda, n)。 |
输入 |
lda |
整型数 |
矩阵A中每列存储的地址距离，要求lda大于或等于max(1, n)。 |
输入 |
X |
|
向量X，向量规模至少是(1+(N-1)*abs(incX))。 |
输入 |
incX |
整型数 |
表示X中向量增长步长，不能为零。 |
输入 |
beta |
|
乘法系数。 |
输入 |
Y |
|
向量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