Category

Development Platform
Home > SourceCode/Document > Communication-Mobile > Voice Compress > View Code
gainncd2.f
Package: celp_3.2a.tar.gz　[view]
Upload User: szhypcb168
Upload Date: 2007-01-06
Package Size: 2187k
Code Size: 2k
Category:

Voice Compress

Development Platform:

Unix_Linux

gainncd2.f：Code Content
  1. c==========================================================================
  2. c
  3. c ROUTINE
  4. c gainencode2
  5. c
  6. c FUNCTION
  7. c  
  8. c encode and quantize code book gain
  9. c
  10. c SYNOPSIS
  11. c subroutine gainencode2(numer, denom, bits, type, index)
  12. c
  13. c   formal 
  14. c
  15. c                       data    I/O
  16. c       name            type    type    function
  17. c       -------------------------------------------------------------------
  18. c numer i i code book gain numerator
  19. c denom i i code book gain denominator
  20. c bits i i # bits for encode
  21. c type c i type of quantization
  22. c index r o encoded code book gain ZERO BASED index
  23. c gainencode2 r f encoded code book gain
  24. c
  25. c==========================================================================
  26. c
  27. c DESCRIPTION
  28. c
  29. c   Fast nonuniform division for code book gain quantization to allow
  30. c practical quantization inside the code book search loop.  A binary
  31. c tree search with cross multiply quantization is implemented below.
  32. c
  33. C==========================================================================
  34. C 
  35. C INPUT FILES
  36. C cbgain.tbl code book gain coding
  37. C
  38. c**************************************************************************
  39. c
  40. real function gainencode2(numer, denom, bits, type, index)
  41. implicit undefined(a-z)
  42. c
  43. real numer, denom
  44. integer bits, index
  45. character*10 type
  46. c
  47. c Hard coded for 5 bit quantization to achieve high speed
  48. c
  49. integer bit, bitm1
  50. parameter (bit   = 5)
  51. parameter (bitm1 = bit-1)
  52. real gainlog5(2**bit), midpoints(2**bit-1)
  53. integer i, shift(bitm1)
  54. data shift/8, 4, 2, 1/
  55. c
  56. include './cbgain.tbl'
  57. convex #include "./cbgain.tbl"
  58. c
  59. ccc if (bits .ne. bit) stop ' gainencode2:  bad bits'
  60. c
  61. c *Binary tree search for closest gain
  62. c
  63. index = 2**bitm1
  64. do 69 i = 1, bitm1
  65.    if (numer .gt. denom*midpoints(index)) then
  66.       index = index + shift(i)
  67.    else
  68.       index = index - shift(i)
  69.    end if
  70. 69 continue
  71. if (numer .gt. denom*midpoints(index)) index = index + 1
  72. c
  73. c *Return quantized gain and ZERO based index
  74. c
  75. gainencode2= gainlog5(index)
  76. index      = index - 1
  78. return
  79. end