diff --git a/BMM算子优化分析_Release/BMM算子优化分析_v0.4.html b/BMM算子优化分析_Release/BMM算子优化分析_v0.4.html new file mode 100644 index 0000000..73d91aa --- /dev/null +++ b/BMM算子优化分析_Release/BMM算子优化分析_v0.4.html @@ -0,0 +1,659 @@ + + + + + +BMM 算子优化分析 v0.4 — 昇腾 950PR + + + + + +
+

BMM 算子优化分析

+
版本:v0.4 | 目标芯片:昇腾 950PR(DAV_3510,32 AIC / 64 AIV,GM 1.6TB/s,L2 128MB/5.2TB/s)
基于 issue#3 (v0.3) 完善:补全 StreamK/ASW_Basic 实现方案,精化全部分支的进入条件与实现步骤,加入源码对照与硬件依据
+
+

一、算子功能与接口说明

+

1.1 算子功能

+

Batch Matmul(BMM)完成带 batch 维的矩阵乘:对 batch 维的每个索引独立执行一次矩阵乘,再加可选偏置:

+
$$ +C[b, m, n] = \sum_{k=0}^{K-1} A[b, m, k] \cdot B[b, k, n] + bias[b, 1, n] +$$
+

C = A @ B + bias。A、B 输入维度典型为 3 维(支持最多 4 级 batch 维 $b = b_0 b_1 b_2 b_3$ 展平),最后两维做矩阵乘。

+

1.2 接口参数

+ + + + +
参数形状数据类型layout说明
左矩阵 A[BatchA, M, K]dtype(FP16/BF16/FP8/…)典型 ND可带转置标记(isATrans)
右矩阵 B[BatchB, K, N]dtype典型 ND可带转置标记(isBTrans),推理场景可为 weightNz
偏置 bias[B, 1, N]dtype固定 ND可为空
输出 C[BatchC, M, N]dtype(可随路量化)典型 NDBatchC = broadcast(BatchA, BatchB)
+

广播语义:batch 维兼容广播——每一级 batch 维上,两侧取值要么相等、要么为 1;输出该维取两者的最大值。两类广播的形态差异:

+ +

1.3 case 空间

+

一个 BMM case 的实现方式完全由以下特征量决定:

+
$$ +(\; B,\; M,\; N,\; K,\; \text{dtype},\; \text{layout/转置},\; \text{广播形态} \;) +$$
+

其中 $B$ = BatchC(展平后的输出 batch 数)。后续全部分析就是回答:给定这组特征量,最优实现是什么。

+
+

二、性能模型:什么叫"最优"

+

2.1 总时延 = 流水线最慢的一级

+

NPU 上 BMM 的执行是核内多级硬件流水的并行——Cube 计算(MMAD)、GM/L2→L1 搬移(MTE2)、L1→L0 搬移(MTE1)、L0C 写出(Fixpipe)。多核并行时各流水级时延可被双缓冲(double buffer)等机制相互掩盖,最终:

+
$$ +\boxed{\;T_{total} = \max\big(T_{MMAD},\; T_{MTE2},\; T_{MTE1},\; T_{Fixpipe}\;[,\;T_{Reduce}]\big)\;} +$$
+

总时延 = 流水线最慢的一级。 算子优化的关键就是对瓶颈流水级的优化。由此得出一个重要的设计自由度——瓶颈交换

+
当 MTE2(搬移)是瓶颈、MMAD(计算)不是瓶颈时,可以牺牲一定 MMAD 时延(例如冗余计算)换取 MTE2 性能提升;反之,当 MMAD 是瓶颈时,可以牺牲一定 MTE2 时延(例如重复搬移)换取计算效率提升。只要瓶颈级时延下降,总时延就下降。
+

后文会看到:MergeBatch 就是"牺牲算力换搬移效率"的典型,ASW_Basic 切 M/N 就是"牺牲搬移(重复读)换并行度"的典型——它们的存在正当性都来自这个 max 模型。

+

2.2 目标芯片关键规格与经验常数

+ + + + + + +
规格项数值对设计的意义
AIC / AIV 核数32 / 64(1:2)核间并行度上限 $C = 32$;StreamK 归约可用 AIV
Cube 算力 FP16/BF16≈486 TFLOPS(单核 ≈ 16×16×16 MAC/拍 × 1.65GHz)算存比分子
GM 带宽1.6 TB/s(读写共享总线)访存 Bound 的分母;读+写互相挤占
L2 Cache128MB,5.2 TB/s(读写各独享)重复读取的吸收层;Cube 输出写 L2 即算完成
L1 / L0A / L0B / L0C512KB / 64KB / 64KB / 256KB(每 AIC)核内分块的容量约束;L0C 按 FP32 累加计
UB512KB(每 AIV)向量通路/归约
+

GM→L1 搬移效率(ND2NZ)达到 90%+ 带宽利用率的四条经验约束(重要性递减):

+
    +
  1. 参与核数:建议 32 核并行搬移,至少 3/4 核(24 核)并发——核数不足时带宽利用率上限被压低;
  2. +
  3. 单核搬移数据总量 ≥ 480KB(min_DatamountPerCore);
  4. +
  5. 单次搬移 tile 大小 ≥ 16KB(min_TileSize);
  6. +
  7. dValue(单次搬移的连续内轴字节数,如非转置 A 的 K 向、非转置 B 的 N 向)≥ 128B,建议 256B,最好 512B。
  8. +
+
这些常数(480KB、16KB、256B/512B,以及后文的 $b_{thr}=4$、$k_{thr}$=32B/dtype、$minCoreNum \approx 0.8 \times 32$)都是该档芯片的实测经验值;换芯片时所有逻辑结构不变,只需替换常数表。
+

2.3 算存比与 Bound 判定

+

约定:以"元素"为访存单位(1 元素 = dtype 字节),算存比单位为 FLOP/元素。

+ +
$$ +AI = \frac{2MNK}{(MK + KN)} = \frac{2MN}{M + N} \quad [\text{FLOP/元素}] +$$
+ +
$$ +AI_{full} = \frac{2MNK}{MK + KN + MN} \quad [\text{FLOP/元素}] +$$
+ +
$$ +R = \frac{\text{Cube 峰值算力}}{\text{GM 带宽} / \text{dtype 字节}} = \frac{486 \times 2}{1.6} \approx 607.5 \quad [\text{FLOP/元素,BF16/FP16}] +$$
+

判定规则:$AI < R$ → 访存 Bound(瓶颈在 MTE2,优化重心是搬移效率与复用);$AI > R$ → 计算 Bound(瓶颈在 MMAD,优化重心是 Cube 利用率与流水掩盖)。

+
大白话:矩阵乘就像搬砖盖楼——Cube 是起重机(算力),GM 带宽是供货卡车(带宽)。算存比就是"每块砖要盖多少平米"。每块砖盖的平米数低于起重机的额定配比时,楼永远盖不快,问题出在供货而不在起重机——这时候优化要围绕"怎么少堵车、怎么一次多运点",而不是换更大的起重机。
+
+

三、分支划分的系统推导

+

本章是全文逻辑主线:从分块计算的本质出发,不靠经验罗列,推导出六大分支的划分。

+

3.1 分块计算的本质与四个可切维度

+

