好的，言归正传！一般的FIR算法如下：
<!--[if !supportEmptyParas]--> <!--[endif]-->

那么如何使用SSE指令来优化它呢？

SSE的指令在于它使用了向量计算，极大的提高了运算速度，所以在优化之前，我们首先要对原有的c代码进行向量构造。因为水平有限，构造了一个粗糙向量（如果有那位网友能够指导我进行优化，我将万分感谢！）。

这里假设相关的所有长度都是16个字节的倍数。

<!--[if !supportEmptyParas]--> <!--[endif]-->

上面方法的主要思路是，在保持原有算法的基础上，对代码进行展开和重组，形成若干个向量组。在这个基础上，我们再使用SSE指令进行快速的矢量计算以获取高性能。从上的构建可以看出，我们有两组向量，分别用in和coeff指向。in指向的向量in[j]是16个字节地址对齐，而in[j+1], in[j+2], in[j+2]是不对齐的，故除了第一组能够使用movaps之外，其他都使用movups来装载。coeff仅使用一组向量和in的4组向量进行相乘，所以可以考虑将它放入xmm寄存器中提高运行效率。从上面的c代码可以看出，对于每组向量计算的结果将被放在4个单元中（[127,96], [95,64],[63:32], [31:0]），所以在计算最后要将它们水平相加后再放入out数组中。

<!--[if !supportEmptyParas]--> <!--[endif]-->

相关的汇编代码如下：

<!--[if !supportEmptyParas]--> <!--[endif]-->

为了防止与汇编保留字冲突，我将原有的C代码中的部分变量进行了重新命名。此外，汇编代码假设他们的首地址都是16字节对齐的。16字节地址对齐的方法很简单，只要在变相声明的前面加上“__declspec( align( 16 ) )”修饰，比如：

它告诉编译器，上面声明的变量必须16字节地址对齐J。

<!--[if !supportEmptyParas]--> <!--[endif]-->

小结：

上面的汇编代码有几个较大的缺陷：

1. 没有对水平相加进行优化，SSE3有一个HADDPS的水平相加指令，它对性能的提升有重大作用，因为我现在无法调试（计算机仅仅支持SSE2）所以使用了上面一对难看的指令。

2. 没有对out向量进行一次向量写操作（在这个实现中，用了很傻的4次标量写操作）。主要原因是还没有充分理解所有相关的SSE指令，并将它们融会贯通。

3. 对原有代码的向量构造不够好，我想一定有一种方法可以避免最后的差劲的水平向量加操作和4次标量写操作。我的理想状态是没有水平向量相加和一次向量写操作。

<!--[if !supportEmptyParas]--> <!--[endif]-->

我想，随着学习的深入，这些方面将会被逐一解决。我将在后继的文章中对这些遗憾的方面进行逐步探讨以获取更好的性能。如果有大侠能够出手相救，给我醍醐灌顶的提示，小弟在此叩拜了J！

<!--[if !supportEmptyParas]--> <!--[endif]-->