sad.c
上传用户:sun1608
上传日期:2007-02-02
资源大小:6116k
文件大小:3k
源码类别:

流媒体/Mpeg4/MP4

开发平台:

Visual C++

  1. /**************************************************************************
  2.  *
  3.  * XVID MPEG-4 VIDEO CODEC
  4.  * sum of absolute difference
  5.  *
  6.  * This program is an implementation of a part of one or more MPEG-4
  7.  * Video tools as specified in ISO/IEC 14496-2 standard.  Those intending
  8.  * to use this software module in hardware or software products are
  9.  * advised that its use may infringe existing patents or copyrights, and
  10.  * any such use would be at such party's own risk.  The original
  11.  * developer of this software module and his/her company, and subsequent
  12.  * editors and their companies, will have no liability for use of this
  13.  * software or modifications or derivatives thereof.
  14.  *
  15.  * This program is free software; you can redistribute it and/or modify
  16.  * it under the terms of the GNU General Public License as published by
  17.  * the Free Software Foundation; either version 2 of the License, or
  18.  * (at your option) any later version.
  19.  *
  20.  * This program is distributed in the hope that it will be useful,
  21.  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  22.  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  23.  * GNU General Public License for more details.
  24.  *
  25.  * You should have received a copy of the GNU General Public License
  26.  * along with this program; if not, write to the Free Software
  27.  * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  28.  *
  29.  *************************************************************************/
  30. /**************************************************************************
  31.  *
  32.  * History:
  33.  *
  34.  * 10.11.2001 initial version; (c)2001 peter ross <pross@cs.rmit.edu.au>
  35.  *
  36.  *************************************************************************/
  37. #include "../portab.h"
  38. #include "sad.h"
  39. sad16FuncPtr sad16;
  40. sad8FuncPtr sad8;
  41. dev16FuncPtr dev16;
  42. sadInitFuncPtr sadInit;
  43. #define ABS(X) (((X)>0)?(X):-(X))
  44. uint32_t sad16_c(const uint8_t * const cur,
  45.  const uint8_t * const ref,
  46.  const uint32_t stride,
  47.  const uint32_t best_sad)
  48. {
  49. uint32_t sad = 0;
  50. uint32_t i,j;
  51. uint8_t const *ptr_cur = cur;
  52. uint8_t const *ptr_ref = ref;
  53. for (j = 0; j < 16; j++) {
  54. for (i = 0; i < 16; i++) {
  55. sad += ABS(*(ptr_cur + i) - *(ptr_ref + i));
  56. if (sad >= best_sad) {
  57. return sad;
  58. }
  59. }
  60. ptr_cur += stride;
  61. ptr_ref += stride;
  62. }
  63. return sad;
  64. }
  65. uint32_t sad8_c(const uint8_t * const cur,
  66. const uint8_t * const ref,
  67. const uint32_t stride)
  68. {
  69. uint32_t sad = 0;
  70. uint32_t i,j;
  71. uint8_t const *ptr_cur = cur;
  72. uint8_t const *ptr_ref = ref;
  73. for (j = 0; j < 8; j++) {
  74. for (i = 0; i < 8; i++) {
  75. sad += ABS(*(ptr_cur + i) - *(ptr_ref + i));
  76. }
  77. ptr_cur += stride;
  78. ptr_ref += stride;
  79. }
  80. return sad;
  81. }
  82. /* average deviation from mean */
  83. uint32_t dev16_c(const uint8_t * const cur,
  84. const uint32_t stride)
  85. {
  86. uint32_t mean = 0;
  87. uint32_t dev = 0;
  88. uint32_t i,j;
  89. uint8_t const *ptr_cur = cur;
  90. for (j = 0; j < 16; j++) {
  91. for (i = 0; i < 16; i++)
  92. mean += *(ptr_cur + i);
  93. ptr_cur += stride;
  94. }
  95. mean /= (16 * 16);
  96. ptr_cur = cur;
  97. for (j = 0; j < 16; j++) {
  98. for (i = 0; i < 16; i++)
  99. dev += ABS(*(ptr_cur + i) - (int32_t)mean);
  100. ptr_cur += stride;
  101. }
  102. return dev;
  103. }