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