1 /* 2 * simple math operations 3 * Copyright (c) 2006 Michael Niedermayer <michaelni@gmx.at> et al 4 * 5 * This file is part of FFmpeg. 6 * 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. 11 * 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. 16 * 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 20 */ 21 22 #ifndef FFMPEG_I386_MATHOPS_H 23 #define FFMPEG_I386_MATHOPS_H 24 25 #ifdef FRAC_BITS 26 # define MULL(ra, rb) \ 27 ({ int rt, dummy; asm (\ 28 "imull %3 \n\t"\ 29 "shrdl %4, %%edx, %%eax \n\t"\ 30 : "=a"(rt), "=d"(dummy)\ 31 : "a" (ra), "rm" (rb), "i"(FRAC_BITS));\ 32 rt; }) 33 #endif 34 35 #define MULH(ra, rb) \ 36 ({ int rt, dummy;\ 37 asm ("imull %3\n\t" : "=d"(rt), "=a"(dummy): "a" (ra), "rm" (rb));\ 38 rt; }) 39 40 #define MUL64(ra, rb) \ 41 ({ int64_t rt;\ 42 asm ("imull %2\n\t" : "=A"(rt) : "a" (ra), "g" (rb));\ 43 rt; }) 44 45 #endif /* FFMPEG_I386_MATHOPS_H */ 46