sysdeps/alpha/remq.S

   1       /* This file is generated from divrem.m4; DO NOT EDIT! */
   2 /* For each N divided by D, we do:
   3       result = (double) N / (double) D
   4    Then, for each N mod D, we do:
   5       result = N - (D * divMODE (N, D))
   6
   7    FIXME:
   8    The q and qu versions won't deal with operands > 50 bits.  We also
   9    don't check for divide by zero.  */
  10
  11 #include "DEFS.h"
  12 #if 0
  13 /* We do not handle div by zero yet.  */
  14 #include <machine/pal.h>
  15 #endif
  16 #include <sysdep.h>
  17
  18 /* Avoid the definition of ret that we set in the alpha sysdep.h.  */
  19 #undef ret
  20
  21
  22
  23
  24
  25
  26 FUNC__(remq)
  27         /* First set up the dividend.  */
  28
  29         stq t10,0(sp)
  30         ldt $f10,0(sp)
  31         cvtqt $f10,$f10
  32
  33
  34         /* Then set up the divisor.  */
  35
  36         stq t11,0(sp)
  37         ldt $f1,0(sp)
  38         cvtqt $f1,$f1
  39
  40
  41         /* Do the division.  */
  42         divt $f10,$f1,$f10
  43         cvttqc $f10,$f10
  44
  45         /* Put the result in t12.  */
  46         stt $f10,0(sp)
  47         ldq t12,0(sp)
  48
  49
  50                 /* Compute the remainder.  */
  51         mulq t11, t12, t11
  52         subq t10, t11, t12
  53
  54
  55         lda sp,16(sp)
  56         ret zero,(t9),1
  57         .end NAME__(remq)