Mercurial > libavcodec.hg
comparison ra288.c @ 7913:3551c55a7ccb libavcodec
Merge the 2 lpc loops.
This changes the output by (stddev: 0.21 PSNR:109.51, file:ddinterview28.ra),
the changes are due to float rounding inaccuracies and do not happen if doubles
are used.
| author | michael |
|---|---|
| date | Tue, 23 Sep 2008 23:16:48 +0000 |
| parents | 9fbae87668db |
| children | 2d530c3d155e |
comparison
equal
deleted
inserted
replaced
| 7912:9fbae87668db | 7913:3551c55a7ccb |
|---|---|
| 101 memmove(gain_block, gain_block + 1, 9 * sizeof(*gain_block)); | 101 memmove(gain_block, gain_block + 1, 9 * sizeof(*gain_block)); |
| 102 | 102 |
| 103 gain_block[9] = 10 * log10(sum) - 32; | 103 gain_block[9] = 10 * log10(sum) - 32; |
| 104 | 104 |
| 105 for (i=0; i < 5; i++) { | 105 for (i=0; i < 5; i++) { |
| 106 block[i] = 0; | 106 block[i] = buffer[i]; |
| 107 for (j=0; j < 36; j++) | 107 for (j=0; j < 36; j++) |
| 108 block[i] -= block[i-1-j]*ractx->sp_lpc[j]; | 108 block[i] -= block[i-1-j]*ractx->sp_lpc[j]; |
| 109 for (j=0; j < i; j++) | |
| 110 buffer[i] -= buffer[i-1-j]*ractx->sp_lpc[j]; | |
| 111 } | 109 } |
| 112 | 110 |
| 113 /* output */ | 111 /* output */ |
| 114 for (i=0; i < 5; i++) | 112 for (i=0; i < 5; i++) |
| 115 block[i] = av_clipf(block[i] + buffer[i], -4095, 4095); | 113 block[i] = av_clipf(block[i], -4095, 4095); |
| 116 } | 114 } |
| 117 | 115 |
| 118 static void convolve(float *tgt, const float *src, int len, int n) | 116 static void convolve(float *tgt, const float *src, int len, int n) |
| 119 { | 117 { |
| 120 for (; n >= 0; n--) | 118 for (; n >= 0; n--) |