深度学习编译器之Layerout Transform优化
在本文的描述中，存在一些接口和interface的混用，这两个是一样的都表示mlir的interface。
0x0. 背景
继续深度学习编译器的优化工作解读，本篇文章要介绍的是oneflow系统中如何基于mlir实现layerout transform。在2d卷积神经网络中，除了nchw数据格式之外一般还存在nhwc的数据格式，对于卷积操作来说使用nhwc格式进行计算可能会获得更好的性能。
但深度学习网络的训练一般来说是采用nchw进行的，我们一般只有在推理时才做nchw到nhwc的layerout transform。这里存在两个问题：首先对于一个算子比如conv2d，它以nchw方式训练时保存的权重格式是[out_channels, in_channels, *kernel_size]，但是要以nhwc格式进行推理时我们需要对权重的格式进行转换；然后对于没有权重的算子来说，我们也需要尽量的让算子支持nhwc的运算，来减少因为卷积算子前后插入的transpose操作带来的额外开销。举个例子，假设有如下的一个小网络 x->conv->relu->conv->relu->out，如果我们要以nhwc格式执行那么我们除了对2个卷积的权重进行改动之外，我们还需要在conv前后插入transpose来修改输入到conv算子的数据格式，也就是x->transpose(0, 2, 3, 1)->conv->transpose(0, 3, 1, 2) -> relu -> transpose(0, 2, 3, 1)->conv->transpose(0, 3, 1, 2) -> relu->out。然后细心的读者可以发现，实际上这里存在很多冗余的transpose，因为relu是支持以nhwc格式进行运算的，那么这个网络可以化简为x->transpose(0, 2, 3, 1)->conv->relu->conv->relu->transpose(0, 3, 1, 2)->out。这样可以减少一半的transpose op开销。
之所以要做transpose的化简是因为transpose算子本身也有运行以及调度的开销，如果我们不尽量减少transpose的个数，那么因为改用nhwc带来的计算加速可能会被 transpose 的开销掩盖住。我们基于oneflow实现了上述的layerout transform优化，以下给出测试结果。
在v100上对这个优化进行了测试，测试代码见 https://github.com/oneflow-inc/oneflow/blob/master/oneflow/ir/test/oneflow/auto_nhwc/test_resnet101_benchmark.py ，性能结果如下：
开启nn.graph的amp选项。
网络选取resnet101，对其做前向推理。
batch_size nchw auto nhwc
16 14s 13s
32 24s 22s
64 44s 38s
在batchsize=64时得到了13.6%的加速，随着batchsize减少加速比会减小，但始终会保持一些加速。需要注意的是，这里对权重参数部分提前进行了transpose，所以这部分是没有额外开销的。实际上，我们采用了常量折叠的方式来完成，这个下篇文章再讲。
0x1. 实现解析
在实现上主要需要搞定3个问题，第一个是如何确定哪些算子支持nhwc的运算，第二个是插入transpose算子，第三个是消除多余的transpose对。
0x1.1 基于interface确定哪些算子支持nhwc运算
在oneflow中如果我们要让某个op支持nhwc的计算，只需在op定义时声明一个nchwcompatibleinterface。以卷积为例：
def oneflow_conv2dop : oneflow_convolutionbaseop {}
这里的 declareopinterfacemethods 表示这个 operator 实现了 nchwcompatibleinterface 接口,该接口定义了与 nchw 格式兼容的 operator 需要实现的方法。
我们想让其它的任意 op 支持 nhwc 的运算，只需要定义这个接口并且重写这个接口的成员函数即可，接下来我们看一下nchwcompatibleinterface 的定义。
def nchwcompatibleinterface : opinterface { let description = [{ interface of nchw compatibility }]; let methods = [ interfacemethod, interfacemethod, interfacemethod, interfacemethod, ]; let cppnamespace = ::oneflow;}
这个接口继承自 opinterface 接口, opinterface 是 mlir 框架中描述 operator interface 的基类。nchwcompatibleinterface 表示一个与 nchw 格式兼容的 operator interface。nchwcompatibleinterface定义了几个方法：
isnchw: 返回一个 bool 值, 表示当前的 operator 在什么条件下是处理输入为 nchw 格式的数据。
nchwtonhwc: 接受 transpose 后的输入和重写器 (rewriter), 用于从 nchw 格式转换为 nhwc 格式。
operandstotranspose: 返回需要 transpose 的输入值集合。
resultstotranspose:返回需要 transpose 的输出值集合。
接下来我们看一下conv2d op对应的 nchwcompatibleinterface 接口实现：
bool conv2dop::isnchw() { return this->getdataformat().str() == channels_first; }llvm::denseset conv2dop::operandstotranspose() { if (this->get_addtooutput()) { return {this->getin(), this->getweight(), this->get_addtooutput()}; } else { return {this->getin(), this->getweight()}; }}llvm::denseset conv2dop::resultstotranspose() { return {this->getout()}; }llvm::smallvector conv2dop::nchwtonhwc(llvm::smallvector value, patternrewriter& rewriter) { auto conv_op = *this; smallvector operands; operands.push_back(value[0]); operands.push_back(value[1]); if (conv_op.getbias()) operands.push_back(conv_op.getbias()); if (this->get_addtooutput()) { operands.push_back(value[2]); } namedattrlist attributes = conv_op->getattrs(); attributes.set(conv_op.getdataformatattrname(), rewriter.getstringattr(channels_last)); auto res = rewriter .create(conv_op.getloc(), getnhwcresulttypes(conv_op), operands, attributes) ->getresults(); llvm::smallvector results; results.push_back(res[0]); return results;}
其中，isnchw 方法返回一个 bool 值,表示该 conv2dop operation 是否使用 nchw 格式。它通过检查 operation 的data_format 属性来判断。operandstotranspose 方法返回需要 transpose 的输入值集合。对于 conv2dop 来说,主要输入包括input、weight、bias(可选) 和 addto_output(可选)，其中bias不需要 transpose，并且这个addto_output是oneflow的一个特殊的输出用来做算子融合读者可以忽略。
resultstotranspose 方法返回需要 transpose 的输出值集合。对于 conv2dop 来说,仅有一个输出, 所以返回输出特征图的值。nchwtonhwc 方法接受 nchw 格式的输入值和重写器,并返回 nhwc 格式的结果值。它通过创建一个新的 conv2dop operation, 并将 data_format 属性设置为 channels_last, 来实现从 nchw 到 nhwc 的转换。
0x1.2 插入transpose算子
接下来就是贪心的给网络里的算子插入transpose算子，这里的思路是我们尽可能的对网络里面的所有算子都前后分别插入一个transpose，这样的话在消除transopose对的时候才能获得最优的解。给网络中的算子插入transpose的逻辑如下面的pattern代码所述：
struct autonhwcpattern : public opinterfacerewritepattern { explicit autonhwcpattern(mlir::mlircontext* context) : opinterfacerewritepattern(context, /*benefit=*/1) {} public: logicalresult matchandrewrite(nchwcompatible op, patternrewriter& rewriter) const override { if (op->hastrait()) { for (mlir::value operand : op.operandstotranspose()) { if (operand.gettype().cast().getshape().size() != 4) { return failure(); } } const auto device_name = optrait::isopconfcompatible::getdevicetag(op) .cast() .getvalue() .str(); if (device_name == cpu) { return failure(); } } llvm::smallvector perm = getchannellasttransposeperm(); llvm::smallvector result_perm = getchannelfirsttransposeperm(); namedattrlist transpose_attributes; if (inittransposeattributes(op, transpose_attributes, rewriter).succeeded()) { transpose_attributes.append(llvm::stringref(perm), getsi32arrayattr(rewriter, perm)); } else { return failure(); } // when op op has no sense of data_format and pre op is transpose, we greedily insert transpose // into this op, seeking more opportunities to eliminate transpose pattern. const bool greedily_transpose_flag = !op.isnchw() && isinserttransposeopbefore(op, rewriter); if (op.isnchw() || greedily_transpose_flag) { // create transpose op for input operand smallvector tranposed_operands; llvm::denseset operand_transpose = op.operandstotranspose(); int num_transposed_operand = 0; for (value operand : op->getoperands()) { if (operand_transpose.find(operand) != operand_transpose.end()) { smallvector input_res = getinputoperandtransposeop( op, operand, transpose_attributes, num_transposed_operand, rewriter); tranposed_operands.push_back(input_res[0]); num_transposed_operand += 1; } } // create nhwc op smallvector created_results = op.nchwtonhwc(tranposed_operands, rewriter); // create transpose op for results int num_transposed_result = 0; transpose_attributes.set(llvm::stringref(perm), getsi32arrayattr(rewriter, result_perm)); llvm::denseset transpose_result = op.resultstotranspose(); for (value result : op->getopresults()) { if (transpose_result.find(result) != transpose_result.end()) { if (auto result_transpose_op = getresulttransposeop(op, created_results[num_transposed_result], transpose_attributes, num_transposed_result, rewriter)) { result.replacealluseswith(result_transpose_op); num_transposed_result += 1; } else { return failure(); } } } } return success(); }};
首先 autonhwcpattern 类继承自 opinterfacerewritepattern，opinterfacerewritepattern 是一个用于重写 operation 的基类。autonhwcpattern 针对实现了 nchwcompatible interface 的 operation 进行重写,以实现 nchw 到 nhwc 的格式转换。然后，autonhwcpattern 重写了 matchandrewrite 方法。该方法会在遇到 nchwcompatible interface 的 operation 时被调用,来实现从 nchw 到 nhwc 的转换。接下来，matchandrewrite 方法首先会检查 operation 是否满足转换条件,如是否 4 维、是否在 cpu 设备上等。
如果不满足则返回 failure。如果满足, matchandrewrite 方法会获取 nchw 到nhwc 和 nhwc 到 nchw 的转换顺序。并初始化 transpose operation 的属性。然后对于当前 op 是 nchw 格式或者这个 op 的前一个 op 是transpose op，这里都进行插入 transpose op的操作来获得更多的优化机会。
这里还涉及到几个相关的工具函数，我们也解释一下：
llvm::smallvector getchannellasttransposeperm() { return {0, 2, 3, 1}; }llvm::smallvector getchannelfirsttransposeperm() { return {0, 3, 1, 2}; }llvm::smallvector getinputoperandtransposeop(nchwcompatible op, value val, namedattrlist transpose_attributes, int num_transposed_operand, patternrewriter& rewriter) { std::string transpose_name = optrait::isopconfcompatible::getopname(op).str() + _transpose_input_ + std::to_string(num_transposed_operand); transpose_attributes.set(llvm::isopconfcompatible::getopnameattr()), rewriter.getstringattr(transpose_name)); smallvector input_operands; input_operands.push_back(val); auto res = rewriter .create(op.getloc(), getnhwctype(val.gettype()), input_operands, transpose_attributes) ->getresults(); return res;}transposeop getresulttransposeop(nchwcompatible op, value val, namedattrlist transpose_attributes, int num_transposed_result, patternrewriter& rewriter) { std::string transpose_name = optrait::isopconfcompatible::getopname(op).str() + _transpose_output_ + std::to_string(num_transposed_result); transpose_attributes.set(llvm::isopconfcompatible::getopnameattr()), rewriter.getstringattr(transpose_name)); smallvector operands; operands.push_back(val); transposeop transpose_op = rewriter.create( op.getloc(), getnchwtype(val.gettype()), operands, transpose_attributes); return transpose_op;}bool isinserttransposeopbefore(nchwcompatible op, patternrewriter& rewriter) { bool insert_transpose_op_flag = false; for (mlir::value operand : op->getoperands()) { transposeop transposeinputop = operand.getdefiningop(); if (!transposeinputop) continue; const auto perm = transposeinputop.getpermattr(); if (perm.size() == 4 && perm[0] == rewriter.getsi32integerattr(0) && perm[1] == rewriter.getsi32integerattr(3) && perm[2] == rewriter.getsi32integerattr(1) && perm[3] == rewriter.getsi32integerattr(2)) { insert_transpose_op_flag = true; break; } } return insert_transpose_op_flag;}
其中 getchannellasttransposeperm 和 getchannelfirsttransposeperm 方法分别返回 nhwc 到 nchw 和 nchw 到nhwc 的转换顺序。getinputoperandtransposeop 方法为 operation 的输入创建一个transpose operation。它使用输入值、transpose属性 和 重写器创建一个 transposeop , 并返回其结果。
类似的，getresulttransposeop 方法为 operation 的输出创建一个transpose operation。它使用输出值、transpose属性和重写器创建一个transposeop,并返回该operation。isinserttransposeopbefore方法检查operation的输入是否已有 transpose operation。如果有,并且该 transpose operation 将 nhwc 转为 nchw, 则返回 true, 否则返回false。
0x1.3 消除多余的transpose对
接下来，我们需要把插入transpose op的图中所有相邻的transpose对尽可能的消除，代码实现如下：
bool isredundanttransposematch(arrayattr pre, arrayattr afe, mlir::patternrewriter& rewriter) { const auto preperm = pre.getvalue().vec(); const auto afeperm = afe.getvalue().vec(); if (preperm.size() == 4 && afeperm.size() == 4) { // handle nchw->nhwc->nchw: (0, 2, 3, 1) -> (0, 3, 1, 2) if (preperm[0] == afeperm[0] && preperm[1] == afeperm[3] && preperm[2] == afeperm[1] && preperm[3] == afeperm[2] && preperm[0] == rewriter.getsi32integerattr(0) && preperm[1] == rewriter.getsi32integerattr(2) && preperm[2] == rewriter.getsi32integerattr(3) && preperm[3] == rewriter.getsi32integerattr(1)) return true; // handle nhwc->nchw->nhwc: (0, 3, 1, 2) -> (0, 2, 3, 1) if (preperm[0] == afeperm[0] && preperm[1] == afeperm[2] && preperm[2] == afeperm[3] && preperm[3] == afeperm[1] && preperm[0] == rewriter.getsi32integerattr(0) && preperm[1] == rewriter.getsi32integerattr(3) && preperm[2] == rewriter.getsi32integerattr(1) && preperm[3] == rewriter.getsi32integerattr(2)) return true; } return false;}struct autonhwceliminateredundanttransposepattern : public mlir::oprewritepattern { explicit autonhwceliminateredundanttransposepattern(mlir::mlircontext* context) : oprewritepattern(context, /*benefit=*/1) {} mlir::logicalresult matchandrewrite(transposeop op, mlir::patternrewriter& rewriter) const override { mlir::value transposeinput = op.getoperand(); transposeop transposeinputop = transposeinput.getdefiningop(); if (!transposeinputop || !isredundanttransposematch(op.getpermattr(), transposeinputop.getpermattr(), rewriter)) { return failure(); } rewriter.replaceop(op, {transposeinputop.getoperand()}); return success(); }};
isredundanttransposematch 方法检查两个 transpose operation的顺序是否会导致冗余。它通过比较两个 transpose 的 perm 属性来判断。类似 autonhwcpattern ，autonhwceliminateredundanttransposepattern 类继承自 oprewritepattern 。它对transposeop 进行重写以实现 transpose 消除。如果顺序是 nhwc->nchw->nhwc 或nchw->nhwc->nchw , 则判定为冗余 transpose 。
如果输入也来自transposeop且两个 transpose 顺序导致冗余,matchandrewrite方法会用transposeop的输入替换transposeop。实现 transpose 消除。matchandrewrite 方法首先获取 transposeop 的输入,并检查该输入是否也来自一个 transposeop。如果不是, 或两个 transpose 的顺序不导致冗余, 则返回 failure。最后返回 success 表示成功消除冗余 transpose 。
最终，上面介绍的2个pass都被封装到 autonhwcpass 中作用在 mlir 的计算图上完成全局优化。从下面的代码可以看到这个优化只有在打开 oneflow_mlir_prefer_nhwc 环境变量时才正常生效。
void populateautonhwcpatterns(::rewritepatternset& patterns) { bool enable_nhwc = ::parsebooleanfromenv(oneflow_mlir_prefer_nhwc, false); if (enable_nhwc) { patterns.add(patterns.getcontext()); patterns.add(patterns.getcontext()); }}class autonhwcpass : public autonhwcpassbase { void runonoperation() override { operation* op = getoperation(); rewritepatternset patterns(op->getcontext()); oneflow::populateautonhwcpatterns(patterns); (void)applypatternsandfoldgreedily(op, std::move(patterns)); }};
补充：0x1.4 weight的transpose消除
这里还需要粗略的说明一下对于 weight 的 transpose 是如何处理的。在0x1.2中我们为 weight（常量constant op） 也插入了 transpose op，然后我们知道 weight 是常量，所以针对 weight 的 transpose op 完全可以在编译期折叠起来。这个过程是在 https://github.com/oneflow-inc/oneflow/blob/master/oneflow/ir/oneflow-translate/lib/oneflow/mlironeflowtranslation.cpp#l808-l811 这里完成的，我们后面会单独介绍一下 constant folding 的实现。
0x2. 结论
本文介绍了一下oneflow的编译器中的 layerout transform，这个技术在后来 oneflow 版本的 stable diffusion 中也发挥了重要作用，提升了推理速度。在 tvm 的 ansor 中也有类似的优化，通过将不同的 layerout 设定为 op 的 strategy 进而影响 op 的 schedule，在搜索的时候考虑到 layerout transform 来获得更大的搜索空间和更好的结果。在处理transpose 额外开销的方法并不是唯一的，这里只是采用了一种个人认为比较简单的方式，读者们如果有类似需要可以自由发挥。
基站蓄电池组远程在线监测解决方案
西部数据SN750带散热片版本的优缺点
迎接3D IC时代 半导体供应链加速投入研发
当物联网遇上人工智能时，成为创新的超能力
安卓8.0最新消息：安卓8.0初体验，竟然这么流畅还有这么多黑科技功能
深度学习编译器之Layerout Transform优化
一种基于人工智能的脑图像处理方法
大众CEO：加速业务改革，不要成为诺基亚
魅族pro7什么时候上市?魅族pro7最新消息:魅族pro7终极曝光:价格、配置抢先看,外观设计一流你不期待?
磁阻电机的工作原理 磁阻电机的优缺点
电流互感器和电压互感器的极性及测定方法
罗姆推出业界最小晶体管封装“VML0806”
柠檬豆：把握采购方的真需求就是我们的追求
多速率转换信号频谱研究
无线电发射功率与发射距离有什么关系
西陇科学与宝洁公司在西陇科学总部举行了合作20周年纪念会议
贴片电容的基本原理和功能 贴片电容在直流电路中的应用
中国科技巨头实现弯道超车，发布全球最强5G芯片
通过引导分区破解树莓派的6种方法
星恒电源对锂电行业发展的思考与经验
0x0. 背景
继续深度学习编译器的优化工作解读，本篇文章要介绍的是oneflow系统中如何基于mlir实现layerout transform。在2d卷积神经网络中，除了nchw数据格式之外一般还存在nhwc的数据格式，对于卷积操作来说使用nhwc格式进行计算可能会获得更好的性能。
但深度学习网络的训练一般来说是采用nchw进行的，我们一般只有在推理时才做nchw到nhwc的layerout transform。这里存在两个问题：首先对于一个算子比如conv2d，它以nchw方式训练时保存的权重格式是[out_channels, in_channels, *kernel_size]，但是要以nhwc格式进行推理时我们需要对权重的格式进行转换；然后对于没有权重的算子来说，我们也需要尽量的让算子支持nhwc的运算，来减少因为卷积算子前后插入的transpose操作带来的额外开销。举个例子，假设有如下的一个小网络 x->conv->relu->conv->relu->out，如果我们要以nhwc格式执行那么我们除了对2个卷积的权重进行改动之外，我们还需要在conv前后插入transpose来修改输入到conv算子的数据格式，也就是x->transpose(0, 2, 3, 1)->conv->transpose(0, 3, 1, 2) -> relu -> transpose(0, 2, 3, 1)->conv->transpose(0, 3, 1, 2) -> relu->out。然后细心的读者可以发现，实际上这里存在很多冗余的transpose，因为relu是支持以nhwc格式进行运算的，那么这个网络可以化简为x->transpose(0, 2, 3, 1)->conv->relu->conv->relu->transpose(0, 3, 1, 2)->out。这样可以减少一半的transpose op开销。
之所以要做transpose的化简是因为transpose算子本身也有运行以及调度的开销，如果我们不尽量减少transpose的个数，那么因为改用nhwc带来的计算加速可能会被 transpose 的开销掩盖住。我们基于oneflow实现了上述的layerout transform优化，以下给出测试结果。
在v100上对这个优化进行了测试，测试代码见 https://github.com/oneflow-inc/oneflow/blob/master/oneflow/ir/test/oneflow/auto_nhwc/test_resnet101_benchmark.py ，性能结果如下：
开启nn.graph的amp选项。
网络选取resnet101，对其做前向推理。
batch_size nchw auto nhwc
16 14s 13s
32 24s 22s
64 44s 38s
在batchsize=64时得到了13.6%的加速，随着batchsize减少加速比会减小，但始终会保持一些加速。需要注意的是，这里对权重参数部分提前进行了transpose，所以这部分是没有额外开销的。实际上，我们采用了常量折叠的方式来完成，这个下篇文章再讲。
0x1. 实现解析
在实现上主要需要搞定3个问题，第一个是如何确定哪些算子支持nhwc的运算，第二个是插入transpose算子，第三个是消除多余的transpose对。
0x1.1 基于interface确定哪些算子支持nhwc运算
在oneflow中如果我们要让某个op支持nhwc的计算，只需在op定义时声明一个nchwcompatibleinterface。以卷积为例：
def oneflow_conv2dop : oneflow_convolutionbaseop {}
这里的 declareopinterfacemethods 表示这个 operator 实现了 nchwcompatibleinterface 接口,该接口定义了与 nchw 格式兼容的 operator 需要实现的方法。
我们想让其它的任意 op 支持 nhwc 的运算，只需要定义这个接口并且重写这个接口的成员函数即可，接下来我们看一下nchwcompatibleinterface 的定义。
def nchwcompatibleinterface : opinterface { let description = [{ interface of nchw compatibility }]; let methods = [ interfacemethod, interfacemethod, interfacemethod, interfacemethod, ]; let cppnamespace = ::oneflow;}
这个接口继承自 opinterface 接口, opinterface 是 mlir 框架中描述 operator interface 的基类。nchwcompatibleinterface 表示一个与 nchw 格式兼容的 operator interface。nchwcompatibleinterface定义了几个方法：
isnchw: 返回一个 bool 值, 表示当前的 operator 在什么条件下是处理输入为 nchw 格式的数据。
nchwtonhwc: 接受 transpose 后的输入和重写器 (rewriter), 用于从 nchw 格式转换为 nhwc 格式。
operandstotranspose: 返回需要 transpose 的输入值集合。
resultstotranspose:返回需要 transpose 的输出值集合。
接下来我们看一下conv2d op对应的 nchwcompatibleinterface 接口实现：
bool conv2dop::isnchw() { return this->getdataformat().str() == channels_first; }llvm::denseset conv2dop::operandstotranspose() { if (this->get_addtooutput()) { return {this->getin(), this->getweight(), this->get_addtooutput()}; } else { return {this->getin(), this->getweight()}; }}llvm::denseset conv2dop::resultstotranspose() { return {this->getout()}; }llvm::smallvector conv2dop::nchwtonhwc(llvm::smallvector value, patternrewriter& rewriter) { auto conv_op = *this; smallvector operands; operands.push_back(value[0]); operands.push_back(value[1]); if (conv_op.getbias()) operands.push_back(conv_op.getbias()); if (this->get_addtooutput()) { operands.push_back(value[2]); } namedattrlist attributes = conv_op->getattrs(); attributes.set(conv_op.getdataformatattrname(), rewriter.getstringattr(channels_last)); auto res = rewriter .create(conv_op.getloc(), getnhwcresulttypes(conv_op), operands, attributes) ->getresults(); llvm::smallvector results; results.push_back(res[0]); return results;}
其中，isnchw 方法返回一个 bool 值,表示该 conv2dop operation 是否使用 nchw 格式。它通过检查 operation 的data_format 属性来判断。operandstotranspose 方法返回需要 transpose 的输入值集合。对于 conv2dop 来说,主要输入包括input、weight、bias(可选) 和 addto_output(可选)，其中bias不需要 transpose，并且这个addto_output是oneflow的一个特殊的输出用来做算子融合读者可以忽略。
resultstotranspose 方法返回需要 transpose 的输出值集合。对于 conv2dop 来说,仅有一个输出, 所以返回输出特征图的值。nchwtonhwc 方法接受 nchw 格式的输入值和重写器,并返回 nhwc 格式的结果值。它通过创建一个新的 conv2dop operation, 并将 data_format 属性设置为 channels_last, 来实现从 nchw 到 nhwc 的转换。
0x1.2 插入transpose算子
接下来就是贪心的给网络里的算子插入transpose算子，这里的思路是我们尽可能的对网络里面的所有算子都前后分别插入一个transpose，这样的话在消除transopose对的时候才能获得最优的解。给网络中的算子插入transpose的逻辑如下面的pattern代码所述：
struct autonhwcpattern : public opinterfacerewritepattern { explicit autonhwcpattern(mlir::mlircontext* context) : opinterfacerewritepattern(context, /*benefit=*/1) {} public: logicalresult matchandrewrite(nchwcompatible op, patternrewriter& rewriter) const override { if (op->hastrait()) { for (mlir::value operand : op.operandstotranspose()) { if (operand.gettype().cast().getshape().size() != 4) { return failure(); } } const auto device_name = optrait::isopconfcompatible::getdevicetag(op) .cast() .getvalue() .str(); if (device_name == cpu) { return failure(); } } llvm::smallvector perm = getchannellasttransposeperm(); llvm::smallvector result_perm = getchannelfirsttransposeperm(); namedattrlist transpose_attributes; if (inittransposeattributes(op, transpose_attributes, rewriter).succeeded()) { transpose_attributes.append(llvm::stringref(perm), getsi32arrayattr(rewriter, perm)); } else { return failure(); } // when op op has no sense of data_format and pre op is transpose, we greedily insert transpose // into this op, seeking more opportunities to eliminate transpose pattern. const bool greedily_transpose_flag = !op.isnchw() && isinserttransposeopbefore(op, rewriter); if (op.isnchw() || greedily_transpose_flag) { // create transpose op for input operand smallvector tranposed_operands; llvm::denseset operand_transpose = op.operandstotranspose(); int num_transposed_operand = 0; for (value operand : op->getoperands()) { if (operand_transpose.find(operand) != operand_transpose.end()) { smallvector input_res = getinputoperandtransposeop( op, operand, transpose_attributes, num_transposed_operand, rewriter); tranposed_operands.push_back(input_res[0]); num_transposed_operand += 1; } } // create nhwc op smallvector created_results = op.nchwtonhwc(tranposed_operands, rewriter); // create transpose op for results int num_transposed_result = 0; transpose_attributes.set(llvm::stringref(perm), getsi32arrayattr(rewriter, result_perm)); llvm::denseset transpose_result = op.resultstotranspose(); for (value result : op->getopresults()) { if (transpose_result.find(result) != transpose_result.end()) { if (auto result_transpose_op = getresulttransposeop(op, created_results[num_transposed_result], transpose_attributes, num_transposed_result, rewriter)) { result.replacealluseswith(result_transpose_op); num_transposed_result += 1; } else { return failure(); } } } } return success(); }};
首先 autonhwcpattern 类继承自 opinterfacerewritepattern，opinterfacerewritepattern 是一个用于重写 operation 的基类。autonhwcpattern 针对实现了 nchwcompatible interface 的 operation 进行重写,以实现 nchw 到 nhwc 的格式转换。然后，autonhwcpattern 重写了 matchandrewrite 方法。该方法会在遇到 nchwcompatible interface 的 operation 时被调用,来实现从 nchw 到 nhwc 的转换。接下来，matchandrewrite 方法首先会检查 operation 是否满足转换条件,如是否 4 维、是否在 cpu 设备上等。
如果不满足则返回 failure。如果满足, matchandrewrite 方法会获取 nchw 到nhwc 和 nhwc 到 nchw 的转换顺序。并初始化 transpose operation 的属性。然后对于当前 op 是 nchw 格式或者这个 op 的前一个 op 是transpose op，这里都进行插入 transpose op的操作来获得更多的优化机会。
这里还涉及到几个相关的工具函数，我们也解释一下：
llvm::smallvector getchannellasttransposeperm() { return {0, 2, 3, 1}; }llvm::smallvector getchannelfirsttransposeperm() { return {0, 3, 1, 2}; }llvm::smallvector getinputoperandtransposeop(nchwcompatible op, value val, namedattrlist transpose_attributes, int num_transposed_operand, patternrewriter& rewriter) { std::string transpose_name = optrait::isopconfcompatible::getopname(op).str() + _transpose_input_ + std::to_string(num_transposed_operand); transpose_attributes.set(llvm::isopconfcompatible::getopnameattr()), rewriter.getstringattr(transpose_name)); smallvector input_operands; input_operands.push_back(val); auto res = rewriter .create(op.getloc(), getnhwctype(val.gettype()), input_operands, transpose_attributes) ->getresults(); return res;}transposeop getresulttransposeop(nchwcompatible op, value val, namedattrlist transpose_attributes, int num_transposed_result, patternrewriter& rewriter) { std::string transpose_name = optrait::isopconfcompatible::getopname(op).str() + _transpose_output_ + std::to_string(num_transposed_result); transpose_attributes.set(llvm::isopconfcompatible::getopnameattr()), rewriter.getstringattr(transpose_name)); smallvector operands; operands.push_back(val); transposeop transpose_op = rewriter.create( op.getloc(), getnchwtype(val.gettype()), operands, transpose_attributes); return transpose_op;}bool isinserttransposeopbefore(nchwcompatible op, patternrewriter& rewriter) { bool insert_transpose_op_flag = false; for (mlir::value operand : op->getoperands()) { transposeop transposeinputop = operand.getdefiningop(); if (!transposeinputop) continue; const auto perm = transposeinputop.getpermattr(); if (perm.size() == 4 && perm[0] == rewriter.getsi32integerattr(0) && perm[1] == rewriter.getsi32integerattr(3) && perm[2] == rewriter.getsi32integerattr(1) && perm[3] == rewriter.getsi32integerattr(2)) { insert_transpose_op_flag = true; break; } } return insert_transpose_op_flag;}
其中 getchannellasttransposeperm 和 getchannelfirsttransposeperm 方法分别返回 nhwc 到 nchw 和 nchw 到nhwc 的转换顺序。getinputoperandtransposeop 方法为 operation 的输入创建一个transpose operation。它使用输入值、transpose属性 和 重写器创建一个 transposeop , 并返回其结果。
类似的，getresulttransposeop 方法为 operation 的输出创建一个transpose operation。它使用输出值、transpose属性和重写器创建一个transposeop,并返回该operation。isinserttransposeopbefore方法检查operation的输入是否已有 transpose operation。如果有,并且该 transpose operation 将 nhwc 转为 nchw, 则返回 true, 否则返回false。
0x1.3 消除多余的transpose对
接下来，我们需要把插入transpose op的图中所有相邻的transpose对尽可能的消除，代码实现如下：
bool isredundanttransposematch(arrayattr pre, arrayattr afe, mlir::patternrewriter& rewriter) { const auto preperm = pre.getvalue().vec(); const auto afeperm = afe.getvalue().vec(); if (preperm.size() == 4 && afeperm.size() == 4) { // handle nchw->nhwc->nchw: (0, 2, 3, 1) -> (0, 3, 1, 2) if (preperm[0] == afeperm[0] && preperm[1] == afeperm[3] && preperm[2] == afeperm[1] && preperm[3] == afeperm[2] && preperm[0] == rewriter.getsi32integerattr(0) && preperm[1] == rewriter.getsi32integerattr(2) && preperm[2] == rewriter.getsi32integerattr(3) && preperm[3] == rewriter.getsi32integerattr(1)) return true; // handle nhwc->nchw->nhwc: (0, 3, 1, 2) -> (0, 2, 3, 1) if (preperm[0] == afeperm[0] && preperm[1] == afeperm[2] && preperm[2] == afeperm[3] && preperm[3] == afeperm[1] && preperm[0] == rewriter.getsi32integerattr(0) && preperm[1] == rewriter.getsi32integerattr(3) && preperm[2] == rewriter.getsi32integerattr(1) && preperm[3] == rewriter.getsi32integerattr(2)) return true; } return false;}struct autonhwceliminateredundanttransposepattern : public mlir::oprewritepattern { explicit autonhwceliminateredundanttransposepattern(mlir::mlircontext* context) : oprewritepattern(context, /*benefit=*/1) {} mlir::logicalresult matchandrewrite(transposeop op, mlir::patternrewriter& rewriter) const override { mlir::value transposeinput = op.getoperand(); transposeop transposeinputop = transposeinput.getdefiningop(); if (!transposeinputop || !isredundanttransposematch(op.getpermattr(), transposeinputop.getpermattr(), rewriter)) { return failure(); } rewriter.replaceop(op, {transposeinputop.getoperand()}); return success(); }};
isredundanttransposematch 方法检查两个 transpose operation的顺序是否会导致冗余。它通过比较两个 transpose 的 perm 属性来判断。类似 autonhwcpattern ，autonhwceliminateredundanttransposepattern 类继承自 oprewritepattern 。它对transposeop 进行重写以实现 transpose 消除。如果顺序是 nhwc->nchw->nhwc 或nchw->nhwc->nchw , 则判定为冗余 transpose 。
如果输入也来自transposeop且两个 transpose 顺序导致冗余,matchandrewrite方法会用transposeop的输入替换transposeop。实现 transpose 消除。matchandrewrite 方法首先获取 transposeop 的输入,并检查该输入是否也来自一个 transposeop。如果不是, 或两个 transpose 的顺序不导致冗余, 则返回 failure。最后返回 success 表示成功消除冗余 transpose 。
最终，上面介绍的2个pass都被封装到 autonhwcpass 中作用在 mlir 的计算图上完成全局优化。从下面的代码可以看到这个优化只有在打开 oneflow_mlir_prefer_nhwc 环境变量时才正常生效。
void populateautonhwcpatterns(::rewritepatternset& patterns) { bool enable_nhwc = ::parsebooleanfromenv(oneflow_mlir_prefer_nhwc, false); if (enable_nhwc) { patterns.add(patterns.getcontext()); patterns.add(patterns.getcontext()); }}class autonhwcpass : public autonhwcpassbase { void runonoperation() override { operation* op = getoperation(); rewritepatternset patterns(op->getcontext()); oneflow::populateautonhwcpatterns(patterns); (void)applypatternsandfoldgreedily(op, std::move(patterns)); }};
补充：0x1.4 weight的transpose消除
这里还需要粗略的说明一下对于 weight 的 transpose 是如何处理的。在0x1.2中我们为 weight（常量constant op） 也插入了 transpose op，然后我们知道 weight 是常量，所以针对 weight 的 transpose op 完全可以在编译期折叠起来。这个过程是在 https://github.com/oneflow-inc/oneflow/blob/master/oneflow/ir/oneflow-translate/lib/oneflow/mlironeflowtranslation.cpp#l808-l811 这里完成的，我们后面会单独介绍一下 constant folding 的实现。
0x2. 结论
本文介绍了一下oneflow的编译器中的 layerout transform，这个技术在后来 oneflow 版本的 stable diffusion 中也发挥了重要作用，提升了推理速度。在 tvm 的 ansor 中也有类似的优化，通过将不同的 layerout 设定为 op 的 strategy 进而影响 op 的 schedule，在搜索的时候考虑到 layerout transform 来获得更大的搜索空间和更好的结果。在处理transpose 额外开销的方法并不是唯一的，这里只是采用了一种个人认为比较简单的方式，读者们如果有类似需要可以自由发挥。
基站蓄电池组远程在线监测解决方案
西部数据SN750带散热片版本的优缺点
迎接3D IC时代 半导体供应链加速投入研发
当物联网遇上人工智能时，成为创新的超能力
安卓8.0最新消息：安卓8.0初体验，竟然这么流畅还有这么多黑科技功能
深度学习编译器之Layerout Transform优化
一种基于人工智能的脑图像处理方法
大众CEO：加速业务改革，不要成为诺基亚
魅族pro7什么时候上市?魅族pro7最新消息:魅族pro7终极曝光:价格、配置抢先看,外观设计一流你不期待?
磁阻电机的工作原理 磁阻电机的优缺点
电流互感器和电压互感器的极性及测定方法
罗姆推出业界最小晶体管封装“VML0806”
柠檬豆：把握采购方的真需求就是我们的追求
多速率转换信号频谱研究
无线电发射功率与发射距离有什么关系
西陇科学与宝洁公司在西陇科学总部举行了合作20周年纪念会议
贴片电容的基本原理和功能 贴片电容在直流电路中的应用
中国科技巨头实现弯道超车，发布全球最强5G芯片
通过引导分区破解树莓派的6种方法
星恒电源对锂电行业发展的思考与经验