New single instruction math operation header
[ffmpeg.git] / libavcodec / mathops.h
1 /*
2  * simple math operations
3  * Copyright (c) 2001, 2002 Fabrice Bellard.
4  * Copyright (c) 2006 Michael Niedermayer <michaelni@gmx.at> et al
5  *
6  * This library is free software; you can redistribute it and/or
7  * modify it under the terms of the GNU Lesser General Public
8  * License as published by the Free Software Foundation; either
9  * version 2 of the License, or (at your option) any later version.
10  *
11  * This library is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14  * Lesser General Public License for more details.
15  *
16  * You should have received a copy of the GNU Lesser General Public
17  * License along with this library; if not, write to the Free Software
18  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19  */
20 #ifndef MATHOPS_H
21 #define MATHOPS_H
22
23 #ifdef ARCH_X86
24
25 #include "i386/mathops.h"
26
27 #elif defined(ARCH_ARMV4L)
28
29 #include "armv4l/mathops.h"
30
31 #elif defined(ARCH_PPC)
32
33 #include "ppc/mathops.h"
34
35 #endif
36
37 /* generic implementation */
38
39 #ifndef MULL
40 #   define MULL(a,b) (((int64_t)(a) * (int64_t)(b)) >> FRAC_BITS)
41 #endif
42
43 #ifndef MULH
44 //gcc 3.4 creates an incredibly bloated mess out of this
45 //#    define MULH(a,b) (((int64_t)(a) * (int64_t)(b))>>32)
46
47 static always_inline int MULH(int a, int b){
48     return ((int64_t)(a) * (int64_t)(b))>>32;
49 }
50 #endif
51
52 #ifndef MUL64
53 #   define MUL64(a,b) ((int64_t)(a) * (int64_t)(b))
54 #endif
55
56 /* signed 16x16 -> 32 multiply add accumulate */
57 #ifndef MAC16
58 #   define MAC16(rt, ra, rb) rt += (ra) * (rb)
59 #endif
60
61 /* signed 16x16 -> 32 multiply */
62 #ifndef MUL16
63 #   define MUL16(ra, rb) ((ra) * (rb))
64 #endif
65
66 #endif //MATHOPS_H
67