summaryrefslogtreecommitdiff
path: root/sysdeps/alpha/reml.S
diff options
context:
space:
mode:
Diffstat (limited to 'sysdeps/alpha/reml.S')
-rw-r--r--sysdeps/alpha/reml.S57
1 files changed, 57 insertions, 0 deletions
diff --git a/sysdeps/alpha/reml.S b/sysdeps/alpha/reml.S
new file mode 100644
index 0000000000..2ece297f09
--- /dev/null
+++ b/sysdeps/alpha/reml.S
@@ -0,0 +1,57 @@
+ /* This file is generated from divrem.m4; DO NOT EDIT! */
+/* For each N divided by D, we do:
+ result = (double) N / (double) D
+ Then, for each N mod D, we do:
+ result = N - (D * divMODE (N, D))
+
+ FIXME:
+ The q and qu versions won't deal with operands > 50 bits. We also
+ don't check for divide by zero. */
+
+#include "DEFS.h"
+#if 0
+/* We do not handle div by zero yet. */
+#include <machine/pal.h>
+#endif
+#include <sysdep.h>
+
+
+
+
+
+
+FUNC__(reml)
+ /* First set up the dividend. */
+ sextl t10, t10
+
+ stq t10,0(sp)
+ ldt $f10,0(sp)
+ cvtqt $f10,$f10
+
+
+ /* Then set up the divisor. */
+ sextl t11, t11
+
+ stq t11,0(sp)
+ ldt $f1,0(sp)
+ cvtqt $f1,$f1
+
+
+ /* Do the division. */
+ divt $f10,$f1,$f10
+ cvttqc $f10,$f10
+
+ /* Put the result in t12. */
+ stt $f10,0(sp)
+ ldq t12,0(sp)
+ sextl t12, t12
+
+
+ /* Compute the remainder. */
+ mull t11, t12, t11
+ subl t10, t11, t12
+
+
+ lda sp,16(sp)
+ ret zero,(t9),1
+ .end NAME__(reml)