~vcs-imports/gawk/master

        /* die silently on EPIPE to stdout */

        if (fp == stdout && errno == EPIPE)

                gawk_exit(EXIT_FATAL);

 

        /* otherwise die verbosely */

                d = floor(d);

                d = ceil(d);

                        prec = 1;

                                if (count > 0) {

                                        prec = count;

                                        /* may need to increase fw so that padding happens, see pr_tail code */

                                        if (fw > 0)

                                                fw += count - 1;

                                }

                        else if (gawk_mb_cur_max > 1) {

                                if (cs1 == 's') {

                                        assert(cp == arg->stptr || cp == cpbuf);

                                        copy_count = mbc_byte_count(arg->stptr, prec);

                                }

                                /* prec was set by code for %c */

                                /* else

                                        copy_count = prec; */

        if (do_lint && (f0->flags & NULL_FIELD) != 0)

        double tmprand;

#define RAND_DIVISOR ((double)GAWK_RANDOM_MAX+1.0)

        /*

         * Date: Wed, 28 Aug 2013 17:52:46 -0700

         * From: Bob Jewett <jewett@bill.scs.agilent.com>

         *

         * Call random() twice to fill in more bits in the value

         * of the double.  Also, there is a bug in random() such

         * that when the values of successive values are combined

         * like (rand1*rand2)^2, (rand3*rand4)^2,  ...  the

         * resulting time series is not white noise.  The

         * following also seems to fix that bug. 

         *

         * The add/subtract 0.5 keeps small bits from filling

         * below 2^-53 in the double, not that anyone should be

         * looking down there. 

         *

         * Date: Wed, 25 Sep 2013 10:45:38 -0600 (MDT)

         * From: "Nelson H. F. Beebe" <beebe@math.utah.edu>

         * (4) The code is typical of many published fragments for converting

         *     from integer to floating-point, and I discuss the serious pitfalls

         *     in my book, because it leads to platform-dependent behavior at the

         *     end points of the interval [0,1]

         * 

         * (5) the documentation in the gawk info node says

         * 

         *     `rand()'

         *       Return a random number.  The values of `rand()' are uniformly

         *       distributed between zero and one.  The value could be zero but is

         *       never one.(1)

         * 

         *     The division by RAND_DIVISOR may not guarantee that 1.0 is never

         *     returned: the programmer forgot the platform-dependent issue of

         *     rounding.

         * 

         * For points 4 and 5, the safe way is a loop:

         * 

         *         double 

         *         rand(void)           // return value in [0.0, 1.0)

         *         {

         *          value = internal_rand();

         *

         *          while (value == 1.0) 

         *                 value = internal_rand();

         * 

         *          return (value);

         *         }

         */

 

        do {

                tmprand = 0.5 + ( (random()/RAND_DIVISOR + random())

                                        / RAND_DIVISOR);

                tmprand -= 0.5;

        } while (tmprand == 1.0);

 

        return make_number((AWKNUM) tmprand);

        if ((matches == 0 || (flags & LITERAL) != 0) && buf != NULL) {

                buf = NULL;

        }

                        assert(buf != NULL);