add Berkeley copyright; NOTE: this code is AT&T derived!

[unix-history] / usr / src / lib / libm / common_source / floor.c

/*
* Copyright (c) 1985 Regents of the University of California.
* All rights reserved.  The Berkeley software License Agreement
* specifies the terms and conditions for redistribution.
*/

#ifndef lint
static char sccsid[] = "@(#)floor.c     5.2 (Berkeley) %G%";
#endif /* not lint */

/*
* floor and ceil-- greatest integer <= arg
* (resp least >=)
*/

double  modf();

double
floor(d)
double d;
{
       double fract;

       if (d<0.0) {
               d = -d;
               fract = modf(d, &d);
               if (fract != 0.0)
                       d += 1;
               d = -d;
       } else
               modf(d, &d);
       return(d);
}

double
ceil(d)
double d;
{
       return(-floor(-d));
}

#ifndef national                        /* rint() is in ./NATIONAL/support.s */
/*
* algorithm for rint(x) in pseudo-pascal form ...
*
* real rint(x): real x;
*      ... delivers integer nearest x in direction of prevailing rounding
*      ... mode
* const        L = (last consecutive integer)/2
*        = 2**55; for VAX D
*        = 2**52; for IEEE 754 Double
* real s,t;
* begin
*      if x != x then return x;                ... NaN
*      if |x| >= L then return x;              ... already an integer
*      s := copysign(L,x);
*      t := x + s;                             ... = (x+s) rounded to integer
*      return t - s
* end;
*
* Note: Inexact will be signaled if x is not an integer, as is
*      customary for IEEE 754.  No other signal can be emitted.
*/
#if defined(vax)||defined(tahoe)
#ifdef vax
#define _0x(A,B)        0x/**/A/**/B
#else   /* vax */
#define _0x(A,B)        0x/**/B/**/A
#endif  /* vax */
static long Lx[] = {_0x(0000,5c00),_0x(0000,0000)};     /* 2**55 */
#define L *(double *) Lx
#else   /* defined(vax)||defined(tahoe) */
static double L = 4503599627370496.0E0;         /* 2**52 */
#endif  /* defined(vax)||defined(tahoe) */
double
rint(x)
double x;
{
       double s,t,one = 1.0,copysign();
#if !defined(vax)&&!defined(tahoe)
       if (x != x)                             /* NaN */
               return (x);
#endif  /* !defined(vax)&&!defined(tahoe) */
       if (copysign(x,one) >= L)               /* already an integer */
           return (x);
       s = copysign(L,x);
       t = x + s;                              /* x+s rounded to integer */
       return (t - s);
}
#endif  /* not national */