1 ! SPARC v9 32-bit __mpn_submul_1 -- Multiply a limb vector with a limb
2 ! and subtract the result from a second limb vector.
4 ! Copyright (C) 2013-2015 Free Software Foundation, Inc.
5 ! This file is part of the GNU C Library.
6 ! Contributed by David S. Miller <davem@davemloft.net>
8 ! The GNU C Library 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.
13 ! The GNU C Library 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.
18 ! You should have received a copy of the GNU Lesser General Public
19 ! License along with the GNU C Library; if not, see
20 ! <http://www.gnu.org/licenses/>.
41 srl s2l_arg, 0, s2_limb
43 be,pn %icc, .Lfinal_limb
47 lduw [s1_ptr + 0x00], tmp1
48 lduw [res_ptr + 0x00], tmp3
49 lduw [s1_ptr + 0x04], tmp2
50 lduw [res_ptr + 0x04], tmp4
51 mulx tmp1, s2_limb, tmp64_1
53 mulx tmp2, s2_limb, tmp64_2
55 add res_ptr, 8, res_ptr
56 addx carry, tmp64_1, tmp64_1
57 srlx tmp64_1, 32, carry
58 subcc tmp3, tmp64_1, tmp64_1
59 stw tmp64_1, [res_ptr - 0x08]
60 addx carry, tmp64_2, tmp64_2
61 srlx tmp64_2, 32, carry
62 subcc tmp4, tmp64_2, tmp64_2
64 stw tmp64_2, [res_ptr - 0x04]
70 lduw [s1_ptr + 0x00], tmp1
71 lduw [res_ptr + 0x00], tmp3
72 mulx tmp1, s2_limb, tmp64_1
73 addx carry, tmp64_1, tmp64_1
74 srlx tmp64_1, 32, carry
75 subcc tmp3, tmp64_1, tmp64_1
76 stw tmp64_1, [res_ptr + 0x00]