?hpr

三角展开的埃尔米特矩阵秩1更新， $\text{[math]}$ 。

其中，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中，alpha为单精度复数类型。在zhpr中，alpha为双精度复数类型。	乘法系数。	输入
X	在chpr中，X为单精度复数类型。在zhpr中，X为双精度复数类型。	矩阵X，长度至少 1+(n-1)*abs(incX)。	输入
incX	整型数	向量X的增长步长，不能为零。	输入
A	在chpr中，A为单精度复数类型。在zhpr中，A为双精度复数类型。	矩阵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) |

父主题： KML_BLAS Level 2函数