BMM 在 NPU 上实现的本质是:把参与计算的数据分块(tile),由 32 个 AIC 核并行 + 串行地完成这些分块的计算,再组合成最终结果:

+
$$ +C[B, M, N] = \sum_{\kappa} A[B, M, K_\kappa] \cdot B[B, K_\kappa, N] +$$
+

分块有 4 个维度:B、M、N、K核间怎么分这 4 个维度,就是分支划分的第一性问题(核内分块是第二性问题,属于各分支内部的 tiling)。

+

3.2 四个维度的切分特征:一张决定性的表

+

从"读入 / 计算 / 写出"三个视角考察每个维度的核间切分特征(这是后续一切推导的基石):

+ + + +
切分维度读入特征计算特征写出特征
切 B每个数据块只被固定的 1 个核读取,核间零重复读每个输出块由 1 个核独立完成,无核间依赖只写最终结果,无中间结果
切 M / 切 N切 M 则同一右矩阵块被多核重复读;切 N 则同一左矩阵块被多核重复读每个输出块由 1 个核独立完成,无核间依赖只写最终结果,无中间结果
切 K每个数据块只被固定的 1 个核读取,零重复读每个输出块由多核共同完成,存在核间依赖有中间结果写出,需核间 Reduce 归约
+

为什么特征差异这么大?两条硬件层面的根本原因:

+ +

3.3 完备枚举:15 种核间切分组合

+

4 个维度的任意非空子集都可作为一种核间切分方案,共 $2^4 - 1 = 15$ 种:

+
$$ +\{B\},\{M\},\{N\},\{K\},\{B,M\},\{B,K\},\{B,N\},\{M,K\},\{M,N\},\{K,N\},\{B,M,K\},\{B,M,N\},\{B,K,N\},\{M,K,N\},\{B,M,K,N\} +$$
+

任何实现方案必属于其中一种 ⇒ 这 15 种是完备的。按 §3.2 的特征分组:

+ + + +
组合共同特征优化重心
纯 B{B}零重复读,读写数据量固定访存 Bound:搬移效率高;计算 Bound:计算效率高
含 M/N 不含 K{M},{N},{M,N},{B,M},{B,N},{B,M,N}可能有重复读访存 Bound:重复读尽量少 + 搬移效率高;计算 Bound:计算效率高
含 K{K},{B,K},{M,K},{K,N},{B,M,K},{B,K,N},{M,K,N},{B,M,K,N}有中间结果写出 + 归约计算效率高,且归约引入的额外 MTE2/Fixpipe 时延不能成为新瓶颈
+

3.4 代价不对称性:切分维度的"价格表"

+

15 种组合的代价结构只由两个布尔特征决定——是否含 K、是否含 M/N。三个维度的"核间切分价格"严格排序:

+
$$ +\text{cost}(\text{切}B) = 0 \;<\; \text{cost}(\text{切}M/N) \;\ll\; \text{cost}(\text{切}K) +$$
+ +

"价格表"不是经验,是 BMM 语义 + L0C 累加机制 + L2/GM 带宽结构三条事实的推论。 整条分支决策树就是一句话:按价格从低到高购买并行度,买不够才加价。

+

3.5 推导主链:从价格表到六大分支

+

第 0 层:问题归约与通路选择。

+ +

第 1 层:归约不掉的 case(BatchA = BatchB = B > 1,K ≥ 2),必须在 BMM 框架内组织 4 维分块。先买免费的 B 维。

+ +

第 2 层:免费的 B 买不满 32 核,加价买廉价的 M/N。

+ +

第 3 层:廉价维度也买不满,才买昂贵的 K。

+ +
case (B, M, N, K, dtype, layout, 广播形态)
+   │
+   ▼
+[0] K=0? → AIV 清零;K=1? → AIV 逐元素乘          ← 前置通路层(正交于切分)
+   │ K ≥ 2
+   ▼
+[1] BatchA=1 或 BatchB=1? ──是──→ 转Matmul(问题归约)
+   │ 否(BatchA=BatchB=B>1,或交叉广播)
+   ▼
+[2] 并行度账本:P = B·⌈M/16⌉·⌈N/16⌉ 与 C=32 比较
+   │
+   ├─ B ≥ 32(切 B 可满核)─────────────────────────┐
+   │    M×N 大,Cube 饱 → IterBatch                 │
+   │    M×N 小,Cube 饿,且访存 Bound → MergeBatch  │
+   │                                                │
+   ├─ B < 32 且 P ≥ 32 → ASW_Basic(切 M/N 补并行,  │
+   │                     共享读取由 L2+swizzle 吸收) │
+   │                                                │
+   └─ P < 32(B/M/N 用尽仍缺并行,K 大)→ StreamK    │
+      (切 K 买并行,付归约代价)◀──────────────────┘
+

由此推导出六大 BMM 本体分支(转Matmul、MergeBatch、IterBatch、ASW_Basic、StreamK)加一个前置通路层(K 退化),共七个决策路径。

+

3.6 完备性与极小性:不多不少

+

完备性。对归约不掉的 case,15 种组合按"(是否含 K,是否含 M/N)"坍缩:

+ +

$1 + 1 + 2 = 4$,加上前置的归约分支(转Matmul)与通路层(K 退化)= 7。任何合法 case 必落其一,无空洞。

+

极小性。每个分支都有它"唯一最优"的 shape 区域,去掉任何一个都会有 case 失去最优实现:

+ + + + + + +
分支独占最优的代表 case(BF16)替代方案为何更差
特殊(K=0/1)K=0 或 K=1Cube 通路完全或几乎无用
转MatmulBatchB=1, B=128, M=N=K=2048免费折叠后 Matmul 体系的 L2 切分/全载直接可用;BMM 分支内重做无收益
IterBatchB=32, M=N=K=4096ASW 切 M/N 引入无谓共享读;MergeBatch 引入无谓冗余算力
MergeBatchB=128, M=32, N=128, K=64IterBatch 的 L0C 利用率仅 ~1.6%,Cube 空转,搬移 tile 碎(<16KB)带宽利用率崩
ASW_BasicB=2, M=N=8192, K=1024切 B 仅 2 核干活;StreamK 付无谓归约
StreamKB=1, M=N=64, K=65536不切 K 时 P=16 < 32,近半核闲置,时延差数量级
+

⇒ 七大路径构成极小完备集

+

3.7 大白话总结

+
把 BMM 想成给 32 个工人分一批"矩阵乘订单":
- K=0/1 的订单太简单,不需要起重机:Cube 派不上用场,直接让搬运工(AIV)干完(特殊分支);
- 能整单外包的就别自己干:一侧只有一个 batch,问题其实就是一次普通矩阵乘,直接交给成熟的 Matmul 产线(转Matmul);
- 按订单分(切 B)最省心:每人几单,互不干扰(零重复搬料、零协调)。订单多时,单大就一单接一单干(IterBatch);单太小时一人同时干几单、拼成大活干,虽然多做了点无用功,反正瓶颈在等料不在干活(MergeBatch);
- 订单不够分,就按行/列拆单(切 M/N):大家会重复领同一份料,但料放在近处仓库(L2),多跑几趟近仓库很便宜(ASW_Basic);
- 行/列也不够拆,只好沿深度 K 拆:几个人合做同一块输出,各算一段再汇总——汇总是要额外开会的(Reduce),所以只在实在分不满、且 K 足够长时才这么干(StreamK)。
+
+

