annotate src/opus-1.3/celt/float_cast.h @ 69:7aeed7906520

Add Opus sources and macOS builds
author Chris Cannam
date Wed, 23 Jan 2019 13:48:08 +0000
parents
children
rev   line source
Chris@69 1 /* Copyright (C) 2001 Erik de Castro Lopo <erikd AT mega-nerd DOT com> */
Chris@69 2 /*
Chris@69 3 Redistribution and use in source and binary forms, with or without
Chris@69 4 modification, are permitted provided that the following conditions
Chris@69 5 are met:
Chris@69 6
Chris@69 7 - Redistributions of source code must retain the above copyright
Chris@69 8 notice, this list of conditions and the following disclaimer.
Chris@69 9
Chris@69 10 - Redistributions in binary form must reproduce the above copyright
Chris@69 11 notice, this list of conditions and the following disclaimer in the
Chris@69 12 documentation and/or other materials provided with the distribution.
Chris@69 13
Chris@69 14 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
Chris@69 15 ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
Chris@69 16 LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
Chris@69 17 A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER
Chris@69 18 OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
Chris@69 19 EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
Chris@69 20 PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
Chris@69 21 PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
Chris@69 22 LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
Chris@69 23 NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
Chris@69 24 SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
Chris@69 25 */
Chris@69 26
Chris@69 27 /* Version 1.1 */
Chris@69 28
Chris@69 29 #ifndef FLOAT_CAST_H
Chris@69 30 #define FLOAT_CAST_H
Chris@69 31
Chris@69 32
Chris@69 33 #include "arch.h"
Chris@69 34
Chris@69 35 /*============================================================================
Chris@69 36 ** On Intel Pentium processors (especially PIII and probably P4), converting
Chris@69 37 ** from float to int is very slow. To meet the C specs, the code produced by
Chris@69 38 ** most C compilers targeting Pentium needs to change the FPU rounding mode
Chris@69 39 ** before the float to int conversion is performed.
Chris@69 40 **
Chris@69 41 ** Changing the FPU rounding mode causes the FPU pipeline to be flushed. It
Chris@69 42 ** is this flushing of the pipeline which is so slow.
Chris@69 43 **
Chris@69 44 ** Fortunately the ISO C99 specifications define the functions lrint, lrintf,
Chris@69 45 ** llrint and llrintf which fix this problem as a side effect.
Chris@69 46 **
Chris@69 47 ** On Unix-like systems, the configure process should have detected the
Chris@69 48 ** presence of these functions. If they weren't found we have to replace them
Chris@69 49 ** here with a standard C cast.
Chris@69 50 */
Chris@69 51
Chris@69 52 /*
Chris@69 53 ** The C99 prototypes for lrint and lrintf are as follows:
Chris@69 54 **
Chris@69 55 ** long int lrintf (float x) ;
Chris@69 56 ** long int lrint (double x) ;
Chris@69 57 */
Chris@69 58
Chris@69 59 /* The presence of the required functions are detected during the configure
Chris@69 60 ** process and the values HAVE_LRINT and HAVE_LRINTF are set accordingly in
Chris@69 61 ** the config.h file.
Chris@69 62 */
Chris@69 63
Chris@69 64 /* With GCC, when SSE is available, the fastest conversion is cvtss2si. */
Chris@69 65 #if defined(__GNUC__) && defined(__SSE__)
Chris@69 66
Chris@69 67 #include <xmmintrin.h>
Chris@69 68 static OPUS_INLINE opus_int32 float2int(float x) {return _mm_cvt_ss2si(_mm_set_ss(x));}
Chris@69 69
Chris@69 70 #elif defined(HAVE_LRINTF)
Chris@69 71
Chris@69 72 /* These defines enable functionality introduced with the 1999 ISO C
Chris@69 73 ** standard. They must be defined before the inclusion of math.h to
Chris@69 74 ** engage them. If optimisation is enabled, these functions will be
Chris@69 75 ** inlined. With optimisation switched off, you have to link in the
Chris@69 76 ** maths library using -lm.
Chris@69 77 */
Chris@69 78
Chris@69 79 #define _ISOC9X_SOURCE 1
Chris@69 80 #define _ISOC99_SOURCE 1
Chris@69 81
Chris@69 82 #define __USE_ISOC9X 1
Chris@69 83 #define __USE_ISOC99 1
Chris@69 84
Chris@69 85 #include <math.h>
Chris@69 86 #define float2int(x) lrintf(x)
Chris@69 87
Chris@69 88 #elif (defined(HAVE_LRINT))
Chris@69 89
Chris@69 90 #define _ISOC9X_SOURCE 1
Chris@69 91 #define _ISOC99_SOURCE 1
Chris@69 92
Chris@69 93 #define __USE_ISOC9X 1
Chris@69 94 #define __USE_ISOC99 1
Chris@69 95
Chris@69 96 #include <math.h>
Chris@69 97 #define float2int(x) lrint(x)
Chris@69 98
Chris@69 99 #elif (defined(_MSC_VER) && _MSC_VER >= 1400) && (defined(_M_X64) || (defined(_M_IX86_FP) && _M_IX86_FP >= 1))
Chris@69 100 #include <xmmintrin.h>
Chris@69 101
Chris@69 102 __inline long int float2int(float value)
Chris@69 103 {
Chris@69 104 return _mm_cvtss_si32(_mm_load_ss(&value));
Chris@69 105 }
Chris@69 106 #elif (defined(_MSC_VER) && _MSC_VER >= 1400) && defined (_M_IX86)
Chris@69 107 #include <math.h>
Chris@69 108
Chris@69 109 /* Win32 doesn't seem to have these functions.
Chris@69 110 ** Therefore implement OPUS_INLINE versions of these functions here.
Chris@69 111 */
Chris@69 112
Chris@69 113 __inline long int
Chris@69 114 float2int (float flt)
Chris@69 115 { int intgr;
Chris@69 116
Chris@69 117 _asm
Chris@69 118 { fld flt
Chris@69 119 fistp intgr
Chris@69 120 } ;
Chris@69 121
Chris@69 122 return intgr ;
Chris@69 123 }
Chris@69 124
Chris@69 125 #else
Chris@69 126
Chris@69 127 #if (defined(__GNUC__) && defined(__STDC__) && __STDC__ && __STDC_VERSION__ >= 199901L)
Chris@69 128 /* supported by gcc in C99 mode, but not by all other compilers */
Chris@69 129 #warning "Don't have the functions lrint() and lrintf ()."
Chris@69 130 #warning "Replacing these functions with a standard C cast."
Chris@69 131 #endif /* __STDC_VERSION__ >= 199901L */
Chris@69 132 #include <math.h>
Chris@69 133 #define float2int(flt) ((int)(floor(.5+flt)))
Chris@69 134 #endif
Chris@69 135
Chris@69 136 #ifndef DISABLE_FLOAT_API
Chris@69 137 static OPUS_INLINE opus_int16 FLOAT2INT16(float x)
Chris@69 138 {
Chris@69 139 x = x*CELT_SIG_SCALE;
Chris@69 140 x = MAX32(x, -32768);
Chris@69 141 x = MIN32(x, 32767);
Chris@69 142 return (opus_int16)float2int(x);
Chris@69 143 }
Chris@69 144 #endif /* DISABLE_FLOAT_API */
Chris@69 145
Chris@69 146 #endif /* FLOAT_CAST_H */