mathops.h 2.4 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879
  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_PPC_MATHOPS_H
  23. #define AVCODEC_PPC_MATHOPS_H
  24. #include <stdint.h>
  25. #include "config.h"
  26. #include "libavutil/common.h"
  27. #if HAVE_PPC4XX
  28. /* signed 16x16 -> 32 multiply add accumulate */
  29. #define MAC16(rt, ra, rb) \
  30. __asm__ ("maclhw %0, %2, %3" : "=r" (rt) : "0" (rt), "r" (ra), "r" (rb));
  31. /* signed 16x16 -> 32 multiply */
  32. #define MUL16(ra, rb) \
  33. ({ int __rt; \
  34. __asm__ ("mullhw %0, %1, %2" : "=r" (__rt) : "r" (ra), "r" (rb)); \
  35. __rt; })
  36. #endif
  37. #define MULH MULH
  38. static inline av_const int MULH(int a, int b){
  39. int r;
  40. __asm__ ("mulhw %0, %1, %2" : "=r"(r) : "r"(a), "r"(b));
  41. return r;
  42. }
  43. #if !ARCH_PPC64
  44. static inline av_const int64_t MAC64(int64_t d, int a, int b)
  45. {
  46. union { uint64_t x; unsigned hl[2]; } x = { d };
  47. int h, l;
  48. __asm__ ("mullw %3, %4, %5 \n\t"
  49. "mulhw %2, %4, %5 \n\t"
  50. "addc %1, %1, %3 \n\t"
  51. "adde %0, %0, %2 \n\t"
  52. : "+r"(x.hl[0]), "+r"(x.hl[1]), "=&r"(h), "=&r"(l)
  53. : "r"(a), "r"(b));
  54. return x.x;
  55. }
  56. #define MAC64(d, a, b) ((d) = MAC64(d, a, b))
  57. static inline av_const int64_t MLS64(int64_t d, int a, int b)
  58. {
  59. union { uint64_t x; unsigned hl[2]; } x = { d };
  60. int h, l;
  61. __asm__ ("mullw %3, %4, %5 \n\t"
  62. "mulhw %2, %4, %5 \n\t"
  63. "subfc %1, %3, %1 \n\t"
  64. "subfe %0, %2, %0 \n\t"
  65. : "+r"(x.hl[0]), "+r"(x.hl[1]), "=&r"(h), "=&r"(l)
  66. : "r"(a), "r"(b));
  67. return x.x;
  68. }
  69. #define MLS64(d, a, b) ((d) = MLS64(d, a, b))
  70. #endif
  71. #endif /* AVCODEC_PPC_MATHOPS_H */