四、各分支详解:进入条件与实现方案

+
每节按统一结构展开:做什么(定义)→ 什么时候进(进入条件逐条 + 每条的理由)→ 怎么做(实现方案与参数求解)→ 与源码对照(批判性)
+

4.0 前置通路层:K 退化 case

+

在切分决策之前先做计算通路判断:

+ + +
K 值数学本质计算通路数据通路核类型
K = 0无任何计算,C = bias 或 0无 mmadGM → AIV → GMAIV_ONLY
K = 1退化为逐元素乘 $C = A \odot B$,无累加深度无 mmad,向量 MulGM → UB → Mul → GMAIV_ONLY
+

为什么是 AIV 而非 AIC:Cube 阵列一拍完成 16×16×16(fp16)的 fractal 乘加。K=1 时 K 维只有 1/16 被利用——阵列 15/16 的 MAC 空转,还要付出 GM→L1→L0A/L0B 的分形搬运、L0C 累加与 fixpipe 写出全链路开销。而 AIV 每拍处理 256B 连续数据,逐元素乘 + UB 多 batch pingpong 能把 MTE 带宽吃满;64 个 AIV 的聚合向量算力(FP16 54 TFLOPS)对这个计算密度绰绰有余。

+

源码对照:K_EQUAL_ZERO(策略 0)与 TO_MUL(策略 1)两个独立策略处理,tiling 直接下发 AIV_ONLY kernel。触发条件 TO_MUL 要求 batchC ≥ 128(= aivNum×2,开 UB pingpong)、N 不在 (32B/dtype, 256B/dtype] 区间(中间区让给其它分支)、UB 容量放得下。

+

4.1 转Matmul(问题归约分支)

+

做什么:当 BatchA = 1BatchB = 1 时,把 batch 维折叠进矩阵维度,转化为普通 Matmul,计算完成后按需恢复 batch 维。

+

为什么这么做:单边 batch=1 的 BMM 与普通 Matmul 在数学上只差一个维度标签。Matmul 的优化体系(L2 切分、AL1/BL1 全载、swizzle)比 BMM 各分支在"单样本"情形下更成熟——站在巨人肩膀上,不重复造轮子

+

折叠规则与代价(关键的不对称性)

+ + +
情形折叠方式是否免费代价分析
BatchB = 1左矩阵 [B, M, K] 的 batch 维与 M 维在 ND 布局下内存相邻,直接视图为 [B·M, K];输出 [B·M, N][B, M, N] 的内存布局逐元素一致完全免费零输入重排、零输出 split,直接转 Matmul 无任何代价
BatchA = 1需将右矩阵 [B, K, N] 折叠为 [K, B·N]有代价B 的 batch 维与 N 维在内存中不相邻(中间隔 K),折叠等价于一次 [B,K,N]→[K,B,N] 的转置重排(O(B·K·N) 读写),且输出 [M, B·N] 与目标 [B, M, N] 之间存在置换,需要随路 scatter
+

BatchA = 1 时的决策规则

+ +

进入条件BatchA = 1 || BatchB = 1。其中 BatchB = 1 恒进(免费);BatchA = 1 时按上述代价比较决定。

+

源码对照:源码中的 MergeBatchAndMAxis() 只做了 batchB=1 方向的折叠(args_.mValue = batchA * mValue),与"该方向免费"的判断互为印证。batchA=1 方向源码未做折叠,而是走 AL1_FULL_LOAD 或 ASW_Basic 的广播处理。

+

4.2 MergeBatch(多 batch 合并计算)

+

4.2.1 做什么

+

核间按 B 分核(每核负责 $b_{core}$ 个 batch),核间无同步无通信。核内将多个 batch 合并计算

+
$$ +[b, M, K] @ [b, K, N] \;\Rightarrow\; [bM, K] @ [K, bN] = [bM, bN] \;\xrightarrow{\text{BlockTrace}}\; [b, M, N] +$$
+

其中 BlockTrace 指以 $[M, N]$ 的 block 粒度取结果矩阵的块对角线作为各 batch 的有效输出:$C[i, m, n] = R[iM + m,\; iN + n]$。交叉项(第 i 个 batch 的 A 乘第 j≠i 个 batch 的 B)被算出但丢弃——这就是"算力浪费",浪费比例 $(b-1)/b$。

+

为什么允许浪费:进入该分支的 case 必然是访存 Bound(条件三保证),瓶颈在 MTE2 不在 MMAD,浪费的算力被搬移时延掩盖(§2.1 瓶颈交换)——用本来闲置的 Cube 算力,换 tile 变大后的搬移效率与 Cube 利用率

+

4.2.2 进入条件(逐条 + 理由)

+

设计原则四条:① 硬件时延可流水掩盖(double buffer 乒乓);② 满足容量约束(每次计算 L1/L0A/L0B/L0C 放得下);③ GM→L1 搬移高效(§2.2 四条经验约束);④ 算力有浪费但计算不能成为瓶颈。

+

形式化后,进入 MergeBatch 需同时满足:

+

条件 1(batch 够分且够合并)BatchA = BatchB(无广播)且

+
$$ +b_{core} = \frac{B}{C} \ge b_{thr} +$$
+

$b_{thr}$ 是"多 batch 合并搬移能拿到效率收益"的最小合并数——$b$ 太小时合并的搬移收益抵不过实现复杂度,不如 IterBatch 且完全不浪费算力(功耗)。经验值 $b_{thr} = 4$(源码 MIN_BATCH_L0 = 4,注释 "each aic should process at least 4 batchs")。

+

条件 2(合并数 b 的上下界):设核内单次 Cube 计算合并 $b$ 个 batch($b \le b_{core}$),要求 $b \ge b_{thr}$ 且 $b$ 由以下约束共同封顶:

+ +
$$ +AI(b) < R \;\Longleftrightarrow\; b < b_{AI} = \frac{R \cdot (M + N)}{2MN} +$$
+ +
$$ +2 \cdot (bM)(bN) \cdot 4\text{B} \le 256\text{KB} \;\Rightarrow\; b \le b_{L0C} = \sqrt{\frac{256\text{KB}}{2 \cdot MN \cdot 4\text{B}}} +$$
+ +
乒乓取舍:若 $b_{core}$ 小于"不乒乓时 L0C 允许的合并数上限",说明 batch 余量不足,只能不乒乓(牺牲流水掩盖换合并数);否则乒乓。形式上 $b_{max} = \min(b_{AI},\; b_{L0C}^{(\text{是否乒乓})},\; b_{core})$。
+

条件 3(访存 Bound 性质)

+
$$ +\frac{2MN}{M + N} < \frac{R}{b} +$$
+

即 case 固有算存比显著低于平衡点(低一个合并倍数 $b$ 的量级)——这是"浪费可被掩盖"的定量保证。

+

条件 4(K 向搬移效率):L1 级 K 切分 $k_{L1} = K / StepK$(StepK 为正整数),要求 $k_{L1} \ge k_{thr}$,$k_{thr} = 32\text{B}/\text{dtype}$(BF16 即 16 元素)。理由:$k_{L1}$ 决定 ND2NZ 的 dValue,过碎则搬移指令效率崩(§2.2 第 4 条)。

+

条件 5(L1 驻留与 tile 效率):L1 级 batch 驻留数

