?hemv
向量与埃尔米特矩阵乘积。
即:
。alpha和beta是乘法系数,x和y是含有n个元素的向量,A是n阶埃尔米特矩阵。
接口定义
C interface:
void cblas_chemv(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, const void *alpha, const void *A, const BLASINT lda, const void *X, const BLASINT incX, const void *beta, void *Y, const BLASINT incY);
void cblas_zhemv(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, const void *alpha, const void *A, const BLASINT lda, const void *X, const BLASINT incX, const void *beta, void *Y, const BLASINT incY);
Fortran interface:
CALL CHEMV(UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)
CALL ZHEMV(UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)
参数
参数名  | 
类型  | 
描述  | 
输入/输出  | 
|---|---|---|---|
order  | 
枚举类型CBLAS_ORDER  | 
表示矩阵是行主序或列主序。  | 
输入  | 
Uplo  | 
枚举类型CBLAS_UPLO  | 
表示使用矩阵A的上三角或下三角, 
  | 
输入  | 
N  | 
整型数  | 
表示矩阵A的阶数,要求N大于或等于零。  | 
输入  | 
alpha  | 
  | 
表示系数。  | 
输入  | 
A  | 
  | 
埃尔米特矩阵A(lda, n)。  | 
输入  | 
lda  | 
整型数  | 
矩阵A中主维度长度,要求lda大于或等于(k+1)。  | 
输入  | 
X  | 
  | 
向量X,向量规模至少是(1+(N-1)*abs(incX))。  | 
输入  | 
incX  | 
整型数  | 
表示X中向量增长步长,不能为零。  | 
输入  | 
beta  | 
  | 
乘法系数。  | 
输入  | 
Y  | 
  | 
向量Y,向量规模至少是(1+(N-1)*abs(incY))。  | 
输入/输出  | 
incY  | 
整型数  | 
表示Y中向量增长步长,不能为零。  | 
输入  | 
依赖
#include "kblas.h"
示例
C interface:
    int n = 3, lda = 3; 
    float alpha[2] = {1.0, 0}, beta[2] = {1.0, 0}; 
    int incx = 1, incy = 1; 
    /*            |  (21.0,  * )      *           *      | 
     *      A =   | (1.0, -1.0)  (3.0,  * )       *      | 
     *            |  (1.0, 3.0)  (2.0, 2.0)   (2.0, *  ) | 
     */ 
    float a[18] = {21.0, 0, 1.0, -1.0, 1.0, 3.0, 0, 0, 3.0, 0, 2.0, 2.0, 0, 0, 0, 0, 2.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_chemv(CblasColMajor,CblasLower, n, alpha, a, lda, x, incx, beta, y, incy); 
    /* 
     * Output y = | 65.0, 51.0, 11.0, 11.0, -8.0, 20.0| 
     */
Fortran interface:
      INTEGER :: N=3 
      INTEGER :: LDA=3 
      COMPLEX(4) :: ALPHA=(1.0, 0) 
      COMPLEX(4) :: BETA=(1.0, 0) 
      INTEGER :: INCX=1 
      INTEGER :: INCY=1 
      COMPLEX(4) :: A(3, 3) 
      DATA A/(21.0, 0), (1.0, -1.0), (1.0, 3.0), 
     $       (0, 0), (3.0, 0), (2.0, 2.0), 
     $       (0, 0), (0, 0), (2.0, 0)/ 
 
      COMPLEX(4) :: X(3) 
      DATA X/(3.0, 2.0), (1.0, 2.0), (-1.0, 1.0)/ 
      COMPLEX(4) :: Y(3) 
      DATA Y/(1.0, 2.0), (3.0, 2.0), (-1.0, 1.0)/ 
      EXTERNAL CHEMV 
      CALL CHEMV('L', N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY) 
 
*     Output Y : (65.0, 51.0), (11.0, 11.0), (-8.0, 20.0)