zoukankan      html  css  js  c++  java
  • 内存拷贝的优化方法

    在复杂的底层网络程序中,内存拷 贝、字符串比较和搜索操作很容易成为性能瓶颈所在。编译器自带的此类函数虽然做了一些通用性的优化工作,但因为在使用指令集方面受到兼容性的约束,远远没 有达到最大限度利用硬件能力的地步。而通过针对特定硬件平台的优化,可以大大提高此类操作的性能。下面我将以P4平台下内存拷贝操作为例,根据AMD提供 的一份优化文档中的例子,简要介绍一下如何通过特定指令集,优化内存带宽的使用。虽然因为硬件限制没有达到AMD文档中所说memcpy函数300%的性 能提升,但在我机器上实测也有%175-%200的明显性能提升(此数据可能根据机器情况不同)。

    Optimizing Memory Bandwidth from AMD

    按照众所周 知的“摩尔”定律,CPU的运算速度每18个月翻一翻,但与此同时内存和外存(硬盘)的速度并无法达到同步增长。这就造成高速CPU与相对低速的内存和外 设之间的不同步发展,成为很多程序的瓶颈所在。而如何最大限度提升对现有硬件的利用程度,是算法以下层面优化的主要途径。对内存拷贝操作来说,了解和合理 使用Cache是最关键的一点。为追求性能,我们将以牺牲兼容性为代价,因此以下讨论和代码都以P4及以上级别CPU为主,AMD芯片虽然实现上有所区 别,但在指令集和整体结构上相同。

    首先我们来看一个最简单的memcpy的汇编实现:

    以下为引用:

    ;
    ; Flier Lu (flier@nsfocus.com)
    ;
    ; nasmw.exe -f win32 fastmemcpy.asm -o fastmemcpy.obj
    ;
    ; extern "C" {
    ; extern void fast_memcpy1(void *dst, const void *src, size_t size);
    ; }
    ;
    cpu p4

    segment .text use32

    global _fast_memcpy1

    %define param esp+8+4
    %define src param+0
    %define dst param+4
    %define len param+8

    _fast_memcpy1:
    push esi
    push edi

    mov esi, [src] ; source array
    mov edi, [dst] ; destination array
    mov ecx, [len]

    rep movsb

    pop edi
    pop esi
    ret



    这里我为了代码可移植性,使用的是NASM格式的汇编代码。NASM是一个非常出色的开源汇编编译器,支持各种平台和中间格式,被开源项目广泛使用,这样可 以避免同时使用 VC 的嵌入式汇编和 GCC 中麻烦的 unix 风格 AT&T 格式汇编 :P

    代码初始的cpu p4定义使用p4指令集,因为后面的很多优化工作使用了P4指令集和相关特性;接着的segment .text use32定义此代码在32位代码段;然后global定义标签_fast_memcpy1为全局符号,使得C++代码中可以LINK其.obj后访问此 代码;最后%define定义多个宏,用于访问函数参数。

    在C++中只需要定义fast_memcpy1函数格式并链接nasm编译生成 的.obj文件即可。NASM编译时 -f 参数指定生成中间文件格式为 MS 的 32 位 COFF 格式,-o 参数指定输出文件名。

    上 面这段代码非常简单,适合小内存块的快速拷贝。实际上VC编译器在处理小内存拷贝时,会自动根据情况使用 rep movsb 直接替换 memcpy 函数,通过忽略函数调用和堆栈操作,优化代码长度和性能。

    不过在 32 位的 x86 架构下,完全没有必要逐字节进行操作,使用 movsd 替换 movsb 是必然的选择。

    以下为引用:

    global _fast_memcpy2

    %define param esp+8+4
    %define src param+0
    %define dst param+4
    %define len param+8

    _fast_memcpy2:
    push esi
    push edi

    mov esi, [src] ; source array
    mov edi, [dst] ; destination array
    mov ecx, [len]
    shr ecx, 2 ; convert to DWORD count

    rep movsd

    pop edi
    pop esi
    ret



    为了展示方便,这里假设源和目标内存块本身长度都是64字节的整数倍,并且已经4K页对齐。前者保证单条指令不会出现跨CACHE 行访问的情况;后者保证测试速度时不会因为跨页操作影响测试结果。等会分析CACHE时再详细解释为什么要做这种假设。

    不过因为现代 CPU大多使用了很长的指令流水线,多条指令并行工作往往比一条指令效率更高,因此 AMD 文档中给出了这样的优化:

    以下为引用:

    global _fast_memcpy3

    %define param esp+8+4
    %define src param+0
    %define dst param+4
    %define len param+8

    _fast_memcpy3:
    push esi
    push edi

    mov esi, [src] ; source array
    mov edi, [dst] ; destination array
    mov ecx, [len]
    shr ecx, 2 ; convert to DWORD count

    .copyloop:
    mov eax, dword [esi]
    mov dword [edi], eax

    add esi, 4
    add edi, 4

    dec ecx
    jnz .copyloop

    pop edi
    pop esi
    ret



    标签.copyloop中那段循环实际上完成跟rep movsd指令完全相同的工作,但是因为是多条指令,理论上CPU指令流水线可以并行处理之。故而在AMD的文档中指出能有1.5%的性能提高,不过就我 实测效果不太明显。相对而言,当年从486向pentium架构迁移时,这两种方式的区别非常明显。记得Delphi 3还是4中就只是通过做这一种优化,其字符串处理性能就有较大提升。而目前主流CPU厂商,实际上都是通过微代码技术,内核中使用RISC微指令模拟 CISC指令集,因此现在效果并不明显。

    然后,可以通过循环展开的优化策略,增加每次处理数据量并减少循环次数,达到性能提升目的。

    以下为引用:

    global _fast_memcpy4

    %define param esp+8+4
    %define src param+0
    %define dst param+4
    %define len param+8

    _fast_memcpy4:
    push esi
    push edi

    mov esi, [src] ; source array
    mov edi, [dst] ; destination array
    mov ecx, [len]
    shr ecx, 4 ; convert to 16-byte size count

    .copyloop:
    mov eax, dword [esi]
    mov dword [edi], eax

    mov ebx, dword [esi+4]
    mov dword [edi+4], ebx

    mov eax, dword [esi+8]
    mov dword [edi+8], eax

    mov ebx, dword [esi+12]
    mov dword [edi+12], ebx

    add esi, 16
    add edi, 16

    dec ecx
    jnz .copyloop

    pop edi
    pop esi
    ret



    但这种操作就 AMD 文档上评测反而有 %1.5 性能降低,呵呵。其自己的说法是需要将读取内存和写入内存的操作分组,以使CPU可以一次性搞定。改称以下分组操作就可以比_fast_memcpy3提 高3% -_-b

    以下为引用:

    global _fast_memcpy5

    %define param esp+8+4
    %define src param+0
    %define dst param+4
    %define len param+8

    _fast_memcpy5:
    push esi
    push edi

    mov esi, [src] ; source array
    mov edi, [dst] ; destination array
    mov ecx, [len]
    shr ecx, 4 ; convert to 16-byte size count

    .copyloop:
    mov eax, dword [esi]
    mov ebx, dword [esi+4]
    mov dword [edi], eax
    mov dword [edi+4], ebx

    mov eax, dword [esi+8]
    mov ebx, dword [esi+12]
    mov dword [edi+8], eax
    mov dword [edi+12], ebx

    add esi, 16
    add edi, 16

    dec ecx
    jnz .copyloop

    pop edi
    pop esi
    ret



    可惜我在P4上实在测不出什么区别,呵呵,大概P4和AMD实现流水线的思路有细微的出入吧 :D

    既然进行循环展 开,为什么不干脆多展开一些呢?虽然x86下面通用寄存器只有那么几个,但是现在有MMX啊,呵呵,大把的寄存器啊 :D 改称使用MMX寄存器后,一次载入/写入操作可以处理64字节的数据,呵呵,比_fast_memcpy5可以再有7%的性能提升。

    以下为引用:

    global _fast_memcpy6

    %define param esp+8+4
    %define src param+0
    %define dst param+4
    %define len param+8

    _fast_memcpy6:
    push esi
    push edi

    mov esi, [src] ; source array
    mov edi, [dst] ; destination array
    mov ecx, [len] ; number of QWORDS (8 bytes) assumes len / CACHEBLOCK is an integer
    shr ecx, 3

    lea esi, [esi+ecx*8] ; end of source
    lea edi, [edi+ecx*8] ; end of destination
    neg ecx ; use a negative offset as a combo pointer-and-loop-counter

    .copyloop:
    movq mm0, qword [esi+ecx*8]
    movq mm1, qword [esi+ecx*8+8]
    movq mm2, qword [esi+ecx*8+16]
    movq mm3, qword [esi+ecx*8+24]
    movq mm4, qword [esi+ecx*8+32]
    movq mm5, qword [esi+ecx*8+40]
    movq mm6, qword [esi+ecx*8+48]
    movq mm7, qword [esi+ecx*8+56]

    movq qword [edi+ecx*8], mm0
    movq qword [edi+ecx*8+8], mm1
    movq qword [edi+ecx*8+16], mm2
    movq qword [edi+ecx*8+24], mm3
    movq qword [edi+ecx*8+32], mm4
    movq qword [edi+ecx*8+40], mm5
    movq qword [edi+ecx*8+48], mm6
    movq qword [edi+ecx*8+56], mm7

    add ecx, 8
    jnz .copyloop

    emms

    pop edi
    pop esi

    ret



    优化到这个份上,常规的优化手段基本上已经用尽,需要动用非常手段了,呵呵。
    让我们回过头来看看P4架构下的Cache结 构。

    The IA-32 Intel Architecture Software Developer's Manual, Volume 3: System Programming Guide

    Intel的系统变成手册中第十章介绍 了IA32架构下的内存缓存控制。因为CPU速度和内存速度的巨大差距,CPU厂商通过在CPU中内置和外置多级缓存提高频繁使用数据的访问速度。一般来 说,在CPU和内存之间存在L1, L2和L3三级缓存(还有几种TLB缓存在此不涉及),每级缓存的速度有一个数量级左右的差别,容量也有较大差别(实际上跟$有关,呵呵),而L1缓存更 是细分为指令缓存和数据缓存,用于不同的目的。就P4和Xeon的处理器来说,L1指令缓存由Trace Cache取代,内置在NetBust微架构中;L1数据缓存和L2缓存则封装在CPU中,根据CPU档次不同,分别在8-16K和256-512K之 间;而L3缓存只在Xeon处理器中实现,也是封装在CPU中,512K-1M左右。
    可以通过查看CPU信息的软件如CPUInfo查看当前机器的缓存信息,如我的系统为:
    P4 1.7G, 8K L1 Code Cache, 12K L1 Data Cache, 256K L2 Cache。

    而缓存在实现上是若干行(slot or line)组成的,每行对应内存中的一个地址上的连续数据,由高速缓存管理器控制读写中的数据载入和命中。其原理这里不多罗嗦,有兴趣的朋友可以自行查看 Intel手册。需要知道的就是每个slot的长度在P4以前是32字节,P4开始改成64字节。而对缓存行的操作都是完整进行的,哪怕只读一个字节也需 要将整个缓存行(64字节)全部载入,后面的优化很大程度上基于这些原理。

    就缓存的工作模式来说,P4支持的有六种之多,这里就不一一介 绍了。对我们优化有影响的,实际上就是写内存时缓存的表现。最常见的WT(Write-through)写通模式在写数据到内存的同时更新数据到缓存中; 而WB(Write-back)写回模式,则直接写到缓存中,暂不进行较慢的内存读写。这两种模式在操作频繁操作(每秒百万次这个级别)的内存变量处理上 有较大性能差别。例如通过编写驱动模块操作MTRR强行打开WB模式,在Linux的网卡驱动中曾收到不错的效果,但对内存复制的优化帮助不大,因为我们 需要的是完全跳过对缓存的操作,无论是缓存定位、载入还是写入。

    好在P4提供了MOVNTQ指令,使用WC(Write- combining)模式,跳过缓存直接写内存。因为我们的写内存操作是纯粹的写,写入的数据一定时间内根本不会被使用,无论使用WT还是WB模式,都会 有冗余的缓存操作。优化代码如下:

    以下为引用:

    global _fast_memcpy7

    %define param esp+8+4
    %define src param+0
    %define dst param+4
    %define len param+8

    _fast_memcpy7:
    push esi
    push edi

    mov esi, [src] ; source array
    mov edi, [dst] ; destination array
    mov ecx, [len] ; number of QWORDS (8 bytes) assumes len / CACHEBLOCK is an integer
    shr ecx, 3

    lea esi, [esi+ecx*8] ; end of source
    lea edi, [edi+ecx*8] ; end of destination
    neg ecx ; use a negative offset as a combo pointer-and-loop-counter

    .copyloop:
    movq mm0, qword [esi+ecx*8]
    movq mm1, qword [esi+ecx*8+8]
    movq mm2, qword [esi+ecx*8+16]
    movq mm3, qword [esi+ecx*8+24]
    movq mm4, qword [esi+ecx*8+32]
    movq mm5, qword [esi+ecx*8+40]
    movq mm6, qword [esi+ecx*8+48]
    movq mm7, qword [esi+ecx*8+56]

    movntq qword [edi+ecx*8], mm0
    movntq qword [edi+ecx*8+8], mm1
    movntq qword [edi+ecx*8+16], mm2
    movntq qword [edi+ecx*8+24], mm3
    movntq qword [edi+ecx*8+32], mm4
    movntq qword [edi+ecx*8+40], mm5
    movntq qword [edi+ecx*8+48], mm6
    movntq qword [edi+ecx*8+56], mm7

    add ecx, 8
    jnz .copyloop

    sfence ; flush write buffer
    emms

    pop edi
    pop esi

    ret



    写内存的movq指令全部改为movntq指令,并在复制操作完成后,调用sfence刷新写缓存,因为缓存中内容可能已经失效 了。这样一来在写内存外的载入缓存操作,以及缓存本身的操作都被省去,大大减少了冗余内存操作。按AMD的说法能有60%的性能提升,我实测也有50%左 右明显的性能提升。
    movntq和sfence等指令可以参考Intel的指令手册:

    The IA-32 Intel Architecture Software Developer's Manual, Volume 2A: Instruction Set Reference, A-M
    The IA-32 Intel Architecture Software Developer's Manual, Volume 2B: Instruction Set Reference, N-Z


    在优化完写内存 后,同样可以通过对读内存的操作进行优化提升性能。虽然CPU在读取数据时,会有一个自动的预读优化,但在操作连续内存区域时显式要求CPU预读数据,还 是可以明显地优化性能。

    以下为引用:

    global _fast_memcpy8

    %define param esp+8+4
    %define src param+0
    %define dst param+4
    %define len param+8

    _fast_memcpy8:
    push esi
    push edi

    mov esi, [src] ; source array
    mov edi, [dst] ; destination array
    mov ecx, [len] ; number of QWORDS (8 bytes) assumes len / CACHEBLOCK is an integer
    shr ecx, 3

    lea esi, [esi+ecx*8] ; end of source
    lea edi, [edi+ecx*8] ; end of destination
    neg ecx ; use a negative offset as a combo pointer-and-loop-counter

    .writeloop:
    prefetchnta [esi+ecx*8 + 512] ; fetch ahead by 512 bytes

    movq mm0, qword [esi+ecx*8]
    movq mm1, qword [esi+ecx*8+8]
    movq mm2, qword [esi+ecx*8+16]
    movq mm3, qword [esi+ecx*8+24]
    movq mm4, qword [esi+ecx*8+32]
    movq mm5, qword [esi+ecx*8+40]
    movq mm6, qword [esi+ecx*8+48]
    movq mm7, qword [esi+ecx*8+56]

    movntq qword [edi+ecx*8], mm0
    movntq qword [edi+ecx*8+8], mm1
    movntq qword [edi+ecx*8+16], mm2
    movntq qword [edi+ecx*8+24], mm3
    movntq qword [edi+ecx*8+32], mm4
    movntq qword [edi+ecx*8+40], mm5
    movntq qword [edi+ecx*8+48], mm6
    movntq qword [edi+ecx*8+56], mm7

    add ecx, 8
    jnz .writeloop

    sfence ; flush write buffer
    emms

    pop edi
    pop esi

    ret



    增加一个简单的prefetchnta指令,提示CPU在处理当前读取内存操作的同时,预读前面512字节处的一个缓存行64字节 内容。这样一来又可以有10%左右的性能提升。
    最后,对正在处理的内存,可以通过显式的内存读取操作,强制性要求其载入到缓存中,因为 prefetchnta指令还只是一个提示,可以被CPU忽略。这样可以再次获得60%左右的性能提示,我实测没有这么高,但是也比较明显。

    以下为引用:

    global _fast_memcpy9

    %define param esp+12+4
    %define src param+0
    %define dst param+4
    %define len param+8

    %define CACHEBLOCK 400h

    _fast_memcpy9:
    push esi
    push edi
    push ebx

    mov esi, [src] ; source array
    mov edi, [dst] ; destination array
    mov ecx, [len] ; number of QWORDS (8 bytes) assumes len / CACHEBLOCK is an integer
    shr ecx, 3

    lea esi, [esi+ecx*8] ; end of source
    lea edi, [edi+ecx*8] ; end of destination
    neg ecx ; use a negative offset as a combo pointer-and-loop-counter

    .mainloop:
    mov eax, CACHEBLOCK / 16 ; note: .prefetchloop is unrolled 2X
    add ecx, CACHEBLOCK ; move up to end of block

    .prefetchloop:
    mov ebx, [esi+ecx*8-64] ; read one address in this cache line...
    mov ebx, [esi+ecx*8-128] ; ... and one in the previous line
    sub ecx, 16 ; 16 QWORDS = 2 64-byte cache lines
    dec eax
    jnz .prefetchloop

    mov eax, CACHEBLOCK / 8

    .writeloop:
    prefetchnta [esi+ecx*8 + 512] ; fetch ahead by 512 bytes

    movq mm0, qword [esi+ecx*8]
    movq mm1, qword [esi+ecx*8+8]
    movq mm2, qword [esi+ecx*8+16]
    movq mm3, qword [esi+ecx*8+24]
    movq mm4, qword [esi+ecx*8+32]
    movq mm5, qword [esi+ecx*8+40]
    movq mm6, qword [esi+ecx*8+48]
    movq mm7, qword [esi+ecx*8+56]

    movntq qword [edi+ecx*8], mm0
    movntq qword [edi+ecx*8+8], mm1
    movntq qword [edi+ecx*8+16], mm2
    movntq qword [edi+ecx*8+24], mm3
    movntq qword [edi+ecx*8+32], mm4
    movntq qword [edi+ecx*8+40], mm5
    movntq qword [edi+ecx*8+48], mm6
    movntq qword [edi+ecx*8+56], mm7

    add ecx, 8
    dec eax
    jnz .writeloop

    or ecx, ecx ; assumes integer number of cacheblocks
    jnz .mainloop

    sfence ; flush write buffer
    emms

    pop ebx
    pop edi
    pop esi

    ret



    至此,一个完整的内存复制函数的优化流程就结束了,通过对缓存的了解和使用,一次又一次地超越自己,最终获得一个较为令人满意地结 果。(号称300%性能提示,实测175%-200%,也算相当不错了)

    在编写测试代码的时候需要注意两点:

    一是计时精 度的问题,需要使用高精度的物理计数器,避免误差。推荐使用rdtsc指令,然后根据CPU主频计算时间。CPU主频可以通过高精度计时器动态计算,我这 儿偷懒直接从注册表里面读取了 :P
    代码如下:

    以下为引用:

    #ifdef WIN32
    typedef __int64 uint64_t;
    #else
    #include <stdint.h>
    #endif

    bool GetPentiumClockEstimateFromRegistry(uint64_t& frequency)
    {
    HKEY hKey;

    frequency = 0;

    LONG rc = ::RegOpenKeyEx(HKEY_LOCAL_MACHINE, "Hardware\\Description\\System\\CentralProcessor\\0", 0, KEY_READ, &hKey);

    if(rc == ERROR_SUCCESS)
    {
    DWORD cbBuffer = sizeof (DWORD);
    DWORD freq_mhz;

    rc = ::RegQueryValueEx(hKey, "~MHz", NULL, NULL, (LPBYTE)(&freq_mhz), &cbBuffer);

    if (rc == ERROR_SUCCESS)
    frequency = freq_mhz * MEGA;

    RegCloseKey (hKey);
    }

    return frequency > 0;
    }

    void getTimeStamp(uint64_t& timeStamp)
    {
    #ifdef WIN32
    __asm
    {
    push edx
    push ecx
    mov ecx, timeStamp
    //_emit 0Fh // RDTSC
    //_emit 31h
    rdtsc
    mov [ecx], eax
    mov [ecx+4], edx
    pop ecx
    pop edx
    }
    #else
    __asm__ __volatile__ ("rdtsc" : "=A" (timeStamp));
    #endif
    }




    二是测试内存复制的缓冲区的大小,如果缓冲区过小,第一次拷贝两个缓冲区时就会导致所有数据都被载入L2缓存中,得出比普 通内存操作高一个数量级的数值。例如我的L2缓冲为256K,如果我用两个128K的缓冲区对着拷贝,无论循环多少次,速度都在普通内存复制的10倍左 右。因此设置一个较大的值是必要的。

  • 相关阅读:
    新闻发布项目——接口类(newsTbDao)
    新闻发布项目——接口类(commentDao)
    新闻发布项目——接口类(newsTbDao)
    新闻发布项目——接口类(UserDao)
    新闻发布项目——接口类(UserDao)
    新闻发布项目——接口类(UserDao)
    新闻发布项目——实体类(categoryTB)
    新闻发布项目——实体类(categoryTB)
    新闻发布项目——实体类(categoryTB)
    新闻发布项目——实体类(comment)
  • 原文地址:https://www.cnblogs.com/nsnow/p/1724629.html
Copyright © 2011-2022 走看看