+
$$ +b_{L1} = \frac{L1_{size}}{(M k_{L1} + k_{L1} N) \cdot \text{dtype}}, \qquad b_{L1} > b +$$
+

(L1 驻留的 batch 组必须大于单次计算的合并数,否则合并无从谈起),且搬移 tile 大小满足效率:左矩阵非转置时 $\max(b_{L1} M k_{L1} \cdot \text{dtype},\; k_{L1} N \cdot \text{dtype}) > min\_TileSize$(16KB);左矩阵转置时相应调整。

+

条件 6(单核搬移总量)

+
$$ +b_{core} \cdot (M k_{L1} + k_{L1} N) \cdot \text{dtype} \ge min\_DatamountPerCore = 480\text{KB} +$$
+

对应 §2.2 第 2 条:单核搬移数据总量不足时带宽利用率上限被压低。

+

条件汇总逻辑:条件 1 定资格(batch 够多),条件 2/3 定上限(别算出瓶颈、别撑爆 L0),条件 4/5/6 定下限(搬移效率不能崩)。上限与下限之间必须有交集,交集为空则该 case 与 MergeBatch 无缘。

+

4.2.3 核内参数求解

+

进入分支后,$b$ 与 baseK 的具体取值有一个重要性质:**在访存 Bound 前提下,$b < b_{AI}$ 时 b 的取值不影响性能,baseK 只需满足分形约束(16 倍数)也不影响算存比**——因为瓶颈是搬移,算力余量内怎么切都一样。据此:

+
    +
  1. 先由条件 2 求 $b_{max}$,实际 $b$ 在 $[b_{thr}, b_{max}]$ 内取,且尽量均匀($b_{core}$ 整除 $b$)——每次计算的 $b$ 均匀一致对功耗更有利;
  2. +
  3. 再由 $b$ 反查 L0AB 允许的最大 baseK:
  4. +
+
$$ +baseK_{max} = \min\left(\frac{64\text{KB}}{2 \cdot bM \cdot \text{dtype}},\; \frac{64\text{KB}}{2 \cdot bN \cdot \text{dtype}},\; \frac{C0_{size}}{\text{dtype}}\right) +$$
+

向下取 16 倍数。

+
    +
  1. L1 级参数:$k_{L1} \ge \min(k_{L0\_max},\; 128\text{B}/\text{dtype})$,$b_{L1} = \min(b_{L1\_max},\; b_{core})$。
  2. +
+

4.2.4 数值例子

+

case:B=128, M=32, K=64, N=128, BF16,32 核。

+
    +
  1. 资格:$b_{core} = 128/32 = 4 \ge b_{thr} = 4$ ✓;
  2. +
  3. 算存比:$AI = 2 \times 32 \times 128 / (32 + 128) = 51.2$ FLOP/元素 $< R = 607.5$ ✓ 访存 Bound;$b_{AI} = 607.5 / 51.2 \approx 11.86$,即 $b \le 11$ 都不会进入算力 Bound——算存比约束很宽;
  4. +
  5. L0C 约束(乒乓):$2(b \cdot 32)(b \cdot 128) \times 4\text{B} \le 256\text{KB} \Rightarrow b^2 \le 8 \Rightarrow b \le 2$(不乒乓则 $b \le 4$)——L0C 才是真正的紧约束
  6. +
  7. 取 b = 2(= $b_{core}$ 的因子,均匀):L0B 给出 $k_{L0} \le 64\text{KB} / (2 \times 2 \times 128 \times 2\text{B}) = 64$,$K = 64$ 一步到位,baseK = 64;
  8. +
  9. 搬移效率:$k_{L1} = 64 \ge k_{thr} = 16$ ✓;单核搬移量 $4 \times (32 \times 64 + 64 \times 128) \times 2\text{B} = 80\text{KB}$ —— 不足 480KB,条件 6 不满足!
  10. +
+

结论:该 case 在"单核搬移总量"一项上不达标——这正是 MergeBatch 条件体系的用处:它提前告诉你瓶颈不在算力而在搬移效率

+

4.2.5 执行流程

+
核间:32 核,每核 b_core 个 batch
+核内(每核):
+  for k_l1 in range(0, K, kL1):                    # L1 级 K 切分
+    MTE2: A[bL1组, M, k_l1:k_l1+kL1]、B[bL1组, k_l1:k_l1+kL1, N] → L1(乒乓)
+    for b_start in range(0, bL1, b):               # 按合并数 b 分组
+      for k_l0 in range(0, kL1, kL0):              # L0 级 K 切分
+        MTE1: A[bM, kL0] → L0A;B[kL0, bN] → L0B
+        Cube: [bM, kL0] @ [kL0, bN] → L0C 累加
+      Fixpipe: BlockTrace 取 b 个 [M, N] 对角块 → L2/GM
+

4.2.6 与源码对照

+

源码 batch_matmul_v3_mergebatch_basicapi_tiling.cpp 的 IsCapable 条件:各级 batchA_i == batchB_i、batchC ≥ 4 × aicNumalignK ≥ 64M ≤ N、无 bias、非 NZ、拒绝非连续转置等。对照分析:

+ + + + +
源码条件本文对应差异
batchC ≥ 4 × aicNum条件 1:$b_{core} \ge b_{thr} = 4$一致(每核至少 4 batch)——合理
alignK ≥ 64条件 4:$k_{L1} \ge k_{thr}$=32B/dtype(BF16 为 16 元素)源码保守 4 倍,可能误杀 K∈[16,64) 的可获益 case
M ≤ N无此限制(M > N 的镜像 case 原理上同样可合并,交换合并方向即可)源码放弃镜像 case,属于覆盖缺口
无显式算存比/搬移效率判定条件 3/4/6源码把"是否最优"推给优先级顺序,本文条件体系是补全
+

4.3 IterBatch(逐 batch 计算)

+

4.3.1 做什么

+

核间按 B 分核(每核 1 个或多个 batch),核间无同步无通信;核内逐个 batch 分别执行标准 Matmul 分块(L1→L0A/L0B→Cube→L0C→Fixpipe)并输出。无算力浪费、无跨 batch 依赖——是"切 B"最朴素的形态。

+

4.3.2 进入条件(三大类 + 理由)

+

满足以下任一

+

类 1(计算 Bound 型)

+
$$ +AI_{full} = \frac{2MKN}{MK + KN + MN} \ge R \quad \text{且} \quad BatchA = BatchB,\; b_{core} = \frac{B}{C} \in \mathbb{Z}^+ +$$
+

理由:连输出写出都计入仍是计算 Bound,则瓶颈恒在 Cube——只要负载均衡($b_{core}$ 整除,各核同量),逐 batch 计算就是满算力实现,无需任何花哨。

+

类 2(输出驻留 L2 的均衡型)

+
$$ +\frac{2MN}{M+N} \ge R_{读GM} \;\;\text{且}\;\; 2K \ge R_{写L2} \;\;\text{且}\;\; MN \cdot \text{dtype} \le L2_{size} \;\;\text{且}\;\; BatchA = BatchB,\; b_{core} \in \mathbb{Z}^+ +$$
+

理由:读侧访存不弱($AI \ge R_{读GM}$);单 batch 输出 $MN \cdot \text{dtype} \le 128$MB 可驻留 L2——Cube 输出写 L2 即算完成,写 GM 的流量被省掉;$2K$ 是"每写出一个元素对应的计算量",$2K \ge R_{写L2}$ 保证写 L2 的 5.2TB/s 也不是瓶颈。三个条件合起来 = 读、算、写三条路都不堵

