?syr
对称矩阵秩1更新
即：
。alpha为乘法系数，x为包含n个元素的向量，A为n*n的对称矩阵。
接口定义
C interface：
void cblas_ssyr(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, const float alpha, const float *X, const BLASINT incX, float *A, const BLASINT lda);
void cblas_dsyr(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, const double alpha, const double *X, const BLASINT incX, double *A, const BLASINT lda);
Fortran interface：
CALL SSYR(UPLO, N, ALPHA, X, INCX, A, LDA)
CALL DSYR(UPLO, N, ALPHA, X, INCX, A, LDA)
参数
参数名 |
类型 |
描述 |
输入/输出 |
|---|---|---|---|
order |
枚举类型CBLAS_ORDER |
表示矩阵是行主序或列主序。 |
输入 |
Uplo |
枚举类型CBLAS_UPLO |
使用矩阵A的上三角或下三角。
|
输入 |
N |
整型数 |
向量X的元素个数。 |
输入 |
alpha |
|
乘法系数。 |
输入 |
X |
|
矩阵X，长度至少1+(n-1)*abs(incX)。 |
输入 |
incX |
整型数 |
向量X的增长步长，不能为零。 |
输入 |
A |
|
矩阵A(lda, n)。 |
输入/输出 |
lda |
整型数 |
矩阵A的主维度长度，要求lda大于或等于max(1, n)。 |
输入 |
依赖
#include "kblas.h"
示例
C interface：
int n = 3, lda = 4, incx = 1;
float alpha = 1.0;
float x[3] = {1.0, 3.0, 2.0};
/**
* | 3.0 14.0 12.0 |
* A = | . 16.0 17.0 |
* | . . 13.0 |
* | . . . |
*/
float a[12] = {3.0, 0, 0, 0, 14.0, 16.0, 0, 0, 12, 17.0, 13.0, 0};
cblas_ssyr(CblasColMajor,CblasUpper, n, alpha, x, incx, a, lda);
/**
* | 4.0 17.0 14.0 |
* Output A = | . 25.0 23.0 |
* | . . 17.0 |
* | . . . |
*/
Fortran interface：
INTEGER :: N=3, LDA=4
INTEGER :: INCX=1
REAL(4) :: ALPHA=1.0
REAL(4) :: X(3), A(12)
DATA X/1.0, 3.0, 2.0/
DATA A/3.0, 0, 0, 0, 14.0, 16.0, 0, 0, 12, 17.0, 13.0, 0/
EXTERNAL SSYR
CALL SSYR('U', N, ALPHA, X, INCX, A, LDA)
* | 4.0 17.0 14.0 |
* Output A = | . 25.0 23.0 |
* | . . 17.0 |
* | . . . |