Update.
[glibc.git] / sysdeps / vax / submul_1.s
blobc473937caf95bd70045cf22a17000bbcd7723851
1 # VAX __mpn_submul_1 -- Multiply a limb vector with a limb and subtract
2 # the result from a second limb vector.
4 # Copyright (C) 1992, 1994 Free Software Foundation, Inc.
6 # This file is part of the GNU MP Library.
8 # The GNU MP Library is free software; you can redistribute it and/or modify
9 # it under the terms of the GNU Library General Public License as published by
10 # the Free Software Foundation; either version 2 of the License, or (at your
11 # option) any later version.
13 # The GNU MP Library is distributed in the hope that it will be useful, but
14 # WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
15 # or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Library General Public
16 # License for more details.
18 # You should have received a copy of the GNU Library General Public License
19 # along with the GNU MP Library; see the file COPYING.LIB. If not, write to
20 # the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
21 # MA 02111-1307, USA.
24 # INPUT PARAMETERS
25 # res_ptr (sp + 4)
26 # s1_ptr (sp + 8)
27 # size (sp + 12)
28 # s2_limb (sp + 16)
30 .text
31 .align 1
32 .globl ___mpn_submul_1
33 ___mpn_submul_1:
34 .word 0xfc0
35 movl 12(ap),r4
36 movl 8(ap),r8
37 movl 4(ap),r9
38 movl 16(ap),r6
39 jlss s2_big
41 clrl r3
42 incl r4
43 ashl $-1,r4,r7
44 jlbc r4,L1
45 clrl r11
47 # Loop for S2_LIMB < 0x80000000
48 Loop1: movl (r8)+,r1
49 jlss L1n0
50 emul r1,r6,$0,r2
51 addl2 r11,r2
52 adwc $0,r3
53 subl2 r2,(r9)+
54 adwc $0,r3
55 L1: movl (r8)+,r1
56 jlss L1n1
57 L1p1: emul r1,r6,$0,r10
58 addl2 r3,r10
59 adwc $0,r11
60 subl2 r10,(r9)+
61 adwc $0,r11
63 jsobgtr r7,Loop1
64 movl r11,r0
65 ret
67 L1n0: emul r1,r6,$0,r2
68 addl2 r11,r2
69 adwc r6,r3
70 subl2 r2,(r9)+
71 adwc $0,r3
72 movl (r8)+,r1
73 jgeq L1p1
74 L1n1: emul r1,r6,$0,r10
75 addl2 r3,r10
76 adwc r6,r11
77 subl2 r10,(r9)+
78 adwc $0,r11
80 jsobgtr r7,Loop1
81 movl r11,r0
82 ret
85 s2_big: clrl r3
86 incl r4
87 ashl $-1,r4,r7
88 jlbc r4,L2
89 clrl r11
91 # Loop for S2_LIMB >= 0x80000000
92 Loop2: movl (r8)+,r1
93 jlss L2n0
94 emul r1,r6,$0,r2
95 addl2 r11,r2
96 adwc r1,r3
97 subl2 r2,(r9)+
98 adwc $0,r3
99 L2: movl (r8)+,r1
100 jlss L2n1
101 L2p1: emul r1,r6,$0,r10
102 addl2 r3,r10
103 adwc r1,r11
104 subl2 r10,(r9)+
105 adwc $0,r11
107 jsobgtr r7,Loop2
108 movl r11,r0
111 L2n0: emul r1,r6,$0,r2
112 addl2 r11,r2
113 adwc r6,r3
114 subl2 r2,(r9)+
115 adwc r1,r3
116 movl (r8)+,r1
117 jgeq L2p1
118 L2n1: emul r1,r6,$0,r10
119 addl2 r3,r10
120 adwc r6,r11
121 subl2 r10,(r9)+
122 adwc r1,r11
124 jsobgtr r7,Loop2
125 movl r11,r0