lab10.c

北京理工大学DSP C6000课程的实验程序

 /* lab3字长优化改进程序 */

  short DP(short *m, short *n, short count) ;
  
  short a[30] = { 1, 3, 5, 7, 9, 11, 13, 15, 17, 19,
                             21, 23, 25, 27, 29, 31, 33, 35, 37, 39, 
                             41, 43, 45, 47, 49, 51, 53, 55, 57, 59 } ;

  short x[30] = {  60, 58, 56, 54, 52, 50, 48, 46, 44, 42,
                              40, 38, 36,34, 32, 30, 28, 26, 24, 22,
                              20, 18, 16, 14, 12, 10, 8, 6, 4, 2 } ;

  short y = 0 ;

  main()
    { 
       y = DP(a, x, 30) ;
    }
 
  short DP(short *m, short *n, short count)
    {
       short i ;
	int pro_h, pro_l ;
       int sum_h = 0 ;
       int sum_l = 0 ;
       int sum = 0 ;

       int *data_a = ( int *) m ;
       int *data_x = ( int *) n ;
       count = count >> 1 ;
       for(i=0;i<count;i++)
         {
            pro_l = _mpy( data_a[i], data_x[i]) ;
            pro_h = _mpyh( data_a[i], data_x[i]) ;
            sum_l += pro_l ; 
            sum_h += pro_h ;
         }
       sum = sum_l + sum_h ;
       return( sum) ;
    }