1. 项目概述与核心价值在嵌入式数字信号处理领域尤其是在对功耗、成本和实时性有严苛要求的场景里如何将复杂的算法高效地“翻译”成芯片能理解并快速执行的指令是每个工程师都会面临的硬核挑战。今天我想结合自己过去在通信和音频处理项目中的实战经验深入聊聊两个经典算法——自适应滤波LMS和卷积编码——在TI C55x DSP平台上的实现细节。自适应滤波的核心思想是“动态追踪”它不像传统FIR滤波器那样系数固定而是能根据输入信号和期望信号的误差实时调整自身系数从而在噪声环境中“揪出”我们想要的信号或者把不想要的干扰“抵消”掉。LMS算法则是实现这一思想的“实干家”它用最直接的梯度下降法来更新系数计算量小结构清晰非常适合在DSP上跑。而卷积编码则是通信链路里的“纠错卫士”它通过给原始数据比特加入有规律的冗余让接收端在信号受到干扰时依然有办法把原始信息还原出来其编码过程本质上是一个基于移位寄存器的有限状态机。为什么选C55x DSP来聊因为它的指令集设计非常“聪明”为这类信号处理算法量身定制了专用指令。比如一条LMS指令就能同时完成乘累加和系数更新一条MAXDIFF指令就能搞定维特比解码中最耗时的“加-比-选”操作。理解这些指令的用法不仅仅是写几行汇编代码更是理解如何让算法和硬件特性深度结合榨干芯片的每一分性能。接下来我会从算法原理、C55x的指令级优化再到具体的代码实现和避坑指南带你走完从理论到落地的全过程。2. 自适应滤波LMS算法原理与C55x实现精讲2.1 延迟LMS算法核心思想拆解标准的LMS算法系数更新公式是w(n1) w(n) μ * e(n) * x(n)。这里w是滤波器系数向量μ是步长收敛因子e(n)是当前时刻的误差x(n)是当前的输入向量。这个公式要求在同一时刻n我们既要用系数w(n)算出输出y(n)和误差e(n)又要立刻用这个误差去更新系数用于下一个采样点。但在实际的流水线处理器如C55x中这可能会带来数据冲突或增加流水线停顿。因此工程上更常用的是延迟LMS。它的核心改动在于用上一个时刻的误差e(n-1)来更新当前时刻的系数。即w(n1) w(n) μ * e(n-1) * x(n-1)。这个“延迟”带来了什么好处它解耦了滤波和更新之间的严格时序依赖。在计算y(n)需要w(n)的同时我们可以并行地准备用于更新w(n1)的数据即e(n-1)和x(n-1)使得处理器流水线能被更充分地利用指令调度更灵活从而提升整体吞吐率。当然代价是算法的收敛速度会受微小影响但在大多数实时系统中这点延迟是可接受的。2.2 C55x的LMS专用指令深度剖析C55x DSP为延迟LMS算法量身打造了lms指令其语法为lms(Xmem, Ymem, ACx, ACy)。这条指令在一个周期内并行完成两件大事乘累加MACACy ACy (Xmem * Ymem)。这用于计算FIR滤波器的输出。系数更新带舍入的加法ACx rnd(ACx (Xmem #16))。这用于更新滤波器系数。这里面的门道很多我们逐一拆解操作数角色分配通常Xmem指向系数数组b[k]Ymem指向数据延迟线数组x[n-k]。ACy累加器被清零后用于累加计算滤波器输出y[n]。ACx累加器则预先装载了μ * e[n-1]即步长与上一时刻误差的乘积用于系数更新。更新公式的映射指令中的(Xmem #16)非常关键。在C55x中数据通常是Q15格式1位符号15位小数。Xmem即系数b[k]是16位。左移16位后它被放置到ACx累加器的高32位区域ACx(39:16)。然后与ACx中原有的μ * e[n-1]相加。因为μ通常也是一个Q15格式的小数μ * e[n-1]的结果是Q30格式。b[k]左移16位后也变成了Q30格式两者相加再经过舍入rnd操作结果被写回ACx的高位。之后我们通过mov HI(ACx), *ar_coef这样的指令将ACx的高16位即更新后的Q15格式系数存回系数内存。这个过程完美地实现了定点数下的系数更新b[k] b[k] μ * e * x。“延迟”体现在哪里在代码中误差项e[n-1]是在上一个样本处理周期结束时计算并保存到T3寄存器或某个内存单元的。在本周期开始时我们将其与步长μ相乘得到μ * e[n-1]并放入ACx然后才开始本轮的滤波和系数更新循环。这就是“延迟”的具体实现。实操心得数据对齐与指针初始化使用lms指令时确保系数数组b[]和数据数组x[]在内存中是长字对齐地址是4的倍数的这能保证双存储器操作数Xmem, Ymem访问的最高效率。指针ar_coef和ar_data在初始化时也应指向各自数组的起始地址。一个常见的坑是在循环开始前没有正确预装第一个误差项到ACx或者没有将ACy清零导致第一次滤波输出或系数更新出错。务必在循环外做好初始化。2.3 延迟LMS内核循环代码逐行解读让我们结合文档中的示例代码看看一个高效的延迟LMS自适应滤波器内核是如何运作的。这里我将其转换为更易读的代数指令格式并加上详细注释; 初始化阶段 AC1 #0 ; 初始化误差寄存器实际是μ*e初始为0 || localrepeat { ; 开始外循环处理每个输入样本 *ar_data *ar_input ; 将新输入样本x[n]移入延迟线缓冲区 T3 HI(AC1) ; 将上一轮计算得到的μ*e[n-1]存入T3准备用于系数更新 ; 内循环初始化计算第一个抽头的贡献并准备更新项 AC0 T3 * *ar_data ; AC0 μ*e[n-1] * x[n-1] (第一个更新项) || AC1 #0 ; 同时清零AC1用于累加当前输出y[n] ; 核心内循环完成滤波输出计算和系数更新 LMS(*ar_coef, *ar_data, AC0, AC1) ; 首次LMSAC1开始累加输出AC0更新第一个系数 || localrepeat { ; 内循环处理剩余N-1个抽头 *ar_coef HI(AC0) ; 存储上一个抽头更新后的系数 || AC0 T3 * *ar_data ; 并行计算下一个抽头的更新项 μ*e[n-1] * x[n-k] LMS(*ar_coef, *ar_data, AC0, AC1) ; 核心AC1累加AC0更新 } ; 内循环结束后处理 *ar_coef HI(AC0) ; 存储最后一个抽头更新后的系数 || *ar_output HI(rnd(AC1)) ; 存储当前滤波器的输出y[n]舍入到16位 ; 计算当前时刻的误差 e[n] d[n] - y[n] AC2 (*ar_des #16) - AC1 ; 将期望信号d[n]左移16位转为Q31后减去输出y[n] || mar(*ar_data) ; 调整数据指针可能用于滑动延迟线 ; 为下一个样本计算 μ*e[n]存入AC1供下一轮使用 AC1 rnd(T_step * AC2) ; T_step μ计算 μ * e[n]舍入 } ; 外循环结束代码逻辑流梳理样本输入与误差准备移入新数据取出上一轮的μ*e。滤波与更新循环这是最核心的部分。内循环巧妙地交织了系数存储、下一更新项计算和当前的LMS操作。注意LMS指令中的ACx用于更新和ACy用于输出累加是分开的实现了并行。收尾与误差计算循环结束后存储最后一个系数和当前输出。然后计算当前输出与期望信号的误差并预计算下一轮需要的μ*e。避坑指南指针管理与缓冲区延迟线管理ar_data指针通常指向一个循环缓冲区延迟线。在每次处理完一个样本后需要更新指针以丢弃最旧的数据x[n-N]并为新样本x[n1]腾出位置。文档代码中的mar(*ar_data)可能与此相关但更完整的实现可能需要模寻址circular addressing或手动缓冲区滑动。系数指针复位内循环结束后ar_coef指针指向了系数数组末尾。在下一次外循环开始前必须将其重置回系数数组的起始地址。示例代码中这一点被隐含了在实际编写时绝对不能忘记否则会导致指针越界和错误的数据访问。步长μ的选择T_step即μ的值至关重要。太大可能导致算法发散不稳定太小则收敛缓慢。它需要根据输入信号的功率进行归一化一个经验公式是μ 2 / (N * 输入功率)。在实际系统中常通过实验或自适应调整来确定。3. 卷积编码算法原理与比特流操作3.1 卷积编码器一个移位寄存器的故事卷积编码可以形象地理解为一个带有抽头的移位寄存器。假设约束长度K5寄存器有K-14级。每输入一个信息比特寄存器就右移一位。编码器有多个输出比如2个对应码率R1/2每个输出都是当前输入比特与寄存器中某些特定位置比特的模2加异或运算。文档中的图7-3清晰地展示了这一点。生成多项式G0 1 D^3 D^4和G1 1 D D^3 D^4其中D代表延迟定义了抽头位置。G01D^3D^4意味着第一个输出G0是当前输入比特、3个时刻前的比特和4个时刻前的比特三者异或的结果。这种结构使得当前的输出不仅取决于当前输入还取决于过去的K-1个输入引入了“记忆性”从而具备了纠错能力。编码器状态可以用移位寄存器的内容来表示对于K5有2^(K-1)16种可能状态这引出了后续维特比解码的网格图Trellis。3.2 C55x的比特域操作指令field_expand与field_extract在卷积编码后多个输出比特流如G0和G1需要被复用成一个高速串行流进行传输在接收端则需要从接收到的串行流中解复用出各个分支的比特流。C55x提供的field_expand比特域扩展和field_extract比特域提取指令就是为了高效完成这类比特级交织/解交织操作而生的。这两条指令都使用一个16位的掩码k16来控制操作。掩码的每一位对应目标寄存器对于field_expand或源寄存器对于field_extract的一个比特位。field_expand(ACx, k16)复用/交织。它扫描掩码k16的每一位从LSB到MSB如果该位为1则将源累加器ACx从LSB开始的当前比特复制到目标累加器的当前比特位然后ACx的比特指针和目标的比特指针都加1。如果该位为0则只将目标的比特指针加1ACx的指针不动。结果就是根据掩码k16定义的“空洞”模式将ACx中的连续比特流“喷洒”到目标寄存器的指定位置上。field_extract(ACx, k16)解复用/解交织。操作与field_expand相反扫描掩码k16。如果该位为1则将源累加器ACx的当前比特复制到目标累加器的当前比特位从LSB开始然后ACx和目标的比特指针都加1。如果该位为0则只将ACx的比特指针加1目标指针不动。结果就是从ACx中根据掩码“收集”比特形成连续的目标比特流。示例解析假设我们要将G0和G1两个比特流交织成G0G1模式是G0, G1, G0, G1,...即比特交错。发送端复用G0G1 field_expand(G0, #5555h) ; 掩码0101 0101 ...将G0比特放到目标的位置0,2,4,6... Temp G0G1 ; 暂存 G0G1 field_expand(G1, #AAAAh) ; 掩码1010 1010 ...将G1比特放到目标的位置1,3,5,7... G0G1 G0G1 | Temp ; 合并形成完整的交织流接收端解复用G0 field_extract(G0G1, #5555h) ; 用同样的掩码#5555h提取出G0流 G1 field_extract(G0G1, #AAAAh) ; 用掩码#AAAAh提取出G1流注意事项数据准备与对齐使用这些指令前需要将比特流数据正确地装入40位的累加器如AC0。通常我们会将16个比特对应一个short型数据放入累加器的低16位AC0(15:0)。掩码k16的设计是灵活性的关键它可以实现任意的交织图案不仅限于简单的奇偶交错。务必根据通信标准规定的交织器结构来设计掩码。4. 维特比解码算法与C55x的蝶形运算加速4.1 维特比算法在网格图中寻找最优路径卷积编码在接收端需要用维特比算法进行最大似然序列检测。其核心是“加-比-选”ACS操作在编码器的网格图上进行。对于每个状态在每一个时间点会有两条路径对应输入比特0或1汇聚过来。维特比算法需要加Add计算这两条路径的累积路径度量Path Metric即旧路径度量加上从上一状态转移到当前状态的分支度量Branch Metric通常用接收符号与期望符号之间的“距离”如汉明距或欧氏距离的负值来表示。比Compare比较这两条新路径的度量。选Select选择度量较大对于最大似然或较小对于最小距离的一条作为幸存路径并记录选择结果即“追溯比特”。对于R1/2, K5的GSM编码器有16个状态。由于网格图的对称性这16个状态的ACS操作可以两两配对形成8个“蝶形”运算单元如图7-6所示。每个蝶形运算同时计算两个新状态如j和j8的路径度量。4.2 C55x的蝶形运算专用指令套件C55x为高效实现维特比蝶形运算提供了三条强大的指令ADDSUB、SUBADD和MAXDIFF。ADDSUB和SUBADD这两条指令充分利用了C55x累加器可进行双16位运算的特性。它们在一个周期内同时对累加器的高16位和低16位进行不同的操作。ADDSUBHI(ACy) HI(ACx) Smem,LO(ACy) LO(ACx) - Smem。即高16位加低16位减。SUBADDHI(ACy) HI(ACx) - Smem,LO(ACy) LO(ACx) Smem。即高16位减低16位加。在维特比蝶形中Smem通常存放本地距离LD或-LDACx的高低位分别存放两个旧路径度量。这样一条指令就能同时计算出两条候选新路径的度量。MAXDIFF这是实现“比-选”的关键。指令MAXDIFF AC0, AC1, AC2, AC3比较AC0和AC1的高16位将较大者存入AC2的高16位比较AC0和AC1的低16位将较大者存入AC2的低16位。同时它将比较结果哪个更大以比特形式记录到两个独立的追溯寄存器TRN0和TRN1中TRN0对应高16位比较结果TRN1对应低16位。AC3在比较过程中被用作临时寄存器。4.3 维特比蝶形运算宏解析文档中的BFLY_DIR_MNEM和BFLY_REV_MNEM宏定义展示了完整的蝶形运算。我们以BFLY_DIR_MNEM为例结合GSM的蝶形结构来看BFLY_DIR_MNEM .MACRO ; 假设: T3 本地距离 LD ; AR5指向包含两个旧路径度量的内存位置Old_Met(2*j) 和 Old_Met(2*j1) ; AR3和AR4分别指向存储新路径度量低位和高位结果的内存位置 ADDSUB T3, *AR5, AC0 ; AC0高Old_Met(2*j)LD, AC0低Old_Met(2*j1)-LD SUBADD T3, *AR5, AC1 ; AC1高Old_Met(2*j)-LD, AC1低Old_Met(2*j1)LD MAXDIFF AC0, AC1, AC2, AC3 ; 比较并选择AC2高max(AC0高,AC1高), AC2低max(AC0低,AC1低) MOV AC2, *AR3, *AR4 ; 存储新路径度量低位到*AR3高位到*AR4 .ENDM运算对应关系参考图7-6对于新状态New_Metric(0)其两条候选路径来自旧状态Old_Metric(0)输入比特0和Old_Metric(1)输入比特1。对应的分支度量分别是LD和-LD。因此AC0的高16位Old_Met(2*j)LD和AC1的高16位Old_Met(2*j)-LD正好对应了这两条路径的新度量。MAXDIFF选择较大的一个存入AC2的高16位即New_Metric(0)。对于新状态New_Metric(8)其两条候选路径来自相同的两个旧状态但分支度量符号相反。这正好对应了AC0的低16位Old_Met(2*j1)-LD和AC1的低16位Old_Met(2*j1)LD。MAXDIFF选择较大的一个存入AC2的低16位即New_Metric(8)。一条MAXDIFF指令配合TRN0/TRN1的记录就完成了一个蝶形单元中两个状态的全部ACS操作效率极高。核心技巧软件流水与指令并行文档中的例7-7展示了如何通过软件流水和用户自定义并行来进一步加速维特比解码循环。其核心思想是将当前蝶形运算的加载操作如加载新的本地距离LD与上一个蝶形运算的存储操作和比较操作重叠起来。localrepeat { ADDSUB T3, *AR0, AC0 ; 当前蝶形计算1 || MOV *AR5, AR7 ; 并行为下一个蝶形加载地址指针 SUBADD T3, *AR0, AC1 ; 当前蝶形计算2 || MOV *AR6, T3 ; 并行加载下一个蝶形的本地距离LD到T3 MOV AC2, *AR2, *AR2(T0) ; 存储上一个蝶形的结果 || MAXDIFF AC0, AC1, AC2, AC3 ; 并行比较当前蝶形 ... ; 继续下一个蝶形 }注意看在计算当前蝶形使用T3中的LD的同时我们已经在为下一个蝶形加载新的LD到T3寄存器。同时存储的是上一个蝶形MAXDIFF的结果。这种“提前加载”和“延迟存储”交织在一起形成了高效的流水线消除了数据依赖带来的停顿是DSP性能优化的经典手法。实现时需要精心安排指针和寄存器确保数据流正确无误。5. 从C代码到优化汇编DSPLIB的角色与使用策略5.1 为什么需要DSPLIBTI C55x DSPLIB是一个用高度优化的汇编语言编写的函数库提供了超过50个C语言可调用的通用信号处理函数。对于大多数开发者而言直接使用DSPLIB有三大好处极致的性能这些函数由TI的专家手工汇编优化充分挖掘了C55x的双MAC单元、并行指令、循环缓冲等硬件特性其执行速度远超用C语言编写的同等功能代码。显著的开发效率无需从头实现复杂的算法直接调用API即可大幅缩短开发周期。可靠的正确性库函数经过充分测试并与MATLAB脚本进行过比对保证了算法的正确性。5.2 关键数据类型与调用约定DSPLIB主要使用两种定点数据类型DATA(Q15格式)16位有符号整数表示范围为[-1, 1-2^(-15)]。这是最常用的格式对应C语言的short类型。LDATA(Q31格式)32位有符号整数表示范围更精确对应C语言的long类型。在调用DSPLIB函数时必须注意其数据组织方式向量存储所有向量元素必须连续存储步长为1。复数存储采用交错存储格式即[实部0, 虚部0, 实部1, 虚部1, ...]。原位计算很多函数支持源操作数和目的操作数是同一块内存以节省内存空间。5.3 调用示例与混合编程要点一个典型的C语言调用示例如下#include dsplib.h // 1. 包含头文件 DATA x[3] {12398, 23167, 564}; DATA r[NX]; DATA rexp[NX]; float rf1[NX]; void main() { short i; // 2. 调用DSPLIB函数例如计算倒数 recip16(x, r, rexp, NX); // 3. 可能的数据格式转换 q15tofl(r, rf1, NX); // ... 后续处理 }使用前必须在工程中包含dsplib.h头文件。链接时加入55xdsp.lib库文件。配置好链接器命令文件.cmd正确分配程序段和数据段到DSP的物理内存。关于从汇编调用虽然DSPLIB是C可调用的理论上也可以从汇编调用但强烈不建议。因为C调用约定会带来额外的栈帧管理和参数传递开销在纯汇编环境中直接实现核心算法循环或者将DSPLIB汇编源码集成到自己的汇编模块中是更高效的做法。5.4 查找与参考样例DSPLIB的安装目录下例如c:\ti\C5500\dsplib\examples为每个函数都提供了完整的样例工程包括_t.c测试主函数。test.h由MATLAB脚本生成的测试输入和期望输出数据。test.c/ftest.c/ltest.c比较函数输出与期望值的验证代码。55x.cmd链接器命令文件示例。在实现自己的算法时首先查阅这些样例是快速上手和验证理解的最佳途径。更重要的是可以仔细研读DSPLIB中相关函数的汇编源码位于55x_src目录例如dlms.asm或卷积解码相关的函数这是学习C55x汇编优化技巧的“活教材”。6. 常见问题、调试技巧与性能优化实录6.1 自适应滤波LMS实现中的典型问题算法发散输出或系数溢出症状滤波器输出或系数值变得异常大最终饱和。排查首要检查步长μ这是最常见的原因。使用公式μ 2 / (N * 输入功率)进行估算并在实际系统中通过实验选择一个更保守的值例如估算值的1/10开始尝试。检查输入信号幅度确保输入信号x[n]和期望信号d[n]在Q15格式的表示范围内-1到~0.9999。如果信号来自ADC需要正确进行缩放。启用饱和模式在系数更新指令如lms前后可以临时开启累加器的饱和保护位SATD但这只是防止异常扩散治标不治本。收敛速度慢或稳态误差大排查步长μ太小在保证稳定的前提下适当增大μ。滤波器阶数N不足无法充分建模系统。期望信号d[n]不准确在回声消除中d[n]是近端语音回声如果近端语音太强会影响对回声路径的建模。指针错误导致数据错乱症状输出结果完全无规律或程序跑飞。排查单步调试在CCS中单步运行观察ar_coef和ar_data指针在循环前后的值是否符合预期。确保系数指针在内循环结束后被正确复位。内存查看在算法运行前后导出系数数组和延迟线数组的内容检查是否被意外修改。边界检查确保循环次数N抽头数与数组声明的大小匹配防止指针越界。6.2 卷积编码与维特比解码调试要点编码输出与标准不符验证生成多项式这是根本。确认代码中实现的多项式抽头位置与标准文档完全一致。一个比特的差异会导致整个编码序列错误。检查移位寄存器初始化编码器开始工作前移位寄存器通常需要初始化为全零状态。确保在编码每个帧之前进行了正确的初始化。验证field_expand掩码用一组已知的输入比特和G0/G1输出手动计算交织后的比特流与程序输出对比确认掩码k16使用正确。维特比解码性能差误码率高分支度量计算这是解码器的“尺子”。对于软判决解码接收到的符号是多个比特的量化值如3比特软信息。需要将软信息映射为分支度量通常是相关值或负距离。确保这个映射关系正确且度量值范围适合累加器避免过快溢出。路径度量归一化路径度量在ACS过程中会不断累加可能导致溢出。常见的做法是定期找到所有路径度量中的最小值然后所有度量减去这个值归一化。需要实现一个高效的查找最小值并归一化的子程序。追溯深度追溯深度Traceback Depth通常取约束长度K的5到10倍。深度太浅解码性能下降深度太深增加延迟和存储开销。需要通过仿真确定一个合适的值。TRN0/TRN1管理MAXDIFF指令会更新TRN寄存器。你需要将每个时间点的TRN值保存到追溯记忆数组中。注意TRN0和TRN1分别对应蝶形运算中两个输出的选择比特。存储时需要按照状态顺序正确组织这些比特。6.3 C55x汇编级性能优化关键技巧充分利用双MAC和并行指令C55x有两个乘累加单元。在FIR滤波、复数乘法等计算密集型循环中通过安排数据确保每个周期能发射两条MAC指令或LMS指令。仔细研究指令集将不冲突的指令安排在同一行用||并行执行。例如一个存储操作完全可以和一个算术逻辑运算并行。软件流水消除循环依赖这是将循环性能推向极限的技术。如维特比示例所示核心思想是将当前迭代的加载、上一个迭代的计算、上上一个迭代的存储重叠起来。实现步骤a) 创建循环体核b) 添加循环前代码填充流水线c) 添加循环后代码排空流水线。编译器在高级别优化-o3 -pm时会自动进行软件流水但对于最核心的循环手动汇编优化往往能做得更好。明智使用循环缓冲与模寻址对于延迟线、滤波器系数等需要循环访问的数据结构务必使用C55x的循环缓冲Circular Buffer功能。正确设置缓冲区起始地址BSAxx、大小BKxx和指针的循环配置位ARxLC。这能省去大量手动检查指针边界和回绕的指令开销。在汇编中使用.arms_off伪指令告知汇编器你正在使用DSP模式而非ARM兼容模式以便使用模寻址。数据对齐与内存访问优化长字对齐对于32位数据long或两个连续的short访问确保其地址是4的倍数。这能使dbl()内存访问指令在一个周期内完成。双操作数访问*ARx和*ARy这类双操作数访问要求数据对象在内存中按特定方式对齐通常是长字对齐并且两个辅助寄存器ARx和ARy属于不同的辅助寄存器组。合理规划数据布局可以最大化利用这种高效访问方式。合理分配寄存器与减少内存访问将最内层循环频繁使用的变量如循环计数器、步长μ、当前误差等分配给寄存器如T0-T3,AC0-AC3,AR0-AR7。避免在循环内部进行不必要的内存加载/存储。例如在LMS循环中误差项μ*e应始终保存在寄存器如T3中。在我实际调试一个语音降噪项目的LMS滤波器时就曾因为指针复位问题折腾了大半天。算法看似在工作但噪声抑制效果总是不稳定。后来单步跟踪发现在处理完一帧数据后系数指针没有回到起始位置导致下一帧数据使用了错误的系数开头。这个教训让我养成了一个习惯在编写任何涉及指针循环的汇编代码后立刻在循环入口和出口处打上断点验证指针的初始值和终值是否符合预期。对于维特比解码初期最大的挑战是追溯记忆的管理和路径度量的归一化。我当时的做法是先用C语言实现一个功能正确的、但速度慢的版本作为“黄金参考”然后逐步将核心的ACS蝶形运算用汇编替换并逐阶段比对中间结果路径度量、追溯比特确保汇编优化没有引入逻辑错误。这个过程虽然繁琐但能从根本上保证最终优化代码的正确性。