+

类 3(访存 Bound 型):$\dfrac{2MN}{M+N} < R$ 且同时满足:

+
    +
  1. $BatchA = BatchB$ 且 $b_{core} = B / C \ge 1$;
  2. +
  3. 负载均衡:访存 Bound 时核负载利用率建议 $b_{Avg}/b_{Max} > 0.8$(计算 Bound 时建议 100%)。操作化表述:B 整除核数,或尾波(B mod C)活跃的核数 $\ge minCoreNum$(取 $0.8 \times 32 \approx 26$)。注意:若只写 B mod C > minCoreNum 会误杀整除 case(余数 0 恰恰是完全均衡),判据应以均衡比为准、余数规则为其近似
  4. +
  5. 单核搬移不重复读(访存 Bound 下重复读就是纯损失),按 L1 容量分五种形态之一:
  6. + +
+

4.3.3 设计原理

+
    +
  1. 核间切 B 零共享零依赖,唯一的系统性风险是负载不均——所以均衡是第一条件;
  2. +
  3. 输出须满足 L0C 容量:$MN \cdot 4\text{B} \le L0C$(否则核内还要切 M/N,那就不是纯 IterBatch 而是 ASW 行为——但 IterBatch 核内允许对单 batch 做标准 M/N/K tiling,此处的准确含义是"核内 tiling 不构成跨 batch 的耦合");
  4. +
  5. 访存 Bound 时单核数据不得重复读——L1 装得下才不重复,装不下就按 Step 切分且切分后仍满足搬移效率下限;
  6. +
  7. 与 MergeBatch 的分工:IterBatch 不浪费算力,但需要"单 batch 足够大"撑搬移效率与 Cube 利用率;MergeBatch 用浪费换效率,专治小 M×N。两者在 $M \times N$ 的中段重叠,由时延模型仲裁(第五章)。
  8. +
+

4.3.4 核内参数求解

+

**情形 (a):$b_{core} = 1$ 且 L1 可放完整单 batch**

+
if (L0C_Size >= M * N * L0C_ElementSize):
+    BaseM = M;  BaseN = N
+    BaseK = min(align(L0A / BaseM, 16), align(L0B / BaseN, 16))
+else:
+    if M < N:
+        BaseM = align(M, 16)
+        BaseN = floor(L0C / BaseM)
+        BaseK = min(floor_align(L0A / BaseM, 16), floor_align(L0B / BaseN, 16))
+    else:
+        BaseN = align(N, 16)
+        BaseM = floor(L0C / BaseN)
+        BaseK = min(floor_align(L0A / BaseM, 16), floor_align(L0B / BaseN, 16))
+

**情形 (b):$b_{core} > 1$ 且 L1 可放 2 batch 乒乓**

+
if (L0C_Size >= 2 * M * N * L0C_ElementSize):
+    BaseM = M;  BaseN = N
+    BaseK = min(floor_align(L0A / BaseM, 16), floor_align(L0B / BaseN, 16))
+else:
+    if M < N:
+        BaseM = align(M, 16)
+        BaseK = min(floor_align(L0A / 2 / BaseM, 16), K)
+        BaseN = max(floor_align(L0C / 2 / BaseM, 16), floor_align(L0B / 2 / BaseK, 16))
+    else:
+        BaseN = align(N, 16)
+        BaseK = min(floor_align(L0B / 2 / BaseN, 16), K)
+        BaseM = max(floor_align(L0C / 2 / BaseN, 16), floor_align(L0A / 2 / BaseK, 16))
+

情形 (c)/(d)/(e):L1 放不下完整 batch

+

一矩阵不切、另一切分:假设 L1 放完整左矩阵和部分右矩阵,则右矩阵应搬入的 K 向长度:

+
$$ +k_{L1\_b} = \min\!\left(\frac{L1_{size} - M \cdot K \cdot \text{dtype}}{N \cdot \text{dtype}},\; K\right) +$$
+

且 $k_{L1\_b} \ge k_{thr}$,切分后 tile ≥ 16KB。

+

4.3.5 执行流程

+
核间:32 核分 batch(尽量整除,尾波核数 ≥ minCoreNum)
+核内(每核):
+  for batch in 本核的 b_core 个 batch:           # 逐个 batch
+    for m_tile / n_tile(核内标准 tiling):
+      for k_tile in range(0, K, baseK):
+        MTE2 预取下一 k_tile → L1(双缓冲)
+        MTE1: L1 → L0A/L0B
+        Cube: mmad → L0C 原地累加               # K 循环不出核
+      Fixpipe: L0C → L2(写 L2 即完成,GM 回写可异步)
+

4.3.6 与源码对照

+

源码有三个 iterbatch 变体:ITER_BATCH_BROADCAST_BASICAPI(单边广播)、ITER_BATCH_BASICAPI(基础 API)、ITER_BATCH(高阶 API)。关键差异:

+ + + + +
维度源码 basicapi源码高阶 API本文
L1 容量计算$(sizeA + sizeB + bias) \times 2 \le l1Size$(除 DB)$iterBatch = l1Size / inputSizeOneBatch$(不除 DB)按 L1 形态 (a)~(e) 分类
截断mmadCount=8(issue queue)、fullCopySize=64KBiterBatch ≤ 4 时 singleCoreK 减半统一为搬移效率 + 均衡率约束
均衡率0.80.80.8(访存 Bound)/ 1.0(计算 Bound)
广播独立分支(单边单轴)不支持由转Matmul 或 ASW_Basic 吸收
+

4.4 StreamK(K 维核间切分)

+

4.4.1 做什么

+

当 B、M、N 三个维度切到最碎仍填不满 32 核时,把 K 维切到核间:多核各算一段 K 的部分和,再归约:

+
$$ +C_{(\beta,\mu,\nu)} = \sum_{c \in \text{group}} C^{(c)}_{(\beta,\mu,\nu)} \quad \text{(部分和写 workspace,AIV 归约或原子加)} +$$
+

4.4.2 进入条件(两条缺一不可)

+

条件 1(并行缺口存在):不切 K 时的独立输出块数

+
$$ +P = B \times \lceil M / 16 \rceil \times \lceil N / 16 \rceil < C = 32 +$$
+

(源码取更保守的 $B \cdot mCnt \cdot nCnt \le aicNum/2$,并附加 ND-only、无交叉广播、确定性等级 ≤ 1 等工程限制)。注意严格的 StreamK 不要求核间完全不切 B/M/N——它是一般框架 grid_K × grid_B × grid_M × grid_N ≤ C,纯切 K 只是 grid_B=grid_M=grid_N=1 的特例;当 B/M/N 能提供部分并行度时,应该用组合 grid 把归约组 $grid_K$ 压到最小(grid_K 小一档,K 的门槛降一档平方级)。

+

条件 2(归约代价可接受):每核计算时延须远大于归约时延(安全系数 $\alpha = 10$):

+
$$ +T_{MMAD/core} \ge \alpha \cdot T_{Reduce} +\;\Longleftrightarrow\; +\frac{K}{grid_K} \;\gtrsim\; grid_K \times 1690 +$$
+

