?geru
一般共轭矩阵秩1更新。
即：
。alpha是乘法系数，A是m*n的一般矩阵，x是含有m个元素的向量，y是含有n个元素的向量。
接口定义
C interface：
void cblas_cgeru(const enum CBLAS_ORDER order, const BLASINT M, const BLASINT N, const void *alpha, const void *X, const BLASINT incX, const void *Y, const BLASINT incY, void *A, const BLASINT lda);
void cblas_zgeru(const enum CBLAS_ORDER order, const BLASINT M, const BLASINT N, const void *alpha, const void *X, const BLASINT incX, const void *Y, const BLASINT incY, void *A, const BLASINT lda);
Fortran interface：
CALL CGERU(M, N, ALPHA, X, INCX, Y, INCY, A, LDA)
CALL ZGERU(M, N, ALPHA, X, INCX, Y, INCY, A, LDA)
参数
参数名 |
类型 |
描述 |
输入/输出 |
|---|---|---|---|
order |
枚举类型CBLAS_ORDER |
表示矩阵是行主序或列主序。 |
输入 |
M |
整型数 |
矩阵A的行数。 |
输入 |
N |
整型数 |
矩阵A的列数。 |
输入 |
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：
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 | int m = 4, n = 4, lda = 4; float alpha[2] = {1.0, 2.0}; int incx = 1, incy = 1; /** * X: * 0.340188, -0.105617, 0.283099, 0.298440, 0.411647, -0.302449, -0.164777, 0.268230 * Y: * -0.222225, 0.053970, -0.022603, 0.128871, -0.135216, 0.013401, 0.452230, 0.416195 * A: * 0.135712, 0.217297, -0.343321, -0.099056, 0.112640, -0.203968, 0.026745, 0.269914 * -0.358397, 0.106969, -0.370210, -0.391191, 0.137552, 0.024287, -0.099771, 0.391529 * -0.483699, -0.257113, 0.498924, -0.281743, -0.006417, 0.472775, -0.216685, -0.147542 * -0.362768, 0.304177, 0.012932, 0.339112, -0.207483, 0.271358, 0.307725, 0.419026 */ float x[8] = {0.340188, -0.105617, 0.283099, 0.298440, 0.411647, -0.302449, -0.164777, 0.268230}; float y[8] = {-0.222225, 0.053970, -0.022603, 0.128871, -0.135216, 0.013401, 0.452230, 0.416195}; float a[32] = {0.135712, 0.217297, -0.343321, -0.099056, 0.112640, -0.203968, 0.026745, 0.269914, -0.358397, 0.106969, -0.370210, -0.391191, 0.137552, 0.024287, -0.099771, 0.391529, -0.483699, -0.257113, 0.498924, -0.281743, -0.006417, 0.472775, -0.216685, -0.147542, -0.362768, 0.304177, 0.012932, 0.339112, -0.207483, 0.271358, 0.307725, 0.419026}; cblas_cgeru(CblasColMajor, m, n, alpha, x, incx, y, incy, a, 4); /** * Output A: * -0.017847 0.119331 -0.444931 0.165040 -0.565962 -0.327440 -0.352610 0.793599 * -0.320256 -0.308135 -0.474544 -0.451172 0.529765 -0.402861 -0.488827 0.599533 * -0.141372 -0.264850 0.047453 0.143518 -0.150850 0.415971 0.035456 0.929981 * 0.185887 0.245696 -0.076018 0.302546 -0.121045 -0.148647 0.016127 0.099442 */ |
Fortran interface：
INTEGER :: M=4
INTEGER :: N=4
INTEGER :: LDA=4
INTEGER :: INCX=1
INTEGER :: INCY=1
COMPLEX(4) :: ALPHA=(1.0, 2.0)
COMPLEX(4) :: A(4, 4)
DATA A/(0.135712, 0.217297), (-0.343321, -0.099056),
$ (0.112640, -0.203968), (0.026745, 0.269914),
$ (-0.358397, 0.106969), (-0.370210, -0.391191),
$ (0.137552, 0.024287), (-0.099771, 0.391529),
$ (-0.483699, -0.257113), (0.498924, -0.281743),
$ (-0.006417, 0.472775), (-0.216685, -0.147542),
$ (-0.362768, 0.304177), (0.012932, 0.339112),
$ (-0.207483, 0.271358), (0.307725, 0.419026)/
COMPLEX(4) :: X(4)
DATA X/(0.340188, -0.105617), (0.283099, 0.298440),
$ (0.411647, -0.302449), (-0.164777, 0.268230)/
COMPLEX(4) :: Y(4)
DATA Y/(-0.222225, 0.053970), (-0.022603, 0.128871),
$ (-0.135216, 0.013401), (0.452230, 0.416195)/
EXTERNAL CGERU
CALL CGERU(M, N, ALPHA, X, INCX, Y, INCY, A, LDA)
* Output A:
* -0.017847 0.119331 -0.444931 0.165040 -0.565962 -0.327440 -0.352610 0.793599
* -0.320256 -0.308135 -0.474544 -0.451172 0.529765 -0.402861 -0.488827 0.599533
* -0.141372 -0.264850 0.047453 0.143518 -0.150850 0.415971 0.035456 0.929981
* 0.185887 0.245696 -0.076018 0.302546 -0.121045 -0.148647 0.016127 0.099442