Return to modf.c CVS log

Up to [MW Coherent from dump] / coherent / b / lib / libc / crt / i386

Annotation of coherent/b/lib/libc/crt/i386/modf.c, revision 1.1

1.1     ! root        1: /*
        !             2:  * libc/crt/modf.c
        !             3:  */
        !             4: 
        !             5: /*
        !             6:  * double
        !             7:  * modf(real, dp)
        !             8:  * double real, *dp;
        !             9:  *
        !            10:  * modf() finds g = greatest integer <= real,
        !            11:  * then stores g through dp and returns real - g.
        !            12:  * This version uses a machine dependent subroutine
        !            13:  *     double _modf(real, dp, e)
        !            14:  * which assumes real >= 0, 1 <= e = exponent(real) <= DMBITS+1,
        !            15:  * where DMBITS is the number of mantissa bits in a double.
        !            16:  */
        !            17: 
        !            18: #if    _IEEE
        !            19: #define        DEBITS  11
        !            20: #define        DMBITS  52
        !            21: #define        FEBITS   8
        !            22: #define        FMBITS  23
        !            23: #endif
        !            24: 
        !            25: #if    _DECVAX
        !            26: #define        DEBITS   8
        !            27: #define        DMBITS  55
        !            28: #define        FEBITS   8
        !            29: #define        FMBITS  23
        !            30: #endif
        !            31: 
        !            32: extern double  _modf();
        !            33: extern double  frexp();
        !            34: 
        !            35: double
        !            36: modf(d, dp) double d; register double *dp;
        !            37: {
        !            38:        int e;
        !            39: 
        !            40:        frexp(d, &e);
        !            41:        if (e >= DMBITS+1) {                    /* very large d */
        !            42:                *dp = d;
        !            43:                return 0.0;
        !            44:        } else if (e <= 0) {                    /* -1. < d < 1. */
        !            45:                *dp = 0.0;
        !            46:                return d;
        !            47:        } else if (d >= 0.0)
        !            48:                return _modf(d, dp, e);         /* d >= 0., 1 <= e < DMBITS+1 */
        !            49: 
        !            50:        /* d is negative. */
        !            51:        d = _modf(-d, dp, e);
        !            52:        *dp = -*dp;
        !            53:        return -d;
        !            54: }
        !            55: 
        !            56: /* end of modf.c */

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.