dither.c
上传用户:aoeyumen
上传日期:2007-01-06
资源大小:3329k
文件大小:2k
源码类别:

DVD

开发平台:

Unix_Linux

  1. /* 
  2.  *    dither.c
  3.  *
  4.  * Copyright (C) Aaron Holtzman - May 1999
  5.  *
  6.  *  This file is part of ac3dec, a free Dolby AC-3 stream decoder.
  7.  *
  8.  *  ac3dec is free software; you can redistribute it and/or modify
  9.  *  it under the terms of the GNU General Public License as published by
  10.  *  the Free Software Foundation; either version 2, or (at your option)
  11.  *  any later version.
  12.  *   
  13.  *  ac3dec is distributed in the hope that it will be useful,
  14.  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
  15.  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  16.  *  GNU General Public License for more details.
  17.  *   
  18.  *  You should have received a copy of the GNU General Public License
  19.  *  along with GNU Make; see the file COPYING.  If not, write to
  20.  *  the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. 
  21.  *
  22.  */
  23. #include <stdlib.h>
  24. #include <stdio.h>
  25. #include "ac3.h"
  26. #include "dither.h"
  27. static uint_32 lfsr_state = 1;
  28. /* 
  29.  * Generate eight bits of pseudo-entropy using a 16 bit linear
  30.  * feedback shift register (LFSR). The primitive polynomial used
  31.  * is 1 + x^4 + x^14 + x^16.
  32.  *
  33.  * The distribution is uniform, over the range [-0.707,0.707]
  34.  *
  35.  */
  36. uint_16 dither_gen(void)
  37. {
  38. int i;
  39. uint_32 state;
  40. //explicitly bring the state into a local var as gcc > 3.0?
  41. //doesn't know how to optimize out the stores
  42. state = lfsr_state;
  43. //Generate eight pseudo random bits
  44. for(i=0;i<8;i++)
  45. {
  46. state <<= 1;
  47. if(state & 0x10000)
  48. state ^= 0xa011;
  49. }
  50. lfsr_state = state;
  51. return (((((sint_32)state<<8)>>8) * (sint_32) (0.707106 * 256.0))>>16);
  52. }