aboutsummaryrefslogtreecommitdiff
path: root/arch/c6x/lib/mpyll.S
blob: d07c13ec4fd4c1a54c522b4d5be4b77b201d7946 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
;; SPDX-License-Identifier: GPL-2.0-or-later
;;  Copyright (C) 2010 Texas Instruments Incorporated
;;  Contributed by Mark Salter <msalter@redhat.com>.
;;

#include <linux/linkage.h>

	;; uint64_t __c6xabi_mpyll(uint64_t x, uint64_t y)
	;;
	;; 64x64 multiply
	;; First compute partial results using 32-bit parts of x and y:
	;;
	;;   b63	 b32 b31	  b0
	;;    -----------------------------
	;;    |      1	    |	   0	  |
	;;    -----------------------------
	;;
	;;   P0 = X0*Y0
	;;   P1 = X0*Y1 + X1*Y0
	;;   P2 = X1*Y1
	;;
	;;   result = (P2 << 64) + (P1 << 32) + P0
	;;
	;; Since the result is also 64-bit, we can skip the P2 term.

	.text
ENTRY(__c6xabi_mpyll)
	mpy32u	.m1x	A4,B4,A1:A0	; X0*Y0
	b	.s2	B3
 ||	mpy32u	.m2x	B5,A4,B1:B0	; X0*Y1 (don't need upper 32-bits)
 ||	mpy32u	.m1x	A5,B4,A3:A2	; X1*Y0 (don't need upper 32-bits)
	nop
	nop
	mv	.s1	A0,A4
	add	.l1x	A2,B0,A5
	add	.s1	A1,A5,A5
ENDPROC(__c6xabi_mpyll)