summaryrefslogtreecommitdiff
path: root/libm/s_remquo.c
blob: 18853e211d0fd03a29b95d0ac1e39412914a7f9c (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
/* Copyright (C) 2002 by  Red Hat, Incorporated. All rights reserved.
 *
 * Permission to use, copy, modify, and distribute this software
 * is freely granted, provided that this notice is preserved.
 */

#include "math.h"
#include "math_private.h"

libm_hidden_proto(remquo)
#ifdef __STDC__
	double remquo(double x, double y, int *quo)	/* wrapper remquo */
#else
	double remquo(x,y,quo)			/* wrapper remquo */
	double x,y;
        int *quo;
#endif
{
        int signx, signy, signres;
        int mswx;
        int mswy;
        double x_over_y;

        GET_HIGH_WORD(mswx, x);
        GET_HIGH_WORD(mswy, y);

        signx = (mswx & 0x80000000) >> 31;
        signy = (mswy & 0x80000000) >> 31;

        signres = (signx ^ signy) ? -1 : 1;

        x_over_y = fabs(x / y);

        *quo = signres * (lrint(x_over_y) & 0x7f);

        return remainder(x,y);
}
libm_hidden_def(remquo)