aclnnQuantConvolution
产品支持情况
| 产品 | 是否支持 |
|---|---|
| Ascend 950PR/Ascend 950DT | √ |
| Atlas A3 训练系列产品/Atlas A3 推理系列产品 | √ |
| Atlas A2 训练系列产品/Atlas A2 推理系列产品 | √ |
| Atlas 200I/500 A2 推理产品 | × |
| Atlas 推理系列产品 | × |
| Atlas 训练系列产品 | × |
功能说明
-
接口功能:完成 per-channel 量化的 2D、3D 卷积计算,其中卷积计算过程与 aclnnConvolution 接口一致。
-
计算公式: 我们假定输入(input)的 shape 是 (N,Cin,D,H,W)(N, C_{\text{in}}, D, H, W),weight 的 shape 是 (Cout,Cin,Kd,Kh,Kw)(C_{\text{out}}, C_{\text{in}}, K_d, K_h, K_w),scale 的 shape 是 (Cout)(C_{\text{out}}),bias 的 shape 是 CoutC_{\text{out}},输出(output)的 shape 是 (N,Cout,Dout,Hout,Wout)(N, C_{\text{out}}, D_{\text{out}}, H_{\text{out}}, W_{\text{out}}),其中 NN 表示批次大小(batch size),CC 是通道数,DD、HH 和 WW 分别是样本的深度、高度和宽度,KdK_d、KhK_h 和 KwK_w 分别是卷积核的深度、高度和宽度,那输出将被表示为:
output(Ni,Coutj,Dout,Hout,Wout)=[∑k=0Cin−1weight(Coutj,k)⋆input(Ni,k)]×scale(Coutj)+bias(Coutj)\text{output}(N_i, C_{\text{out}_j}, D_{\text{out}}, H_{\text{out}}, W_{\text{out}}) = \left[\sum_{k = 0}^{C_{\text{in}} - 1} \text{weight}(C_{\text{out}_j}, k) \star \text{input}(N_i, k)\right] \times \text{scale}(C_{\text{out}_j}) + \text{bias}(C_{\text{out}_j})
其中,⋆\star 表示卷积计算,根据卷积输入的维度,卷积的类型(空洞卷积、分组卷积)而定。NN 代表批次大小(batch size),CC 代表通道数,DD、HH 和 WW 分别代表深度、高度和宽度,相应输出维度的计算公式如下:
Dout=[(D+2×padding[0]−dilation[0]×(Kd−1)−1)/stride[0]]+1Hout=[(H+2×padding[1]−dilation[1]×(Kh−1)−1)/stride[1]]+1Wout=[(W+2×padding[2]−dilation[2]×(Kw−1)−1)/stride[2]]+1D_{\text{out}}=[(D + 2 \times padding[0] - dilation[0] \times (K_d - 1) - 1 ) / stride[0]] + 1 \\ H_{\text{out}}=[(H + 2 \times padding[1] - dilation[1] \times (K_h - 1) - 1 ) / stride[1]] + 1 \\ W_{\text{out}}=[(W + 2 \times padding[2] - dilation[2] \times (K_w - 1) - 1 ) / stride[2]] + 1
函数原型
每个算子分为两段式接口,必须先调用 aclnnQuantConvolutionGetWorkspaceSize 接口获取计算所需 workspace 大小以及包含了算子计算流程的执行器,再调用 aclnnQuantConvolution 接口执行计算。
aclnnStatus aclnnQuantConvolutionGetWorkspaceSize(
const aclTensor *input,
const aclTensor *weight,
const aclTensor *bias,
const aclTensor *scale,
const aclTensor *offset,
const aclIntArray *stride,
const aclIntArray *padding,
const aclIntArray *dilation,
bool transposed,
const aclIntArray *outputPadding,
int64_t groups,
int32_t offsetx,
const char *roundMode,
aclTensor *output,
uint64_t *workspaceSize,
aclOpExecutor **executor)
aclnnStatus aclnnQuantConvolution(
void *workspace,
const uint64_t workspaceSize,
aclOpExecutor *executor,
aclrtStream stream)
aclnnQuantConvolutionGetWorkspaceSize
-
参数说明
参数名 输入/输出 描述 使用说明 数据类型 数据格式 维度(shape) 非连续 Tensor input(aclTensor*) 输入 公式中的 input,表示卷积输入。 - 不支持空 Tensor。
- 数据类型与 weight 的数据类型需满足数据类型推导规则(参见互推导关系)。
- input、weight、output 的维度需要相同。
- N≥0,C≥1,D≥0,H≥0,W≥0。
INT8、FLOAT8_E4M3FN、HIFLOAT8 NCHW、NCDHW 4-5 √ weight(aclTensor*) 输入 公式中的 weight,表示卷积权重。 - 不支持空 Tensor。
- 数据类型与 input 的数据类型需满足数据类型推导规则(参见互推导关系)。
- 其 shape 的 C 维度需要与 input 的 C 维度保持一致。
- 所有维度≥1。
INT8、FLOAT8_E4M3FN、HIFLOAT8 NCHW、NCDHW 4-5 √ bias(aclTensor*) 输入 公式中的 bias,表示卷积偏置。 一维且与 weight 第一维相等。 BFLOAT16、FLOAT16、FLOAT、INT32 ND 1 √ scale(aclTensor*) 输入 公式中的 scale,表示量化参数。 一维且与 weight 第一维相等。 FLOAT、INT64、UINT64 ND 1 √ offset(aclTensor*) 输入 预留量化参数。 目前暂不支持,传入空指针 nullptr 即可。 - - - - stride(aclIntArray*) 输入 卷积扫描步长。 - 2d 场景下数组长度=2,3d 场景下数组长度=3。
- strideH 和 strideW 应在 [1,63] 范围内。
- conv3d 场景下 strideD 应在 [1,1000000] 范围内。
INT32 - - - padding(aclIntArray*) 输入 对 input 的填充。 - 值应≥0。
- paddingH 和 paddingW 应在 [0,255] 范围内。
- conv3d 场景下 paddingD 应在 [0,1000000] 范围内。
INT32 - - - dilation(aclIntArray*) 输入 卷积核中元素的间隔。 - 2d 场景下数组长度=2,3d 场景下数组长度=3。
- 值应>0。
- dilationH 和 dilationW 应在 [1,255] 范围内。
- conv3d 场景下 dilationD 应在 [1,1000000] 范围内。
INT32 - - - transposed(bool) 输入 预留参数。表示是否为转置量化卷积。 目前暂不支持,传入 false 即可。 BOOL - - - outputPadding(aclIntArray*) 输入 预留参数。表示转置卷积情况下,对输出所有边的填充。 非转置卷积情况下,忽略该属性配置。目前暂不支持,传入空指针 nullptr 即可。 INT32 - - - groups(int64_t) 输入 表示从输入通道到输出通道的块链接个数。 值≥1,且满足 groups*weight 的 C 维度=input 的 C 维度。 INT64 - - - offsetx(int32_t) 输入 表示量化因子。 [-128,127] 或 0。 INT32 - - - roundMode(char*) 输入 表示取整模式。 rint、round 或 nullptr。 CHAR* - - - output(aclTensor*) 输出 公式中的 output,表示卷积输出。 - 不支持空 Tensor 输出。
- 其 shape 满足卷积的推导规则。
- 通道数等于 weight 第一维,其他维度≥0。
BFLOAT16、FLOAT16、FLOAT、FLOAT8_E4M3FN、HIFLOAT8 NCHW、NCDHW 4-5 √ workspaceSize(uint64_t*) 输出 返回需要在 Device 侧申请的 workspace 大小。 - - - - - executor(aclOpExecutor**) 输出 返回 op 执行器,包含算子计算流程。 - - - - - -
返回值
aclnnStatus:返回状态码,具体参见 aclnn 返回码。
一段接口完成入参校验,出现以下场景时报错:
返回值 错误码 描述 ACLNN_ERR_PARAM_NULLPTR 161001 传入的指针类型入参是空指针。 ACLNN_ERR_PARAM_INVALID 161002 input、weight、bias、scale、offset、output 数据类型和数据格式不在支持的范围之内。 stride、padding、dilation 输入 shape 不对。 groups 输入不对的情况。 scale 和 bias 输入 shape 不对。 output 的 shape 不满足 infershape 结果。 传入 tensor 中任意维度为零的均不满足要求。 input 空间尺度在 padding 操作后小于 weight(经过 dilation 扩张(如存在 dilation>1 的情况))的空间尺度。 weight 和 input 通道数不满足要求。 stride、dilation 小于 0 情况下不满足要求。 当前处理器不支持卷积。 ACLNN_ERR_INNER_NULLPTR 561103 API 内部校验错误,通常由于输入数据或属性的规格不在支持的范围之内导致。 ACLNN_ERR_RUNTIME_ERROR 361001 API 调用 npu runtime 的接口异常,如 SocVersion 不支持。
aclnnQuantConvolution
-
参数说明
参数名 输入/输出 描述 workspace 输入 在 Device 侧申请的 workspace 内存地址。 workspaceSize 输入 在 Device 侧申请的 workspace 大小,由第一段接口 aclnnQuantConvolutionGetWorkspaceSize 获取。 executor 输入 op 执行器,包含了算子计算流程。 stream 输入 指定执行任务的 Stream。 -
返回值
aclnnStatus:返回状态码,具体参见 aclnn 返回码。
约束说明
- 确定性计算
- aclnnQuantConvolution默认确定性实现。
| 约束类型 | Atlas A2 训练系列产品/Atlas A2 推理系列产品、Atlas A3 训练系列产品/Atlas A3 推理系列产品 |
|---|---|
| input、weight | input、weight 数据类型不支持 FLOAT8_E4M3FN、HIFLOAT8,数据格式不支持 NCHW。 |
| bias | bias 数据类型不支持 INT32,会转成 FLOAT 参与计算。 |
| scale | scale 数据类型不支持INT64、UINT64。 |
| padding | padding 的数组长度需要等于 3。 |
| groups | groups 数值必须为 1。 |
| offsetx | offsetx 暂不支持,传入 0 值即可。 |
| roundMode | roundMode 暂不支持,传入空指针 nullptr。 |
| output | output 数据类型支持 BFLOAT16、FLOAT16,数据格式仅支持 NCDHW。 |
| 其他约束 |
|
调用示例
示例代码如下,仅供参考,具体编译和执行过程请参考编译与运行样例。
不同产品型号请参考使用不同的 main 函数。
#include <iostream>
#include <memory>
#include <vector>
#include "acl/acl.h"
#include "aclnnop/aclnn_quant_convolution.h"
#define CHECK_RET(cond, return_expr) \
do { \
if (!(cond)) { \
return_expr; \
} \
} while (0)
#define CHECK_FREE_RET(cond, return_expr) \
do { \
if (!(cond)) { \
Finalize(deviceId, stream); \
return_expr; \
} \
} while (0)
#define LOG_PRINT(message, ...) \
do { \
printf(message, ##__VA_ARGS__); \
} while (0)
int64_t GetShapeSize(const std::vector<int64_t>& shape) {
int64_t shape_size = 1;
for (auto i: shape) {
shape_size *= i;
}
return shape_size;
}
int Init(int32_t deviceId, aclrtStream* stream) {
// 固定写法,资源初始化
auto ret = aclInit(nullptr);
CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("aclInit failed. ERROR: %d\n", ret); return ret);
ret = aclrtSetDevice(deviceId);
CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("aclrtSetDevice failed. ERROR: %d\n", ret); return ret);
ret = aclrtCreateStream(stream);
CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("aclrtCreateStream failed. ERROR: %d\n", ret); return ret);
return 0;
}
template <typename T>
int CreateAclTensor(const std::vector<T>& hostData, const std::vector<int64_t>& shape, void** deviceAddr,
aclDataType dataType, aclTensor** tensor) {
auto size = GetShapeSize(shape) * sizeof(T);
// 调用 aclrtMalloc 申请 device 侧内存
auto ret = aclrtMalloc(deviceAddr, size, ACL_MEM_MALLOC_HUGE_FIRST);
CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("aclrtMalloc failed. ERROR: %d\n", ret); return ret);
// 调用 aclrtMemcpy 将 host 侧数据拷贝到 device 侧内存上
ret = aclrtMemcpy(*deviceAddr, size, hostData.data(), size, ACL_MEMCPY_HOST_TO_DEVICE);
CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("aclrtMemcpy failed. ERROR: %d\n", ret); return ret);
// 计算连续tensor的strides
std::vector<int64_t> strides(shape.size(), 1);
for (int64_t i = shape.size() - 2; i >= 0; i--) {
strides[i] = shape[i + 1] * strides[i + 1];
}
// 调用aclCreateTensor接口创建aclTensor
*tensor = aclCreateTensor(shape.data(), shape.size(), dataType, strides.data(), 0, aclFormat::ACL_FORMAT_NCDHW,
shape.data(), shape.size(), *deviceAddr);
return 0;
}
template <typename T>
int CreateAclTensorND(const std::vector<T>& hostData, const std::vector<int64_t>& shape, void** deviceAddr,
aclDataType dataType, aclTensor** tensor) {
auto size = GetShapeSize(shape) * sizeof(T);
// 调用 aclrtMalloc 申请 device 侧内存
auto ret = aclrtMalloc(deviceAddr, size, ACL_MEM_MALLOC_HUGE_FIRST);
CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("aclrtMalloc failed. ERROR: %d\n", ret); return ret);
// 调用 aclrtMemcpy 将 host 侧数据拷贝到 device 侧内存上
ret = aclrtMemcpy(*deviceAddr, size, hostData.data(), size, ACL_MEMCPY_HOST_TO_DEVICE);
CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("aclrtMemcpy failed. ERROR: %d\n", ret); return ret);
// 计算连续tensor的strides
std::vector<int64_t> strides(shape.size(), 1);
for (int64_t i = shape.size() - 2; i >= 0; i--) {
strides[i] = shape[i + 1] * strides[i + 1];
}
// 调用aclCreateTensor接口创建aclTensor
*tensor = aclCreateTensor(shape.data(), shape.size(), dataType, strides.data(), 0, aclFormat::ACL_FORMAT_ND,
shape.data(), shape.size(), *deviceAddr);
return 0;
}
void Finalize(int32_t deviceId, aclrtStream& stream)
{
aclrtDestroyStream(stream);
aclrtResetDevice(deviceId);
aclFinalize();
}
int aclnnQuantConvolutionTest(int32_t deviceId, aclrtStream& stream, std::vector<aclDataType> dtypesInfo)
{
auto ret = Init(deviceId, &stream);
// check根据自己的需要处理
CHECK_FREE_RET(ret == ACL_SUCCESS, LOG_PRINT("Init acl failed. ERROR: %d\n", ret); return ret);
// 2. 构造输入与输出,需要根据 API 的接口自定义构造
std::vector<int64_t> shapeInput = {2, 2, 32, 32, 32};
std::vector<int64_t> shapeWeight = {2, 2, 3, 3, 3};
std::vector<int64_t> shapeScale = {2};
std::vector<int64_t> shapeBias = {2};
std::vector<int64_t> shapeResult = {2, 2, 32, 32, 32};
std::vector<int64_t> convStrides;
std::vector<int64_t> convPads;
std::vector<int64_t> convOutPads;
std::vector<int64_t> convDilations;
void* deviceDataA = nullptr;
void* deviceDataB = nullptr;
void* deviceDataScale = nullptr;
void* deviceDataBias = nullptr;
void* deviceDataResult = nullptr;
aclTensor* input = nullptr;
aclTensor* weight = nullptr;
aclTensor* scale= nullptr;
aclTensor* bias= nullptr;
aclTensor* result = nullptr;
std::vector<int8_t> inputData(GetShapeSize(shapeInput), 1);
std::vector<int8_t> weightData(GetShapeSize(shapeWeight), 1);
std::vector<float> biasData(GetShapeSize(shapeBias), 1);
std::vector<float> scaleData(GetShapeSize(shapeScale), 1);
std::vector<uint16_t> outputData(GetShapeSize(shapeResult), 1);
convStrides = {1, 1, 1};
convPads = {1, 1, 1};
convOutPads = {1, 1, 1};
convDilations = {1, 1, 1};
aclDataType inputDtype = dtypesInfo[0];
aclDataType weightDtype = dtypesInfo[1];
aclDataType biasDtype = dtypesInfo[2];
aclDataType scaleDtype = dtypesInfo[3];
aclDataType outputDtype = dtypesInfo[4];
// 创建input aclTensor
ret = CreateAclTensor(inputData, shapeInput, &deviceDataA, inputDtype, &input);
std::unique_ptr<aclTensor, aclnnStatus (*)(const aclTensor *)> inputTensorPtr(input, aclDestroyTensor);
std::unique_ptr<void, aclError (*)(void *)> deviceDataAPtr(deviceDataA, aclrtFree);
CHECK_FREE_RET(ret == ACL_SUCCESS, return ret);
// 创建weight aclTensor
ret = CreateAclTensor(weightData, shapeWeight, &deviceDataB, weightDtype, &weight);
std::unique_ptr<aclTensor, aclnnStatus (*)(const aclTensor *)> weightTensorPtr(weight, aclDestroyTensor);
std::unique_ptr<void, aclError (*)(void *)> deviceDataBPtr(deviceDataB, aclrtFree);
CHECK_FREE_RET(ret == ACL_SUCCESS, return ret);
// 创建scale
ret = CreateAclTensorND(scaleData, shapeScale, &deviceDataScale, scaleDtype, &scale);
std::unique_ptr<aclTensor, aclnnStatus (*)(const aclTensor *)> scaleTensorPtr(scale, aclDestroyTensor);
std::unique_ptr<void, aclError (*)(void *)> deviceDataScalePtr(deviceDataScale, aclrtFree);
CHECK_FREE_RET(ret == ACL_SUCCESS, return ret);
// 创建bias
ret = CreateAclTensorND(biasData, shapeBias, &deviceDataBias, biasDtype, &bias);
std::unique_ptr<aclTensor, aclnnStatus (*)(const aclTensor *)> biasTensorPtr(bias, aclDestroyTensor);
std::unique_ptr<void, aclError (*)(void *)> deviceDataBiasPtr(deviceDataBias, aclrtFree);
CHECK_FREE_RET(ret == ACL_SUCCESS, return ret);
// 创建out aclTensor
ret = CreateAclTensor(outputData, shapeResult, &deviceDataResult, outputDtype, &result);
std::unique_ptr<aclTensor, aclnnStatus (*)(const aclTensor *)> outputTensorPtr(result, aclDestroyTensor);
std::unique_ptr<void, aclError (*)(void *)> deviceDataResultPtr(deviceDataResult, aclrtFree);
CHECK_FREE_RET(ret == ACL_SUCCESS, return ret);
aclIntArray *strides = aclCreateIntArray(convStrides.data(), 3);
std::unique_ptr<aclIntArray, aclnnStatus (*)(const aclIntArray *)> stridesPtr(strides, aclDestroyIntArray);
CHECK_FREE_RET(strides != nullptr, return ACL_ERROR_INTERNAL_ERROR);
aclIntArray *pads = aclCreateIntArray(convPads.data(), 3);
std::unique_ptr<aclIntArray, aclnnStatus (*)(const aclIntArray *)> padsPtr(pads, aclDestroyIntArray);
CHECK_FREE_RET(pads != nullptr, return ACL_ERROR_INTERNAL_ERROR);
aclIntArray *outPads = aclCreateIntArray(convOutPads.data(), 3);
std::unique_ptr<aclIntArray, aclnnStatus (*)(const aclIntArray *)> outPadsPtr(outPads, aclDestroyIntArray);
CHECK_FREE_RET(outPads != nullptr, return ACL_ERROR_INTERNAL_ERROR);
aclIntArray *dilations = aclCreateIntArray(convDilations.data(), 3);
std::unique_ptr<aclIntArray, aclnnStatus (*)(const aclIntArray *)> dilationsPtr(dilations, aclDestroyIntArray);
CHECK_FREE_RET(dilations != nullptr, return ACL_ERROR_INTERNAL_ERROR);
// 3. 调用 CANN 算子库 API,需要修改为具体的 API
uint64_t workspaceSize = 0;
aclOpExecutor* executor;
// 调用aclnnConvolution第一段接口
ret = aclnnQuantConvolutionGetWorkspaceSize(input, weight, bias, scale, nullptr, strides, pads, dilations,
false, outPads, 1, 0, nullptr, result, &workspaceSize, &executor);
CHECK_FREE_RET(ret == ACL_SUCCESS, LOG_PRINT("aclnnQuantConvolutionGetWorkspaceSize failed. ERROR: %d\n", ret); return ret);
// 根据第一段接口计算出的workspaceSize申请device内存
void* workspaceAddr = nullptr;
std::unique_ptr<void, aclError (*)(void *)> workspaceAddrPtr(nullptr, aclrtFree);
if (workspaceSize > 0) {
ret = aclrtMalloc(&workspaceAddr, workspaceSize, ACL_MEM_MALLOC_HUGE_FIRST);
CHECK_FREE_RET(ret == ACL_SUCCESS, LOG_PRINT("allocate workspace failed. ERROR: %d\n", ret); return ret);
workspaceAddrPtr.reset(workspaceAddr);
}
// 调用aclnnConvolution第二段接口
ret = aclnnQuantConvolution(workspaceAddr, workspaceSize, executor, stream);
CHECK_FREE_RET(ret == ACL_SUCCESS, LOG_PRINT("aclnnQuantConvolution failed. ERROR: %d\n", ret); return ret);
// 4. (固定写法)同步等待任务执行结束
ret = aclrtSynchronizeStream(stream);
CHECK_FREE_RET(ret == ACL_SUCCESS, LOG_PRINT("aclrtSynchronizeStream failed. ERROR: %d\n", ret); return ret);
// 5. 获取输出的值,将 device 侧内存上的结果拷贝至 host 侧,需要根据具体 API 的接口定义修改
auto size = GetShapeSize(shapeResult);
std::vector<uint16_t> resultData(size, 0);
ret = aclrtMemcpy(resultData.data(), resultData.size() * sizeof(resultData[0]), deviceDataResult,
size * sizeof(uint16_t), ACL_MEMCPY_DEVICE_TO_HOST);
CHECK_FREE_RET(ret == ACL_SUCCESS, LOG_PRINT("copy result from device to host failed. ERROR: %d\n", ret); return ret);
for (int64_t i = 0; i < size; i++) {
LOG_PRINT("result[%ld] is: %d\n", i, resultData[i]);
}
return ACL_SUCCESS;
}
int main() {
// 1. (固定写法)device/stream 初始化,参考 acl API 手册
// 根据自己的实际device填写deviceId
int32_t deviceId = 0;
aclrtStream stream;
std::vector<aclDataType> dtypesInfo = {aclDataType::ACL_INT8, aclDataType::ACL_INT8, aclDataType::ACL_FLOAT,
aclDataType::ACL_FLOAT, aclDataType::ACL_BF16}; // 分别是input/weight/bias/scale/output的datatype
auto ret = aclnnQuantConvolutionTest(deviceId, stream, dtypesInfo);
CHECK_FREE_RET(ret == ACL_SUCCESS, LOG_PRINT("aclnnQuantConvolutionTest failed. ERROR: %d\n", ret); return ret);
Finalize(deviceId, stream);
return 0;
}