?hpr2
三角展开的埃尔米特矩阵秩2更新。
即：
。alpha是乘法系数，x和y是含有n个元素的向量，A是n*n的埃尔米特矩阵。
接口定义
C interface：
void cblas_chpr2(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, const void *alpha, const void *X, const BLASINT incX, const void *Y, const BLASINT incY, void *Ap);
void cblas_zhpr2(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, const void *alpha, const void *X, const BLASINT incX, const void *Y, const BLASINT incY, void *Ap);
Fortran interface：
CALL CHPR2(UPLO, N, ALPHA, X, INCX, Y, INCY, AP)
CALL ZHPR2(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的增长步长，不能为零。 |
输入 |
Ap |
|
矩阵A三角存储，矩阵规模至少是(n*(n+1)/2)。 |
输入/输出 |
依赖
#include "kblas.h"
示例
C interface：
int n = 3;
float alpha[2] = {1.0, 2.9};
int incx = 1, incy = 1;
/**
* | (2.0, 0.0) (13.0, 5.0) (12.0, -13.0) |
* A = | (13.0, -5.0) (2.0, 0.0) (4.0, -18.0) |
* | (12.0, 13.0) (4.0, 18.0) (1.0, 0.0) |
*/
float ap[12] = {2.0, 0, 13.0, -5.0, 12.0, 13.0, 2.0, 0, 4.0, 18.0, 1.0, 0};
float x[6] = {2.0, -2.0, 14.0, 1, 3.0, 4.0};
float y[6] = {2.0, -2.0, 14.0, 1, 3.0, 4.0};
cblas_chpr2(CblasColMajor, CblasLower, n, alpha, x, incx, y, incy, ap);
/**
* Output ap = |(18.0, 0.0), (65.0, 55.0), (8.0, 41.0), (396.0, 0.0), (96.0, 124.0), (51.0, 0.0)|
*/
Fortran interface：
INTEGER :: N=3
COMPLEX(4) :: ALPHA=(1.0, 2.9)
INTEGER :: INCX=1, INCY=1
COMPLEX(4) :: AP(6), X(3), Y(3)
DATA AP/(2.0, 0), (13.0, -5.0), (12.0, 13.0), (2.0, 0),
$ (4.0, 18.0), (1.0, 0)/
DATA X/(2.0, -2.0), (14.0, 1), (3.0, 4.0)/
DATA Y/(2.0, -2.0), (14.0, 1), (3.0, 4.0)/
EXTERNAL CHPR2
CALL CHPR2('L', N, ALPHA, X, INCX, Y, INCY, AP)
* Output AP = |(18.0, 0.0), (65.0, 55.0), (8.0, 41.0), (396.0, 0.0), (96.0, 124.0), (51.0, 0.0)|