DSP

从别人的视频优化中看dsp算法优化

2019-07-13 18:53发布

根据C6000的软件编程流程,对Viterbi算法的编程及其优化可分为三个阶段来进行。这三个阶段分别为:开发C代码、优化C代码、编写线性汇编代码。在代码编写和优化过程中,这三个阶段不是必须都要经过的,只要在某一阶段已经满足了算法代码的功能和性能要求,就不必继续进行下面的阶段。 ①开发C代码。这一阶段完全是根据任务要求来完成算法的代码编写工作。在C6000的集成开发环境CCS(Code Composer Studio)下进行代码的编译和功能验证,然后可用CCS的调试工具(如Profiler),利用在程序中设置断点的方法可找出程序中耗时最多、最影响整体性能的代码段。为改进代码性能,可进入下一阶段。如下是针对(2,1,3)码的Viterbi算法代码中完成算法功能的核心循环,也是最耗时、最影响代码整体性能的低效率段。 for(c=0;caccum_err_metric[j][0]+branch_metric] {accum_err_metric[nextstate[j][i][1]=accum_err_metric[j][0]+branch_metri; state_history[nextstate[j][i]][sh_ptr]=j; } }*/end of i<2*/ }/*end of j=0;i--) hamm+=(target_vector>>i)&0x01; return hamm; } 在验证了算法代码实现功能并以设置断点的方法测试代码的性能,这段循环运行耗时(时钟周期)为1790。显然,性能不能达到要求,就要进入代码优化的第二阶段了。 ②一般在代码调试中,最影响性能的是其中的循环代码段。而软件流水是一种用于安排循环内的指令运行方式,尽可能充分利用CPU的功能单元等资源,使循环的多次迭代能够并行执行的一种技术。在C6000的C/C++编译器里,采用软件流水使编译出来的程序代码优化是一项核心技术。所以在进一步优化之前,需要调整并尽可能简化代码的结构并去除影响软件流水的因素使其能够被编译器充分流水,这对大幅提高整个代码的性能非常重要。 所以,在考虑影响因素同时对Viterbi算法的循环代码进行如下调整; *使用内联函数(intrinsics)替代复杂的C语言程序。C6000编译器提供了许多intrinsics,可以快速优化C代码。Intrinsics是直接参与C6000汇编指令映射的内联函数。在这里使用了_extu(x,y,z),以简化其中hamm代码部分。 *尽管软件流水循环可包含intrinsics,但不能包含函数调用。所以需要把调用函数hamm在循环中展开实现。 *由于编译器仅对最内部的循环执行流水,所以为了提高性能应尽可能创造一比较大的内循环。在代码中可以看到,在最内循环是i的两次循环,仅对它进行流水,对整个代码的性能提高不大。所以一个想法是,将i和j循环全部展开,使编译器直接面对最大的C循环以最大发挥软件流水的作用。 *另外,展开循环后代码中的变量如果可以确定其运行中的值,就尽量以实值代入,这样减少了变量个数,也就是减少了所需分配的寄存器个数(C62xxCPU中有32个寄存器)。 在进行上述调整后运行代码,进行测试发展,性能没有太大改善;用编译器反馈表(feedback)进行观察发现,循环并没有发生流水。这是为什么呢?原来在展开内部循环后导致C循环内代码尺寸太大,需要的寄存器数目大于C62XX的32个寄存器,所以不能进行软件流水。为了解决这问题,需要简化循环或将循环拆成几个小循环。在这里先将C循环内部的小循环展开,然后将其拆成分别完成度量计算和累计度量比较的两个循环,这样就减小了每个循环中的代码尺寸。限于篇幅这里只写出累计度量比较的循环代码。 /*完成累计度量比较的循环*/ accum00=accum_err_metric[0];accum10=accum_err_metric[1]; accum20=accum_err_metric[2];accum30=accum_err_metric[3]; for(c=0;cadd2){accum00=add2;state_history[0][sh_ptr]=0;} else{accum00=add1;state_history[0][sh_ptr]=1;} if(add3>add4){accum20=add4;state_history[2][sh_ptr]=0;} else{accum20=add3;state_history[2][sh_ptr]=1;} if(add5>add6){accum30=add6;state_history[3][sh_ptr]=2;} else{accum30=add5;state_history[3][sh_ptr]=3;} if(add7>add8){accum10=add8;state_history[1][sh_ptr]=2;} else{accum10=add7;state_history[1][sh_ptr]=3;} } accum_err_metric[0]=accum00;accum_err_metic[1]=accum10; accum_err_metric[2]=accum20;accum_err_metric[3]=accum30; 其中accum_err_metric[i]为状态i的累计度量值,branch_metric_array[][]为计算得到的各时刻量值,原来代码中的二维数码mextstate[j][i]被以实值代入。另外在编程考虑时要注意一点:程序中对数据的取命令(load)是非常耗时的,所以应考虑尽量减少对数据数组的操作。在上面程序的改进中,先从数组中取出要进行循环处理的累计度量值,再使用accumXX及addX作为各次迭代的中间变量,在循环后将最后的结果放入数据。这样就大大减少了对数组的操作,从而使优化进一步提高。 *编译器优化选项的选择。C6000 C/C++编译器提供了大量的编译选项,供用户在编译时选择使用。这些选项中的部分会直接影响或控制编译器优化过程,因而会影响编译输出的代码优化性能。选择适合的选项,能极大地提高优化性能。在这里使用的优化选项有: -03——表示可得到最高程度的优化,编译器将执行各种优化循环的方法,如软件流水、循环展开等等。 -pm——在使用-o3选项进行优化时尽量联合使用-pm选项,-pm是程序级优化,使优化器访问整个程序,了解循环次数。 -op1——使用了外部变量,但未使用外部函数调用。 -g——使能符号调试和汇编源语句调试。 另外,还有不少考虑因素和优化调试方法,如消除存储器相关性、对短字长的数据使用宽离长度的存储器访问等。由于篇幅所限不能在这里一一列出,详细资料可参考TMS320C6000 Code Composer Studio Manuals中的TMS320C6000 Optimizing C Compiler User's Guide。 测试结果:在经过上述优化后运行耗时(时钟周期)已降为406个,代码的性能大为提高,已经满足系统要求。 ③由上述可知,在程序中影响性能的主要代码通常是循环。优化一个循环较好的方法是抽出这个循环,使之成为一个单独文件,对其进行重新编写、重新编译和单独运行。为了提高代码性能,对影响速度的关键C代码段可以用线性汇编重新编写,使用汇编优化器进行优化后效率是非常高的。若代码性能仍未满足要求,则可进行第三阶段,将其抽出,全部用线性汇编来编写,在代码中以函数的形式将改写的部分调用。将循环代码段改写为线性汇编调用函数的格式如下: .global_KernelLoop ;函数名定义前加_ _KernelLoop:.cproc channel_data,branch_metric_array,depunc ;定义入口形参变量 .reg c,q0,q1,y1,y2,x1,x2,cc,temp,temp0,temp1,temp2,temp3 .reg counter,valuel,value2,value3定义中间变量 no_mdep ;表明存储器地址不相关 zer0 c ;初始化变量 zero cc loop: .trip 24 ;声明循环24次 · · · 运行语句 · · · [counter]sub counter,23,counter;循环计数 [counter]b loop ;循环跳转 .return ;完成返回 .endproc ;结束 编写线性汇编的工作量大,开发周期长且不能像C语言程序一样移植到其它类型DSP上,所以尽量在第一、二阶段完成工作。若仍满足不了性能要求,则再对关键代码段进行线性汇编的改写。     这个算法的优化过程分为3个部分。 1.使用内联函数,提高函数的读取,一般的作用不是很大。 2.使用流水作业。将大循环拆分成多个小循环,这样实现流水作业,大大提高了效率。(问题的关键是怎样实现流水作业。) 3.关键的位置使用线性汇编。