| rev |
line source |
|
nengel@2
|
1 /*
|
|
nengel@2
|
2 * simple math operations
|
|
nengel@2
|
3 * Copyright (c) 2001, 2002 Fabrice Bellard
|
|
nengel@2
|
4 * Copyright (c) 2006 Michael Niedermayer <michaelni@gmx.at> et al
|
|
nengel@2
|
5 *
|
|
nengel@2
|
6 * This file is part of FFmpeg.
|
|
nengel@2
|
7 *
|
|
nengel@2
|
8 * FFmpeg is free software; you can redistribute it and/or
|
|
nengel@2
|
9 * modify it under the terms of the GNU Lesser General Public
|
|
nengel@2
|
10 * License as published by the Free Software Foundation; either
|
|
nengel@2
|
11 * version 2.1 of the License, or (at your option) any later version.
|
|
nengel@2
|
12 *
|
|
nengel@2
|
13 * FFmpeg is distributed in the hope that it will be useful,
|
|
nengel@2
|
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
nengel@2
|
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
nengel@2
|
16 * Lesser General Public License for more details.
|
|
nengel@2
|
17 *
|
|
nengel@2
|
18 * You should have received a copy of the GNU Lesser General Public
|
|
nengel@2
|
19 * License along with FFmpeg; if not, write to the Free Software
|
|
nengel@2
|
20 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
|
|
nengel@2
|
21 */
|
|
nengel@2
|
22 #ifndef AVCODEC_MATHOPS_H
|
|
nengel@2
|
23 #define AVCODEC_MATHOPS_H
|
|
nengel@2
|
24
|
|
nengel@2
|
25 // #include "libavutil/common.h"
|
|
nengel@2
|
26 // #include "libavutil/internal.h"
|
|
nengel@2
|
27 //
|
|
nengel@2
|
28 // /* generic implementation */
|
|
nengel@2
|
29 //
|
|
nengel@2
|
30 // #ifndef MULL
|
|
nengel@2
|
31 // # define MULL(a,b,s) (((int64_t)(a) * (int64_t)(b)) >> (s))
|
|
nengel@2
|
32 // #endif
|
|
nengel@2
|
33 //
|
|
nengel@2
|
34 // #ifndef MULH
|
|
nengel@2
|
35 // //gcc 3.4 creates an incredibly bloated mess out of this
|
|
nengel@2
|
36 // //# define MULH(a,b) (((int64_t)(a) * (int64_t)(b))>>32)
|
|
nengel@2
|
37 //
|
|
nengel@2
|
38 // static av_always_inline int MULH(int a, int b){
|
|
nengel@2
|
39 // return ((int64_t)(a) * (int64_t)(b))>>32;
|
|
nengel@2
|
40 // }
|
|
nengel@2
|
41 // #endif
|
|
nengel@2
|
42 //
|
|
nengel@2
|
43 // #ifndef UMULH
|
|
nengel@2
|
44 // static av_always_inline unsigned UMULH(unsigned a, unsigned b){
|
|
nengel@2
|
45 // return ((uint64_t)(a) * (uint64_t)(b))>>32;
|
|
nengel@2
|
46 // }
|
|
nengel@2
|
47 // #endif
|
|
nengel@2
|
48 //
|
|
nengel@2
|
49 // #ifndef MUL64
|
|
nengel@2
|
50 // # define MUL64(a,b) ((int64_t)(a) * (int64_t)(b))
|
|
nengel@2
|
51 // #endif
|
|
nengel@2
|
52 //
|
|
nengel@2
|
53 // #ifndef MAC64
|
|
nengel@2
|
54 // # define MAC64(d, a, b) ((d) += MUL64(a, b))
|
|
nengel@2
|
55 // #endif
|
|
nengel@2
|
56 //
|
|
nengel@2
|
57 // #ifndef MLS64
|
|
nengel@2
|
58 // # define MLS64(d, a, b) ((d) -= MUL64(a, b))
|
|
nengel@2
|
59 // #endif
|
|
nengel@2
|
60 //
|
|
nengel@2
|
61 // /* signed 16x16 -> 32 multiply add accumulate */
|
|
nengel@2
|
62 // #ifndef MAC16
|
|
nengel@2
|
63 // # define MAC16(rt, ra, rb) rt += (ra) * (rb)
|
|
nengel@2
|
64 // #endif
|
|
nengel@2
|
65 //
|
|
nengel@2
|
66 // /* signed 16x16 -> 32 multiply */
|
|
nengel@2
|
67 // #ifndef MUL16
|
|
nengel@2
|
68 // # define MUL16(ra, rb) ((ra) * (rb))
|
|
nengel@2
|
69 // #endif
|
|
nengel@2
|
70 //
|
|
nengel@2
|
71 // #ifndef MLS16
|
|
nengel@2
|
72 // # define MLS16(rt, ra, rb) ((rt) -= (ra) * (rb))
|
|
nengel@2
|
73 // #endif
|
|
nengel@2
|
74
|
|
nengel@2
|
75 /* median of 3 */
|
|
nengel@2
|
76 #ifndef mid_pred
|
|
nengel@2
|
77 #define mid_pred mid_pred
|
|
nengel@2
|
78 static inline av_const int mid_pred(int a, int b, int c)
|
|
nengel@2
|
79 {
|
|
nengel@2
|
80 #if 0
|
|
nengel@2
|
81 int t= (a-b)&((a-b)>>31);
|
|
nengel@2
|
82 a-=t;
|
|
nengel@2
|
83 b+=t;
|
|
nengel@2
|
84 b-= (b-c)&((b-c)>>31);
|
|
nengel@2
|
85 b+= (a-b)&((a-b)>>31);
|
|
nengel@2
|
86
|
|
nengel@2
|
87 return b;
|
|
nengel@2
|
88 #else
|
|
nengel@2
|
89 if(a>b){
|
|
nengel@2
|
90 if(c>b){
|
|
nengel@2
|
91 if(c>a) b=a;
|
|
nengel@2
|
92 else b=c;
|
|
nengel@2
|
93 }
|
|
nengel@2
|
94 }else{
|
|
nengel@2
|
95 if(b>c){
|
|
nengel@2
|
96 if(c>a) b=c;
|
|
nengel@2
|
97 else b=a;
|
|
nengel@2
|
98 }
|
|
nengel@2
|
99 }
|
|
nengel@2
|
100 return b;
|
|
nengel@2
|
101 #endif
|
|
nengel@2
|
102 }
|
|
nengel@2
|
103 #endif
|
|
nengel@2
|
104
|
|
nengel@2
|
105 // #ifndef sign_extend
|
|
nengel@2
|
106 // static inline av_const int sign_extend(int val, unsigned bits)
|
|
nengel@2
|
107 // {
|
|
nengel@2
|
108 // return (val << (INT_BIT - bits)) >> (INT_BIT - bits);
|
|
nengel@2
|
109 // }
|
|
nengel@2
|
110 // #endif
|
|
nengel@2
|
111 //
|
|
nengel@2
|
112 // #ifndef zero_extend
|
|
nengel@2
|
113 // static inline av_const unsigned zero_extend(unsigned val, unsigned bits)
|
|
nengel@2
|
114 // {
|
|
nengel@2
|
115 // return (val << (INT_BIT - bits)) >> (INT_BIT - bits);
|
|
nengel@2
|
116 // }
|
|
nengel@2
|
117 // #endif
|
|
nengel@2
|
118 //
|
|
nengel@2
|
119 // #ifndef COPY3_IF_LT
|
|
nengel@2
|
120 // #define COPY3_IF_LT(x, y, a, b, c, d)\
|
|
nengel@2
|
121 // if ((y) < (x)) {\
|
|
nengel@2
|
122 // (x) = (y);\
|
|
nengel@2
|
123 // (a) = (b);\
|
|
nengel@2
|
124 // (c) = (d);\
|
|
nengel@2
|
125 // }
|
|
nengel@2
|
126 // #endif
|
|
nengel@2
|
127 //
|
|
nengel@2
|
128 // #ifndef NEG_SSR32
|
|
nengel@2
|
129 // # define NEG_SSR32(a,s) ((( int32_t)(a))>>(32-(s)))
|
|
nengel@2
|
130 // #endif
|
|
nengel@2
|
131 //
|
|
nengel@2
|
132 // #ifndef NEG_USR32
|
|
nengel@2
|
133 // # define NEG_USR32(a,s) (((uint32_t)(a))>>(32-(s)))
|
|
nengel@2
|
134 // #endif
|
|
nengel@2
|
135
|
|
nengel@2
|
136 #endif /* AVCODEC_MATHOPS_H */
|
|
nengel@2
|
137
|
