开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 通讯/手机编程 > 语音压缩 > 查看源码
bugfix.notes
资源名称：celp_3.2a.tar.gz　[点击查看]
上传用户：szhypcb168
上传日期：2007-01-06
资源大小：2187k
文件大小：17k
源码类别：

语音压缩

开发平台：

Unix_Linux

bugfix.notes：源码内容
  1. Date: Wed, 24 Mar 93 12:30:44 -0500
  2. From: kabal@aldebaran.ee.mcgill.ca (Peter Kabal)
  3. Message-Id: <9303241730.AA07509@aldebaran.EE.McGill.CA>
  4. To: jpcampb@afterlife.ncsc.mil (Joe Campbell)
  5. Subject: celp distribution
  6. Content-Length: 17181
  7. X-Lines: 539
  8. Status: RO
  10. Joe, I grabbed the latest version from furmint.nectar.cs.cmu.edu.  Vince
  11. Cale has now put everything in a tar.Z file.  I had grabbed the stuff earlier
  12. when it was just loose in some directories.  Some of the stuff has been badly
  13. mangled.  I ran a diff with the latest version, and there is an apology
  14. in the comments to the effect that he has messed up the ifdefs.  Indeed he
  15. has, that is why the C-version does not run.
  17. I suggest you send him (and me), the correct celp.c and makefile.  That
  18. should put the distribution back in order.  There are other bits and
  19. pieces lying around which probably were not in the original distribution.
  21. Back to the Fortran version.  The version I have now runs on both Sun and
  22. DEC without further modification.  The DEC version complains about speech
  23. differences, but after processing the DEC and Sun output files measure
  24. up to within 80 dB.  The DEC output file is slightly longer than the Sun,
  25. but that is presumably due to differences in how they handle partial
  26. records at the end.
  28. Fortran routines for 1016 CELP
  29. The following renames were required:
  30. bitperm.h  ->  bitpermute.h
  31. bitprot.h  -> bitprotect.h
  32. cbgain.tlb -> cbgain.table
  33. cbsearch.f -> cbsearch.F
  34. csub.f     -> csub.F
  35. dcodham.f  -> decodeham.f
  36. disto.f    -> distortion.f
  37. encodham.f -> encodeham.f
  38. gaindcod.f -> gaindecode.f
  39. gainncd2.f -> gainencode2.f
  40. gainncod.f -> gainencode.f
  41. ldelaynr.f -> ldelay_nr.f
  42. lsp34.tbl  -> lsp34.table
  43. lspd34.f   -> lspdecode34.f
  44. mexcite.f  -> mexcite.F
  45. mtrxgen.f  -> matrixgen.f
  46. pgain.tbl  -> pagain.table
  47. ptdcode.f  -> pitchdecode.f
  48. ptncod.f   -> pitchencode.f
  49. postfilt.f -> postfilt2.f
  50. read_err.f -> read_error.f
  51. readme.mke -> readme.make
  52. savesg.f   -> save_sg.f
  53. smcbgn.f   -> smoothcbgain.f
  54. smpgn.f    -> smoothpgain.f
  55. smt.f      -> smoothtau.f
  56. wrterr.f   -> write_error.f
  58. In addition, the following routines were more extensively modified and
  59. renamed:
  60. celp.f     -> celp.F     changes + ultrix changes
  61. cli.f      -> cli.F      changes + ultrix changes
  62. iodisk.f   -> iodisk.F   changes + ultrix changes
  63. makefile.mke -> Makefile
  64. psearch.f  -> psearch.F  changes + ultrix changes
  66. Sun refers to Fortran under SunOS 4.1.3.
  67. DEC refers to DEC Fortran 3.2 running under Ultrix 4.2.
  69. Changes:
  70. (1) "implicit undefined" is a Sun'ism.  The equivalent for DEC is
  71.     "implicit none".  However, the best strategy is not to have such a
  72.     statement at all.  Both Sun & DEC Fortran accept the -u option which
  73.     turns on this option without having to modify the source.  This is
  74.     the preferred way for portable code.  This affects ALL fortran routines.
  75. (2) DEC Fortran (and the F77 standard) do not allow null character strings.
  76.     Since strings are padded with blanks anyway, the portable solution is
  77.     to use ' ' instead of ''. (cli.F)
  78. (3) Both Sun and DEC use the extension .F to signify that the file contains
  79.     conditionals which need the C-preprocessor.
  80. (4) Consider a file which is read only.  For Sun an OPEN succeeds on this
  81.     file, and only if one attempts to write will an error be signalled.
  82.     For DEC, an open fails unless the "readonly" option is specified in
  83.     the open.  The readonly option is ignored by Sun, but with a warning.
  84.     I have put conditionals (ifdef ultrix) around such changes.
  85. (5) The interpretation of the record size for unformatted reads is according
  86.     to the F77 standard in terms of "processor defined units".  Sun uses
  87.     bytes and DEC uses longwords (4 bytes).  I have put conditionals around
  88.     such changes. (iodisk.F)
  89. (6) The routine celp.f (renamed to celp.F) had a number of unexplained
  90.     errors.  First the header had been mangled so that comment characters
  91.     were not in place on some statements.  Second, the reference to
  92.     bitprotect treated it as a character, when in fact it is an integer.
  93.     Third the only routine available is postfilt2, yet the code calls out
  94.     postfilt.
  95. (7) There are some incompatibilities for an unformatted read of a character
  96.     string in routine cli.F.  The F77 standard indicates that on the
  97.     external record, the string should be in quotes.  I changed this to a
  98.     formatted read.
  99. (8) The code in psearch.F had some consistency checks on parameter values.
  100.     With DEC Fortran, these were optimized right out of the code.  This
  101.     made part of the code unreachable, causing a compilation error.
  103. Here are my diffs (not including the "implicit undefined" stuff which applies
  104. practically to all files).
  106. *** celp32/celp.f Mon Oct 26 17:19:36 1992
  107. --- ../celp-1016/celp32/celp.F Wed Mar 10 17:32:26 1993
  108. ***************
  109. *** 7,13 ****
  110.   * resulting from this code.  Further reproduction or distribution   *
  111.   * of this code without prior written permission of the U.S.         *
  112.   * Government is prohibited.  See NCS Technical Information Bulletin *
  113. ! * on FS 1016 CELP Speech Coding.   *                                                                         *   *
  114.   ***************************************************************************
  115.   C
  116.   C ROUTINE
  117. --- 7,14 ----
  118.   * resulting from this code.  Further reproduction or distribution   *
  119.   * of this code without prior written permission of the U.S.         *
  120.   * Government is prohibited.  See NCS Technical Information Bulletin *
  121. ! * on FS 1016 CELP Speech Coding.                                    *
  122. ! *                                                                         *
  123.   ***************************************************************************
  124.   C
  125.   C ROUTINE
  126. ***************
  127. *** 146,152 ****
  128.   C**************************************************************************
  129.   C*-
  130.           program celp
  131. !  implicit undefined (a-z)
  132.    include 'ccsub.com'
  133.   convex #include "ccsub.com"
  134.           real sold(maxll), snew(maxll), ssub(maxll)
  135. --- 147,153 ----
  136.   C**************************************************************************
  137.   C*-
  138.           program celp
  139. ! 
  140.    include 'ccsub.com'
  141.   convex #include "ccsub.com"
  142.           real sold(maxll), snew(maxll), ssub(maxll)
  143. ***************
  144. *** 285,292 ****
  145. --- 286,295 ----
  146.   c This issue can be reviewed by reading the UNIX man pages
  147.   c on ieee floating point (man -k ieee).
  148.   c
  149. + #ifdef sun
  150.    call abrupt_underflow()
  151.   #include <f77/f77_floatingpoint.h>
  152. + #endif
  153.   cSUN? call ieee_handler("set", "underflow", SIGFPE_DEFAULT)
  154.   c
  155.   c *** parse the command line
  156. ***************
  157. *** 377,385 ****
  158.   c
  159.   c *bit stream channel file
  160.    if (analy) then
  161. !     open(unit=25, file=ofile(1:oflen)//'.chan')
  162.    else
  163.       open(unit=25, file=ifile(1:iflen)//'.chan', status='old')
  164.    end if
  165.   c
  166.   c *quantizer design (data collection) files
  167. --- 380,393 ----
  168.   c
  169.   c *bit stream channel file
  170.    if (analy) then
  171. !     open(unit=25, file=ofile(1:oflen)//'.chan', status='new')
  172.    else
  173. + #ifdef ultrix
  174. +     open(unit=25, file=ifile(1:iflen)//'.chan', status='old',
  175. +      -         readonly)
  176. + #else
  177.       open(unit=25, file=ifile(1:iflen)//'.chan', status='old')
  178. + #endif
  179.    end if
  180.   c
  181.   c *quantizer design (data collection) files
  182. ***************
  183. *** 387,393 ****
  184. --- 395,406 ----
  185.   copt open(unit=21, file='pgain.data')
  186.   c
  187.   c *read stochastic code book vector file (by 20)
  188. + #ifdef ultrix
  189. +  open(unit=11, file='./codebook.h', status='old', err=112,
  190. +      -      readonly)
  191. + #else
  192.    open(unit=11, file='./codebook.h', status='old', err=112)
  193. + #endif
  194.    i = 1
  195.    do 110 j = 1, (2*(maxncsize-1)+l)/20
  196.       read(11, *, err=112) (x(k), k=i,i-1+20)
  197. ***************
  198. *** 400,406 ****
  199. --- 413,424 ----
  200.   115 continue
  201.   c
  202.   c *read adaptive code book index (pitch delay) file
  203. + #ifdef ultrix
  204. +  open(unit=12, file='./pdelay.h', status='old', err=122,
  205. +      -      readonly)
  206. + #else
  207.    open(unit=12, file='./pdelay.h', status='old', err=122)
  208. + #endif
  209.    do 120 i = 0, maxpd-1
  210.       read(12, *, err=122) pdelay(i)
  211.   120 continue
  212. ***************
  213. *** 412,419 ****
  214. --- 430,444 ----
  215.   c *read pitch delay coding tables for bit assignment
  216.   c *pdencode.h for encoding, pddecode.h for decoding
  217.   c *generate pdtabi for delta delay coding
  218. + #ifdef ultrix
  219. +  open(unit=13, file='./pdencode.h', status='old', err=132,
  220. +      -      readonly)
  221. +  open(unit=14, file='./pddecode.h', status='old', err=142,
  222. +      -      readonly)
  223. + #else
  224.    open(unit=13, file='./pdencode.h', status='old', err=132)
  225.    open(unit=14, file='./pddecode.h', status='old', err=142)
  226. + #endif
  227.    do 130 i = 0, maxpd-1
  228.       read(13, '(z)', err=132) pdencode(i)
  229.       read(14, *, err=142) pddecode(i)
  230. ***************
  231. *** 427,441 ****
  232.   145 continue
  233.   c
  234.   c *read bit protection file
  235.    open(unit=15, file='./bitprotect.h', status='old', err=152)
  236. !     read(15, *, err=152) bitprotect(1:codelength2)
  237.    close(15)
  238.    go to 155
  239.   152    stop ' celp:  Problem with file "bitprotect.h"'
  240.   155 continue
  241.   c
  242.   c *read bit permutation file (by 12)
  243.    open(unit=16, file='./bitpermute.h', status='old', err=162)
  244.    i = 1
  245.    do 160 j = 1, streambits/12
  246.       read(16, *, err=162) (bitpermute(k), k=i,i+11)
  247. --- 452,476 ----
  248.   145 continue
  249.   c
  250.   c *read bit protection file
  251. + #ifdef ultrix
  252. +  open(unit=15, file='./bitprotect.h', status='old', err=152,
  253. +      -      readonly)
  254. + #else
  255.    open(unit=15, file='./bitprotect.h', status='old', err=152)
  256. ! #endif
  257. !     read(15, *, err=152) (bitprotect(i),i=1,codelength2)
  258.    close(15)
  259.    go to 155
  260.   152    stop ' celp:  Problem with file "bitprotect.h"'
  261.   155 continue
  262.   c
  263.   c *read bit permutation file (by 12)
  264. + #ifdef ultrix
  265. +  open(unit=16, file='./bitpermute.h', status='old', err=162,
  266. +      -      readonly)
  267. + #else
  268.    open(unit=16, file='./bitpermute.h', status='old', err=162)
  269. + #endif
  270.    i = 1
  271.    do 160 j = 1, streambits/12
  272.       read(16, *, err=162) (bitpermute(k), k=i,i+11)
  273. ***************
  274. *** 916,922 ****
  275.   c
  276.   c       *** post filtering
  277.   c
  278. !    call postfilt(vdecoded(k),l,alpha,beta,ip,op,dp1,dp2,dp3,dp4)
  279.   canaly   call postfilt(v(k),l,alpha,beta,ip,op,dp1,dp2,dp3,dp4)
  280.   clocomp   call postfilt2(vdecoded(k),l,alpha,beta,ip,op,dp1,dp2,dp3,dp4)
  281.   c
  282. --- 951,957 ----
  283.   c
  284.   c       *** post filtering
  285.   c
  286. !    call postfilt2(vdecoded(k),l,alpha,beta,ip,op,dp1,dp2,dp3,dp4)
  287.   canaly   call postfilt(v(k),l,alpha,beta,ip,op,dp1,dp2,dp3,dp4)
  288.   clocomp   call postfilt2(vdecoded(k),l,alpha,beta,ip,op,dp1,dp2,dp3,dp4)
  289.   c
  290. *** celp32/cli.f Mon Oct 26 17:19:37 1992
  291. --- ../celp-1016/celp32/cli.F Wed Mar 10 20:28:50 1993
  292. ***************
  293. *** 152,158 ****
  294.   c
  295.   c note:  ncsize, no, and gamma are passed through common
  296.   c
  297. -  implicit undefined(a-z)
  298.    include 'ccsub.com'
  299.   convex  #include "ccsub.com"
  300.    character*80 ifile, ofile, stype*10
  301. --- 152,157 ----
  302. ***************
  303. *** 174,180 ****
  304.    logical flag, qflag
  305.   c
  306.    data switch /'-i','-o', '-p', '-q', '-m', '-l', '-c'/
  307. !  data pfile /''/, qfile /''/, mfile /''/, lfile /'celp.log'/
  308.    data type /'cbgain', 'pitch', 'spectrum'/
  309.    data qtype /'kang','log','max2','arcsin','none'/
  310.   c
  311. --- 173,179 ----
  312.    logical flag, qflag
  313.   c
  314.    data switch /'-i','-o', '-p', '-q', '-m', '-l', '-c'/
  315. !  data pfile /' '/, qfile /' '/, mfile /' '/, lfile /'celp.log'/
  316.    data type /'cbgain', 'pitch', 'spectrum'/
  317.    data qtype /'kang','log','max2','arcsin','none'/
  318.   c
  319. ***************
  320. *** 225,232 ****
  321.   c *** read input parameter file if requested
  322.   c           (all parameters are required)
  323.   c
  324. !  if (pfile .ne. '') then
  325.       open(unit=69, file=pfile)
  326.       read(69, *) ncsize 
  327.       read(69, *) l
  328.       read(69, *) ll
  329. --- 224,235 ----
  330.   c *** read input parameter file if requested
  331.   c           (all parameters are required)
  332.   c
  333. !  if (pfile .ne. ' ') then
  334. ! #ifdef ultrix
  335. !     open(unit=69, file=pfile, status='OLD', readonly)
  336. ! #else
  337.       open(unit=69, file=pfile)
  338. + #endif
  339.       read(69, *) ncsize 
  340.       read(69, *) l
  341.       read(69, *) ll
  342. ***************
  343. *** 243,257 ****
  344.       if (mxsw2 .ne. 1) mxsw = .false.
  345.   c
  346.       read(69, *) prewt
  347. !     read(69, *) pstype
  348.       close(unit=69)
  349.    end if
  350.   c
  351.   c *** read quantization characteristics file if requested
  352.   c           (not all parameters are required)
  353.   c
  354. !  if (qfile .ne. '') then
  355.       open(unit=70, file=qfile)
  356.       do 300 i = 1, maxtype
  357.          read(70, 69) temp
  358.   69       format(a)
  359. --- 246,264 ----
  360.       if (mxsw2 .ne. 1) mxsw = .false.
  361.   c
  362.       read(69, *) prewt
  363. !     read(69, '(A)') pstype
  364.       close(unit=69)
  365.    end if
  366.   c
  367.   c *** read quantization characteristics file if requested
  368.   c           (not all parameters are required)
  369.   c
  370. !  if (qfile .ne. ' ') then
  371. ! #ifdef ultrix
  372. !     open(unit=70, file=qfile, status='OLD', readonly)
  373. ! #else
  374.       open(unit=70, file=qfile)
  375. + #endif
  376.       do 300 i = 1, maxtype
  377.          read(70, 69) temp
  378.   69       format(a)
  379. ***************
  380. *** 299,306 ****
  381.   c
  382.   c *** read bit error characteristics file if requested
  383.   c
  384. !  if (mfile .ne. '') then
  385.       open(unit=72, file=mfile)
  386.       do 74 i = 1, bits
  387.          read(72, 73) mask(i)
  388.   73       format(i1)
  389. --- 306,317 ----
  390.   c
  391.   c *** read bit error characteristics file if requested
  392.   c
  393. !  if (mfile .ne. ' ') then
  394. ! #ifdef ultrix
  395. !     open(unit=72, file=mfile, status='OLD', readonly)
  396. ! #else
  397.       open(unit=72, file=mfile)
  398. + #endif
  399.       do 74 i = 1, bits
  400.          read(72, 73) mask(i)
  401.   73       format(i1)
  402. ***************
  403. *** 397,408 ****
  404. --- 408,427 ----
  405.    if (lfile .ne. 'none') then  
  406.       if (lfile .eq. 'celp.log') then
  407.   cvax       open(unit=71, file=lfile, access='append')
  408. + #ifdef ultrix
  409. +        open(unit=71, file=lfile, status='UNKNOWN')
  410. + #else
  411.          open(unit=71, file=lfile)
  412. + #endif
  413.   601       read(71, 69, end=602) temp
  414.          goto 601
  415.   602       continue
  416.       else
  417. + #ifdef ultrix
  418. +        open(unit=71, file=lfile, status='NEW')
  419. + #else
  420.          open(unit=71, file=lfile)
  421. + #endif
  422.       end if
  423.       write(71, *) getlog(), host, fdate()
  424.       write(71, *)'          Program: '//prog(1:lnblnk(prog))
  425. *** celp32/iodisk.f Mon Oct 26 17:19:42 1992
  426. --- ../celp-1016/celp32/iodisk.F Fri Mar 12 11:11:33 1993
  427. ***************
  428. *** 41,47 ****
  429.   C**************************************************************************
  430.   C*-
  431.    function iodisk(mode, lun, fname, nrec, iar, size)
  432. !  implicit undefined(a-z)
  433.    integer mode, lun, nrec, size, iodisk
  434.    integer*2 iar(size)    
  435.    character*(*) fname
  436. --- 41,47 ----
  437.   C**************************************************************************
  438.   C*-
  439.    function iodisk(mode, lun, fname, nrec, iar, size)
  440. ! 
  441.    integer mode, lun, nrec, size, iodisk
  442.    integer*2 iar(size)    
  443.    character*(*) fname
  444. ***************
  445. *** 68,81 ****
  446. --- 68,93 ----
  447.       nrec = nrec+1
  448.    else if (mode .eq. 3) then
  449.   c *open file for read
  450. + #ifdef ultrix
  451. +     if (mod(size,2).ne.0) stop 'frame size must be even'
  452.       open(lun, file=fname, access='direct', status='old',
  453. +      +         recl=size/2, readonly, err=333)
  454. + #else
  455. +     open(lun, file=fname, access='direct', status='old',
  456.        +         recl=2*size, err=333)
  457. + #endif
  458.    else if (mode .eq. 4) then
  459.   c *open file for write
  460.   c *unlink = UNIX file delete
  461.       call unlink(fname)
  462. + #ifdef ultrix
  463. +     if (mod(size,2).ne.0) stop 'frame size must be even'
  464. +     open(lun, file=fname, access='direct', status='new',
  465. +      +         recl=size/2, err=444)
  466. + #else
  467.       open(lun, file=fname, access='direct', status='new',
  468.        +         recl=2*size, err=444)
  469. + #endif
  470.    else if (mode .eq. -1) then
  471.   c *close file
  472.       close(lun, err=555)
  473. *** celp32/psearch.f Mon Oct 26 17:19:50 1992
  474. --- ../celp-1016/celp32/psearch.F Wed Mar 10 20:33:25 1993
  475. ***************
  476. *** 83,89 ****
  477.   C**************************************************************************
  478.   C*-
  479.    subroutine psearch(s, l)
  480. !  implicit undefined (a-z)
  481.    integer l
  482.    real s(l)
  483.    include 'ccsub.com'
  484. --- 83,89 ----
  485.   C**************************************************************************
  486.   C*-
  487.    subroutine psearch(s, l)
  488. ! 
  489.    integer l
  490.    real s(l)
  491.    include 'ccsub.com'
  492. ***************
  493. *** 102,107 ****
  494. --- 102,110 ----
  495.    real v0(maxbufptr), v0shift(maxlp), frac
  496.   c *length of truncated impulse response
  497.    parameter (len = 30)
  498. +  save oldptr, submult
  499. + c
  500. +  data oldptr/1/
  501.   c
  502.   c Choose type of pitch delay search:
  503.   c
  504. ***************
  505. *** 132,146 ****
  506.       stop ' psearch: incorrect pitch search type (pstype)'
  507.    end if
  508.   c
  509. -  save oldptr, submult
  510. - c
  511. -  data oldptr/1/
  512.   c
  513.   c
  514.    if (len .gt. l) stop ' psearch:  impulse response too long'
  515.   c
  516.    bufptr = mmax+no+2*l+maxnp-1
  517. !  if (maxlp .lt. maxl) stop ' psearch:  maxlp < maxl'
  518.   #ifdef SUNGRAPH
  519.   c *save impulse response in file 'pitch'
  520.    call save_sg(pitch_ir_vid,h,l,'save pitch_ir_vid')
  521. --- 135,146 ----
  522.       stop ' psearch: incorrect pitch search type (pstype)'
  523.    end if
  524.   c
  525.   c
  526.   c
  527.    if (len .gt. l) stop ' psearch:  impulse response too long'
  528.   c
  529.    bufptr = mmax+no+2*l+maxnp-1
  530. ! 
  531.   #ifdef SUNGRAPH
  532.   c *save impulse response in file 'pitch'
  533.    call save_sg(pitch_ir_vid,h,l,'save pitch_ir_vid')
  534. ***************
  535. *** 159,165 ****
  536. --- 159,170 ----
  537.       bb(3) = 0.0
  538.       bb(1) = mmin
  539.   c *load pitch submultiple delay table
  540. + #ifdef ultrix
  541. +     open(unit=8, file='./submult.h', status='old', err=22,
  542. +      -         readonly)
  543. + #else
  544.       open(unit=8, file='./submult.h', status='old', err=22)
  545. + #endif
  546.       do 20 i = 0, maxpd-1
  547.          read(8, *, err=22) (submult(i,j), j=0,3)
  548.   20    continue