toolchain/gcc/newlib/libm/common/s_remquo.c

   1 /* Copyright (C) 2002 by  Red Hat, Incorporated. All rights reserved.
   2  *
   3  * Permission to use, copy, modify, and distribute this software
   4  * is freely granted, provided that this notice is preserved.
   5  */
   6
   7 #include "fdlibm.h"
   8
   9 #ifndef _DOUBLE_IS_32BITS
  10
  11 #ifdef __STDC__
  12         double remquo(double x, double y, int *quo)     /* wrapper remquo */
  13 #else
  14         double remquo(x,y,quo)                  /* wrapper remquo */
  15         double x,y;
  16         int *quo;
  17 #endif
  18 {
  19         int signx, signy, signres;
  20         int mswx;
  21         int mswy;
  22         double x_over_y;
  23
  24         GET_HIGH_WORD(mswx, x);
  25         GET_HIGH_WORD(mswy, y);
  26
  27         signx = (mswx & 0x80000000) >> 31;
  28         signy = (mswy & 0x80000000) >> 31;
  29
  30         signres = (signx ^ signy) ? -1 : 1;
  31
  32         x_over_y = fabs(x / y);
  33
  34         *quo = signres * (lrint(x_over_y) & 0x7f);
  35
  36         return remainder(x,y);
  37 }
  38
  39 #endif /* defined(_DOUBLE_IS_32BITS) */