| OLD | NEW |
| 1 /* libFLAC - Free Lossless Audio Codec library | 1 /* libFLAC - Free Lossless Audio Codec library |
| 2 * Copyright (C) 2004,2005,2006,2007 Josh Coalson | 2 * Copyright (C) 2004-2009 Josh Coalson |
| 3 * Copyright (C) 2011-2014 Xiph.Org Foundation |
| 3 * | 4 * |
| 4 * Redistribution and use in source and binary forms, with or without | 5 * Redistribution and use in source and binary forms, with or without |
| 5 * modification, are permitted provided that the following conditions | 6 * modification, are permitted provided that the following conditions |
| 6 * are met: | 7 * are met: |
| 7 * | 8 * |
| 8 * - Redistributions of source code must retain the above copyright | 9 * - Redistributions of source code must retain the above copyright |
| 9 * notice, this list of conditions and the following disclaimer. | 10 * notice, this list of conditions and the following disclaimer. |
| 10 * | 11 * |
| 11 * - Redistributions in binary form must reproduce the above copyright | 12 * - Redistributions in binary form must reproduce the above copyright |
| 12 * notice, this list of conditions and the following disclaimer in the | 13 * notice, this list of conditions and the following disclaimer in the |
| 13 * documentation and/or other materials provided with the distribution. | 14 * documentation and/or other materials provided with the distribution. |
| 14 * | 15 * |
| 15 * - Neither the name of the Xiph.org Foundation nor the names of its | 16 * - Neither the name of the Xiph.org Foundation nor the names of its |
| 16 * contributors may be used to endorse or promote products derived from | 17 * contributors may be used to endorse or promote products derived from |
| 17 * this software without specific prior written permission. | 18 * this software without specific prior written permission. |
| 18 * | 19 * |
| 19 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS | 20 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS |
| 20 * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT | 21 * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT |
| 21 * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR | 22 * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR |
| 22 * A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR | 23 * A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR |
| 23 * CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, | 24 * CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, |
| 24 * EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, | 25 * EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, |
| 25 * PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR | 26 * PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR |
| 26 * PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF | 27 * PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF |
| 27 * LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING | 28 * LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING |
| 28 * NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS | 29 * NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS |
| 29 * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. | 30 * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. |
| 30 */ | 31 */ |
| 31 | 32 |
| 32 #if HAVE_CONFIG_H | 33 #ifdef HAVE_CONFIG_H |
| 33 # include <config.h> | 34 # include <config.h> |
| 34 #endif | 35 #endif |
| 35 | 36 |
| 36 #include "FLAC/assert.h" | 37 #include "FLAC/assert.h" |
| 37 | 38 #include "share/compat.h" |
| 38 #include "private/float.h" | 39 #include "private/float.h" |
| 39 | 40 |
| 40 #ifdef FLAC__INTEGER_ONLY_LIBRARY | 41 #ifdef FLAC__INTEGER_ONLY_LIBRARY |
| 41 | 42 |
| 42 /* adjust for compilers that can't understand using LLU suffix for uint64_t lite
rals */ | |
| 43 #ifdef _MSC_VER | |
| 44 #define FLAC__U64L(x) x | |
| 45 #else | |
| 46 #define FLAC__U64L(x) x##LLU | |
| 47 #endif | |
| 48 | |
| 49 const FLAC__fixedpoint FLAC__FP_ZERO = 0; | 43 const FLAC__fixedpoint FLAC__FP_ZERO = 0; |
| 50 const FLAC__fixedpoint FLAC__FP_ONE_HALF = 0x00008000; | 44 const FLAC__fixedpoint FLAC__FP_ONE_HALF = 0x00008000; |
| 51 const FLAC__fixedpoint FLAC__FP_ONE = 0x00010000; | 45 const FLAC__fixedpoint FLAC__FP_ONE = 0x00010000; |
| 52 const FLAC__fixedpoint FLAC__FP_LN2 = 45426; | 46 const FLAC__fixedpoint FLAC__FP_LN2 = 45426; |
| 53 const FLAC__fixedpoint FLAC__FP_E = 178145; | 47 const FLAC__fixedpoint FLAC__FP_E = 178145; |
| 54 | 48 |
| 55 /* Lookup tables for Knuth's logarithm algorithm */ | 49 /* Lookup tables for Knuth's logarithm algorithm */ |
| 56 #define LOG2_LOOKUP_PRECISION 16 | 50 #define LOG2_LOOKUP_PRECISION 16 |
| 57 static const FLAC__uint32 log2_lookup[][LOG2_LOOKUP_PRECISION] = { | 51 static const FLAC__uint32 log2_lookup[][LOG2_LOOKUP_PRECISION] = { |
| 58 { | 52 { |
| (...skipping 216 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
| 275 FLAC__uint32 FLAC__fixedpoint_log2(FLAC__uint32 x, unsigned fracbits, unsigned p
recision) | 269 FLAC__uint32 FLAC__fixedpoint_log2(FLAC__uint32 x, unsigned fracbits, unsigned p
recision) |
| 276 { | 270 { |
| 277 const FLAC__uint32 ONE = (1u << fracbits); | 271 const FLAC__uint32 ONE = (1u << fracbits); |
| 278 const FLAC__uint32 *table = log2_lookup[fracbits >> 2]; | 272 const FLAC__uint32 *table = log2_lookup[fracbits >> 2]; |
| 279 | 273 |
| 280 FLAC__ASSERT(fracbits < 32); | 274 FLAC__ASSERT(fracbits < 32); |
| 281 FLAC__ASSERT((fracbits & 0x3) == 0); | 275 FLAC__ASSERT((fracbits & 0x3) == 0); |
| 282 | 276 |
| 283 if(x < ONE) | 277 if(x < ONE) |
| 284 return 0; | 278 return 0; |
| 285 » | 279 |
| 286 if(precision > LOG2_LOOKUP_PRECISION) | 280 if(precision > LOG2_LOOKUP_PRECISION) |
| 287 precision = LOG2_LOOKUP_PRECISION; | 281 precision = LOG2_LOOKUP_PRECISION; |
| 288 | 282 |
| 289 /* Knuth's algorithm for computing logarithms, optimized for base-2 with
lookup tables */ | 283 /* Knuth's algorithm for computing logarithms, optimized for base-2 with
lookup tables */ |
| 290 { | 284 { |
| 291 FLAC__uint32 y = 0; | 285 FLAC__uint32 y = 0; |
| 292 FLAC__uint32 z = x >> 1, k = 1; | 286 FLAC__uint32 z = x >> 1, k = 1; |
| 293 while (x > ONE && k < precision) { | 287 while (x > ONE && k < precision) { |
| 294 if (x - z >= ONE) { | 288 if (x - z >= ONE) { |
| 295 x -= z; | 289 x -= z; |
| 296 z = x >> k; | 290 z = x >> k; |
| 297 y += table[k]; | 291 y += table[k]; |
| 298 } | 292 } |
| 299 else { | 293 else { |
| 300 z >>= 1; | 294 z >>= 1; |
| 301 k++; | 295 k++; |
| 302 } | 296 } |
| 303 } | 297 } |
| 304 return y; | 298 return y; |
| 305 } | 299 } |
| 306 } | 300 } |
| 307 | 301 |
| 308 #endif /* defined FLAC__INTEGER_ONLY_LIBRARY */ | 302 #endif /* defined FLAC__INTEGER_ONLY_LIBRARY */ |
| OLD | NEW |