?spr2
对称展开矩阵秩2更新。
即：
。alpha为乘法系数，x和y为包含n个元素的向量，A为n*n的三角展开对称矩阵。
接口定义
C interface：
void cblas_sspr2(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, const float alpha, const float *X, const BLASINT incX, const float *Y, const BLASINT incY, float *A);
void cblas_dspr2(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, const double alpha, const double *X, const BLASINT incX, const double *Y, const BLASINT incY, double *A);
Fortran interface：
CALL SSPR2(UPLO, N, ALPHA, X, INCX, Y, INCY, AP)
CALL DSPR2(UPLO, N, ALPHA, X, INCX, Y, INCY, AP)
参数
参数名 |
类型 |
描述 |
输入/输出 |
|---|---|---|---|
order |
枚举类型CBLAS_ORDER |
表示矩阵是行主序或列主序。 |
输入 |
Uplo |
枚举类型CBLAS_UPLO |
对称矩阵展开存储方式（上三角或下三角展开）。
|
输入 |
N |
整型数 |
向量X的元素个数。 |
输入 |
alpha |
|
乘法系数。 |
输入 |
X |
|
矩阵X，长度至少 1+(n-1)*abs(incX)。 |
输入 |
incX |
整型数 |
向量X的增长步长，不能为零。 |
输入 |
Y |
|
矩阵Y，长度至少 1+(n-1)*abs(incY)。 |
输入 |
incY |
整型数 |
向量Y的增长步长，不能为零。 |
输入 |
A |
|
矩阵A三角存储。 |
输入/输出 |
依赖
#include "kblas.h"
示例
C interface：
float alpha = 1.0;
int n = 3, incx = 1, incy = 1;
float x[3] = {2.0, 2.0, 1.0};
float y[3] = {1.0, 2.0, 1.0};
/**
* | 3.0 1.0 2.0 |
* A = | 1.0 6.0 3.0 |
* | 2.0 3.0 3.0 |
*/
float a[6] = {3.0, 1.0, 2.0, 6.0, 3.0, 3.0};
cblas_sspr2(CblasColMajor,CblasLower, n, alpha, x, incx, y, incy, a);
/**
* Output a = |7.0, 7.0, 5.0, 14.0, 7.0, 5.0|
*/
Fortran interface：
REAL(4) :: ALPHA=1.0
INTEGER :: N=3
INTEGER :: INCX=1
INTEGER :: INCY=1
REAL(4) :: X(3), Y(3), A(6)
DATA X/2.0, 2.0, 1.0/
DATA Y/1.0, 2.0, 1.0/
DATA A/3.0, 1.0, 2.0, 6.0, 3.0, 3.0/
EXTERNAL SSPR2
CALL SSPR2('L', N, ALPHA, X, INCX, Y, INCX, A)
* Output A = |7.0, 7.0, 5.0, 14.0, 7.0, 5.0|