当前位置: 首页 > news >正文

自定义算子开发系列:TilingKey模板化编程介绍

基础知识准备

本文内容基于Ascend C算子开发衍生而来,对于算子开发还不了解的读者可以通过以下资源进行学习:
《Ascend C算子开发文档手册》:

https://www.hiascend.com/document/detail/zh/canncommercial/82RC1/opdevg/Ascendcopdevg/atlas_ascendc_10_0001.html

《Ascend C算子开发系列课程:入门》

https://www.hiascend.com/developer/courses/detail/1691696509765107713

《Ascend C算子开发系列课程:进阶》:

https://www.hiascend.com/developer/courses/detail/1696414606799486977

《Ascend C算子开发系列课程:高级》:

https://www.hiascend.com/developer/courses/detail/1696690858236694530

背景介绍

随着越来越多的用户基于Ascend C进行算子开发,算子开发过程中的易用性相关问题陆续被反馈出来,其中多场景算子的高效开发实现以及维护管理一直是用户的痛点。本次在CANN全面开源开放中做了很多算子开发易用性提升的工作,针对多场景算子的开发CANN推出了TilingKey模板化编程,有效解决了多场景算子开发管理的问题。

亮点介绍

TilingKey模板化编程具备以下特点:

  • 减少icache miss和scalar耗时,提升算子调用性能
  • 通过模板化的开发方式提升算子开发效率
  • 算子代码的可读性和可维护性提升

技术详细解析

在算子开发过程中,会遇到同一个算子有多种不同的实现的情形;为了能够在运行时能够正确加载到对应的算子实现,在Tiling的实现中会引入TilingKey的概念。TilingKey是一个算子内为了区分不同场景的实现而将Kernel代码进行区分的方法,该方法类似于C++的Template模板机制,可减少不必要的icache miss以及scalar耗时,有助于优化单次调用Kernel的性能。不同的Kernel实现分支可以通过TilingKey来标识,Host侧设置TilingKey后,可以选择对应的分支,例如:一个算子在不同的shape下,有不同的算法逻辑,Kernel侧可以通过TilingKey来选择不同的算法逻辑,在Host侧Tiling算法也有差异,Host/Kernel侧通过相同的TilingKey进行关联。

假如有如下Kernel代码:

if (condition) { ProcessA(); } else { ProcessB(); }

如果函数ProcessA、ProcessB两个函数是个非常大的函数,那么上述代码在编译后会变得更大,而每次Kernel运行只会选择1个分支,条件的判断和跳转在代码大到一定程度(16-32K,不同芯片存在差异)后会出现icache miss。通过TilingKey可以对这种情况进行优化,给2个Kernel的处理函数设置不同的TilingKey 1和2:

if (TILING_KEY_IS(1)) { ProcessA(); } elseif (TILING_KEY_IS(2)) { ProcessB(); }

这样Device Kernel编译时会自动识别到2个TilingKey并编译2个Kernel入口函数,将条件判断进行常量折叠。同时需要和Host Tiling函数配合,判断走ProcessA的场景设置TilingKey为1,走ProcessB的场景设置TilingKey为2:

