| wdenk | c609719 | 2002-11-03 00:24:07 +0000 | [diff] [blame] | 1 | /* # 1 "libgcc1.S" */ |
| 2 | @ libgcc1 routines for ARM cpu. |
| 3 | @ Division routines, written by Richard Earnshaw, (rearnsha@armltd.co.uk) |
| 4 | dividend .req r0 |
| 5 | divisor .req r1 |
| 6 | result .req r2 |
| 7 | curbit .req r3 |
| wdenk | ac6dbb8 | 2003-03-26 11:42:53 +0000 | [diff] [blame] | 8 | /* ip .req r12 */ |
| 9 | /* sp .req r13 */ |
| 10 | /* lr .req r14 */ |
| 11 | /* pc .req r15 */ |
| wdenk | c609719 | 2002-11-03 00:24:07 +0000 | [diff] [blame] | 12 | .text |
| 13 | .globl __udivsi3 |
| 14 | .type __udivsi3 ,function |
| 15 | .align 0 |
| 16 | __udivsi3 : |
| 17 | cmp divisor, #0 |
| 18 | beq Ldiv0 |
| 19 | mov curbit, #1 |
| 20 | mov result, #0 |
| 21 | cmp dividend, divisor |
| 22 | bcc Lgot_result |
| 23 | Loop1: |
| 24 | @ Unless the divisor is very big, shift it up in multiples of |
| 25 | @ four bits, since this is the amount of unwinding in the main |
| 26 | @ division loop. Continue shifting until the divisor is |
| 27 | @ larger than the dividend. |
| 28 | cmp divisor, #0x10000000 |
| 29 | cmpcc divisor, dividend |
| 30 | movcc divisor, divisor, lsl #4 |
| 31 | movcc curbit, curbit, lsl #4 |
| 32 | bcc Loop1 |
| 33 | Lbignum: |
| 34 | @ For very big divisors, we must shift it a bit at a time, or |
| 35 | @ we will be in danger of overflowing. |
| 36 | cmp divisor, #0x80000000 |
| 37 | cmpcc divisor, dividend |
| 38 | movcc divisor, divisor, lsl #1 |
| 39 | movcc curbit, curbit, lsl #1 |
| 40 | bcc Lbignum |
| 41 | Loop3: |
| 42 | @ Test for possible subtractions, and note which bits |
| 43 | @ are done in the result. On the final pass, this may subtract |
| 44 | @ too much from the dividend, but the result will be ok, since the |
| 45 | @ "bit" will have been shifted out at the bottom. |
| 46 | cmp dividend, divisor |
| 47 | subcs dividend, dividend, divisor |
| 48 | orrcs result, result, curbit |
| 49 | cmp dividend, divisor, lsr #1 |
| 50 | subcs dividend, dividend, divisor, lsr #1 |
| 51 | orrcs result, result, curbit, lsr #1 |
| 52 | cmp dividend, divisor, lsr #2 |
| 53 | subcs dividend, dividend, divisor, lsr #2 |
| 54 | orrcs result, result, curbit, lsr #2 |
| 55 | cmp dividend, divisor, lsr #3 |
| 56 | subcs dividend, dividend, divisor, lsr #3 |
| 57 | orrcs result, result, curbit, lsr #3 |
| 58 | cmp dividend, #0 @ Early termination? |
| 59 | movnes curbit, curbit, lsr #4 @ No, any more bits to do? |
| 60 | movne divisor, divisor, lsr #4 |
| 61 | bne Loop3 |
| 62 | Lgot_result: |
| 63 | mov r0, result |
| 64 | mov pc, lr |
| 65 | Ldiv0: |
| 66 | str lr, [sp, #-4]! |
| 67 | bl __div0 (PLT) |
| 68 | mov r0, #0 @ about as wrong as it could be |
| 69 | ldmia sp!, {pc} |
| 70 | .size __udivsi3 , . - __udivsi3 |
| 71 | /* # 235 "libgcc1.S" */ |
| 72 | /* # 320 "libgcc1.S" */ |
| 73 | /* # 421 "libgcc1.S" */ |
| 74 | /* # 433 "libgcc1.S" */ |
| 75 | /* # 456 "libgcc1.S" */ |
| 76 | /* # 500 "libgcc1.S" */ |
| 77 | /* # 580 "libgcc1.S" */ |