中文
注册
我要评分
文档获取效率
文档正确性
内容完整性
文档易理解
在线提单
论坛求助

?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的左边还是右边。

  • 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中是单精度浮点类型。
  • 在dsymm中是双精度浮点类型。
  • 在csymm中是单精度复数类型。
  • 在zsymm中是双精度复数类型。

乘法系数。

输入

A

  • 在ssymm中是单精度浮点类型。
  • 在dsymm中是双精度浮点类型。
  • 在csymm中是单精度复数类型。
  • 在zsymm中是双精度复数类型。

矩阵A(lda, ka)。

如果Side = CblasLeft,ka = m,否则ka = n。

输入

lda

整型数

如果Side = CblasLeft,lda至少max(1, m),否则至少max(1, n)。

输入

B

  • 在ssymm中是单精度浮点类型。
  • 在dsymm中是双精度浮点类型。
  • 在csymm中是单精度复数类型。
  • 在zsymm中是双精度复数类型。

矩阵B(ldb, n)。

输入

ldb

整型数

  • 矩阵为列存,ldb至少max(1, m)。
  • 矩阵为行存,ldb至少max(1, n)。

输入

beta

  • 在ssymm中是单精度浮点类型。
  • 在dsymm中是双精度浮点类型。
  • 在csymm中是单精度复数类型。
  • 在zsymm中是双精度复数类型。

乘法系数。

输入

C

  • 在ssymm中是单精度浮点类型。
  • 在dsymm中是双精度浮点类型。
  • 在csymm中是单精度复数类型。
  • 在zsymm中是双精度复数类型。

矩阵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| 
*                       |   .      .      .      .     .      .   | 
*                       |   .      .      .      .     .      .   |