static ge::graphStatus TilingFunc(gert::TilingContext* context) { // some code if (condition) { context->SetTilingKey(1); } else { context->SetTilingKey(2); } return ge::GRAPH_SUCCESS; }

但是,上述编程方式会让其他阅读这份代码的用户产生疑问:假如用户从函数调用入口查看代码时,看到了两个分支分别是TILING_KEY_IS(1)和TILING_KEY_IS(2),那么用户会无法理解“1”和“2”的含义。尤其是在涉及多个TilingKey的场景中,开发者依赖TilingKey来管理Kernel的实现,无论是在管理还是使用上都会遇到相当大的复杂性。因此为了简化这一过程,消除TilingKey不易于记忆和理解的问题,我们引入了一种采用模板编程的方法来替代传统的TilingKey编程的方案,从而减少对TilingKey数值标识的依赖,使Kernel的管理更加直观和高效。

具体的开发步骤如下:

下面我们以add自定义算子开发为样例,详细介绍一下具体的使用方法:

1.在自定义算子工程的op_kernel目录下,新增定义模板参数和模板参数组合的头文件,本示例中头文件命名为tiling_key_add_custom.h

该头文件中需要包含模板头文件ascendc/host_api/tiling/template_argument.h。

定义模板参数ASCENDC_TPL_ARGS_DECL和模板参数组合ASCENDC_TPL_ARGS_SEL(即可使用的模板)。具体API参考见模板参数定义。

模板参数定义:

https://www.hiascend.com/document/detail/zh/canncommercial/82RC1/API/ascendcopapi/atlasascendc_api_07_00011.html

#include"ascendc/host_api/tiling/template_argument.h" #define ADD_TPL_FP16 1 // 数据类型定义 #define ADD_TPL_FP32 0 #define ADD_TPL_ND 2 // 数据格式定义 #define ADD_TPL_NZ 29 // 模板参数 ASCENDC_TPL_ARGS_DECL(AddTemplateCustom, // 算子OpType ASCENDC_TPL_DTYPE_DECL(D_T_X, ADD_TPL_FP16, ADD_TPL_FP32), // DataType类型的模板参数定义:输入参数x的数据类型,取值范围为float16/float32 ASCENDC_TPL_DTYPE_DECL(D_T_Y, ADD_TPL_FP16, ADD_TPL_FP32), // DataType类型的模板参数定义:输入参数y的数据类型,取值范围为float16/float32 ASCENDC_TPL_DTYPE_DECL(D_T_Z, ADD_TPL_FP16, ADD_TPL_FP32), // DataType类型的模板参数定义:输入参数z的数据类型,取值范围为float16/float32 ASCENDC_TPL_UINT_DECL(TILE_NUM, ASCENDC_TPL_8_BW, ASCENDC_TPL_UI_MIX, 2, 0, 2, 3, 5, 10, 12, 13, 9, 8),// 自定义UINT类型(无符号整形)的模板参数定义:模板参数为切分的块数,编码位宽为ASCENDC_TPL_8_BW即8比特,表示该模板参数的个数不超过8比特能表达的范围;ASCENDC_TPL_UI_MIX表示通过混合模式表达取值范围,有2组的数据{0-2}、{3-5}和穷举值10、12、13、9、8,最后结果为{0, 1, 2, 3, 4, 5, 10, 12, 13, 9, 8} ASCENDC_TPL_BOOL_DECL(IS_SPLIT, 0, 1), // 自定义bool类型的模板参数定义:模板参数为是否切分标志位,取值范围为0和1,1表示切分,0表示不切分 ); // 模板参数组合 // 用于调用GET_TPL_TILING_KEY获取TilingKey时,接口内部校验TilingKey是否合法 ASCENDC_TPL_SEL( ASCENDC_TPL_ARGS_SEL( ASCENDC_TPL_DTYPE_SEL(D_T_X, ADD_TPL_FP16), ASCENDC_TPL_DTYPE_SEL(D_T_Y, ADD_TPL_FP16), ASCENDC_TPL_DTYPE_SEL(D_T_Z, ADD_TPL_FP16), ASCENDC_TPL_UINT_SEL(TILE_NUM, ASCENDC_TPL_UI_LIST, 1, 8), ASCENDC_TPL_BOOL_SEL(IS_SPLIT, 0, 1), ), ASCENDC_TPL_ARGS_SEL( ASCENDC_TPL_DTYPE_SEL(D_T_X, ADD_TPL_FP32), ASCENDC_TPL_DTYPE_SEL(D_T_Y, ADD_TPL_FP32), ASCENDC_TPL_DTYPE_SEL(D_T_Z, ADD_TPL_FP32), ASCENDC_TPL_UINT_SEL(TILE_NUM, ASCENDC_TPL_UI_LIST, 1, 8), ASCENDC_TPL_BOOL_SEL(IS_SPLIT, 0, 1), ), );

2.Host侧调用GET_TPL_TILING_KEY接口生成TilingKey

GET_TPL_TILING_KEY:

https://www.hiascend.com/document/detail/zh/canncommercial/82RC1/API/ascendcopapi/atlas_ascendc_10_00008.html

Host实现文件中包含步骤1中定义模板参数和模板参数组合的头文件。

调用GET_TPL_TILING_KEY接口生成TilingKey,GET_TPL_TILING_KEY输入参数为模板参数的具体值,传入时需要与定义模板参数和模板参数组合的头文件中的模板参数顺序保持一致。

#include "tiling_key_add_custom.h" static ge::graphStatus TilingFunc(gert::TilingContext *context) { TilingData tiling; uint32_t totalLength = context->GetInputShape(0)->GetOriginShape().GetShapeSize(); ge::DataType dtype_x = context->GetInputDesc(0)->GetDataType(); ge::DataType dtype_y = context->GetInputDesc(1)->GetDataType(); ge::DataType dtype_z = context->GetOutputDesc(1)->GetDataType(); uint32_t D_T_X = ADD_TPL_FP32, D_T_Y=ADD_TPL_FP32, D_T_Z=ADD_TPL_FP32, TILE_NUM=1, IS_SPLIT=0; if(dtype_x == ge::DataType::DT_FLOAT){ D_T_X = ADD_TPL_FP32; }elseif(dtype_x == ge::DataType::DT_FLOAT16){ D_T_X = ADD_TPL_FP16; } if(dtype_y == ge::DataType::DT_FLOAT){ D_T_Y = ADD_TPL_FP32; }elseif(dtype_y == ge::DataType::DT_FLOAT16){ D_T_Y = ADD_TPL_FP16; } if(dtype_z == ge::DataType::DT_FLOAT){ D_T_Z = ADD_TPL_FP32; }elseif(dtype_z == ge::DataType::DT_FLOAT16){ D_T_Z = ADD_TPL_FP16; } if(totalLength< MIN_LENGTH_FOR_SPLIT){ IS_SPLIT = 0; TILE_NUM = 1; }else{ IS_SPLIT = 1; TILE_NUM = DEFAULT_TILE_NUM; } context->SetBlockDim(BLOCK_DIM); tiling.set_totalLength(totalLength); tiling.SaveToBuffer(context->GetRawTilingData()->GetData(), context->GetRawTilingData()->GetCapacity()); context->GetRawTilingData()->SetDataSize(tiling.GetDataSize()); constuint64_t tilingKey = GET_TPL_TILING_KEY(D_T_X, D_T_Y, D_T_Z, TILE_NUM, IS_SPLIT); context->SetTilingKey(tilingKey); size_t *currentWorkspace = context->GetWorkspaceSizes(1); currentWorkspace[0] = 0; return ge::GRAPH_SUCCESS; }
  • 3.Kernel侧实现
  • Kernel实现文件中包含步骤1中定义模板参数和模板参数组合的头文件。
  • 核函数添加template模板,以便支持模板参数的传入,参数顺序需要与定义模板参数和模板参数组合的头文件中的模板参数顺序保持一致。
  • 通过对模板参数的分支判断,选择不同的Kernel侧实现。
#include "tiling_key_add_custom.h" ... ... template<int D_T_X, int D_T_Y, int D_T_Z, int TILE_NUM, int IS_SPLIT> __global__ __aicore__ void add_custom(GM_ADDR x, GM_ADDR y, GM_ADDR z, GM_ADDR workspace, GM_ADDR tiling) { GET_TILING_DATA(tiling_data, tiling); if(D_T_X == ADD_TPL_FP32 && D_T_Y == ADD_TPL_FP32 && D_T_Z == ADD_TPL_FP32){ KernelAdd<float, float, float> op; op.Init(x, y, z, tiling_data.totalLength, TILE_NUM, IS_SPLIT); op.Process1(); }else if(D_T_X == ADD_TPL_FP16 && D_T_Y == ADD_TPL_FP16 && D_T_Z == ADD_TPL_FP16){ KernelAdd<half, half, half> op; if(IS_SPLIT == 0){ op.Init(x, y, z, tiling_data.totalLength, TILE_NUM, IS_SPLIT); op.Process1(); }else if(IS_SPLIT==1){ op.Init(x, y, z, tiling_data.totalLength, TILE_NUM, IS_SPLIT); op.Process2(); } }}

当前已开源的ops-nn仓和ops-transfomer仓中的部分核心算子已经基于该模板化编程方式完成优化,感兴趣的同学可以访问ops-nn仓 BatMatMulV3算子和ops-transformer仓 FA算子。

ops-nn仓 BatMatMulV3算子:

https://gitcode.com/cann/ops-nn/tree/master/matmul/batch_mat_mul_v3

ops-transformer仓 FA算子:

https://gitcode.com/cann/ops-transformer/tree/master/attention/flash_attention_score

总结

TilingKey模板化编程既提供了一种简单的多场景算子开发编程范式,简化了算子的开发难度,也实际提升了算子的执行效率,帮助开发者更便捷的开发高性能多场景算子。

​​​​

http://www.cnnetsun.cn/news/1590.html

相关文章:

  • 双馈风机并网在Matlab/Simulink中的高效仿真实现
  • 模糊控制SOC主动均衡:提升电池性能的新途径
  • Flac3d 渗流模拟:探索建筑物对地基孔隙水压力及固结沉降的影响
  • 探索自动驾驶量产核心:车道居中 LKA/LCA 控制算法及 Matlab Simulink 模型
  • 探索基于肺性能的优化(LPO)算法:2024年的肺功能优化新突破
  • ACL配置实验
  • 并联混合动力汽车ECMS能量管理策略:基于SIMULINK的实现与探索
  • Comsol电磁超声导波检测之旅
  • Comsol多孔介质流固耦合:探索孔压与位移的时空奥秘
  • 交换机端口异常的排查流程:从物理层到配置层的完整分析
  • 你真的融会贯通了 javascript 中的异步编程了吗?
  • 基于EKF扩展卡尔曼滤波的PMSM无传感器控制探索
  • 用改进 CPO 优化 CNN - BiGRU - Attention 算法模型实现多变量时间序列分类预测(Matlab 版)
  • 基于滑模控制的异步电机直接转矩控制:提升转速跟踪精度
  • 探索AB胶刮胶机的技术实现:基于C#上位机与研华运动控制卡
  • 无人船目标跟踪控制:NMPC 与 PPO 强化学习的对决
  • 探索三相并网仿真中的准PR控制与电容电流反馈
  • 基于MATLAB与Simulink的单相重合闸方式仿真探索
  • COMSOL超声相控阵聚焦仿真:开启瞬态参数探索之旅
  • 城市更新新范式:蓝色星球为城市安装“数字操作系统”
  • 电池防反接电路
  • 【路径规划】基于RRT快速探索随机树算法在含连续曲面障碍物空间中,从起点到目标点生成无碰撞路径附Matlab代码
  • 【图像增强】基于频域处理、对数变换、色彩恢复、多尺度 Retinex(MSRCR)思想、直流分量优化和对比度增强的视频增强算法附Matlab代码
  • 【能源系统】温室效应与光谱吸收附matlab仿真
  • 【卫星信号】基于信噪比SNR的全球导航卫星系统反射测量GPSGNSS-R附matlab代码
  • 【翼型】确定波音787飞机和F-16战斗猎鹰的着陆速度、性能和稳定性特征研究附Matlab代码
  • 【数据聚类】灰狼算法优化Kmeans亚洲足球水平聚类分析附Matlab代码
  • 基于遗传算法GA、模拟退火SA、粒子群PSO求解港口泊位分配与岸桥调度一体化优化调度,最小化船舶在港总停留时间附Matlab代码
  • 【SLAM】基于扩展卡尔曼滤波(EKF)的移动机器人同时定位与地图构建(SLAM)Matlab仿真程序
  • 【Linux C/C++开发】Docker在嵌入式Linux交叉编译中的完整应用方案