dsputil_ppc.c
上传用户:wstnjxml
上传日期:2014-04-03
资源大小:7248k
文件大小:10k
源码类别:

Windows CE

开发平台:

C/C++

  1. /*
  2.  * Copyright (c) 2002 Brian Foley
  3.  * Copyright (c) 2002 Dieter Shirley
  4.  * Copyright (c) 2003-2004 Romain Dolbeau <romain@dolbeau.org>
  5.  *
  6.  * This library is free software; you can redistribute it and/or
  7.  * modify it under the terms of the GNU Lesser General Public
  8.  * License as published by the Free Software Foundation; either
  9.  * version 2 of the License, or (at your option) any later version.
  10.  *
  11.  * This library is distributed in the hope that it will be useful,
  12.  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  13.  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14.  * Lesser General Public License for more details.
  15.  *
  16.  * You should have received a copy of the GNU Lesser General Public
  17.  * License along with this library; if not, write to the Free Software
  18.  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
  19.  */
  20. #include "../dsputil.h"
  21. #include "dsputil_ppc.h"
  22. #ifdef HAVE_ALTIVEC
  23. #include "dsputil_altivec.h"
  24. #endif
  25. extern void fdct_altivec(int16_t *block);
  26. extern void idct_put_altivec(uint8_t *dest, int line_size, int16_t *block);
  27. extern void idct_add_altivec(uint8_t *dest, int line_size, int16_t *block);
  28. int mm_flags = 0;
  29. int mm_support(void)
  30. {
  31.     int result = 0;
  32. #ifdef HAVE_ALTIVEC
  33.     if (has_altivec()) {
  34.         result |= MM_ALTIVEC;
  35.     }
  36. #endif /* result */
  37.     return result;
  38. }
  39. #ifdef POWERPC_PERFORMANCE_REPORT
  40. unsigned long long perfdata[POWERPC_NUM_PMC_ENABLED][powerpc_perf_total][powerpc_data_total];
  41. /* list below must match enum in dsputil_ppc.h */
  42. static unsigned char* perfname[] = {
  43.   "ff_fft_calc_altivec",
  44.   "gmc1_altivec",
  45.   "dct_unquantize_h263_altivec",
  46.   "fdct_altivec",
  47.   "idct_add_altivec",
  48.   "idct_put_altivec",
  49.   "put_pixels16_altivec",
  50.   "avg_pixels16_altivec",
  51.   "avg_pixels8_altivec",
  52.   "put_pixels8_xy2_altivec",
  53.   "put_no_rnd_pixels8_xy2_altivec",
  54.   "put_pixels16_xy2_altivec",
  55.   "put_no_rnd_pixels16_xy2_altivec",
  56.   "hadamard8_diff8x8_altivec",
  57.   "hadamard8_diff16_altivec",
  58.   "avg_pixels8_xy2_altivec",
  59.   "clear_blocks_dcbz32_ppc",
  60.   "clear_blocks_dcbz128_ppc",
  61.   "put_h264_chroma_mc8_altivec",
  62.   "avg_h264_chroma_mc8_altivec",
  63.   "put_h264_qpel16_h_lowpass_altivec",
  64.   "avg_h264_qpel16_h_lowpass_altivec",
  65.   "put_h264_qpel16_v_lowpass_altivec",
  66.   "avg_h264_qpel16_v_lowpass_altivec",
  67.   "put_h264_qpel16_hv_lowpass_altivec",
  68.   "avg_h264_qpel16_hv_lowpass_altivec",
  69.   ""
  70. };
  71. #include <stdio.h>
  72. #endif
  73. #ifdef POWERPC_PERFORMANCE_REPORT
  74. void powerpc_display_perf_report(void)
  75. {
  76.   int i, j;
  77.   av_log(NULL, AV_LOG_INFO, "PowerPC performance reportn Values are from the PMC registers, and represent whatever the registers are set to record.n");
  78.   for(i = 0 ; i < powerpc_perf_total ; i++)
  79.   {
  80.     for (j = 0; j < POWERPC_NUM_PMC_ENABLED ; j++)
  81.       {
  82. if (perfdata[j][i][powerpc_data_num] != (unsigned long long)0)
  83.   av_log(NULL, AV_LOG_INFO,
  84.   " Function "%s" (pmc%d):ntmin: %lluntmax: %lluntavg: %1.2lf (%llu)n",
  85.   perfname[i],
  86.   j+1,
  87.   perfdata[j][i][powerpc_data_min],
  88.   perfdata[j][i][powerpc_data_max],
  89.   (double)perfdata[j][i][powerpc_data_sum] /
  90.   (double)perfdata[j][i][powerpc_data_num],
  91.   perfdata[j][i][powerpc_data_num]);
  92.       }
  93.   }
  94. }
  95. #endif /* POWERPC_PERFORMANCE_REPORT */
  96. /* ***** WARNING ***** WARNING ***** WARNING ***** */
  97. /*
  98.   clear_blocks_dcbz32_ppc will not work properly
  99.   on PowerPC processors with a cache line size
  100.   not equal to 32 bytes.
  101.   Fortunately all processor used by Apple up to
  102.   at least the 7450 (aka second generation G4)
  103.   use 32 bytes cache line.
  104.   This is due to the use of the 'dcbz' instruction.
  105.   It simply clear to zero a single cache line,
  106.   so you need to know the cache line size to use it !
  107.   It's absurd, but it's fast...
  108.   update 24/06/2003 : Apple released yesterday the G5,
  109.   with a PPC970. cache line size : 128 bytes. Oups.
  110.   The semantic of dcbz was changed, it always clear
  111.   32 bytes. so the function below will work, but will
  112.   be slow. So I fixed check_dcbz_effect to use dcbzl,
  113.   which is defined to clear a cache line (as dcbz before).
  114.   So we still can distinguish, and use dcbz (32 bytes)
  115.   or dcbzl (one cache line) as required.
  116.   see <http://developer.apple.com/technotes/tn/tn2087.html>
  117.   and <http://developer.apple.com/technotes/tn/tn2086.html>
  118. */
  119. void clear_blocks_dcbz32_ppc(DCTELEM *blocks)
  120. {
  121. POWERPC_PERF_DECLARE(powerpc_clear_blocks_dcbz32, 1);
  122.     register int misal = ((unsigned long)blocks & 0x00000010);
  123.     register int i = 0;
  124. POWERPC_PERF_START_COUNT(powerpc_clear_blocks_dcbz32, 1);
  125. #if 1
  126.     if (misal) {
  127.       ((unsigned long*)blocks)[0] = 0L;
  128.       ((unsigned long*)blocks)[1] = 0L;
  129.       ((unsigned long*)blocks)[2] = 0L;
  130.       ((unsigned long*)blocks)[3] = 0L;
  131.       i += 16;
  132.     }
  133.     for ( ; i < sizeof(DCTELEM)*6*64-31 ; i += 32) {
  134. #ifndef __MWERKS__
  135.       asm volatile("dcbz %0,%1" : : "b" (blocks), "r" (i) : "memory");
  136. #else
  137.       __dcbz( blocks, i );
  138. #endif
  139.     }
  140.     if (misal) {
  141.       ((unsigned long*)blocks)[188] = 0L;
  142.       ((unsigned long*)blocks)[189] = 0L;
  143.       ((unsigned long*)blocks)[190] = 0L;
  144.       ((unsigned long*)blocks)[191] = 0L;
  145.       i += 16;
  146.     }
  147. #else
  148.     memset(blocks, 0, sizeof(DCTELEM)*6*64);
  149. #endif
  150. POWERPC_PERF_STOP_COUNT(powerpc_clear_blocks_dcbz32, 1);
  151. }
  152. /* same as above, when dcbzl clear a whole 128B cache line
  153.    i.e. the PPC970 aka G5 */
  154. #ifndef NO_DCBZL
  155. void clear_blocks_dcbz128_ppc(DCTELEM *blocks)
  156. {
  157. POWERPC_PERF_DECLARE(powerpc_clear_blocks_dcbz128, 1);
  158.     register int misal = ((unsigned long)blocks & 0x0000007f);
  159.     register int i = 0;
  160. POWERPC_PERF_START_COUNT(powerpc_clear_blocks_dcbz128, 1);
  161. #if 1
  162.  if (misal) {
  163.    // we could probably also optimize this case,
  164.    // but there's not much point as the machines
  165.    // aren't available yet (2003-06-26)
  166.       memset(blocks, 0, sizeof(DCTELEM)*6*64);
  167.     }
  168.     else
  169.       for ( ; i < sizeof(DCTELEM)*6*64 ; i += 128) {
  170. asm volatile("dcbzl %0,%1" : : "b" (blocks), "r" (i) : "memory");
  171.       }
  172. #else
  173.     memset(blocks, 0, sizeof(DCTELEM)*6*64);
  174. #endif
  175. POWERPC_PERF_STOP_COUNT(powerpc_clear_blocks_dcbz128, 1);
  176. }
  177. #else
  178. void clear_blocks_dcbz128_ppc(DCTELEM *blocks)
  179. {
  180.   memset(blocks, 0, sizeof(DCTELEM)*6*64);
  181. }
  182. #endif
  183. #ifndef NO_DCBZL
  184. /* check dcbz report how many bytes are set to 0 by dcbz */
  185. /* update 24/06/2003 : replace dcbz by dcbzl to get
  186.    the intended effect (Apple "fixed" dcbz)
  187.    unfortunately this cannot be used unless the assembler
  188.    knows about dcbzl ... */
  189. long check_dcbzl_effect(void)
  190. {
  191.   register char *fakedata = (char*)av_malloc(1024);
  192.   register char *fakedata_middle;
  193.   register long zero = 0;
  194.   register long i = 0;
  195.   long count = 0;
  196.   if (!fakedata)
  197.   {
  198.     return 0L;
  199.   }
  200.   fakedata_middle = (fakedata + 512);
  201.   memset(fakedata, 0xFF, 1024);
  202.   /* below the constraint "b" seems to mean "Address base register"
  203.      in gcc-3.3 / RS/6000 speaks. seems to avoid using r0, so.... */
  204.   asm volatile("dcbzl %0, %1" : : "b" (fakedata_middle), "r" (zero));
  205.   for (i = 0; i < 1024 ; i ++)
  206.   {
  207.     if (fakedata[i] == (char)0)
  208.       count++;
  209.   }
  210.   av_free(fakedata);
  211.   
  212.   return count;
  213. }
  214. #else
  215. long check_dcbzl_effect(void)
  216. {
  217.   return 0;
  218. }
  219. #endif
  220. void dsputil_h264_init_ppc(DSPContext* c, AVCodecContext *avctx);
  221. void dsputil_init_ppc(DSPContext* c, AVCodecContext *avctx)
  222. {
  223.     // Common optimizations whether Altivec is available or not
  224.   switch (check_dcbzl_effect()) {
  225.   case 32:
  226.     c->clear_blocks = clear_blocks_dcbz32_ppc;
  227.     break;
  228.   case 128:
  229.     c->clear_blocks = clear_blocks_dcbz128_ppc;
  230.     break;
  231.   default:
  232.     break;
  233.   }
  234. #ifdef HAVE_ALTIVEC
  235.   dsputil_h264_init_ppc(c, avctx);
  236.   
  237.     if (has_altivec()) {
  238.         mm_flags |= MM_ALTIVEC;
  239.         
  240.         // Altivec specific optimisations
  241.         c->pix_abs[0][1] = sad16_x2_altivec;
  242.         c->pix_abs[0][2] = sad16_y2_altivec;
  243.         c->pix_abs[0][3] = sad16_xy2_altivec;
  244.         c->pix_abs[0][0] = sad16_altivec;
  245.         c->pix_abs[1][0] = sad8_altivec;
  246.         c->sad[0]= sad16_altivec;
  247.         c->sad[1]= sad8_altivec;
  248.         c->pix_norm1 = pix_norm1_altivec;
  249.         c->sse[1]= sse8_altivec;
  250.         c->sse[0]= sse16_altivec;
  251.         c->pix_sum = pix_sum_altivec;
  252.         c->diff_pixels = diff_pixels_altivec;
  253.         c->get_pixels = get_pixels_altivec;
  254. // next one disabled as it's untested.
  255. #if 0
  256.         c->add_bytes= add_bytes_altivec;
  257. #endif /* 0 */
  258.         c->put_pixels_tab[0][0] = put_pixels16_altivec;
  259.         /* the two functions do the same thing, so use the same code */
  260.         c->put_no_rnd_pixels_tab[0][0] = put_pixels16_altivec;
  261.         c->avg_pixels_tab[0][0] = avg_pixels16_altivec;
  262.         c->avg_pixels_tab[1][0] = avg_pixels8_altivec;
  263. c->avg_pixels_tab[1][3] = avg_pixels8_xy2_altivec;
  264.         c->put_pixels_tab[1][3] = put_pixels8_xy2_altivec;
  265.         c->put_no_rnd_pixels_tab[1][3] = put_no_rnd_pixels8_xy2_altivec;
  266.         c->put_pixels_tab[0][3] = put_pixels16_xy2_altivec;
  267.         c->put_no_rnd_pixels_tab[0][3] = put_no_rnd_pixels16_xy2_altivec;
  268.         
  269. c->gmc1 = gmc1_altivec;
  270. #ifdef CONFIG_DARWIN // ATM gcc-3.3 and gcc-3.4 fail to compile these in linux...
  271. c->hadamard8_diff[0] = hadamard8_diff16_altivec;
  272. c->hadamard8_diff[1] = hadamard8_diff8x8_altivec;
  273. #endif
  274. #ifdef CONFIG_ENCODERS
  275. if (avctx->dct_algo == FF_DCT_AUTO ||
  276.     avctx->dct_algo == FF_DCT_ALTIVEC)
  277. {
  278.     c->fdct = fdct_altivec;
  279. }
  280. #endif //CONFIG_ENCODERS
  281.       if (avctx->lowres==0)
  282.       {
  283.         if ((avctx->idct_algo == FF_IDCT_AUTO) ||
  284.                 (avctx->idct_algo == FF_IDCT_ALTIVEC))
  285.         {
  286.             c->idct_put = idct_put_altivec;
  287.             c->idct_add = idct_add_altivec;
  288. #ifndef ALTIVEC_USE_REFERENCE_C_CODE
  289.             c->idct_permutation_type = FF_TRANSPOSE_IDCT_PERM;
  290. #else /* ALTIVEC_USE_REFERENCE_C_CODE */
  291.             c->idct_permutation_type = FF_NO_IDCT_PERM;
  292. #endif /* ALTIVEC_USE_REFERENCE_C_CODE */
  293.         }
  294.       }
  295.         
  296. #ifdef POWERPC_PERFORMANCE_REPORT
  297.         {
  298.           int i, j;
  299.           for (i = 0 ; i < powerpc_perf_total ; i++)
  300.           {
  301.     for (j = 0; j < POWERPC_NUM_PMC_ENABLED ; j++)
  302.       {
  303. perfdata[j][i][powerpc_data_min] = 0xFFFFFFFFFFFFFFFFULL;
  304. perfdata[j][i][powerpc_data_max] = 0x0000000000000000ULL;
  305. perfdata[j][i][powerpc_data_sum] = 0x0000000000000000ULL;
  306. perfdata[j][i][powerpc_data_num] = 0x0000000000000000ULL;
  307.       }
  308.   }
  309.         }
  310. #endif /* POWERPC_PERFORMANCE_REPORT */
  311.     } else
  312. #endif /* HAVE_ALTIVEC */
  313.     {
  314.         // Non-AltiVec PPC optimisations
  315.         // ... pending ...
  316.     }
  317. }