即 $K \gtrsim grid_K^2 \times 1690$:grid_K=2 → K ≥ 6.8K;grid_K=4 → K ≥ 27K;grid_K=8 → K ≥ 108K;grid_K=32 → K ≥ 1.7M(仅极端 case)。同时 $K / grid_K \ge 256$(单核 K 段过碎则 tiling 效率崩)。grid_K 越大对 K 的要求越苛刻——StreamK 内部的 grid 搜索自然淘汰归约过重的配置。

+

4.4.3 实现方案

+

核间组织:grid_K 个核组成一个归约组,共享同一个输出块 $(\beta, \mu, \nu)$ 的计算。核间分配:

+
$$ +\text{核 } c \text{ 的 K 段:} \left[\frac{c \cdot K}{grid_K},\; \frac{(c+1) \cdot K}{grid_K}\right) +$$
+

核内流水(每核):

+
输入:K 段 [k_start, k_end),输出块 (β, μ, ν)
+  for m_tile / n_tile(核内标准 tiling):
+    for k_tile in range(k_start, k_end, baseK):
+      MTE2: A[β, m_tile, k_tile:k_tile+baseK] → L1
+      MTE1: L1 → L0A/L0B
+      Cube: mmad → L0C 原地累加              # 本核的 K 段内累加
+    Fixpipe: L0C → workspace(GM 或 L2)     # 写出部分和
+

归约阶段

+ +

fixpipe 优化:当 N 不对齐且输出块足够大时($n > 64$ 且 $n \% 16 \ne 0$ 且 $m > 2$ 且 $m \times n \ge 256$),启用 ND_FIXPIPE_1_2(1 AIC : 2 AIV 的 ND fixpipe 通路),由 AIV 分担搬出。

+

grid 搜索:$grid_K$ 从 2 开始递增(2, 4, 8, ...),每档检查条件 2 是否满足;同时检查 $K / grid_K \ge 256$。取满足条件的最小 $grid_K$(归约代价最小)。

+

核间同步:复用 Matmul 高阶 API 的 StreamK 模板(避免手写 CrossCore flagId 冲突,087篇/054篇)。

+

4.4.4 与不切 K 分支的关系

+

StreamK 从不"硬切换"进入:当不切 K 的分支候选已足够快(如已 Cube Bound),StreamK 候选的归约开销使其自然落败;只有当 B/M/N 并行度不足导致大量核闲置时,StreamK 才以数量级优势胜出(例:B=1、M=N=64、K=65536 时,不切 K 仅 16 核可用,StreamK 32 核满负荷)。

+

4.4.5 与源码对照

+

源码 batch_matmul_v3_basic_streamk_tiling.cpp 的条件:

+ + + + + +
源码条件本文对应差异
batchC × mCnt × nCnt ≤ aicNum/2条件 1:$P < C$源码取一半核数为阈值,更保守;本文用精确并行缺口
CeilAlign(k,256) ≥ max(8192, aicNum×256B/dtype)条件 2:$K \gtrsim grid_K^2 \times 1690$ 且 $K/grid_K \ge 256$源码用固定阈值,本文用 grid 搜索动态判定
fp32 非 hf32 时 K ≤ 200 万无(本文不限)源码因 binary accumulation 精度限制 fp32 超长 K;本文不限 dtype,通过 grid 搜索控制归约深度
aivNum == 2 × aicNum归约方式 A 的硬件前提一致
workspace = aicNum×256×256×4B + 20MB归约方式 A 的 workspace一致
+

4.5 ASW_Basic(通用切分框架)

+

4.5.1 做什么

+

不切 K,允许切 B/M/N 的任意组合——它是"不含 K 且含共享读取"的 6 种切分组合({M},{N},{M,N},{B,M},{B,N},{B,M,N})的统一实现框架,也是 B < 32 但 P ≥ 32 时的最优归宿,同时兜住 B ≥ 32 但 IterBatch/MergeBatch 条件不满足的剩余 case。B 可以大、可以小、可以等于 1;ASW_Basic 是实践中最常命中的分支。

+

4.5.2 核心机制:swizzle + L2 管理

+

切 M/N 的固有代价是共享矩阵的重复读,ASW_Basic 用两件武器把代价压到最低:

+

武器 1:ASW 滑窗蛇形 swizzle

+

把 M 向按窗口 $W$ 分组,窗口内蛇形遍历 N 向:

+
$$ +W = \max\{\,d \mid d \mid C,\; d \le \lfloor\sqrt{C}\rfloor \,\} +$$
+

32 核取 $W = 4$。数学效果:同一时刻 32 个核活跃的工作集被压缩到"$W$ 个 A 行块 + 一条 B 列块带",L2 足迹最小 ⇒ 共享读取基本命中 5.2TB/s 的 L2 而非 1.6TB/s 的 GM。

+

**为什么窗口取 $\lfloor\sqrt{C}\rfloor$ 的最大因子**:窗口越接近正方形,A 行块 + B 列块的 L2 足迹越小(方形窗的二维足迹是 $W \cdot K + K \cdot N$,随 W 偏离 $\sqrt{C}$ 而增大),且因子性保证整窗被核数均分、窗口边界不碎。

+

窗口内的遍历顺序:

+
$$ +\sigma(\mu, \nu) = \big(\mu_{row} \cdot W + \mu_{col}\big) \cdot \tilde{N} + \nu' +$$
+

其中 $\mu_{row} = \lfloor \mu / W \rfloor$,$\mu_{col} = \mu \bmod W$,$\nu'$ 由蛇形决定:

+
$$ +\nu' = \begin{cases} \nu & \text{若 } \mu_{row} \text{ 为偶数(正向)} \\ \tilde{N} - 1 - \nu & \text{若 } \mu_{row} \text{ 为奇数(反向)} \end{cases} +$$
+

核号与轮次:

+
$$ +c = \mathrm{idx}(\beta, \mu, \nu) \bmod C,\qquad r = \lfloor \mathrm{idx}(\beta, \mu, \nu) / C \rfloor +$$
+

武器 2:L2 切分 + 错位分核

+

工作集超过 128MB 时,按 mL2TileNum × nL2TileNum 切分。每个 L2 块错位分核(对角线分配):

+
$$ +\sigma_{diag}(\mu, \nu) = \mu \cdot \tilde{N} + \left(\nu + \left\lfloor \frac{\mu \cdot C}{\mathrm{lcm}(\tilde{M}, \tilde{N})} \right\rfloor \right) \bmod \tilde{N} +$$
+

避免多核同时抢同一地址的读读冲突,并优先选拖尾小的方案。

+

4.5.3 核间切分维度的选择顺序

+

按共享代价从低到高:

+ + + + + +
优先级切分方式共享矩阵条件代价
1切 B(B ≥ 核数)$B \ge C$零共享,永远先试
2切 M(B 不够)右矩阵 [K, N]$KN \cdot \text{dtype} \le 128$MB 则驻留 L2右矩阵重复读,L2 吸收
3切 N左矩阵 [M, K]$MK \cdot \text{dtype} \le 128$MB 则驻留 L2左矩阵重复读,L2 吸收
4混合切(B×M、M×N、…)双向共享靠 swizzle + L2 切分管理代价最高,最后试
5降核P 远小于 32 且 K 也不够格走 StreamK宁可部分核闲置(小 case 时延绝对值小)
+

