Return to ecvt.c CVS log

Up to [Research Unix] / researchv9 / libc / gen

Annotation of researchv9/libc/gen/ecvt.c, revision 1.1.1.1

1.1       root        1: /*
                      2:  *     ecvt converts to decimal
                      3:  *     the number of digits is specified by ndigit
                      4:  *     decpt is set to the position of the decimal point
                      5:  *     sign is set to 0 for positive, 1 for negative
                      6:  */
                      7: 
                      8: char   *cvt();
                      9: 
                     10: #define        NDIG    80
                     11: char*
                     12: ecvt(arg, ndigits, decpt, sign)
                     13: double arg;
                     14: int ndigits, *decpt, *sign;
                     15: {
                     16:        return(cvt(arg, ndigits, decpt, sign, 1));
                     17: }
                     18: 
                     19: char*
                     20: fcvt(arg, ndigits, decpt, sign)
                     21: double arg;
                     22: int ndigits, *decpt, *sign;
                     23: {
                     24:        return(cvt(arg, ndigits, decpt, sign, 0));
                     25: }
                     26: 
                     27: static char*
                     28: cvt(arg, ndigits, decpt, sign, eflag)
                     29: double arg;
                     30: int ndigits, *decpt, *sign;
                     31: {
                     32:        register int r2;
                     33:        double fi, fj;
                     34:        register char *p, *p1;
                     35:        static char buf[NDIG];
                     36:        double modf();
                     37: 
                     38:        if (ndigits<0)
                     39:                ndigits = 0;
                     40:        if (ndigits>=NDIG-1)
                     41:                ndigits = NDIG-2;
                     42:        r2 = 0;
                     43:        *sign = 0;
                     44:        p = &buf[0];
                     45:        if (arg<0) {
                     46:                *sign = 1;
                     47:                arg = -arg;
                     48:        }
                     49:        arg = modf(arg, &fi);
                     50:        p1 = &buf[NDIG];
                     51:        /*
                     52:         * Do integer part
                     53:         */
                     54:        if (fi != 0) {
                     55:                p1 = &buf[NDIG];
                     56:                while (fi != 0) {
                     57:                        fj = modf(fi/10, &fi);
                     58:                        *--p1 = (int)((fj+.03)*10) + '0';
                     59:                        r2++;
                     60:                }
                     61:                while (p1 < &buf[NDIG])
                     62:                        *p++ = *p1++;
                     63:        } else if (arg > 0) {
                     64:                while ((fj = arg*10) < 1) {
                     65:                        arg = fj;
                     66:                        r2--;
                     67:                }
                     68:        }
                     69:        p1 = &buf[ndigits];
                     70:        if (eflag==0)
                     71:                p1 += r2;
                     72:        *decpt = r2;
                     73:        if (p1 < &buf[0]) {
                     74:                buf[0] = '\0';
                     75:                return(buf);
                     76:        }
                     77:        while (p<=p1 && p<&buf[NDIG]) {
                     78:                arg *= 10;
                     79:                arg = modf(arg, &fj);
                     80:                *p++ = (int)fj + '0';
                     81:        }
                     82:        if (p1 >= &buf[NDIG]) {
                     83:                buf[NDIG-1] = '\0';
                     84:                return(buf);
                     85:        }
                     86:        p = p1;
                     87:        *p1 += 5;
                     88:        while (*p1 > '9') {
                     89:                *p1 = '0';
                     90:                if (p1>buf)
                     91:                        ++*--p1;
                     92:                else {
                     93:                        *p1 = '1';
                     94:                        (*decpt)++;
                     95:                        if (eflag==0) {
                     96:                                if (p>buf)
                     97:                                        *p = '0';
                     98:                                p++;
                     99:                        }
                    100:                }
                    101:        }
                    102:        *p = '\0';
                    103:        return(buf);
                    104: }

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