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 file is part of FFmpeg.
7 *
8 * FFmpeg is free software; you can redistribute it and/or
9 * modify it under the terms of the GNU Lesser General Public
10 * License as published by the Free Software Foundation; either
11 * version 2.1 of the License, or (at your option) any later version.
12 *
13 * FFmpeg is distributed in the hope that it will be useful,
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
16 * Lesser General Public License for more details.
17 *
18 * You should have received a copy of the GNU Lesser General Public
19 * License along with FFmpeg; if not, write to the Free Software
20 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
21 */
22#ifndef AVCODEC_MATHOPS_H
23#define AVCODEC_MATHOPS_H
24
25#include "libavutil/common.h"
26
27#if   ARCH_X86
28
29#include "x86/mathops.h"
30
31#elif ARCH_ARM
32
33#include "arm/mathops.h"
34
35#elif ARCH_PPC
36
37#include "ppc/mathops.h"
38
39#elif ARCH_BFIN
40
41#include "bfin/mathops.h"
42
43#endif
44
45/* generic implementation */
46
47#ifndef MULL
48#   define MULL(a,b,s) (((int64_t)(a) * (int64_t)(b)) >> (s))
49#endif
50
51#ifndef MULH
52//gcc 3.4 creates an incredibly bloated mess out of this
53//#    define MULH(a,b) (((int64_t)(a) * (int64_t)(b))>>32)
54
55static av_always_inline int MULH(int a, int b){
56    return ((int64_t)(a) * (int64_t)(b))>>32;
57}
58#endif
59
60#ifndef MUL64
61#   define MUL64(a,b) ((int64_t)(a) * (int64_t)(b))
62#endif
63
64#ifndef MAC64
65#   define MAC64(d, a, b) ((d) += MUL64(a, b))
66#endif
67
68#ifndef MLS64
69#   define MLS64(d, a, b) ((d) -= MUL64(a, b))
70#endif
71
72/* signed 16x16 -> 32 multiply add accumulate */
73#ifndef MAC16
74#   define MAC16(rt, ra, rb) rt += (ra) * (rb)
75#endif
76
77/* signed 16x16 -> 32 multiply */
78#ifndef MUL16
79#   define MUL16(ra, rb) ((ra) * (rb))
80#endif
81
82#ifndef MLS16
83#   define MLS16(rt, ra, rb) ((rt) -= (ra) * (rb))
84#endif
85
86/* median of 3 */
87#ifndef mid_pred
88#define mid_pred mid_pred
89static inline av_const int mid_pred(int a, int b, int c)
90{
91#if 0
92    int t= (a-b)&((a-b)>>31);
93    a-=t;
94    b+=t;
95    b-= (b-c)&((b-c)>>31);
96    b+= (a-b)&((a-b)>>31);
97
98    return b;
99#else
100    if(a>b){
101        if(c>b){
102            if(c>a) b=a;
103            else    b=c;
104        }
105    }else{
106        if(b>c){
107            if(c>a) b=c;
108            else    b=a;
109        }
110    }
111    return b;
112#endif
113}
114#endif
115
116#endif /* AVCODEC_MATHOPS_H */
117
118