[BACK]Return to udiv_qrnnd.s CVS log [TXT][DIR] Up to [local] / OpenXM_contrib / gmp / mpn / hppa

Annotation of OpenXM_contrib/gmp/mpn/hppa/udiv_qrnnd.s, Revision 1.1.1.1

1.1       maekawa     1: ; HP-PA  __udiv_qrnnd division support, used from longlong.h.
                      2: ; This version runs fast on pre-PA7000 CPUs.
                      3:
                      4: ; Copyright (C) 1993, 1994 Free Software Foundation, Inc.
                      5:
                      6: ; This file is part of the GNU MP Library.
                      7:
                      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.
                     12:
                     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.
                     17:
                     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.
                     22:
                     23:
                     24: ; INPUT PARAMETERS
                     25: ; rem_ptr      gr26
                     26: ; n1           gr25
                     27: ; n0           gr24
                     28: ; d            gr23
                     29:
                     30: ; The code size is a bit excessive.  We could merge the last two ds;addc
                     31: ; sequences by simply moving the "bb,< Odd" instruction down.  The only
                     32: ; trouble is the FFFFFFFF code that would need some hacking.
                     33:
                     34:        .code
                     35:        .export         __udiv_qrnnd
                     36: __udiv_qrnnd
                     37:        .proc
                     38:        .callinfo       frame=0,no_calls
                     39:        .entry
                     40:
                     41:        comb,<          %r23,0,L$largedivisor
                     42:         sub            %r0,%r23,%r1            ; clear cy as side-effect
                     43:        ds              %r0,%r1,%r0
                     44:        addc            %r24,%r24,%r24
                     45:        ds              %r25,%r23,%r25
                     46:        addc            %r24,%r24,%r24
                     47:        ds              %r25,%r23,%r25
                     48:        addc            %r24,%r24,%r24
                     49:        ds              %r25,%r23,%r25
                     50:        addc            %r24,%r24,%r24
                     51:        ds              %r25,%r23,%r25
                     52:        addc            %r24,%r24,%r24
                     53:        ds              %r25,%r23,%r25
                     54:        addc            %r24,%r24,%r24
                     55:        ds              %r25,%r23,%r25
                     56:        addc            %r24,%r24,%r24
                     57:        ds              %r25,%r23,%r25
                     58:        addc            %r24,%r24,%r24
                     59:        ds              %r25,%r23,%r25
                     60:        addc            %r24,%r24,%r24
                     61:        ds              %r25,%r23,%r25
                     62:        addc            %r24,%r24,%r24
                     63:        ds              %r25,%r23,%r25
                     64:        addc            %r24,%r24,%r24
                     65:        ds              %r25,%r23,%r25
                     66:        addc            %r24,%r24,%r24
                     67:        ds              %r25,%r23,%r25
                     68:        addc            %r24,%r24,%r24
                     69:        ds              %r25,%r23,%r25
                     70:        addc            %r24,%r24,%r24
                     71:        ds              %r25,%r23,%r25
                     72:        addc            %r24,%r24,%r24
                     73:        ds              %r25,%r23,%r25
                     74:        addc            %r24,%r24,%r24
                     75:        ds              %r25,%r23,%r25
                     76:        addc            %r24,%r24,%r24
                     77:        ds              %r25,%r23,%r25
                     78:        addc            %r24,%r24,%r24
                     79:        ds              %r25,%r23,%r25
                     80:        addc            %r24,%r24,%r24
                     81:        ds              %r25,%r23,%r25
                     82:        addc            %r24,%r24,%r24
                     83:        ds              %r25,%r23,%r25
                     84:        addc            %r24,%r24,%r24
                     85:        ds              %r25,%r23,%r25
                     86:        addc            %r24,%r24,%r24
                     87:        ds              %r25,%r23,%r25
                     88:        addc            %r24,%r24,%r24
                     89:        ds              %r25,%r23,%r25
                     90:        addc            %r24,%r24,%r24
                     91:        ds              %r25,%r23,%r25
                     92:        addc            %r24,%r24,%r24
                     93:        ds              %r25,%r23,%r25
                     94:        addc            %r24,%r24,%r24
                     95:        ds              %r25,%r23,%r25
                     96:        addc            %r24,%r24,%r24
                     97:        ds              %r25,%r23,%r25
                     98:        addc            %r24,%r24,%r24
                     99:        ds              %r25,%r23,%r25
                    100:        addc            %r24,%r24,%r24
                    101:        ds              %r25,%r23,%r25
                    102:        addc            %r24,%r24,%r24
                    103:        ds              %r25,%r23,%r25
                    104:        addc            %r24,%r24,%r24
                    105:        ds              %r25,%r23,%r25
                    106:        addc            %r24,%r24,%r28
                    107:        ds              %r25,%r23,%r25
                    108:        comclr,>=       %r25,%r0,%r0
                    109:        addl            %r25,%r23,%r25
                    110:        stws            %r25,0(0,%r26)
                    111:        bv              0(%r2)
                    112:         addc           %r28,%r28,%r28
                    113:
                    114: L$largedivisor
                    115:        extru           %r24,31,1,%r19          ; r19 = n0 & 1
                    116:        bb,<            %r23,31,L$odd
                    117:         extru          %r23,30,31,%r22         ; r22 = d >> 1
                    118:        shd             %r25,%r24,1,%r24        ; r24 = new n0
                    119:        extru           %r25,30,31,%r25         ; r25 = new n1
                    120:        sub             %r0,%r22,%r21
                    121:        ds              %r0,%r21,%r0
                    122:        addc            %r24,%r24,%r24
                    123:        ds              %r25,%r22,%r25
                    124:        addc            %r24,%r24,%r24
                    125:        ds              %r25,%r22,%r25
                    126:        addc            %r24,%r24,%r24
                    127:        ds              %r25,%r22,%r25
                    128:        addc            %r24,%r24,%r24
                    129:        ds              %r25,%r22,%r25
                    130:        addc            %r24,%r24,%r24
                    131:        ds              %r25,%r22,%r25
                    132:        addc            %r24,%r24,%r24
                    133:        ds              %r25,%r22,%r25
                    134:        addc            %r24,%r24,%r24
                    135:        ds              %r25,%r22,%r25
                    136:        addc            %r24,%r24,%r24
                    137:        ds              %r25,%r22,%r25
                    138:        addc            %r24,%r24,%r24
                    139:        ds              %r25,%r22,%r25
                    140:        addc            %r24,%r24,%r24
                    141:        ds              %r25,%r22,%r25
                    142:        addc            %r24,%r24,%r24
                    143:        ds              %r25,%r22,%r25
                    144:        addc            %r24,%r24,%r24
                    145:        ds              %r25,%r22,%r25
                    146:        addc            %r24,%r24,%r24
                    147:        ds              %r25,%r22,%r25
                    148:        addc            %r24,%r24,%r24
                    149:        ds              %r25,%r22,%r25
                    150:        addc            %r24,%r24,%r24
                    151:        ds              %r25,%r22,%r25
                    152:        addc            %r24,%r24,%r24
                    153:        ds              %r25,%r22,%r25
                    154:        addc            %r24,%r24,%r24
                    155:        ds              %r25,%r22,%r25
                    156:        addc            %r24,%r24,%r24
                    157:        ds              %r25,%r22,%r25
                    158:        addc            %r24,%r24,%r24
                    159:        ds              %r25,%r22,%r25
                    160:        addc            %r24,%r24,%r24
                    161:        ds              %r25,%r22,%r25
                    162:        addc            %r24,%r24,%r24
                    163:        ds              %r25,%r22,%r25
                    164:        addc            %r24,%r24,%r24
                    165:        ds              %r25,%r22,%r25
                    166:        addc            %r24,%r24,%r24
                    167:        ds              %r25,%r22,%r25
                    168:        addc            %r24,%r24,%r24
                    169:        ds              %r25,%r22,%r25
                    170:        addc            %r24,%r24,%r24
                    171:        ds              %r25,%r22,%r25
                    172:        addc            %r24,%r24,%r24
                    173:        ds              %r25,%r22,%r25
                    174:        addc            %r24,%r24,%r24
                    175:        ds              %r25,%r22,%r25
                    176:        addc            %r24,%r24,%r24
                    177:        ds              %r25,%r22,%r25
                    178:        addc            %r24,%r24,%r24
                    179:        ds              %r25,%r22,%r25
                    180:        addc            %r24,%r24,%r24
                    181:        ds              %r25,%r22,%r25
                    182:        addc            %r24,%r24,%r24
                    183:        ds              %r25,%r22,%r25
                    184:        addc            %r24,%r24,%r24
                    185:        ds              %r25,%r22,%r25
                    186:        comclr,>=       %r25,%r0,%r0
                    187:        addl            %r25,%r22,%r25
                    188:        sh1addl         %r25,%r19,%r25
                    189:        stws            %r25,0(0,%r26)
                    190:        bv              0(%r2)
                    191:         addc           %r24,%r24,%r28
                    192:
                    193: L$odd  addib,sv,n      1,%r22,L$FF..           ; r22 = (d / 2 + 1)
                    194:        shd             %r25,%r24,1,%r24        ; r24 = new n0
                    195:        extru           %r25,30,31,%r25         ; r25 = new n1
                    196:        sub             %r0,%r22,%r21
                    197:        ds              %r0,%r21,%r0
                    198:        addc            %r24,%r24,%r24
                    199:        ds              %r25,%r22,%r25
                    200:        addc            %r24,%r24,%r24
                    201:        ds              %r25,%r22,%r25
                    202:        addc            %r24,%r24,%r24
                    203:        ds              %r25,%r22,%r25
                    204:        addc            %r24,%r24,%r24
                    205:        ds              %r25,%r22,%r25
                    206:        addc            %r24,%r24,%r24
                    207:        ds              %r25,%r22,%r25
                    208:        addc            %r24,%r24,%r24
                    209:        ds              %r25,%r22,%r25
                    210:        addc            %r24,%r24,%r24
                    211:        ds              %r25,%r22,%r25
                    212:        addc            %r24,%r24,%r24
                    213:        ds              %r25,%r22,%r25
                    214:        addc            %r24,%r24,%r24
                    215:        ds              %r25,%r22,%r25
                    216:        addc            %r24,%r24,%r24
                    217:        ds              %r25,%r22,%r25
                    218:        addc            %r24,%r24,%r24
                    219:        ds              %r25,%r22,%r25
                    220:        addc            %r24,%r24,%r24
                    221:        ds              %r25,%r22,%r25
                    222:        addc            %r24,%r24,%r24
                    223:        ds              %r25,%r22,%r25
                    224:        addc            %r24,%r24,%r24
                    225:        ds              %r25,%r22,%r25
                    226:        addc            %r24,%r24,%r24
                    227:        ds              %r25,%r22,%r25
                    228:        addc            %r24,%r24,%r24
                    229:        ds              %r25,%r22,%r25
                    230:        addc            %r24,%r24,%r24
                    231:        ds              %r25,%r22,%r25
                    232:        addc            %r24,%r24,%r24
                    233:        ds              %r25,%r22,%r25
                    234:        addc            %r24,%r24,%r24
                    235:        ds              %r25,%r22,%r25
                    236:        addc            %r24,%r24,%r24
                    237:        ds              %r25,%r22,%r25
                    238:        addc            %r24,%r24,%r24
                    239:        ds              %r25,%r22,%r25
                    240:        addc            %r24,%r24,%r24
                    241:        ds              %r25,%r22,%r25
                    242:        addc            %r24,%r24,%r24
                    243:        ds              %r25,%r22,%r25
                    244:        addc            %r24,%r24,%r24
                    245:        ds              %r25,%r22,%r25
                    246:        addc            %r24,%r24,%r24
                    247:        ds              %r25,%r22,%r25
                    248:        addc            %r24,%r24,%r24
                    249:        ds              %r25,%r22,%r25
                    250:        addc            %r24,%r24,%r24
                    251:        ds              %r25,%r22,%r25
                    252:        addc            %r24,%r24,%r24
                    253:        ds              %r25,%r22,%r25
                    254:        addc            %r24,%r24,%r24
                    255:        ds              %r25,%r22,%r25
                    256:        addc            %r24,%r24,%r24
                    257:        ds              %r25,%r22,%r25
                    258:        addc            %r24,%r24,%r24
                    259:        ds              %r25,%r22,%r25
                    260:        addc            %r24,%r24,%r24
                    261:        ds              %r25,%r22,%r25
                    262:        addc            %r24,%r24,%r28
                    263:        comclr,>=       %r25,%r0,%r0
                    264:        addl            %r25,%r22,%r25
                    265:        sh1addl         %r25,%r19,%r25
                    266: ; We have computed (n1,,n0) / (d + 1), q' = r28, r' = r25
                    267:        add,nuv         %r28,%r25,%r25
                    268:        addl            %r25,%r1,%r25
                    269:        addc            %r0,%r28,%r28
                    270:        sub,<<          %r25,%r23,%r0
                    271:        addl            %r25,%r1,%r25
                    272:        stws            %r25,0(0,%r26)
                    273:        bv              0(%r2)
                    274:         addc           %r0,%r28,%r28
                    275:
                    276: ; This is just a special case of the code above.
                    277: ; We come here when d == 0xFFFFFFFF
                    278: L$FF.. add,uv          %r25,%r24,%r24
                    279:        sub,<<          %r24,%r23,%r0
                    280:        ldo             1(%r24),%r24
                    281:        stws            %r24,0(0,%r26)
                    282:        bv              0(%r2)
                    283:         addc           %r0,%r25,%r28
                    284:
                    285:        .exit
                    286:        .procend

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>