?ger

一般矩阵秩1更新， $\text{[math]}$ 。

alpha是乘法系数，A是m*n的一般矩阵，x是含有m个元素的向量，y是含有n个元素的向量。

接口定义

C interface:

void cblas_sger (const enum CBLAS_ORDER order, const BLASINT M, 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_dger (const enum CBLAS_ORDER order, const BLASINT M, 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 SGER(M, N, ALPHA, X, INCX, Y, INCY, A, LDA)

CALL DGER(M, N, ALPHA, X, INCX, Y, INCY, A, LDA)

参数

参数名	类型	描述	输入/输出
order	枚举类型CBLAS_ORDER	表示矩阵是行主序或列主序。	输入
M	整型数	矩阵A的行数。	输入
N	整型数	矩阵A的列数。	输入
alpha	在sger中，alpha为单精度浮点类型。在dger中，alpha为双精度浮点类型类型。	乘法系数。	输入
X	在sger中，X为单精度浮点类型。在dger中，X为双精度浮点类型类型。	矩阵X，向量规模至少是(1+(m-1)*abs(incx))。	输入
incX	整型数	向量X的增长步长，不能为零。	输入
Y	在sger中，Y为单精度浮点类型。在dger中，Y为双精度浮点类型类型。	矩阵Y，向量规模至少是(1+(n-1)*abs(incy))。	输入
incY	整型数	向量Y的增长步长，不能为零。	输入
A	在sger中，A为单精度浮点类型。在dger中，A为双精度浮点类型。	矩阵A(lda, n)。	输出
lda	整型数	矩阵A中主维度长度，如果矩阵A为列存，要求lda大于或等于max(1, m)，否则要求大于或等于max(1,n)。	输入

依赖

#include "kblas.h"

示例

C interface:

    int m = 4, n = 4, lda =4; 
    float alpha = 1.0; 
    int incx = 1, incy = 1; 
    /** 
     * X: 
     *    0.340188, -0.105617, 0.283099, 0.298440 
     * Y: 
     *    0.411647, -0.302449, -0.164777, 0.268230 
     * 
     * Input A: 
     *    -0.222225       -0.135216       0.135712        -0.483699 
     *    0.053970        0.013401        0.217297        -0.257113 
     *    -0.022603       0.452230        -0.358397       -0.362768 
     *    0.128871        0.416195        0.106969        0.304177 
     */ 
    float x[4] = {0.340188, -0.105617, 0.283099, 0.298440}; 
    float y[4] = {0.411647, -0.302449, -0.164777, 0.268230}; 
    float a[16] = {-0.222225, -0.135216, 0.135712, -0.483699, 
                   0.053970 , 0.013401, 0.217297 , -0.257113, 
                   -0.022603, 0.452230, -0.358397, -0.362768, 
                   0.128871 , 0.416195, 0.106969, 0.304177}; 
 
    cblas_sger(CblasColMajor, m, n, alpha, x, incx, y, incy, a, lda); 
    /** 
     * Output A: 
     *    -0.082188       -0.048920       -0.078658        0.220120
     *    -0.178693        0.045345        0.469633        0.387865 
     *    0.252249         0.131674       -0.405045        0.182905
     *    -0.360847       -0.347376       -0.411944        0.384228
     */

Fortran interface:

      INTEGER :: M=4 
      INTEGER :: N=4 
      INTEGER :: LDA=4 
      INTEGER :: INCX=1 
      INTEGER :: INCY=1 
      REAL(4) :: ALPHA=1.0 
      REAL(4) :: A(4, 4) 
      DATA A/-0.222225, -0.135216, 0.135712, -0.483699, 
     $       0.053970 , 0.013401, 0.217297 , -0.257113, 
     $       -0.022603, 0.452230, -0.358397, -0.362768, 
     $       0.128871 , 0.416195, 0.106969, 0.304177/ 
      REAL(4) :: X(4) 
      DATA X/0.340188, -0.105617, 0.283099, 0.298440/ 
      REAL(4) :: Y(4) 
      DATA Y/0.411647, -0.302449, -0.164777, 0.268230/ 
      EXTERNAL SGER 
      CALL SGER(M, N, ALPHA, X, INCX, Y, INCY, A, LDA) 
*     Output A: 
*     -0.082188       -0.048920       -0.078658        0.220120
*     -0.178693        0.045345        0.469633        0.387865 
*      0.252249        0.131674       -0.405045        0.182905
*     -0.360847       -0.347376       -0.411944        0.384228

父主题： KML_BLAS Level 2函数