4.5.4 基本块寻优(cubeBound 模型)

+

ASW_Basic 的 baseM/baseN 由 cubeBound 解析模型求解。平台指标:hbmBW/l2BW 由频率×核数×平台速率算出,computePower = 单核算力 × aicNum。

+

cubeBoundEdge 公式:

+
$$ +\text{cubeBoundEdge} = \frac{l2BW}{\text{computePower}} + l2CacheUsage \cdot \left(1 - \frac{l2BW}{hbmBW}\right) \cdot \text{cmr} - \frac{1 + l2BW/hbmBW}{kValue} +$$
+

其中 $\text{cmr} = (m+n)/(m \cdot n)$(每单位输出元素对应的输入行/列搬运代价),$l2CacheUsage = \max(batch \cdot (m+n) \cdot k \cdot \text{dtype} / l2Size, 1.0)$。

+

逐项含义:

+ +

求解流程:在 baseM/baseN 候选解空间内枚举,按"计算访存比 cubeBoundParam = 1/baseM + 1/baseN 与尾块负载均衡率 balanceRate"综合评分:

+ +

进入枚举前 edge 先乘 CUBE_BOUND_RATIO = 0.85 预留 15% 余量。

+

4.5.5 L1 步进与 buffer 管理

+
$$ +\max StepK = \min\left(\left\lceil\frac{K}{baseK}\right\rceil,\; 8\right) +$$
+

stepK ≤ 8 是 issue queue 深度约束(源码注释:"Shape约束 && issue queue约束")。容量约束:

+
$$ +(aL1 + bL1) \times 2 \le L1_{size} - bias - scale +$$
+

且单边不超过 L1 一半:

+
$$ +\max(aL1, bL1) \times 4 \le L1_{size} +$$
+

stepK 枚举优先级:首个可行解 > 256B 对齐 > 单次搬运 ≥ 48KB。

+

4.5.6 内部特化(不是独立分支)

+ +

4.5.7 执行流程

+
核间:32 核,按 swizzle 滑窗分配 (β, μ, ν) 块
+核内(每核):
+  for round in range(R):                          # R = ⌈B̃·M̃·Ñ / C⌉
+    (β, μ, ν) = swizzle_inv(round × C + core_id)
+    for k_tile in range(0, K, baseK):
+      MTE2: A[β, μ, k_tile:k_tile+baseK] → L1(双缓冲/4buffer)
+      MTE1: L1 → L0A/L0B
+      Cube: mmad → L0C 原地累加                  # K 循环不出核
+    Fixpipe: L0C → L2/GM
+

4.5.8 与源码对照

+

源码 batch_matmul_v3_asw_basic_tiling.cpp 的 IsCapable:A、B 的非连续转置状态必须一致;batch 完全相等;batchBias ≤ 1;dtype 限 fp16/bf16 系。DoOpTiling:ResetBase → GetRebalanceBlock → CalL1Tiling。

+ + + + + +
维度源码本文
基本块默认值baseM=256, baseN=256, baseK=128B/dtype由 cubeBound 模型求解,不固定
寻优模型GetRebalanceBlock(cubeBound 解析模型)同(§4.5.4),但源码的经验常数(0.85/0.9/48KB)需注意移植
L1 buffer4 buffer 判定:abL1TensorSize × 4 ≤ l1Size同(§4.5.5)
swizzleASW 滑窗(窗长取 √C 的最大因子)同(§4.5.2)
fp32 splitKK > 1024(或 8192 当 K > 2^28)时启用由 StreamK grid 搜索统一管理
+
+

五、分支仲裁与整体决策流程

+

5.1 为什么需要仲裁层

+

第四章的进入条件给出的是各分支的"主场",但主场之间有重叠区(例如 B ≥ 32 且 M×N 中等时,MergeBatch 与 IterBatch 都合法;B 较大时 IterBatch 与 ASW_Basic 切 B 等效)。重叠区的最优归属没有解析解,必须由统一的时延模型仲裁:

+
$$ +\text{branch}^* = \arg\min_{cand \in \bigcup \text{各分支候选}} \max\big(T_{MMAD},\; T_{MTE2},\; T_{MTE1},\; T_{Fixpipe}\;[,\;T_{Reduce}]\big) +$$
+

分支体系的价值在于候选集完备且无冗余:每个等价类只派一个代表框架生成候选,时延模型在等价类内部和边界上做精细仲裁。两层缺一不可——只有模型没有分支,搜索空间是 15 种组合 × 全部 grid 参数的爆炸;只有分支没有模型,边界 case 被硬阈值误杀。

+

5.2 总决策流程

+
输入:B, M, N, K, dtype, layout, 广播形态, bias
+ │
+ ├─[0] K = 0 → AIV 清零;K = 1 → AIV 逐元素乘(特殊分支)
+ │
+ ├─[1] BatchB = 1 → 转Matmul(免费折叠 [B·M, K],必选)
+ │     BatchA = 1 → 比较"重排 B + Matmul" vs "广播友好 IterBatch/ASW",
+ │                  按时延模型选小者
+ │
+ ├─[2] BatchA = BatchB = B > 1:
+ │     生成四类候选并逐一估算 T_total:
+ │       MergeBatch(条件 1~6 全过才生成,见 4.2.2)
+ │       IterBatch(三大类条件,见 4.3.2)
+ │       ASW_Basic(按 4.5.3 顺序试切 B/M/N/混合,含 swizzle 与 L2 切分)
+ │       StreamK(P < 32 且 K 满足 4.4.2 条件 2 时生成,grid 搜索)
+ │     → argmin T_total
+ │
+ └─[3] 输出:分支 + tiling 参数(baseM/baseN/baseK、b、bL1、kL1…)
+       + swizzle 方案 + L2 切分方案 + 预估端到端时延
+
+

六、完备性审视(漏洞检查)

+

按"转Matmul / MergeBatch / IterBatch / StreamK / ASW_Basic + 前置通路层 cover 全部 BMM case 最优实现"的设计目标,逐条拷问。

+

6.1 覆盖矩阵:无空洞

+ + + + + + + + +
shape 区域主分支兜底
K = 0 / K = 1前置通路层(AIV)
BatchA=1 或 BatchB=1转MatmulASW_Basic(广播友好形态)
BatchA=BatchB>1,B≥32,M×N 大IterBatchASW_Basic 切 B
BatchA=BatchB>1,B≥32,M×N 小,访存 BoundMergeBatchIterBatch
上述但 B<32,P≥32ASW_Basic降核 ASW
P<32,K 大StreamK降核 ASW
P<32,K 也小降核 ASW_Basic(此时时延绝对值小,调度开销主导,分支选择不敏感)
交叉广播(BatchA≠BatchB 且均>1)ASW_Basic
+

6.2 审视发现的问题与修正

+

设计骨架成立,但有五处需要修正或显式声明,否则会被反例拷问:

+

问题 1:转Matmul 的方向不对称(原设想的最大漏洞)。 "一侧 batch=1 就折叠成 Matmul"在 BatchB=1 方向免费(batch 与 M 内存相邻,输出布局一致);但在 BatchA=1 方向,右矩阵 [B,K,N] 折叠成 [K,B·N] 需要一次真实的转置重排,且输出 [M,B·N] 与目标 [B,M,N] 之间存在置换——"只是输出后按 batch split"的设想在 ND 布局下不成立,split 实际是 scatter。修正:BatchA=1 时按"重排代价 vs Matmul 收益"走模型仲裁,小 case 留在 BMM 分支内。

