5942cf447f3080f5722c3d9d9993a39189c6e961
[kopensolaris-gnu/glibc.git] / sysdeps / alpha / divrem.m4
1 /* For each N divided by D, we do:
2       result = (double) N / (double) D
3    Then, for each N mod D, we do:
4       result = N - (D * divMODE (N, D))
5
6    FIXME:
7    The q and qu versions won't deal with operands > 50 bits.  We also
8    don't check for divide by zero.  */
9
10 #include "DEFS.h"
11 #if 0
12 /* We do not handle div by zero yet.  */
13 #include <machine/pal.h>
14 #endif
15 #include <sysdep.h>
16
17 define(path, `SYSDEP_DIR/macros.m4')dnl
18 include(path)
19
20 FUNC__(OP)
21         /* First set up the dividend.  */
22         EXTEND(t10)
23         stq t10,0(sp)
24         ldt $f10,0(sp)
25         cvtqt $f10,$f10
26         ADJQU($f10)
27
28         /* Then set up the divisor.  */
29         EXTEND(t11)
30         stq t11,0(sp)
31         ldt $f1,0(sp)
32         cvtqt $f1,$f1
33         ADJQU($f1)
34
35         /* Do the division.  */
36         divt $f10,$f1,$f10
37         cvttqc $f10,$f10
38
39         /* Put the result in t12.  */
40         stt $f10,0(sp)
41         ldq t12,0(sp)
42         FULLEXTEND(t12)
43
44         DOREM
45
46         lda sp,16(sp)
47         ret zero,(t9),1
48         .end NAME__(OP)