?ger
一般矩阵秩1更新。
即：
。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的行数，要求M大于或等于零。 |
输入 |
N |
整型数 |
表示矩阵A的列数，要求N大于或等于零。 |
输入 |
alpha |
|
乘法系数。 |
输入 |
X |
|
矩阵X，向量规模至少是(1+(M-1)*abs(incx))。 |
输入 |
incX |
整型数 |
向量X的增长步长，不能为零。 |
输入 |
Y |
|
矩阵Y，向量规模至少是(1+(N-1)*abs(incy))。 |
输入 |
incY |
整型数 |
向量Y的增长步长，不能为零。 |
输入 |
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