?syr2

对称矩阵秩2更新

即： $\text{[math]}$ 。alpha为乘法系数，x和y为包含n个元素的向量，A为n*n的对称矩阵。

接口定义

C interface：

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

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

Fortran interface：

CALL SSYR2(UPLO, N, ALPHA, X, INCX, Y, INCY, A, LDA)

CALL DSYR2(UPLO, N, ALPHA, X, INCX, Y, INCY, A, LDA)

参数

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

依赖

#include "kblas.h"

示例

C interface：

    int n = 3, lda = 4; 
    float alpha = 1.0; 
    int incx = 1, incy = 1; 
    float x[3] = {1.0, 3.0, 2.0}; 
    float y[3] = {2.0, 3.0, 2.0}; 
    /** 
     *            | 3.0  1.0   1.0  | 
     *    A    =  |  .   16.0  10.0 | 
     *            |  .    .    1.0  | 
     *            |  .    .      .  | 
     */ 
    float a[12] = {3.0, 0, 0, 0, 1.0, 16.0, 0, 0, 1.0, 10.0, 1.0, 0}; 
 
    cblas_ssyr2(CblasColMajor,CblasUpper, n, alpha, x, incx, y, incy, a, lda); 
    /*** 
     *            | 7.0  10.  7.0 | 
     *            |  .   34.  22.0| 
     * Output A = |  .    .   9.0 | 
     *            |  .    .    .  | 
     */

Fortran interface：

      INTEGER :: N=3, LDA=4 
      INTEGER :: INCX=1, INCY=1 
      REAL(4) :: ALPHA=1.0 
      REAL(4) :: X(3), Y(3), A(12) 
      DATA X/1.0, 3.0, 2.0/ 
      DATA Y/2.0, 3.0, 2.0/ 
      DATA A/3.0, 0, 0, 0, 1.0, 16.0, 0, 0, 1.0, 10.0, 1.0, 0/ 
      EXTERNAL SSYR2 
      CALL SSYR2('U', N, ALPHA, X, INCX, Y, INCY, A, LDA) 
 
*                | 7.0  10.  7.0 | 
*                |  .   34.  22.0| 
*     Output A = |  .    .   9.0 | 
*                |  .    .    .  |

父主题： KML_BLAS Level 2函数