Chromium Code Reviews| OLD | NEW |
|---|---|
| 1 /* | 1 /* |
| 2 * Copyright (c) 2011 The WebRTC project authors. All Rights Reserved. | 2 * Copyright (c) 2011 The WebRTC project authors. All Rights Reserved. |
| 3 * | 3 * |
| 4 * Use of this source code is governed by a BSD-style license | 4 * Use of this source code is governed by a BSD-style license |
| 5 * that can be found in the LICENSE file in the root of the source | 5 * that can be found in the LICENSE file in the root of the source |
| 6 * tree. An additional intellectual property rights grant can be found | 6 * tree. An additional intellectual property rights grant can be found |
| 7 * in the file PATENTS. All contributing project authors may | 7 * in the file PATENTS. All contributing project authors may |
| 8 * be found in the AUTHORS file in the root of the source tree. | 8 * be found in the AUTHORS file in the root of the source tree. |
| 9 */ | 9 */ |
| 10 | 10 |
| (...skipping 47 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... | |
| 58 lspPtr+=2; | 58 lspPtr+=2; |
| 59 | 59 |
| 60 for(i=2; i<=5; i++) | 60 for(i=2; i<=5; i++) |
| 61 { | 61 { |
| 62 (*fPtr) = fPtr[-2]; | 62 (*fPtr) = fPtr[-2]; |
| 63 | 63 |
| 64 for(j=i; j>1; j--) | 64 for(j=i; j>1; j--) |
| 65 { | 65 { |
| 66 /* Compute f[j] = f[j] + tmp*f[j-1] + f[j-2]; */ | 66 /* Compute f[j] = f[j] + tmp*f[j-1] + f[j-2]; */ |
| 67 high = (int16_t)(fPtr[-1] >> 16); | 67 high = (int16_t)(fPtr[-1] >> 16); |
| 68 low = (int16_t)((fPtr[-1] - ((int32_t)high << 16)) >> 1); | 68 low = (int16_t)((fPtr[-1] & 0xffff) >> 1); |
|
tlegrand
2016/05/17 12:24:32
This is so much nicer!
kwiberg-webrtc
2016/05/17 12:32:29
I've found the "x_low = x - x_high << 16" pattern
| |
| 69 | 69 |
| 70 tmpW32 = ((high * *lspPtr) << 2) + (((low * *lspPtr) >> 15) << 2); | 70 tmpW32 = 4 * high * *lspPtr + 4 * ((low * *lspPtr) >> 15); |
| 71 | 71 |
| 72 (*fPtr) += fPtr[-2]; | 72 (*fPtr) += fPtr[-2]; |
| 73 (*fPtr) -= tmpW32; | 73 (*fPtr) -= tmpW32; |
| 74 fPtr--; | 74 fPtr--; |
| 75 } | 75 } |
| 76 *fPtr -= *lspPtr << 10; | 76 *fPtr -= *lspPtr * (1 << 10); |
| 77 | 77 |
| 78 fPtr+=i; | 78 fPtr+=i; |
| 79 lspPtr+=2; | 79 lspPtr+=2; |
| 80 } | 80 } |
| 81 return; | 81 return; |
| 82 } | 82 } |
| OLD | NEW |