1dnl SPARC v7 __udiv_qrnnd division support, used from longlong.h. 2dnl This is for v7 CPUs without a floating-point unit. 3 4dnl Copyright 1993, 1994, 1996, 2000 Free Software Foundation, Inc. 5 6dnl This file is part of the GNU MP Library. 7 8dnl The GNU MP Library is free software; you can redistribute it and/or modify 9dnl it under the terms of the GNU Lesser General Public License as published 10dnl by the Free Software Foundation; either version 3 of the License, or (at 11dnl your option) any later version. 12 13dnl The GNU MP Library is distributed in the hope that it will be useful, but 14dnl WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 15dnl or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public 16dnl License for more details. 17 18dnl You should have received a copy of the GNU Lesser General Public License 19dnl along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. 20 21 22include(`../config.m4') 23 24C INPUT PARAMETERS 25C rem_ptr o0 26C n1 o1 27C n0 o2 28C d o3 29 30ASM_START() 31PROLOGUE(mpn_udiv_qrnnd) 32 tst %o3 33 bneg L(largedivisor) 34 mov 8,%g1 35 36 b L(p1) 37 addxcc %o2,%o2,%o2 38 39L(plop): 40 bcc L(n1) 41 addxcc %o2,%o2,%o2 42L(p1): addx %o1,%o1,%o1 43 subcc %o1,%o3,%o4 44 bcc L(n2) 45 addxcc %o2,%o2,%o2 46L(p2): addx %o1,%o1,%o1 47 subcc %o1,%o3,%o4 48 bcc L(n3) 49 addxcc %o2,%o2,%o2 50L(p3): addx %o1,%o1,%o1 51 subcc %o1,%o3,%o4 52 bcc L(n4) 53 addxcc %o2,%o2,%o2 54L(p4): addx %o1,%o1,%o1 55 addcc %g1,-1,%g1 56 bne L(plop) 57 subcc %o1,%o3,%o4 58 bcc L(n5) 59 addxcc %o2,%o2,%o2 60L(p5): st %o1,[%o0] 61 retl 62 xnor %g0,%o2,%o0 63 64L(nlop): 65 bcc L(p1) 66 addxcc %o2,%o2,%o2 67L(n1): addx %o4,%o4,%o4 68 subcc %o4,%o3,%o1 69 bcc L(p2) 70 addxcc %o2,%o2,%o2 71L(n2): addx %o4,%o4,%o4 72 subcc %o4,%o3,%o1 73 bcc L(p3) 74 addxcc %o2,%o2,%o2 75L(n3): addx %o4,%o4,%o4 76 subcc %o4,%o3,%o1 77 bcc L(p4) 78 addxcc %o2,%o2,%o2 79L(n4): addx %o4,%o4,%o4 80 addcc %g1,-1,%g1 81 bne L(nlop) 82 subcc %o4,%o3,%o1 83 bcc L(p5) 84 addxcc %o2,%o2,%o2 85L(n5): st %o4,[%o0] 86 retl 87 xnor %g0,%o2,%o0 88 89L(largedivisor): 90 and %o2,1,%o5 C %o5 = n0 & 1 91 92 srl %o2,1,%o2 93 sll %o1,31,%g2 94 or %g2,%o2,%o2 C %o2 = lo(n1n0 >> 1) 95 srl %o1,1,%o1 C %o1 = hi(n1n0 >> 1) 96 97 and %o3,1,%g2 98 srl %o3,1,%g3 C %g3 = floor(d / 2) 99 add %g3,%g2,%g3 C %g3 = ceil(d / 2) 100 101 b L(Lp1) 102 addxcc %o2,%o2,%o2 103 104L(Lplop): 105 bcc L(Ln1) 106 addxcc %o2,%o2,%o2 107L(Lp1): addx %o1,%o1,%o1 108 subcc %o1,%g3,%o4 109 bcc L(Ln2) 110 addxcc %o2,%o2,%o2 111L(Lp2): addx %o1,%o1,%o1 112 subcc %o1,%g3,%o4 113 bcc L(Ln3) 114 addxcc %o2,%o2,%o2 115L(Lp3): addx %o1,%o1,%o1 116 subcc %o1,%g3,%o4 117 bcc L(Ln4) 118 addxcc %o2,%o2,%o2 119L(Lp4): addx %o1,%o1,%o1 120 addcc %g1,-1,%g1 121 bne L(Lplop) 122 subcc %o1,%g3,%o4 123 bcc L(Ln5) 124 addxcc %o2,%o2,%o2 125L(Lp5): add %o1,%o1,%o1 C << 1 126 tst %g2 127 bne L(oddp) 128 add %o5,%o1,%o1 129 st %o1,[%o0] 130 retl 131 xnor %g0,%o2,%o0 132 133L(Lnlop): 134 bcc L(Lp1) 135 addxcc %o2,%o2,%o2 136L(Ln1): addx %o4,%o4,%o4 137 subcc %o4,%g3,%o1 138 bcc L(Lp2) 139 addxcc %o2,%o2,%o2 140L(Ln2): addx %o4,%o4,%o4 141 subcc %o4,%g3,%o1 142 bcc L(Lp3) 143 addxcc %o2,%o2,%o2 144L(Ln3): addx %o4,%o4,%o4 145 subcc %o4,%g3,%o1 146 bcc L(Lp4) 147 addxcc %o2,%o2,%o2 148L(Ln4): addx %o4,%o4,%o4 149 addcc %g1,-1,%g1 150 bne L(Lnlop) 151 subcc %o4,%g3,%o1 152 bcc L(Lp5) 153 addxcc %o2,%o2,%o2 154L(Ln5): add %o4,%o4,%o4 C << 1 155 tst %g2 156 bne L(oddn) 157 add %o5,%o4,%o4 158 st %o4,[%o0] 159 retl 160 xnor %g0,%o2,%o0 161 162L(oddp): 163 xnor %g0,%o2,%o2 164 C q' in %o2. r' in %o1 165 addcc %o1,%o2,%o1 166 bcc L(Lp6) 167 addx %o2,0,%o2 168 sub %o1,%o3,%o1 169L(Lp6): subcc %o1,%o3,%g0 170 bcs L(Lp7) 171 subx %o2,-1,%o2 172 sub %o1,%o3,%o1 173L(Lp7): st %o1,[%o0] 174 retl 175 mov %o2,%o0 176 177L(oddn): 178 xnor %g0,%o2,%o2 179 C q' in %o2. r' in %o4 180 addcc %o4,%o2,%o4 181 bcc L(Ln6) 182 addx %o2,0,%o2 183 sub %o4,%o3,%o4 184L(Ln6): subcc %o4,%o3,%g0 185 bcs L(Ln7) 186 subx %o2,-1,%o2 187 sub %o4,%o3,%o4 188L(Ln7): st %o4,[%o0] 189 retl 190 mov %o2,%o0 191EPILOGUE(mpn_udiv_qrnnd) 192