?hpmv

向量与三角展开的埃尔米特矩阵乘积。

即 $\text{[math]}$ 。alpha和beta是缩放系数，x和y是含有n个元素的向量，A是三角展开的埃尔米特矩阵。

接口定义

C interface：

void cblas_chpmv(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, const void *alpha, const void *Ap, const void *X, const BLASINT incX, const void *beta, void *Y, const BLASINT incY);

void cblas_zhpmv(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, const void *alpha, const void *Ap, const void *X, const BLASINT incX, const void *beta, void *Y, const BLASINT incY);

Fortran interface：

CALL CHPMV(UPLO, N, ALPHA, AP, X, INCX, BETA, Y, INCY)

CALL ZHPMV(UPLO, N, ALPHA, AP, X, INCX, BETA, Y, INCY)

参数

参数名	类型	描述	输入/输出
order	枚举类型CBLAS_ORDER	表示矩阵是行主序或列主序。	输入
Uplo	枚举类型CBLAS_UPLO	表示矩阵A的存储展开方式（上三角或下三角）。如果Uplo=CblasUpper，则使用A的上三角展开。如果Uplo=CblasLower，则使用A的下三角展开。	输入
N	整型数	表示矩阵A的阶数，要求N大于或等于零。	输入
alpha	在chpmv中是单精度复数类型。在zhpmv中是双精度复数类型。	表示系数。	输入
Ap	在chpmv中是单精度复数类型。在zhpmv中是双精度复数类型。	压缩的埃尔米特矩阵，矩阵规模至少是(n*(n+1)/2)。	输入
X	在chpmv中是单精度复数类型。在zhpmv中是双精度复数类型。	向量X，向量规模至少是(1+(N-1)*abs(incX))。	输入
incX	整型数	表示X中向量增长步长，不能为零。	输入
beta	在chpmv中是单精度复数类型。在zhpmv中是双精度复数类型。	乘法系数。	输入
Y	在chpmv中是单精度复数类型。在zhpmv中是双精度复数类型。	向量Y，向量规模至少是(1+(N-1)*abs(incY))。	输入/输出
incY	整型数	表示Y中向量增长步长，不能为零。	输入

依赖

#include "kblas.h"

示例

C interface：

    int n = 3; 
    float alpha[2] = {1.0, 0}, beta[2] = {1.0, 0}; 
    int incx = 1, incy = 1; 
   /*             |  (1.0, 0)  (3.0, 3.0)  (3.0, -3.0) | 
    *       A =   | (3.0, -3.0)  (2.0, 0.0)  (3.0, -1.0) | 
    *             |  (3.0, 3.0)  (3.0, 1.0)   (3.0, 0.0) | 
    */ 
    float ap[12] = {1.0, 0, 3.0, -3.0, 3.0, 3.0, 2.0, 0, 3.0, 1.0, 3.0, 0}; 
    float x[6] = {3.0, 2.0, 1.0, 2.0, -1.0, 1.0}; 
    float y[6] = {1.0, 2.0, 3.0, 2.0, -1.0, 1.0}; 
 
    cblas_chpmv(CblasColMajor,CblasLower, n, alpha, ap, x, incx, beta, y, incy); 
   /* 
    * Output y = |1.000000, 19.000000, 18.000000, 7.000000,  0.000000, 26.000000| 
    */

Fortran interface：

      INTEGER :: N=3 
      COMPLEX(4) :: ALPHA=(1.0, 0) 
      COMPLEX(4) :: BETA=(1.0, 0) 
      INTEGER :: INCX=1 
      INTEGER :: INCY=1 
      COMPLEX(4) :: AP(6) 
      DATA AP/(1.0, 0), (3.0, -3.0), (3.0, 3.0), (2.0, 0), (3.0, 1.0), (3.0, 0)/ 
      COMPLEX(4) :: X(3), Y(3) 
      DATA X/(3.0, 2.0), (1.0, 2.0), (-1.0, 1.0)/ 
      DATA Y/(1.0, 2.0), (3.0, 2.0), (-1.0, 1.0)/ 
      EXTERNAL CHPMV 
      CALL CHPMV('L', N, ALPHA, AP, X, INCX, BETA, Y, INCY) 
 
*     Output Y = |(1.000000, 19.000000), (18.000000, 7.000000),  (0.000000, 26.000000)|

父主题： KML_BLAS Level 2函数