2 * simple math operations
3 * Copyright (c) 2006 Michael Niedermayer <michaelni@gmx.at> et al
5 * This file is part of FFmpeg.
7 * FFmpeg is free software; you can redistribute it and/or
8 * modify it under the terms of the GNU Lesser General Public
9 * License as published by the Free Software Foundation; either
10 * version 2.1 of the License, or (at your option) any later version.
12 * FFmpeg is distributed in the hope that it will be useful,
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 * Lesser General Public License for more details.
17 * You should have received a copy of the GNU Lesser General Public
18 * License along with FFmpeg; if not, write to the Free Software
19 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
22 #ifndef AVCODEC_ARM_MATHOPS_H
23 #define AVCODEC_ARM_MATHOPS_H
26 #include "libavutil/common.h"
29 static inline av_const int MULL(int a, int b, unsigned shift)
32 __asm__("smull %0, %1, %2, %3 \n\t"
33 "mov %0, %0, lsr %4 \n\t"
34 "add %1, %0, %1, lsl %5 \n\t"
35 : "=&r"(lo), "=&r"(hi)
36 : "r"(b), "r"(a), "ir"(shift), "ir"(32-shift));
42 static inline av_const int MULH(int a, int b)
45 __asm__ ("smmul %0, %1, %2" : "=r"(r) : "r"(a), "r"(b));
49 static inline av_const int MULH(int a, int b)
52 __asm__ ("smull %0, %1, %2, %3" : "=&r"(lo), "=&r"(hi) : "r"(b), "r"(a));
57 static inline av_const int64_t MUL64(int a, int b)
59 union { uint64_t x; unsigned hl[2]; } x;
60 __asm__ ("smull %0, %1, %2, %3"
61 : "=r"(x.hl[0]), "=r"(x.hl[1]) : "r"(a), "r"(b));
66 static inline av_const int64_t MAC64(int64_t d, int a, int b)
68 union { uint64_t x; unsigned hl[2]; } x = { d };
69 __asm__ ("smlal %0, %1, %2, %3"
70 : "+r"(x.hl[0]), "+r"(x.hl[1]) : "r"(a), "r"(b));
73 #define MAC64(d, a, b) ((d) = MAC64(d, a, b))
74 #define MLS64(d, a, b) MAC64(d, -(a), b)
78 /* signed 16x16 -> 32 multiply add accumulate */
79 # define MAC16(rt, ra, rb) \
80 __asm__ ("smlabb %0, %1, %2, %0" : "+r"(rt) : "r"(ra), "r"(rb));
82 /* signed 16x16 -> 32 multiply */
84 static inline av_const MUL16(int ra, int rb)
87 __asm__ ("smulbb %0, %1, %2" : "=r"(rt) : "r"(ra), "r"(rb));
93 #define mid_pred mid_pred
94 static inline av_const int mid_pred(int a, int b, int c)
111 #endif /* AVCODEC_ARM_MATHOPS_H */