mpyll.S 918 B

12345678910111213141516171819202122232425262728293031323334353637
  1. ;; SPDX-License-Identifier: GPL-2.0-or-later
  2. ;; Copyright (C) 2010 Texas Instruments Incorporated
  3. ;; Contributed by Mark Salter <msalter@redhat.com>.
  4. ;;
  5. #include <linux/linkage.h>
  6. ;; uint64_t __c6xabi_mpyll(uint64_t x, uint64_t y)
  7. ;;
  8. ;; 64x64 multiply
  9. ;; First compute partial results using 32-bit parts of x and y:
  10. ;;
  11. ;; b63 b32 b31 b0
  12. ;; -----------------------------
  13. ;; | 1 | 0 |
  14. ;; -----------------------------
  15. ;;
  16. ;; P0 = X0*Y0
  17. ;; P1 = X0*Y1 + X1*Y0
  18. ;; P2 = X1*Y1
  19. ;;
  20. ;; result = (P2 << 64) + (P1 << 32) + P0
  21. ;;
  22. ;; Since the result is also 64-bit, we can skip the P2 term.
  23. .text
  24. ENTRY(__c6xabi_mpyll)
  25. mpy32u .m1x A4,B4,A1:A0 ; X0*Y0
  26. b .s2 B3
  27. || mpy32u .m2x B5,A4,B1:B0 ; X0*Y1 (don't need upper 32-bits)
  28. || mpy32u .m1x A5,B4,A3:A2 ; X1*Y0 (don't need upper 32-bits)
  29. nop
  30. nop
  31. mv .s1 A0,A4
  32. add .l1x A2,B0,A5
  33. add .s1 A1,A5,A5
  34. ENDPROC(__c6xabi_mpyll)