1/* SPDX-License-Identifier: GPL-2.0 */
2/*---------------------------------------------------------------------------+
3 |  mul_Xsig.S                                                               |
4 |                                                                           |
5 | Multiply a 12 byte fixed point number by another fixed point number.      |
6 |                                                                           |
7 | Copyright (C) 1992,1994,1995                                              |
8 |                       W. Metzenthen, 22 Parker St, Ormond, Vic 3163,      |
9 |                       Australia.  E-mail billm@jacobi.maths.monash.edu.au |
10 |                                                                           |
11 | Call from C as:                                                           |
12 |   void mul32_Xsig(Xsig *x, unsigned b)                                    |
13 |                                                                           |
14 |   void mul64_Xsig(Xsig *x, unsigned long long *b)                         |
15 |                                                                           |
16 |   void mul_Xsig_Xsig(Xsig *x, unsigned *b)                                |
17 |                                                                           |
18 | The result is neither rounded nor normalized, and the ls bit or so may    |
19 | be wrong.                                                                 |
20 |                                                                           |
21 +---------------------------------------------------------------------------*/
22	.file	"mul_Xsig.S"
23
24
25#include "fpu_emu.h"
26
27.text
28SYM_FUNC_START(mul32_Xsig)
29	pushl %ebp
30	movl %esp,%ebp
31	subl $16,%esp
32	pushl %esi
33
34	movl PARAM1,%esi
35	movl PARAM2,%ecx
36
37	xor %eax,%eax
38	movl %eax,-4(%ebp)
39	movl %eax,-8(%ebp)
40
41	movl (%esi),%eax        /* lsl of Xsig */
42	mull %ecx		/* msl of b */
43	movl %edx,-12(%ebp)
44
45	movl 4(%esi),%eax	/* midl of Xsig */
46	mull %ecx		/* msl of b */
47	addl %eax,-12(%ebp)
48	adcl %edx,-8(%ebp)
49	adcl $0,-4(%ebp)
50
51	movl 8(%esi),%eax	/* msl of Xsig */
52	mull %ecx		/* msl of b */
53	addl %eax,-8(%ebp)
54	adcl %edx,-4(%ebp)
55
56	movl -12(%ebp),%eax
57	movl %eax,(%esi)
58	movl -8(%ebp),%eax
59	movl %eax,4(%esi)
60	movl -4(%ebp),%eax
61	movl %eax,8(%esi)
62
63	popl %esi
64	leave
65	ret
66SYM_FUNC_END(mul32_Xsig)
67
68
69SYM_FUNC_START(mul64_Xsig)
70	pushl %ebp
71	movl %esp,%ebp
72	subl $16,%esp
73	pushl %esi
74
75	movl PARAM1,%esi
76	movl PARAM2,%ecx
77
78	xor %eax,%eax
79	movl %eax,-4(%ebp)
80	movl %eax,-8(%ebp)
81
82	movl (%esi),%eax        /* lsl of Xsig */
83	mull 4(%ecx)		/* msl of b */
84	movl %edx,-12(%ebp)
85
86	movl 4(%esi),%eax	/* midl of Xsig */
87	mull (%ecx)		/* lsl of b */
88	addl %edx,-12(%ebp)
89	adcl $0,-8(%ebp)
90	adcl $0,-4(%ebp)
91
92	movl 4(%esi),%eax	/* midl of Xsig */
93	mull 4(%ecx)		/* msl of b */
94	addl %eax,-12(%ebp)
95	adcl %edx,-8(%ebp)
96	adcl $0,-4(%ebp)
97
98	movl 8(%esi),%eax	/* msl of Xsig */
99	mull (%ecx)		/* lsl of b */
100	addl %eax,-12(%ebp)
101	adcl %edx,-8(%ebp)
102	adcl $0,-4(%ebp)
103
104	movl 8(%esi),%eax	/* msl of Xsig */
105	mull 4(%ecx)		/* msl of b */
106	addl %eax,-8(%ebp)
107	adcl %edx,-4(%ebp)
108
109	movl -12(%ebp),%eax
110	movl %eax,(%esi)
111	movl -8(%ebp),%eax
112	movl %eax,4(%esi)
113	movl -4(%ebp),%eax
114	movl %eax,8(%esi)
115
116	popl %esi
117	leave
118	ret
119SYM_FUNC_END(mul64_Xsig)
120
121
122
123SYM_FUNC_START(mul_Xsig_Xsig)
124	pushl %ebp
125	movl %esp,%ebp
126	subl $16,%esp
127	pushl %esi
128
129	movl PARAM1,%esi
130	movl PARAM2,%ecx
131
132	xor %eax,%eax
133	movl %eax,-4(%ebp)
134	movl %eax,-8(%ebp)
135
136	movl (%esi),%eax        /* lsl of Xsig */
137	mull 8(%ecx)		/* msl of b */
138	movl %edx,-12(%ebp)
139
140	movl 4(%esi),%eax	/* midl of Xsig */
141	mull 4(%ecx)		/* midl of b */
142	addl %edx,-12(%ebp)
143	adcl $0,-8(%ebp)
144	adcl $0,-4(%ebp)
145
146	movl 8(%esi),%eax	/* msl of Xsig */
147	mull (%ecx)		/* lsl of b */
148	addl %edx,-12(%ebp)
149	adcl $0,-8(%ebp)
150	adcl $0,-4(%ebp)
151
152	movl 4(%esi),%eax	/* midl of Xsig */
153	mull 8(%ecx)		/* msl of b */
154	addl %eax,-12(%ebp)
155	adcl %edx,-8(%ebp)
156	adcl $0,-4(%ebp)
157
158	movl 8(%esi),%eax	/* msl of Xsig */
159	mull 4(%ecx)		/* midl of b */
160	addl %eax,-12(%ebp)
161	adcl %edx,-8(%ebp)
162	adcl $0,-4(%ebp)
163
164	movl 8(%esi),%eax	/* msl of Xsig */
165	mull 8(%ecx)		/* msl of b */
166	addl %eax,-8(%ebp)
167	adcl %edx,-4(%ebp)
168
169	movl -12(%ebp),%edx
170	movl %edx,(%esi)
171	movl -8(%ebp),%edx
172	movl %edx,4(%esi)
173	movl -4(%ebp),%edx
174	movl %edx,8(%esi)
175
176	popl %esi
177	leave
178	ret
179SYM_FUNC_END(mul_Xsig_Xsig)
180