---
title: ?hpr
description: "三角展开的埃尔米特矩阵秩1更新。"
url: https://www.hikunpeng.com/document/detail/zh/kunpenghpcs/hpckit/devg/kunpengaccel_kml_0630.html
sourcePath: /source/zh/kunpenghpcs/hpckit/devg/kunpengaccel_kml_0630.html
indexId: a44d4e1cabb5b875d8c8b18b3b7cb2b5a570ad835fcc89f6ec94e3e7b78d377a61
---
# ?hpr

三角展开的埃尔米特矩阵秩1更新。

即：。alpha是缩放系数，x是含有n个元素的向量，A是三角展开的埃尔米特矩阵。

#### 接口定义

C interface：

void cblas_chpr(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, const float alpha, const void *X, const BLASINT incX, void *A);

void cblas_zhpr(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, const double alpha, const void *X,const BLASINT incX, void *A);

Fortran interface：

CALL CHPR(UPLO, N, ALPHA, X, INCX, AP)

CALL ZHPR(UPLO, N, ALPHA, X, INCX, AP)


#### 参数

| 参数名 | 类型 | 描述 | 输入/输出 |
| --- | --- | --- | --- |
| order | 枚举类型CBLAS\_ORDER | 表示矩阵是行主序或列主序。 | 输入 |
| Uplo | 枚举类型CBLAS\_UPLO | 埃尔米特矩阵展开存储方式（上三角或下三角展开）。 如果Uplo=CblasUpper，则使用A的上三角展开。 如果Uplo=CblasLower，则使用A的下三角展开。 | 输入 |
| N | 整型数 | 向量X的元素个数。 | 输入 |
| alpha | 在chpr中是单精度复数类型。 在zhpr中是双精度复数类型。 | 乘法系数。 | 输入 |
| X | 在chpr中是单精度复数类型。 在zhpr中是双精度复数类型。 | 矩阵X，长度至少 1+(n\-1)\*abs(incX)。 | 输入 |
| incX | 整型数 | 向量X的增长步长，不能为零。 | 输入 |
| A | 在chpr中是单精度复数类型。 在zhpr中是双精度复数类型。 | 矩阵A，矩阵规模至少是(n\*(n+1)/2)。 | 输入/输出 |


#### 依赖

#include "kblas.h"


#### 示例

C interface：

```
int n = 3; 
    float alpha = 1.0; 
    int incx = 1; 
    /** 
     *              |  (2.0, 0.0)  (13.0, 5.0)  (12.0, -13.0) | 
     *         A =  | (13.0, -5.0)  (3.0, 0.0)  (4.0, -18.0)  | 
     *              |  (12.0, 13.0)  (4.0, 18.0)   (2.0, 0.0) | 
     */ 
    float ap[12] = {2.0, 0, 13.0, -5.0, 12.0, 13.0, 3.0, 0, 4.0, 18.0, 2.0, 0}; 
    float x[6] = {2.0, -2.0, 14.0, 1, 3.0, 4.0}; 
 
    cblas_chpr(CblasColMajor, CblasLower, n, alpha, x, incx, ap); 
    /** 
     * Output ap = |(10.0, 0.0), (39.0, 25.0), (10.0, 27.0), (200.0, 0.0), (50.0, 71.0), (27.0, 0.0) | 
     */
```

Fortran interface：

```
INTEGER :: N=3
REAL(4) :: ALPHA=1.0
INTEGER :: INCX=1
COMPLEX(4) :: AP(6), X(3)
DATA AP/(2.0, 0), (13.0, -5.0), (12.0, 13.0), (3.0, 0),
$        (4.0, 18.0), (2.0, 0)/
DATA X/(2.0, -2.0), (14.0, 1), (3.0, 4.0)/
EXTERNAL CHPR
CALL CHPR('L', N, ALPHA, X, INCX, AP)
*     Output AP = | (10.0, 0.0), (39.0, 25.0), (10.0, 27.0), (200.0, 0.0), (50.0, 71.0), (27.0, 0.0) |
```