---
title: ?symm
description: "计算矩阵A为对称矩阵的矩阵乘法。"
url: https://www.hikunpeng.com/document/detail/zh/kunpengboostkithistory/2200/accel/kunpengaccel_kml_16_0060.html
sourcePath: /source/zh/kunpengboostkithistory/2200/accel/kunpengaccel_kml_16_0060.html
indexId: eb495fb3db0bc50f7142faf12cff081588c84555764495df293ff6bd86883cc574
---
# ?symm

计算矩阵A为对称矩阵的矩阵乘法。

或

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的左边还是右边。 Side = CblasLeft，。 Side = CblasRight，。 | 输入 |
| Uplo | 枚举类型CBLAS\_UPLO | 表示使用矩阵A的上三角或下三角。 如果Uplo=CblasUpper，则使用A的上三角部分。 如果Uplo=CblasLower，则使用A的下三角部分。 | 输入 |
| M | 整型数 | 矩阵C的行 A\*B为矩阵A的行。 B\*A为矩阵B的行。 | 输入 |
| N | 整型数 | 矩阵C的列 A\*B为矩阵B的列。 B\*A为矩阵A的列。 | 输入 |
| alpha | 在ssymm中，alpha为单精度浮点类型。 在dsymm中，alpha为双精度浮点类型。 在csymm中，alpha为单精度复数类型。 在zsymm中，alpha为双精度复数类型。 | 乘法系数。 | 输入 |
| A | 在ssymm中，A为单精度浮点类型。 在dsymm中，A为双精度浮点类型。 在csymm中，A为单精度复数类型。 在zsymm中，A为双精度复数类型。 | 矩阵A(lda, ka)。 如果Side = CblasLeft，ka = m，否则ka = n。 | 输入 |
| lda | 整型数 | 如果Side = CblasLeft，lda至少max(1, m)，否则至少max(1, n)。 | 输入 |
| B | 在ssymm中，B为单精度浮点类型。 在dsymm中，B为双精度浮点类型。 在csymm中，B为单精度复数类型。 在zsymm中，B为双精度复数类型。 | 矩阵B(ldb, n)。 | 输入 |
| ldb | 整型数 | 矩阵为列存，ldb至少max(1, m)。 矩阵为行存，ldb至少max(1, n)。 | 输入 |
| beta | 在ssymm中，beta为单精度浮点类型。 在dsymm中，beta为双精度浮点类型。 在csymm中，beta为单精度复数类型。 在zsymm中，beta为双精度复数类型。 | 乘法系数。 | 输入 |
| C | 在ssymm中，C为单精度浮点类型。 在dsymm中，C为双精度浮点类型。 在csymm中，C为单精度复数类型。 在zsymm中，C为双精度复数类型。 | 矩阵C。 | 输入/输出 |
| ldc | 整型数 | 矩阵为列存，ldc至少max(1, m)。 矩阵为行存，ldc至少max(1, n)。 | 输入 |


#### 依赖

#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|
*                       |   .      .      .      .     .      .   |
*                       |   .      .      .      .     .      .   |
```