?symm
实数对称矩阵乘矩阵。
即：
或 
。
alpha，beta为乘法系数，A为对称矩阵，B，C为m*n一般矩阵。
接口定义
C interface：
void cblas_ssymm(const enum CBLAS_ORDER Order, const enum CBLAS_SIDE Side, const enum CBLAS_UPLO Uplo, const BLASINT M, const BLASINT N, const float alpha, const float *A, const BLASINT lda, const float *B, const BLASINT ldb, const float beta, float *C, const BLASINT ldc);
void cblas_dsymm(const enum CBLAS_ORDER Order, const enum CBLAS_SIDE Side, const enum CBLAS_UPLO Uplo, const BLASINT M, const BLASINT N, const double alpha, const double *A, const BLASINT lda, const double *B, const BLASINT ldb, const double beta, double *C, const BLASINT ldc);
void cblas_csymm(const enum CBLAS_ORDER Order, const enum CBLAS_SIDE Side, const enum CBLAS_UPLO Uplo, const BLASINT M, const BLASINT N, const void *alpha, const void *A, const BLASINT lda, const void *B, const BLASINT ldb, const void *beta, void *C, const BLASINT ldc);
void cblas_zsymm(const enum CBLAS_ORDER Order, const enum CBLAS_SIDE Side, const enum CBLAS_UPLO Uplo, const BLASINT M, const BLASINT N, const void *alpha, const void *A, const BLASINT lda, const void *B, const BLASINT ldb, const void *beta, void *C, const BLASINT ldc);
Fortran interface：
CALL SSYMM(SIDE, UPLO, M, N, ALPHA, A, LDA, B, LDB, BETA, C, LDC)
CALL DSYMM(SIDE, UPLO, M, N, ALPHA, A, LDA, B, LDB, BETA, C, LDC)
CALL CSYMM(SIDE, UPLO, M, N, ALPHA, A, LDA, B, LDB, BETA, C, LDC)
CALL ZSYMM(SIDE, UPLO, M, N, ALPHA, A, LDA, B, LDB, BETA, C, LDC)
参数
参数名 |
类型 |
描述 |
输入/输出 |
|---|---|---|---|
order |
枚举类型CBLAS_ORDER |
表示矩阵是行主序或列主序。 |
输入 |
Side |
枚举类型CBLAS_SIDE |
对称矩阵A在矩阵B的左边还是右边。
|
输入 |
Uplo |
枚举类型CBLAS_UPLO |
表示使用矩阵A的上三角或下三角。
|
输入 |
M |
整型数 |
矩阵C的行
|
输入 |
N |
整型数 |
矩阵C的列
|
输入 |
alpha |
|
乘法系数。 |
输入 |
A |
|
矩阵A(lda, ka)。 如果Side = CblasLeft，ka = m，否则ka = n。 |
输入 |
lda |
整型数 |
如果Side = CblasLeft，lda至少max(1, m)，否则至少max(1, n)。 |
输入 |
B |
|
矩阵B(ldb, n)。 |
输入 |
ldb |
整型数 |
|
输入 |
beta |
|
乘法系数。 |
输入 |
C |
|
矩阵C。 |
输入/输出 |
ldc |
整型数 |
|
输入 |
依赖
#include "kblas.h"
C interface：
int m = 3, n = 6, lda = 4, ldb =3, ldc = 5;
float alpha = 2.0, beta = 2.0;
/**
* | 4.0 . . |
* A(4 * 3) = | 3.0 4.0 . |
* | 6.0 -7.0 -1.0 |
* | . . . |
*
* | -1.0 -3.0 2.0 2.0 -1.0 8.0 |
* B(3 * 6) = | 2.0 3.0 0.0 4.0 2.0 -2.0 |
* | 11.0 -2.0 -21.0 -1.0 -1.0 5.0 |
*
* | 3.0 2.0 -1.0 3.0 -1.0 -11.0|
* | 9.0 13.0 5.0 5.0 3.0 -5.0 |
* C(5 * 6) = | -1.0 -6.0 -3.0 3.0 -1.0 32.0|
* | . . . . . . |
* | . . . . . . |
*/
float a[12] = {4.0, 3.0, 6.0, 0,
3.0, 4.0, -7.0, 0,
6.0, -7.0, -1.0, 0};
float b[18] = {-1.0, 2.0, 11.0,
-3.0, 3.0, -2.0,
2.0, 0, -21.0,
2.0, 4.0, -1.0,
-1.0, 2.0, -1.0,
8.0, -2.0, 5.0};
float c[30] = {6.0, 9.0, -1.0, 0, 0,
2.0, 13.0, -6.0, 0, 0,
-1.0, 5.0, -3.0, 0, 0,
3.0, 5.0, 3.0, 0, 0,
-1.0, 3.0, -1.0, 0, 0,
-11.0, -5.0, 32.0, 0, 0};
cblas_ssymm(CblasColMajor,CblasLeft,CblasLower, m, n, alpha, a, lda, b, ldb, beta, c, ldc);
/**
* Output C
* | 148.0 -26.0 -238.0 34.0 -10.0 90.0 |
* | -126.0 60.0 316.0 68.0 30.0 -48.0 |
* C(5 * 6) = | -64.0 -86.0 60.0 -24.0 -40.0 178.0|
* | . . . . . . |
* | . . . . . . |
*
*/
Fortran interface：
INTEGER :: M=3, N=6
INTEGER :: LDA=4, LDB=3, LDC=5
REAL(4) :: ALPHA=2.0, BETA=2.0
REAL(4) :: A(4, 3), B(3,6), C(5, 6)
DATA A/4.0, 3.0, 6.0, 0,
$ 3.0, 4.0, -7.0, 0,
$ 6.0, -7.0, -1.0, 0/
DATA B/-1.0, 2.0, 11.0,
$ -3.0, 3.0, -2.0,
$ 2.0, 0, -21.0,
$ 2.0, 4.0, -1.0,
$ -1.0, 2.0, -1.0,
$ 8.0, -2.0, 5.0/
DATA C/6.0, 9.0, -1.0, 0, 0,
$ 2.0, 13.0, -6.0, 0, 0,
$ -1.0, 5.0, -3.0, 0, 0,
$ 3.0, 5.0, 3.0, 0, 0,
$ -1.0, 3.0, -1.0, 0, 0,
$ -11.0, -5.0, 32.0, 0, 0/
EXTERNAL SSYMM
CALL SSYMM('L', 'L', M, N, ALPHA, A, LDA, B, LDB, BETA, C, LDC)
* Output C
* | 148.0 -26.0 -238.0 34.0 -10.0 90.0 |
* | -126.0 60.0 316.0 68.0 30.0 -48.0 |
* C(5 * 6) = | -64.0 -86.0 60.0 -24.0 -40.0 178.0|
* | . . . . . . |
* | . . . . . . |