+

问题 2:交叉广播 case 的归属必须显式声明。 BatchA=(6,1)、BatchB=(1,6) 这类 case 既不进转Matmul(两侧都>1),也不满足 MergeBatch/IterBatch 的 BatchA=BatchB 前提——它们落入 ASW_Basic(对广播侧做驻留)。不声明就会显得"五分支有洞"。声明后完备性闭合。

+

问题 3:K 退化不在五分支内,需前置通路层。 K=0(无计算)、K=1(无累加深度)走 AIV 向量通路;这是"计算通路选择"维度,与"切分维度"正交,放在五分支之前判断。源码用 K_EQUAL_ZERO / TO_MUL 两个独立策略处理,印证了这一层的必要性。

+

问题 4:IterBatch 草稿条件的字面表述有误杀。 B mod C > minCoreNum 会拒绝 B 整除核数的完全均衡 case(余数 0)。正确表述应以负载均衡比为准:访存 Bound 要求 $b_{Avg}/b_{Max} > 0.8$,计算 Bound 要求 100%;余数规则只是它的近似操作化。

+

问题 5:经验常数必须可标定,且源码硬编码条件偏粗。 $b_{thr}=4$、$k_{thr}$=32B/dtype、480KB、16KB、dValue≥256B 都是该档芯片的实测经验值,文档中全部显式参数化(换芯片只换常数表)。对照之下源码把 batchC ≥ 4×aicNumalignK ≥ 64M ≤ N 等直接写死,其中 alignK ≥ 64 比草稿阈值保守 4 倍、M ≤ N 放弃了镜像 case——做最优实现时不应继承这些硬编码,而应按第四章的条件体系重新判定

+

6.3 与源码策略的关系

+

源码 arch35 的 11 个策略(K_EQUAL_ZERO / TO_MUL / STREAM_K / MERGE_BATCH / ITER_BATCH_BROADCAST / ITER_BATCH_BASICAPI / ITER_BATCH / AL1_FULL_LOAD / BL1_FULL_LOAD / ASW_BASIC / BASE)与本文分支不是同层概念:

+ + + + + + +
源码策略本文归属
K_EQUAL_ZERO、TO_MUL前置通路层(§4.0),正交于切分
STREAM_KStreamK(§4.4)
MERGE_BATCHMergeBatch(§4.2)
ITER_BATCH、ITER_BATCH_BASICAPI、ITER_BATCH_BROADCASTIterBatch(§4.3;broadcast 变体是广播输入下的数据复用特化;单边广播在本文更多由转Matmul 吸收)
AL1_FULL_LOAD、BL1_FULL_LOADASW_Basic 的 L1 驻留参数极限(§4.5.6)
ASW_BASIC、BASEASW_Basic(§4.5)
+

结论:源码策略集 = 本文分支 × 计算通路 × 驻留/广播特化的展开。本文的分支是"切分等价类"的最小完备骨架,源码的冗余策略可在新实现中收敛,源码遗漏的(条件 3/4/6 的显式判定、M>N 的 MergeBatch 镜像、转Matmul 的 BatchA=1 方向)正是优化的增量空间。

+
+

七、总结

+
    +
  1. 最优的定义:$T_{total} = \max(T_{MMAD}, T_{MTE2}, T_{MTE1}, T_{Fixpipe}[, T_{Reduce}])$ 最小;瓶颈交换是合法且必要的手段。
  2. +
  3. 分支的推导:4 维可切 → 15 种组合完备 → 切分特征表(切 B 免费 / 切 M/N 廉价被 L2 吸收 / 切 K 昂贵需归约)→ 按价格从低到高购买并行度 → 坍缩为 MergeBatch、IterBatch、ASW_Basic、StreamK 四个等价类,加前置的转Matmul 问题归约层与 K 退化通路层,共七大路径,完备且极小
  4. +
  5. 各分支的条件都不是孤立经验,而是"资格(并行度够不够)→ 上限(容量/算存比封顶)→ 下限(搬移效率托底)"三层逻辑的实例化;所有经验常数可标定、可移植。
  6. +
  7. 边界不靠硬切:重叠区由端到端时延模型统一仲裁,分支体系负责候选集的完备无冗余。
  8. +
  9. 设计经受了漏洞拷问:修正了转Matmul 的方向不对称与 IterBatch 均衡条件的字面误杀,显式声明了交叉广播与 K 退化的归属;相对源码实现,本文条件体系更细、更真、覆盖更全。
  10. +
+
+

附录 A:经验常数表(可标定)

+ + + + + + + + + + +
常数数值出处/依据说明
$b_{thr}$4源码 MIN_BATCH_L0合并搬移的最小有效合并数
$k_{thr}$32B/dtype§2.2 第 4 条ND2NZ dValue 下限(BF16 为 16 元素)
$min\_DatamountPerCore$480KB§2.2 第 2 条单核搬移总量下限
$min\_TileSize$16KB§2.2 第 3 条单次搬移 tile 下限(168篇官方为 16KB)
$minCoreNum$26(≈0.8×32)§2.2 第 1 条尾波活跃核数下限
$CUBE\_BOUND\_RATIO$0.85源码 GetRebalanceBlockcubeBound edge 余量
$balanceRateEdge$0.9源码 GetRebalanceBlock负载均衡剪枝阈值
均衡率(访存 Bound)0.8源码 iterbatch$b_{Avg}/b_{Max}$ 下限
$\alpha$(StreamK 安全系数)10本文设定$T_{MMAD} \ge \alpha \cdot T_{Reduce}$
$R$(BF16/FP16)607.5 FLOP/元素§2.3芯片算存比平衡点
+

附录 B:与 issue#3 (v0.3) 的差异

+ + + + + + + + +
章节v0.3 状态v0.4 补全内容
StreamK 实现方案§4.4.3 完整实现:核间组织、核内流水、归约两方式、fixpipe 优化、grid 搜索、核间同步
ASW_Basic 实现方案§4.5.2~4.5.7 完整实现:swizzle 公式、L2 切分、切分顺序、cubeBound 模型、L1 步进、执行流程
MergeBatch 进入条件5 条(简略)6 条(含公式推导),补充算存比约束 $b_{AI}$、L0C 解析式 $b_{L0C}$、乒乓取舍
IterBatch 实现方案部分伪代码补全 (a)~(e) 五种情形的完整伪代码 + 设计原理 + 源码对照
转Matmul 实现方案简略补充 BatchA=1 的决策规则(A 大小判断 + 广播扩展 + 时延比较)
特殊分支简述补充 AIV vs AIC 的定量分析、源码触发条件
源码对照每个分支新增"与源码对照"小节,指出源码的保守/缺口/一致项
完备性论证§3.6 + 第六章:15 种组合 → 7 路径的完备性与极小性论证
+
+

*参考文档:《BMM分块计算数学公式》《BMM最优软件实现方案设计》《BatchMatMulV3算子分支实现分析》《BatchMatmul算子特性分析》;源码:ops-nn/matmul/batch_mat_mul_v3(arch35);芯片资料:昇腾 950PR 架构白皮书与 CANN 9.0.0 性能建模文档。*

+
+ + \ No newline at end of file