aclnnAllGatherMatmulV2
产品支持情况
| 产品 | 是否支持 |
|---|---|
| Ascend 950PR/Ascend 950DT | √ |
| Atlas A3 训练系列产品/Atlas A3 推理系列产品 | √ |
| Atlas A2 训练系列产品/Atlas A2 推理系列产品 | √ |
| Atlas 200I/500 A2 推理产品 | × |
| Atlas 推理系列产品 | × |
| Atlas 训练系列产品 | × |
功能说明
-
接口功能:
aclnnAllGatherMatmulV2接口是对aclnnAllGatherMatmul接口的功能拓展,在支持x1和x2输入类型为FLOAT16/BFLOAT16的基础上,新增功能如下: -
计算公式:
- 情形1:如果x1和x2数据类型为FLOAT16/BFLOAT16时,入参x1进行AllGather后,对x1、x2进行MatMul计算。
output=AllGather(x1)@x2+biasoutput=AllGather(x1)@x2 + bias
gatherOut=AllGather(x1)gatherOut=AllGather(x1)
- 情形2:如果x1和x2数据类型为FLOAT8_E4M3FN/FLOAT8_E5M2/HIFLOAT8的pertensor场景,或者x1和x2数据类型为INT8/INT4的perchannel、pertoken场景,且不输出amaxOut,入参x1进行AllGather后,对x1、x2进行MatMul计算,然后进行dequant操作。
output=(x1Scale∗x2Scale)∗(AllGather(x1)@x2+bias)output=(x1Scale*x2Scale)*(AllGather(x1)@x2 + bias)
gatherOut=AllGather(x1)gatherOut=AllGather(x1)
-
情形3:如果x1和x2数据类型为FLOAT8_E4M3FN/FLOAT8_E5M2/HIFLOAT8的perblock场景,且不输出amaxOut, 当x1为(m, k)、x2为(k, n)时, x1Scale为(ceilDiv(m, 128), ceilDiv(k, 128))、x2Scale为(ceilDiv(k, 128), ceilDiv(n, 128))时,入参x1和x1Scale进行AllGather后,对x1、x2进行perblock量化MatMul计算,然后进行dequant操作。
output=∑0⌊kblockSize=128⌋(AllGather(x1)pr@x2rq∗(AllGather(x1Scale)pr∗x2Scalerq))output=\sum_{0}^{\left \lfloor \frac{k}{blockSize=128} \right \rfloor} (AllGather(x1)_{pr}@x2_{rq}*(AllGather(x1Scale)_{pr}*x2Scale_{rq}))
gatherOut=AllGather(x1)gatherOut=AllGather(x1)
-
情形4:如果x1和x2数据类型为FLOAT8_E4M3FN/FLOAT8_E5M2的mx量化场景,x1为(m, k)、x2 为(n, k),且x1Scale为(m, ceilDiv(k, 64), 2)、x2Scale为(n, ceilDiv(k, 64), 2),入参x1和x1Scale进行AllGather后,对x1、x2进行MatMul计算,然后进行dequant操作;
output=∑0⌊kblockSize=32⌋(AllGather(x1)pr@x2rq∗(AllGather(x1Scale)pr∗x2Scalerq))output=\sum_{0}^{\left \lfloor \frac{k}{blockSize=32} \right \rfloor} (AllGather(x1)_{pr}@x2_{rq}*(AllGather(x1Scale)_{pr}*x2Scale_{rq}))
gatherOut=AllGather(x1)gatherOut=AllGather(x1)
函数原型
每个算子分为两段式接口,必须先调用aclnnAllGatherMatmulV2GetWorkspaceSize接口获取计算所需workspace大小以及包含了算子计算流程的执行器,再调用aclnnAllGatherMatmulV2接口执行计算。
aclnnStatus aclnnAllGatherMatmulV2GetWorkspaceSize(
const aclTensor* x1,
const aclTensor* x2,
const aclTensor* bias,
const aclTensor* x1Scale,
const aclTensor* x2Scale,
const aclTensor* quantScale,
int64_t blockSize,
const char* group,
int64_t gatherIndex,
int64_t commTurn,
int64_t streamMode,
int64_t groupSize,
const char* commMode,
aclTensor* output,
aclTensor* gatherOut,
aclTensor* amaxOut,
uint64_t* workspaceSize,
aclOpExecutor** executor)
aclnnStatus aclnnAllGatherMatmulV2(
void* workspace,
uint64_t workspaceSize,
aclOpExecutor* executor,
aclrtStream stream)
aclnnAllGatherMatmulV2GetWorkspaceSize
-
参数说明:
参数名 输入/输出 描述 使用说明 数据类型 数据格式 维度(shape) 非连续Tensor x1 (aclTensor*) 输入 MM左矩阵,即计算公式中的x1。 当前版本仅支持两维输入,shape为[m, k],且仅支持不转置场景。 FLOAT16、BFLOAT16、FLOAT8_E4M3FN、FLOAT8_E5M2、HIFLOAT8、INT8、INT4 ND 2 - x2 (aclTensor*) 输入 MM右矩阵,即计算公式中的x2。 - 当前版本仅支持二维输入,shape为[k, n],支持转置/不转置场景。
- 仅支持两根轴转置情况下的[非连续的Tensor]。
FLOAT16、BFLOAT16、FLOAT8_E4M3FN、FLOAT8_E5M2、HIFLOAT8、INT8、INT4 ND 2 √(仅适用转置场景) bias (aclTensor*) 输入 即计算公式中的bias。 - Ascend 950PR/Ascend 950DT:支持传入一维输入或者空指针。
- Atlas A2 训练系列产品/Atlas A2 推理系列产品:当前版本仅支持传入空指针。
FLOAT16、BFLOAT16、FLOAT ND 1 - x1Scale (aclTensor*) 输入 mm左矩阵反量化参数。 支持传入空指针场景。 FLOAT16、BFLOAT16、FLOAT ND 1-3 - x2Scale (aclTensor*) 输入 mm右矩阵反量化参数。 支持传入空指针场景。 FLOAT16、BFLOAT16、FLOAT ND 1-3 √(仅适用转置场景) quantScale (aclTensor*) 输入 即计算公式中的bias。 当前仅支持传入空指针场景。 FLOAT ND 1 - blockSize (int64_t) 输入 用于表示mm输出矩阵在M轴方向上和N轴方向上可以用于对应方向上的多少个数的量化。 blockSize由blockSizeM、blockSizeN、blockSizeK三个值拼接而成,每个值占16位,计算公式为blockSize = blockSizeK | blockSizeN << 16 | blockSizeM << 32,mm输出矩阵不涉及K轴,blockSizeK固定为0, 当前版本只支持blockSizeM=blockSizeN=0。 - - - - group (char*) 输入 通信域名称。 通过Hccl提供的接口“extern HcclResult HcclGetCommName(HcclComm comm, char* commName);”获取,其中commName即为group。 - - - - gatherIndex (int64_t) 输入 标识gather目标。 - 0表示目标为x1,1表示目标为x2。
- 当前版本仅支持输入0。
- - - - commTurn (int64_t) 输入 通信数据切分数,即总数据量/单次通信量。 当前版本仅支持输入0。 - - - - streamMode (int64_t) 输入 流模式的枚举。 当前只支持枚举值1。 - - - - groupSize (int64_t) 输入 用于表示反量化中x1Scale/x2Scale输入的一个数在其所在的对应维度方向上可以用于该方向x1/x2输入的多少个数的反量化。 groupSize输入由3个方向的groupSizeM、groupSizeN、groupSizeK三个值拼接组成,每个值占16位,计算公式为groupSize = groupSizeK | groupSizeN << 16 | groupSizeM << 32。 - - - - commMode (char*) 输入 通信模式。 - - - - - output (aclTensor*) 输出 AllGather通信与MatMul计算的结果,即计算公式中的output。 - Ascend 950PR/Ascend 950DT:支持空Tensor。
- Atlas A2 训练系列产品/Atlas A2 推理系列产品:不支持空Tensor。
FLOAT16、BFLOAT16、FLOAT ND 2 - gatherOut (aclTensor*) 输出 仅输出all_gather通信后的结果。即公式中的gatherOut。 - Ascend 950PR/Ascend 950DT:支持空Tensor。
- Atlas A2 训练系列产品/Atlas A2 推理系列产品:不支持空Tensor。
- 数据类型与x1的数据类型保持一致。
FLOAT16、BFLOAT16、FLOAT8_E4M3FN、FLOAT8_E5M2、HIFLOAT8、INT8、INT4 ND 2 - amaxOut (aclTensor*) 输出 MM计算的最大值结果,即公式中的amaxOut。 当前版本仅支持nullptr或空tensor。 FLOAT ND 1 - workspaceSize (uint64_t*) 输出 返回需要在Device侧申请的workspace大小。 - - - - - executor (aclOpExecutor**) 输出 返回op执行器,包含了算子计算流程。 - - - - - - Atlas A2 训练系列产品/Atlas A2 推理系列产品、Atlas A3 训练系列产品/Atlas A3 推理系列产品:
- x1、x2:数据类型支持FLOAT16、BFLOAT16、INT8、INT4。
- bias:在commMode为aiv时,当前版本仅支持输入nullptr。
- x1Scale:数据类型支持FLOAT。当x1和x2数据类型为FLOAT16/BFLOAT16时,仅支持输入为nullptr。在pertoken场景,shape为(m, 1)。
- x2Scale:数据类型支持FLOAT、INT64。INT64数据类型仅在x1和x2数据类型为INT4或者output数据类型为FLOAT16场景支持。当x1和x2数据类型为FLOAT16/BFLOAT16时,仅支持输入为nullptr。在perchannel场景,shape为(1, n)。
- groupSize:当前版本仅支持输入为0。
- commMode:当前仅支持aiv模式。aiv模式下使用AI VECTOR核完成通信任务。当前版本仅支持输入“aiv”。
- output:数据类型支持FLOAT16、BFLOAT16。 如果x1类型为FLOAT16、BFLOAT16,则output类型与x1保持一致。
- gatherOut:数据类型支持FLOAT16、BFLOAT16、INT8、INT4。
- Ascend 950PR/Ascend 950DT:
- x1、x2:的数据类型支持FLOAT16、BFLOAT16、FLOAT8_E4M3FN、FLOAT8_E5M2、HIFLOAT8。
- bias:如果x1的数据类型是FLOAT16、BFLOAT16,则bias的数据类型必须为FLOAT16、BFLOAT16。如果x1的数据类型是FLOAT8_E4M3FN、FLOAT8_E5M2、HIFLOAT8时,在pertensor和mx量化场景下,bias的数据类型必须为FLOAT。在perblock场景下,仅支持输入为nullptr。
- x1Scale:当x1和x2数据类型为FLOAT16、BFLOAT16时,仅支持输入为nullptr。在pertensor场景下,shape为[1]。在perblock场景下,shape为[ceilDiv(m, 128), ceilDiv(k, 128)]。在pertensor和perblock场景下,数据类型支持FLOAT。在mx量化场景下,数据类型为FLOAT8_E8M0,shape为(m, ceilDiv(k, 64), 2)。
- x2Scale:当x1和x2数据类型为FLOAT16、BFLOAT16时,仅支持输入为nullptr。在pertensor场景下,shape为[1]。在perblock场景下,shape为[ceilDiv(k, 128), ceilDiv(n, 128)]。在pertensor和perblock场景下,数据类型支持FLOAT。在mx场景下,数据类型为FLOAT8_E8M0,shape为(n, ceilDiv(k, 64), 2),仅支持转置场景。
- commMode:当前版本仅支持输入“ccu”。
- output:如果x1类型为FLOAT16、BFLOAT16,则output类型与x1保持一致。如果x1类型为FLOAT8_E4M3FN、FLOAT8_E5M2、HIFLOAT8,则数据类型支持FLOAT16、BFLOAT16、FLOAT。
- gatherOut:数据类型支持FLOAT16、BFLOAT16、FLOAT8_E4M3FN、FLOAT8_E5M2、HIFLOAT8。
- groupSize:
- 仅当x1Scale和x2Scale输入都是2维及以上数据时,groupSize取值有效,其他场景需传入0。
- groupSize值支持公式推导:传入的groupSize内部会按如下公式分解得到groupSizeM、groupSizeN、groupSizeK,当其中有1个或多个为0,会根据x1/x2/x1Scale/x2Scale输入shape重新设置groupSizeM、groupSizeN、groupSizeK用于计算。设置原理:如果groupSizeM=0,表示m方向量化分组值由接口推导,推导公式为groupSizeM = m / scaleM(需保证m能被scaleM整除),其中m与x1 shape中的m一致,scaleM与x1Scale shape中的m一致;如果groupSizeK=0,表示k方向量化分组值由接口推导,推导公式为groupSizeK = k / scaleK(需保证k能被scaleK整除),其中k与x1 shape中的k一致,scaleK与x1Scale shape中的k一致;如果groupSizeN=0,表示n方向量化分组值由接口推导,推导公式为groupSizeN = n / scaleN(需保证n能被scaleN整除),其中n与x2 shape中的n一致,scaleN与x2Scale shape中的n一致。
groupSize=groupSizeK∣groupSizeN<<16∣groupSizeM<<32groupSize = groupSizeK | groupSizeN << 16 | groupSizeM << 32
- 如果满足重新设置条件,一般情况下,当x1Scale、x2Scale输入都是2维,且数据类型都为FLOAT时,[groupSizeM,groupSizeN,groupSizeK]取值组合会推导为[128, 128, 128],对应groupSize的值为549764202624;当x1Scale、x2Scale输入都是3维,且数据类型都为FLOAT8_E8M0时,[groupSizeM, groupSizeN, groupSizeK]取值组合会推导为[1, 1, 32],对应groupSize的值为4295032864。
-
返回值:
返回aclnnStatus状态码,具体参见aclnn返回码。
第一段接口完成入参校验,出现以下场景时报错:
返回值 错误码 描述 ACLNN_ERR_PARAM_NULLPTR 161001 传入的x1、x2或output是空指针。 ACLNN_ERR_PARAM_INVALID 161002 传入的x1、x2、x1Scale、x2Scale、bias、quantScale、output、gatherOut或amaxOut的数据类型和维度不在支持的范围内。
aclnnAllGatherMatmulV2
-
参数说明:
参数名 输入/输出 描述 workspace 输入 在Device侧申请的workspace内存地址。 workspaceSize 输入 在Device侧申请的workspace大小,由第一段接口 aclnnAllGatherMatmulV2GetWorkspaceSize获取。executor 输入 op执行器,包含了算子计算流程。 stream 输入 指定执行任务的stream。 -
返回值
返回aclnnStatus状态码,具体参见aclnn返回码。
约束说明
-
确定性计算:
aclnnAllGatherMatmulV2默认确定性实现。
-
Ascend 950PR/Ascend 950DT:
- 输入x1为2维,其维度为(m, k)。x2必须是2维,其维度为(k, n),轴满足mm算子入参要求,k轴相等,且k轴取值范围为[256, 65535)。
- bias为1维,shape为(n,)。
- 输出output为2维,其维度为(m*rank_size, n),rank_size为卡数。
- 输出gatherout为2维,其维度为(m*rank_size, k),rank_size为卡数。
- 当x1、x2的数据类型为FLOAT16/BFLOAT16时,output计算输出数据类型和x1、x2保持一致。
- 当x1、x2的数据类型为FLOAT8_E4M3FN/FLOAT8_E5M2/HIFLOAT8时,output输出数据类型支持FLOAT16、BFLOAT16、FLOAT。
- 当x1、x2的数据类型为FLOAT16/BFLOAT16/HIFLOAT8时,x1和x2数据类型需要保持一致。
- 当x1、x2数据类型为FLOAT8_E4M3FN/FLOAT8_E5M2时,x1和x2数据类型可以为其中一种。
- 当x1、x2数据类型为FLOAT16/BFLOAT16/HIFLOAT8/FLOAT8_E4M3FN/FLOAT8_E5M2时,x2矩阵支持转置/不转置场景,x1矩阵只支持不转置场景。
- 当groupSize取值为549764202624,bias必须为空。
- 支持2、4、8、16、32、64卡。
- allgather(x1)集合通信数据总量不能超过16*256MB,集合通信数据总量计算方式为:m * k * sizeof(x1_dtype) * 卡数。由于shape不同,算子内部实现可能存在差异,实际支持的总通信量可能略小于该值。
-
Atlas A3 训练系列产品/Atlas A3 推理系列产品、Atlas A2 训练系列产品/Atlas A2 推理系列产品:
- 只支持x2矩阵转置/不转置,x1矩阵仅支持不转置场景。
- 输入x1必须是2维,其shape为(m, k)。
- 输入x2必须是2维,其shape为(k, n),轴满足mm算子入参要求,k轴相等,且k轴取值范围为[256, 65535)。
- bias仅支持输入nullptr。
- 输出为2维,其shape为(m*rank_size, n), rank_size为卡数。
- 不支持空tensor。
- x1和x2的数据类型需要保持一致。
- x1和x2数据类型为INT4时,k与n必须为偶数。
- 支持2、4、8卡。
调用示例
说明:本示例代码调用了部分HCCL集合通信库接口:HcclGetCommName、HcclCommInitAll、HcclCommDestroy, 请参考 <<HCCL API (C)>>。
示例代码如下,仅供参考,具体编译和执行过程请参考编译与运行样例。
-
Atlas A2 训练系列产品/Atlas A2 推理系列产品
#include <iostream> #include <vector> #include <thread> #include "hccl/hccl.h" #include "aclnnop/aclnn_all_gather_matmul_v2.h" #define CHECK_RET(cond, return_expr) \ do { \ if (!(cond)) { \ return_expr; \ } \ } while (0) #define LOG_PRINT(message, ...) \ do { \ printf(message, ##__VA_ARGS__); \ } while(0) constexpr int DEV_NUM = 2; int64_t GetShapeSize(const std::vector<int64_t> &shape) { int64_t shape_size = 1; for (auto i : shape) { shape_size *= i; } return shape_size; } 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); auto ret = aclrtMalloc(deviceAddr, size, ACL_MEM_MALLOC_HUGE_FIRST); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtMalloc failed. ret: %d\n", ret); return ret); ret = aclrtMemcpy(*deviceAddr, size, hostData.data(), size, ACL_MEMCPY_HOST_TO_DEVICE); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtMemcpy failed. ret: %d\n", ret); return ret); 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]; } *tensor = aclCreateTensor(shape.data(), shape.size(), dataType, strides.data(), 0, aclFormat::ACL_FORMAT_ND, shape.data(), shape.size(), *deviceAddr); return 0; } struct Args { int rankId; HcclComm hcclComm; aclrtStream stream; aclrtContext context; }; int LaunchOneThreadAllGatherMmV2(Args &args) { int ret = aclrtSetCurrentContext(args.context); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtSetCurrentContext failed. ret: %d\n", ret); return ret); char hcomName[128] = {0}; ret = HcclGetCommName(args.hcclComm, hcomName); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] HcclGetCommName failed. ret: %d\n", ret); return -1); LOG_PRINT("[INFO] rank = %d, hcomName = %s, stream = %p, context = %p\n", args.rankId, hcomName, args.stream, args.context); std::vector<int64_t> x1Shape = {32, 256}; std::vector<int64_t> x2Shape = {256, 128}; std::vector<int64_t> x1ScaleShape = {32, 1}; std::vector<int64_t> x2ScaleShape = {1, 128}; std::vector<int64_t> biasShape = {128}; std::vector<int64_t> outShape = {32 * DEV_NUM, 128}; std::vector<int64_t> gatherOutShape = {32 * DEV_NUM, 256}; void *x1DeviceAddr = nullptr; void *x2DeviceAddr = nullptr; void *x1ScaleDeviceAddr = nullptr; void *x2ScaleDeviceAddr = nullptr; void *biasDeviceAddr = nullptr; void *outDeviceAddr = nullptr; void *gatherOutDeviceAddr = nullptr; aclTensor *x1 = nullptr; aclTensor *x2 = nullptr; aclTensor *x1Scale = nullptr; aclTensor *x2Scale = nullptr; aclTensor *bias = nullptr; aclTensor *quantScale = nullptr; aclTensor *out = nullptr; aclTensor *gatherOut = nullptr; aclTensor *amax = nullptr; int64_t gatherIndex = 0; int64_t commTurn = 0; int64_t streamMode = 1; int64_t blockSize = 0; int64_t groupSize = 0; uint64_t workspaceSize = 0; aclOpExecutor *executor = nullptr; void *workspaceAddr = nullptr; long long x1ShapeSize = GetShapeSize(x1Shape); long long x2ShapeSize = GetShapeSize(x2Shape); long long x1ScaleShapeSize = GetShapeSize(x1ScaleShape); long long x2ScaleShapeSize = GetShapeSize(x2ScaleShape); long long biasShapeSize = GetShapeSize(biasShape); long long outShapeSize = GetShapeSize(outShape); long long gatherOutShapeSize = GetShapeSize(gatherOutShape); std::vector<int8_t> x1HostData(x1ShapeSize, 0); std::vector<int8_t> x2HostData(x2ShapeSize, 0); std::vector<int32_t> x1ScaleHostData(x1ScaleShapeSize, 0); std::vector<int32_t> x2ScaleHostData(x2ScaleShapeSize, 0); std::vector<int32_t> biasHostData(biasShapeSize, 0); std::vector<int16_t> outHostData(outShapeSize, 0); std::vector<int8_t> gatherOutHostData(gatherOutShapeSize, 0); ret = CreateAclTensor(x1HostData, x1Shape, &x1DeviceAddr, aclDataType::ACL_INT8, &x1); CHECK_RET(ret == ACL_SUCCESS, return ret); ret = CreateAclTensor(x2HostData, x2Shape, &x2DeviceAddr, aclDataType::ACL_INT8, &x2); CHECK_RET(ret == ACL_SUCCESS, return ret); ret = CreateAclTensor(x1ScaleHostData, x1ScaleShape, &x1ScaleDeviceAddr, aclDataType::ACL_FLOAT, &x1Scale); CHECK_RET(ret == ACL_SUCCESS, return ret); ret = CreateAclTensor(x2ScaleHostData, x2ScaleShape, &x2ScaleDeviceAddr, aclDataType::ACL_FLOAT, &x2Scale); CHECK_RET(ret == ACL_SUCCESS, return ret); ret = CreateAclTensor(outHostData, outShape, &outDeviceAddr, aclDataType::ACL_FLOAT16, &out); CHECK_RET(ret == ACL_SUCCESS, return ret); ret = CreateAclTensor(gatherOutHostData, gatherOutShape, &gatherOutDeviceAddr, aclDataType::ACL_INT8, &gatherOut); CHECK_RET(ret == ACL_SUCCESS, return ret); // 调用第一阶段接口 ret = aclnnAllGatherMatmulV2GetWorkspaceSize( x1, x2, bias, x1Scale, x2Scale, quantScale, blockSize, hcomName, gatherIndex, commTurn, streamMode, groupSize, "aiv", out, gatherOut, amax, &workspaceSize, &executor); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclnnAllGatherMatmulV2GetWorkspaceSize failed. ret = %d \n", ret); return ret); // 根据第一阶段接口计算出的workspaceSize申请device内存 if (workspaceSize > 0) { ret = aclrtMalloc(&workspaceAddr, workspaceSize, ACL_MEM_MALLOC_HUGE_FIRST); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtMalloc workspace failed. ret = %d \n", ret); return ret); } // 调用第二阶段接口 ret = aclnnAllGatherMatmulV2(workspaceAddr, workspaceSize, executor, args.stream); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclnnAllGatherMatmulV2 failed. ret = %d \n", ret); return ret); // (固定写法)同步等待任务执行结束 ret = aclrtSynchronizeStreamWithTimeout(args.stream, 10000); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtSynchronizeStreamWithTimeout failed. ret = %d \n", ret); return ret); LOG_PRINT("[INFO] device_%d aclnnAllGatherMatmulV2 execute successfully.\n", args.rankId); // 释放device资源,需要根据具体API的接口定义修改 if (x1 != nullptr) { aclDestroyTensor(x1); } if (x2 != nullptr) { aclDestroyTensor(x2); } if (x1Scale != nullptr) { aclDestroyTensor(x1Scale); } if (x2Scale != nullptr) { aclDestroyTensor(x2Scale); } if (bias != nullptr) { aclDestroyTensor(bias); } if (quantScale != nullptr) { aclDestroyTensor(quantScale); } if (out != nullptr) { aclDestroyTensor(out); } if (gatherOut != nullptr) { aclDestroyTensor(gatherOut); } if (amax != nullptr) { aclDestroyTensor(amax); } if (x1DeviceAddr != nullptr) { aclrtFree(x1DeviceAddr); } if (x2DeviceAddr != nullptr) { aclrtFree(x2DeviceAddr); } if (x1ScaleDeviceAddr != nullptr) { aclrtFree(x1ScaleDeviceAddr); } if (x2ScaleDeviceAddr != nullptr) { aclrtFree(x2ScaleDeviceAddr); } if (biasDeviceAddr != nullptr) { aclrtFree(biasDeviceAddr); } if (outDeviceAddr != nullptr) { aclrtFree(outDeviceAddr); } if (gatherOutDeviceAddr != nullptr) { aclrtFree(gatherOutDeviceAddr); } if (workspaceSize > 0) { aclrtFree(workspaceAddr); } ret = HcclCommDestroy(args.hcclComm); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] HcclCommDestroy failed. ret = %d \n", ret); return ret); ret = aclrtDestroyStream(args.stream); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtDestroyStream failed. ret = %d \n", ret); return ret); ret = aclrtResetDevice(args.rankId); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtResetDevice failed. ret = %d \n", ret); return ret); ret = aclrtDestroyContext(args.context); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtDestroyContext failed. ret = %d \n", ret); return ret); return 0; } int main(int argc, char *argv[]) { int ret = aclInit(nullptr); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclInit failed. ret = %d \n", ret); return ret); aclrtStream stream[DEV_NUM]; aclrtContext context[DEV_NUM]; for (uint32_t rankId = 0; rankId < DEV_NUM; rankId++) { ret = aclrtSetDevice(rankId); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtSetDevice failed. ret = %d \n", ret); return ret); ret = aclrtCreateContext(&context[rankId], rankId); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtCreateContext failed. ret = %d \n", ret); return ret); ret = aclrtCreateStream(&stream[rankId]); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtCreateStream failed. ret = %d \n", ret); return ret); } int32_t devices[DEV_NUM]; for (int i = 0; i < DEV_NUM; i++) { devices[i] = i; } // 初始化集合通信域 HcclComm comms[DEV_NUM]; ret = HcclCommInitAll(DEV_NUM, devices, comms); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] HcclCommInitAll failed. ret = %d \n", ret); return ret); Args args[DEV_NUM]; // 启动多线程 std::vector<std::unique_ptr<std::thread>> threads(DEV_NUM); for (uint32_t rankId = 0; rankId < DEV_NUM; rankId++) { args[rankId].rankId = rankId; args[rankId].hcclComm = comms[rankId]; args[rankId].context = context[rankId]; args[rankId].stream = stream[rankId]; threads[rankId].reset(new(std::nothrow) std::thread(&LaunchOneThreadAllGatherMmV2, std::ref(args[rankId]))); } for (uint32_t rankId = 0; rankId < DEV_NUM; rankId++) { threads[rankId]->join(); } aclFinalize(); return 0; } -
Ascend 950PR/Ascend 950DT:
#include <iostream> #include <vector> #include <thread> #include "hccl/hccl.h" #include "aclnnop/aclnn_all_gather_matmul_v2.h" #define CHECK_RET(cond, return_expr) \ do { \ if (!(cond)) { \ return_expr; \ } \ } while (0) #define LOG_PRINT(message, ...) \ do { \ printf(message, ##__VA_ARGS__); \ } while(0) constexpr int DEV_NUM = 2; int64_t GetShapeSize(const std::vector<int64_t> &shape) { int64_t shape_size = 1; for (auto i : shape) { shape_size *= i; } return shape_size; } 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); auto ret = aclrtMalloc(deviceAddr, size, ACL_MEM_MALLOC_HUGE_FIRST); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtMalloc failed. ret: %d\n", ret); return ret); ret = aclrtMemcpy(*deviceAddr, size, hostData.data(), size, ACL_MEMCPY_HOST_TO_DEVICE); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtMemcpy failed. ret: %d\n", ret); return ret); 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]; } *tensor = aclCreateTensor(shape.data(), shape.size(), dataType, strides.data(), 0, aclFormat::ACL_FORMAT_ND, shape.data(), shape.size(), *deviceAddr); return 0; } struct Args { int rankId; HcclComm hcclComm; aclrtStream stream; aclrtContext context; }; int LaunchOneThreadAllGatherMmV2(Args &args) { int ret = aclrtSetCurrentContext(args.context); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtSetCurrentContext failed. ret: %d\n", ret); return ret); char hcomName[128] = {0}; ret = HcclGetCommName(args.hcclComm, hcomName); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] HcclGetCommName failed. ret: %d\n", ret); return -1); LOG_PRINT("[INFO] rank = %d, hcomName = %s, stream = %p, context = %p\n", args.rankId, hcomName, args.stream, args.context); std::vector<int64_t> x1Shape = {32, 256}; std::vector<int64_t> x2Shape = {256, 128}; std::vector<int64_t> x1ScaleShape = {1}; std::vector<int64_t> x2ScaleShape = {1}; std::vector<int64_t> biasShape = {128}; std::vector<int64_t> outShape = {32 * DEV_NUM, 128}; std::vector<int64_t> gatherOutShape = {32 * DEV_NUM, 256}; void *x1DeviceAddr = nullptr; void *x2DeviceAddr = nullptr; void *x1ScaleDeviceAddr = nullptr; void *x2ScaleDeviceAddr = nullptr; void *biasDeviceAddr = nullptr; void *outDeviceAddr = nullptr; void *gatherOutDeviceAddr = nullptr; aclTensor *x1 = nullptr; aclTensor *x2 = nullptr; aclTensor *x1Scale = nullptr; aclTensor *x2Scale = nullptr; aclTensor *bias = nullptr; aclTensor *quantScale = nullptr; aclTensor *out = nullptr; aclTensor *gatherOut = nullptr; aclTensor *amax = nullptr; int64_t gatherIndex = 0; int64_t commTurn = 0; int64_t streamMode = 1; int64_t blockSize = 0; int64_t groupSize = 0; uint64_t workspaceSize = 0; aclOpExecutor *executor = nullptr; void *workspaceAddr = nullptr; long long x1ShapeSize = GetShapeSize(x1Shape); long long x2ShapeSize = GetShapeSize(x2Shape); long long x1ScaleShapeSize = GetShapeSize(x1ScaleShape); long long x2ScaleShapeSize = GetShapeSize(x2ScaleShape); long long biasShapeSize = GetShapeSize(biasShape); long long outShapeSize = GetShapeSize(outShape); long long gatherOutShapeSize = GetShapeSize(gatherOutShape); std::vector<int8_t> x1HostData(x1ShapeSize, 0); std::vector<int8_t> x2HostData(x2ShapeSize, 0); std::vector<int32_t> x1ScaleHostData(x1ScaleShapeSize, 0); std::vector<int32_t> x2ScaleHostData(x2ScaleShapeSize, 0); std::vector<int32_t> biasHostData(biasShapeSize, 0); std::vector<int16_t> outHostData(outShapeSize, 0); std::vector<int8_t> gatherOutHostData(gatherOutShapeSize, 0); ret = CreateAclTensor(x1HostData, x1Shape, &x1DeviceAddr, aclDataType::ACL_FLOAT8_E4M3FN, &x1); CHECK_RET(ret == ACL_SUCCESS, return ret); ret = CreateAclTensor(x2HostData, x2Shape, &x2DeviceAddr, aclDataType::ACL_FLOAT8_E4M3FN, &x2); CHECK_RET(ret == ACL_SUCCESS, return ret); ret = CreateAclTensor(x1ScaleHostData, x1ScaleShape, &x1ScaleDeviceAddr, aclDataType::ACL_FLOAT, &x1Scale); CHECK_RET(ret == ACL_SUCCESS, return ret); ret = CreateAclTensor(x2ScaleHostData, x2ScaleShape, &x2ScaleDeviceAddr, aclDataType::ACL_FLOAT, &x2Scale); CHECK_RET(ret == ACL_SUCCESS, return ret); ret = CreateAclTensor(outHostData, outShape, &outDeviceAddr, aclDataType::ACL_FLOAT16, &out); CHECK_RET(ret == ACL_SUCCESS, return ret); ret = CreateAclTensor(gatherOutHostData, gatherOutShape, &gatherOutDeviceAddr, aclDataType::ACL_FLOAT8_E4M3FN, &gatherOut); CHECK_RET(ret == ACL_SUCCESS, return ret); // 调用第一阶段接口 ret = aclnnAllGatherMatmulV2GetWorkspaceSize( x1, x2, bias, x1Scale, x2Scale, quantScale, blockSize, hcomName, gatherIndex, commTurn, streamMode, groupSize, "ccu", out, gatherOut, amax, &workspaceSize, &executor); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclnnAllGatherMatmulV2GetWorkspaceSize failed. ret = %d \n", ret); return ret); // 根据第一阶段接口计算出的workspaceSize申请device内存 if (workspaceSize > 0) { ret = aclrtMalloc(&workspaceAddr, workspaceSize, ACL_MEM_MALLOC_HUGE_FIRST); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtMalloc workspace failed. ret = %d \n", ret); return ret); } // 调用第二阶段接口 ret = aclnnAllGatherMatmulV2(workspaceAddr, workspaceSize, executor, args.stream); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclnnAllGatherMatmulV2 failed. ret = %d \n", ret); return ret); // (固定写法)同步等待任务执行结束 ret = aclrtSynchronizeStreamWithTimeout(args.stream, 10000); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtSynchronizeStreamWithTimeout failed. ret = %d \n", ret); return ret); LOG_PRINT("[INFO] device_%d aclnnAllGatherMatmulV2 execute successfully.\n", args.rankId); // 释放device资源,需要根据具体API的接口定义修改 if (x1 != nullptr) { aclDestroyTensor(x1); } if (x2 != nullptr) { aclDestroyTensor(x2); } if (x1Scale != nullptr) { aclDestroyTensor(x1Scale); } if (x2Scale != nullptr) { aclDestroyTensor(x2Scale); } if (bias != nullptr) { aclDestroyTensor(bias); } if (quantScale != nullptr) { aclDestroyTensor(quantScale); } if (out != nullptr) { aclDestroyTensor(out); } if (gatherOut != nullptr) { aclDestroyTensor(gatherOut); } if (amax != nullptr) { aclDestroyTensor(amax); } if (x1DeviceAddr != nullptr) { aclrtFree(x1DeviceAddr); } if (x2DeviceAddr != nullptr) { aclrtFree(x2DeviceAddr); } if (x1ScaleDeviceAddr != nullptr) { aclrtFree(x1ScaleDeviceAddr); } if (x2ScaleDeviceAddr != nullptr) { aclrtFree(x2ScaleDeviceAddr); } if (biasDeviceAddr != nullptr) { aclrtFree(biasDeviceAddr); } if (outDeviceAddr != nullptr) { aclrtFree(outDeviceAddr); } if (gatherOutDeviceAddr != nullptr) { aclrtFree(gatherOutDeviceAddr); } if (workspaceSize > 0) { aclrtFree(workspaceAddr); } ret = HcclCommDestroy(args.hcclComm); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] HcclCommDestroy failed. ret = %d \n", ret); return ret); ret = aclrtDestroyStream(args.stream); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtDestroyStream failed. ret = %d \n", ret); return ret); ret = aclrtResetDevice(args.rankId); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtResetDevice failed. ret = %d \n", ret); return ret); ret = aclrtDestroyContext(args.context); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtDestroyContext failed. ret = %d \n", ret); return ret); return 0; } int main(int argc, char *argv[]) { int ret = aclInit(nullptr); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclInit failed. ret = %d \n", ret); return ret); aclrtStream stream[DEV_NUM]; aclrtContext context[DEV_NUM]; for (uint32_t rankId = 0; rankId < DEV_NUM; rankId++) { ret = aclrtSetDevice(rankId); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtSetDevice failed. ret = %d \n", ret); return ret); ret = aclrtCreateContext(&context[rankId], rankId); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtCreateContext failed. ret = %d \n", ret); return ret); ret = aclrtCreateStream(&stream[rankId]); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] aclrtCreateStream failed. ret = %d \n", ret); return ret); } int32_t devices[DEV_NUM]; for (int i = 0; i < DEV_NUM; i++) { devices[i] = i; } // 初始化集合通信域 HcclComm comms[DEV_NUM]; ret = HcclCommInitAll(DEV_NUM, devices, comms); CHECK_RET(ret == ACL_SUCCESS, LOG_PRINT("[ERROR] HcclCommInitAll failed. ret = %d \n", ret); return ret); Args args[DEV_NUM]; // 启动多线程 std::vector<std::unique_ptr<std::thread>> threads(DEV_NUM); for (uint32_t rankId = 0; rankId < DEV_NUM; rankId++) { args[rankId].rankId = rankId; args[rankId].hcclComm = comms[rankId]; args[rankId].context = context[rankId]; args[rankId].stream = stream[rankId]; threads[rankId].reset(new(std::nothrow) std::thread(&LaunchOneThreadAllGatherMmV2, std::ref(args[rankId]))); } for (uint32_t rankId = 0; rankId < DEV_NUM; rankId++) { threads[rankId]->join(); } aclFinalize(); return 0; }