?spr2

对称展开矩阵秩2更新， $\text{[math]}$ 。

alpha为乘法系数，x和y为包含n个元素的向量，A为n*n的三角展开对称矩阵。

接口定义

C interface:

void cblas_sspr2(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);

void cblas_dspr2(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);

Fortran interface:

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

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

参数

参数名	类型	描述	输入/输出
order	枚举类型CBLAS_ORDER	表示矩阵是行主序或列主序。	输入
Uplo	枚举类型CBLAS_UPLO	对称矩阵展开存储方式（上三角或下三角展开）。如果Uplo=CblasUpper，则使用A的上三角展开。如果Uplo=CblasLower，则使用A的下三角展开。	输入
N	整型数	向量X的元素个数。	输入
alpha	在sspr2中，alpha为单精度复数类型。在dspr2中，alpha为双精度复数类型。	乘法系数。	输入
X	在sspr2中，X为单精度复数类型。在dspr2中，X为双精度复数类型。	矩阵X，长度至少1+(n-1)*abs(incX)。	输入
incX	整型数	向量X的增长步长，不能为零。	输入
Y	在sspr2中，Y为单精度复数类型。在dspr2中，Y为双精度复数类型。	矩阵Y，长度至少1+(n-1)*abs(incY)。	输入
incY	整型数	向量Y的增长步长，不能为零。	输入
A	在sspr2中，A为单精度复数类型。在dspr2中，A为双精度复数类型。	矩阵A三角存储。	输入/输出

依赖

#include "kblas.h"

示例

C interface:

    float alpha = 1.0; 
    int n = 3, incx = 1, incy = 1; 
    float x[3] = {2.0, 2.0, 1.0}; 
    float y[3] = {1.0, 2.0, 1.0}; 
    /** 
     *             | 3.0  1.0  2.0 | 
     *         A = | 1.0  6.0  3.0 | 
     *             | 2.0  3.0  3.0 | 
     */ 
    float a[6] = {3.0, 1.0, 2.0, 6.0, 3.0, 3.0}; 
 
    cblas_sspr2(CblasColMajor,CblasLower, n, alpha, x, incx, y, incy, a); 
    /** 
     *  Output a = |7.0, 7.0, 5.0, 14.0, 7.0, 5.0| 
     */

Fortran interface:

      REAL(4) :: ALPHA=1.0 
      INTEGER :: N=3 
      INTEGER :: INCX=1 
      INTEGER :: INCY=1 
      REAL(4) :: X(3), Y(3), A(6) 
      DATA X/2.0, 2.0, 1.0/ 
      DATA Y/1.0, 2.0, 1.0/ 
      DATA A/3.0, 1.0, 2.0, 6.0, 3.0, 3.0/ 
      EXTERNAL SSPR2 
      CALL SSPR2('L', N, ALPHA, X, INCX, Y, INCX, A) 
 
*     Output A = |7.0, 7.0, 5.0, 14.0, 7.0, 5.0|

父主题： KML_BLAS Level 2函数