1/*
2 * This file is part of FFmpeg.
3 *
4 * FFmpeg is free software; you can redistribute it and/or
5 * modify it under the terms of the GNU Lesser General Public
6 * License as published by the Free Software Foundation; either
7 * version 2.1 of the License, or (at your option) any later version.
8 *
9 * FFmpeg is distributed in the hope that it will be useful,
10 * but WITHOUT ANY WARRANTY; without even the implied warranty of
11 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
12 * Lesser General Public License for more details.
13 *
14 * You should have received a copy of the GNU Lesser General Public
15 * License along with FFmpeg; if not, write to the Free Software
16 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
17 */
18
19/**
20 * @file
21 * Replacements for frequently missing libm functions
22 */
23
24#ifndef AVUTIL_LIBM_H
25#define AVUTIL_LIBM_H
26
27#include <math.h>
28#include "config.h"
29#include "attributes.h"
30#include "intfloat.h"
31
32#if HAVE_MIPSFPU && HAVE_INLINE_ASM
33#include "libavutil/mips/libm_mips.h"
34#endif /* HAVE_MIPSFPU && HAVE_INLINE_ASM*/
35
36#if !HAVE_ATANF
37#undef atanf
38#define atanf(x) ((float)atan(x))
39#endif
40
41#if !HAVE_ATAN2F
42#undef atan2f
43#define atan2f(y, x) ((float)atan2(y, x))
44#endif
45
46#if !HAVE_POWF
47#undef powf
48#define powf(x, y) ((float)pow(x, y))
49#endif
50
51#if !HAVE_CBRT
52static av_always_inline double cbrt(double x)
53{
54    return x < 0 ? -pow(-x, 1.0 / 3.0) : pow(x, 1.0 / 3.0);
55}
56#endif
57
58#if !HAVE_CBRTF
59static av_always_inline float cbrtf(float x)
60{
61    return x < 0 ? -powf(-x, 1.0 / 3.0) : powf(x, 1.0 / 3.0);
62}
63#endif
64
65#if !HAVE_COSF
66#undef cosf
67#define cosf(x) ((float)cos(x))
68#endif
69
70#if !HAVE_EXPF
71#undef expf
72#define expf(x) ((float)exp(x))
73#endif
74
75#if !HAVE_EXP2
76#undef exp2
77#define exp2(x) exp((x) * 0.693147180559945)
78#endif /* HAVE_EXP2 */
79
80#if !HAVE_EXP2F
81#undef exp2f
82#define exp2f(x) ((float)exp2(x))
83#endif /* HAVE_EXP2F */
84
85
86#if !HAVE_ISINF
87static av_always_inline av_const int isinf(float x)
88{
89    uint32_t v = av_float2int(x);
90    if ((v & 0x7f800000) != 0x7f800000)
91        return 0;
92    return !(v & 0x007fffff);
93}
94#endif /* HAVE_ISINF */
95
96#if !HAVE_ISNAN
97static av_always_inline av_const int isnan(float x)
98{
99    uint32_t v = av_float2int(x);
100    if ((v & 0x7f800000) != 0x7f800000)
101        return 0;
102    return v & 0x007fffff;
103}
104#endif /* HAVE_ISNAN */
105
106#if !HAVE_LDEXPF
107#undef ldexpf
108#define ldexpf(x, exp) ((float)ldexp(x, exp))
109#endif
110
111#if !HAVE_LLRINT
112#undef llrint
113#define llrint(x) ((long long)rint(x))
114#endif /* HAVE_LLRINT */
115
116#if !HAVE_LLRINTF
117#undef llrintf
118#define llrintf(x) ((long long)rint(x))
119#endif /* HAVE_LLRINT */
120
121#if !HAVE_LOG2
122#undef log2
123#define log2(x) (log(x) * 1.44269504088896340736)
124#endif /* HAVE_LOG2 */
125
126#if !HAVE_LOG2F
127#undef log2f
128#define log2f(x) ((float)log2(x))
129#endif /* HAVE_LOG2F */
130
131#if !HAVE_LOG10F
132#undef log10f
133#define log10f(x) ((float)log10(x))
134#endif
135
136#if !HAVE_SINF
137#undef sinf
138#define sinf(x) ((float)sin(x))
139#endif
140
141#if !HAVE_RINT
142static inline double rint(double x)
143{
144    return x >= 0 ? floor(x + 0.5) : ceil(x - 0.5);
145}
146#endif /* HAVE_RINT */
147
148#if !HAVE_LRINT
149static av_always_inline av_const long int lrint(double x)
150{
151    return rint(x);
152}
153#endif /* HAVE_LRINT */
154
155#if !HAVE_LRINTF
156static av_always_inline av_const long int lrintf(float x)
157{
158    return (int)(rint(x));
159}
160#endif /* HAVE_LRINTF */
161
162#if !HAVE_ROUND
163static av_always_inline av_const double round(double x)
164{
165    return (x > 0) ? floor(x + 0.5) : ceil(x - 0.5);
166}
167#endif /* HAVE_ROUND */
168
169#if !HAVE_ROUNDF
170static av_always_inline av_const float roundf(float x)
171{
172    return (x > 0) ? floor(x + 0.5) : ceil(x - 0.5);
173}
174#endif /* HAVE_ROUNDF */
175
176#if !HAVE_TRUNC
177static av_always_inline av_const double trunc(double x)
178{
179    return (x > 0) ? floor(x) : ceil(x);
180}
181#endif /* HAVE_TRUNC */
182
183#if !HAVE_TRUNCF
184static av_always_inline av_const float truncf(float x)
185{
186    return (x > 0) ? floor(x) : ceil(x);
187}
188#endif /* HAVE_TRUNCF */
189
190#endif /* AVUTIL_LIBM_H */
191