?her

复数埃尔米特矩阵秩1更新。

即： $\text{[math]}$ 。A是n*n的埃尔米特矩阵，alpha是乘法系数，x是有n个元素的向量。

接口定义

C interface：

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

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

Fortran interface：

CALL CHER(UPLO, N, ALPHA, X, INCX, A, LDA)

CALL ZHER(UPLO, N, ALPHA, X, INCX, A, LDA)

参数

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

依赖

#include "kblas.h"

示例

C interface：

    int n = 3; 
    float alpha = 2.9; 
    int incx = 1; 
    /** 
     *              |  (1.0, 0.0)  (13.0, 5.0)  (12.0, 0) | 
     *         A =  | (13.0, -5.0)  (2.0, 0.0)  (4.0, -18.0) | 
     *              |  (12.0, 0)  (4.0, 18.0)   (3.0, 0.0) | 
     */ 
    float ap[18] = {1.0, 0, 13.0, -5.0, 12.0, 0, 13.0, 5.0, 2.0, 0, 4.0, 18.0, 12.0, 0, 4.0, -18.0, 3.0, 0}; 
    float x[6] = {2.0, -2.0, 14.0, 1, 3.0, 4.0}; 
 
    cblas_cher(CblasColMajor, CblasLower, n, alpha, x, incx, ap, 3); 
    /** 
     *             | (24.200001, 0.000000)  (13.000000   5.000000)     (12.000000  0.000000)  | 
     * Output  A = | (88.400002, 82.000000) (573.300049  0.000000)     (4.000000   -18.000000)| 
     *             | (6.200001, 40.600002)  (137.400009  171.700012)   (75.500000   0.000000) |

     */

Fortran interface：

      INTEGER :: N=3 
      REAL(4) :: ALPHA=2.9 
      INTEGER :: INCX=1 
      COMPLEX(4) :: AP(3, 3) 
      DATA AP/(1.0, 0), (13.0, -5.0), (12.0, 0), (13.0, 5.0), (2.0, 0), (4.0, 18.0), (12.0, 0), (4.0, -18.0), (3.0, 0)/ 
      COMPLEX(4) :: X(3) 
      DATA X/(2.0, -2.0), (14.0, 1), (3.0, 4.0)/ 
      EXTERNAL CHER 
      CALL CHER('L', N, ALPHA, X, INCX, AP, LDA) 
 
*                 | (24.200001, 0.000000)  (13.000000   5.000000)     (12.000000  0.000000)  | 
*     Output  A = | (88.400002, 82.000000) (573.300049  0.000000)     (4.000000   -18.000000)| 
*                 | (6.200001, 40.600002)  (137.400009  171.700012)   (75.500000   0.000000) |

父主题： KML_BLAS Level 2函数