---
title: ?imatsub
description: "矩阵缩放后in-place减计算。"
url: https://www.hikunpeng.com/document/detail/zh/kunpengboostkithistory/240RC2/accel/kunpengaccel_kml_16_0291.html
sourcePath: /source/zh/kunpengboostkithistory/240RC2/accel/kunpengaccel_kml_16_0291.html
indexId: 9258bbc84d984d03d0aed1b6c985947aedb2f28c4dd4f0f33b7ddf595ff83c3876
---
# ?imatsub

矩阵缩放后in-place减计算。

即。

op可取值：。

#### 接口定义

C interface：

void cblas_simatsub(const enum CBLAS_ORDER ordering, const enum CBLAS_TRANSPOSE TransA, const BLASINT m, const BLASINT n, const float alpha, float *A, const BLASINT lda, float *BC, const BLASINT ldb);

void cblas_dimatsub(const enum CBLAS_ORDER ordering, const enum CBLAS_TRANSPOSE TransA, const BLASINT m, const BLASINT n, const double alpha, double *A, const BLASINT lda, double *BC, const BLASINT ldb);

void cblas_cimatsub(const enum CBLAS_ORDER ordering, const enum CBLAS_TRANSPOSE TransA, const BLASINT m, const BLASINT n, const float *alpha, float *A, const BLASINT lda, float *BC, const BLASINT ldb);

void cblas_zimatsub(const enum CBLAS_ORDER ordering, const enum CBLAS_TRANSPOSE TransA, const BLASINT m, const BLASINT n, const double *alpha, double *A, const BLASINT lda, double *BC, const BLASINT ldb);

Fortran interface：

CALL SIMATADD(ORDER, TRANSA, M, N, ALPHA, A, LDA, BC, LDB)

CALL DIMATADD(ORDER, TRANSA, M, N, ALPHA, A, LDA, BC, LDB)

CALL CIMATADD(ORDER, TRANSA, M, N, ALPHA, A, LDA, BC, LDB)

CALL ZIMATADD(ORDER, TRANSA, M, N, ALPHA, A, LDA, BC, LDB)


#### 参数

| 参数名 | 类型 | 描述 | 输入/输出 |
| --- | --- | --- | --- |
| order | 枚举类型CBLAS\_ORDER | 表示矩阵是行主序或列主序。 | 输入 |
| TransA | 枚举类型CBLAS\_TRANSPOSE | 矩阵A为常规矩阵，转置矩阵，共轭矩阵或转置共轭矩阵。 如果TransA = CblasNoTrans，op(A) = A。 如果TransA = CblasTrans，op(A) = A'。 如果TransA = CblasConjNoTrans，op(A) = conj(A)。 如果TransA = CblasConjTrans，op(A) = conj(A')。 | 输入 |
| M | 整型数 | 矩阵op(A)、BC的行。 | 输入 |
| N | 整型数 | 矩阵op(A)、BC的列。 | 输入 |
| alpha | 在simatsub中是单精度浮点类型。 在dimatsub中是双精度浮点类型。 在cimatsub中是单精度复数类型。 在zimatsub中是双精度复数类型。 | 表示乘法系数。 | 输入 |
| A | 在simatsub中是单精度浮点类型。 在dimatsub中是双精度浮点类型。 在cimatsub中是单精度复数类型。 在zimatsub中是双精度复数类型。 | 矩阵A。 | 输入 |
| lda | 整型数 | 矩阵A为列存，TransA = CblasNoTrans，lda至少max(1, m)，否则max(1, n)。 矩阵A为行存，TransA = CblasNoTrans，lda至少max(1, n)，否则max(1, m)。 | 输入 |
| BC | 在simatsub中是单精度浮点类型。 在dimatsub中是双精度浮点类型。 在cimatsub中是单精度复数类型。 在zimatsub中是双精度复数类型。 | 矩阵BC。 | 输入 |
| ldb | 整型数 | 矩阵为列存，TransB = CblasNoTrans，ldb至少max(1, m)，否则max(1, n )。 矩阵为行存，TransB = CblasNoTrans，ldb至少max(1, n)，否则max(1, m)。 | 输入 |


#### 依赖

#include "kblas.h"


#### 示例

C interface：

```
int m = 4;
int n = 3;
int lda = 4;
int ldb = 4;
float alpha = 2.0;
/**
*    A:
*         0.878666, 0.587949, 0.911689,
*         0.537057, 0.738171, 0.775693,
*         0.933167, 0.700325, 0.563810,
*         0.712020, 0.606780, 0.538903
*    BC:
*         0.780532, 0.649196, 0.703486,
*         0.854725, 0.511782, 0.570734,
*         0.550857, 0.832913, 0.937889,
*         0.838400, 0.885429, 0.970482
*/
float a[12] = {0.878666, 0.587949, 0.911689,
0.537057, 0.738171, 0.775693,
0.933167, 0.700325, 0.563810,
0.712020, 0.606780, 0.538903
};
float bc[12] = {0.780532, 0.649196, 0.703486,
0.854725, 0.511782, 0.570734,
0.550857, 0.832913, 0.937889,
0.838400, 0.885429, 0.970482
};
float c[12] = {0};
cblas_simatsub(CblasColMajor, CblasNoTrans, m, n, alpha, a, lda, bc, ldb);
/**
* Output C:
*        -0.976800, -0.526701, -1.119892,
*        -0.219388, -0.964561, -0.980653,
*        -1.315476, -0.567737, -0.189731,
*        -0.585641, -0.328131, -0.107325
*/
```

Fortran interface：

```
PROGRAM IMATSUB
INTEGER :: M=4, N=3
INTEGER :: LDA=4, LDB=4
REAL(4) :: ALPHA=2.0
REAL(4) :: A(4, 3), BC(4, 3), C(4, 3)
DATA A/0.878666, 0.587949, 0.911689,
0.537057, 0.738171, 0.775693,
0.933167, 0.700325, 0.563810,
0.712020, 0.606780, 0.538903/
DATA BC/0.780532, 0.649196, 0.703486,
0.854725, 0.511782, 0.570734,
0.550857, 0.832913, 0.937889,
0.838400, 0.885429, 0.970482/
EXTERNAL SIMATSUB
CALL SIMATSUB('C', 'N', M, N, ALPHA, A, LDA, BC, LDB)
END
* Output BC:
*        0.694343, 0.882727, 0.822963,
*        0.859743, 0.543875, 0.927105,
*        0.516767, 0.659971, 0.950618,
*        0.849312, 0.576374, 0.557569
```