---
title: ?hpr2
description: "三角展开的埃尔米特矩阵秩2更新。"
url: https://www.hikunpeng.com/document/detail/zh/kunpenghpcs/hpckit/devg/kunpengaccel_kml_0631.html
sourcePath: /source/zh/kunpenghpcs/hpckit/devg/kunpengaccel_kml_0631.html
indexId: c195604af2b30588d323183c7481ceed40487617a22b6ee7821bf88008867b4761
---
# ?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 | 埃尔米特矩阵展开存储方式（上三角或下三角展开） 如果Uplo=CblasUpper，则使用A的上三角展开。 如果Uplo=CblasLower，则使用A的下三角展开。 | 输入 |
| N | 整型数 | 向量X的元素个数。 | 输入 |
| alpha | 在chpr2中是单精度复数类型。 在zhpr2中是双精度复数类型。 | 乘法系数。 | 输入 |
| X | 在chpr2中是单精度复数类型。 在zhpr2中是双精度复数类型。 | 矩阵X，长度至少 1+(n\-1)\*abs(incX)。 | 输入 |
| incX | 整型数 | 向量X的增长步长，不能为零。 | 输入 |
| Y | 在cher2中是单精度复数类型。 在zher2中是双精度复数类型。 | 矩阵Y，长度至少1+(n\-1)\*abs(incY)。 | 输入 |
| incY | 整型数 | 向量Y的增长步长，不能为零。 | 输入 |
| Ap | 在chpr2中是单精度复数类型。 在zhpr2中是双精度复数类型。 | 矩阵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)|
```