Return to rhftoa.c CVS log

Up to [Research Unix] / researchv10no / cmd / PDP11 / fpp / aux

Annotation of researchv10no/cmd/PDP11/fpp/aux/rhftoa.c, revision 1.1

1.1     ! root        1: /* ftoa from portable c library, modified to use faster multiply 
        !             2:    up and down. */
        !             3: 
        !             4: ftoa (x, str, prec, format)
        !             5: float x;
        !             6: char *str;
        !             7: {
        !             8:        /* converts a floating point number to an ascii string */
        !             9:        /* x is stored into str, which should be at least 30 chars long */
        !            10:        int ie, i, k, ndig, fstyle;
        !            11:        double y;
        !            12:        ndig =  ( prec<0) ? 7 : (prec > 22 ? 23 : prec+1);
        !            13:        if  (format == 'f' || format == 'F')
        !            14:                fstyle = 1;
        !            15:        else
        !            16:                fstyle = 0;
        !            17:        /* print in e format unless last arg is 'f' */
        !            18:        ie = 0;
        !            19:        /* if x negative, write minus and reverse */
        !            20:        if ( x < 0.0)
        !            21:        {
        !            22:                *str++ = '-';
        !            23:                x = -x;
        !            24:        }
        !            25: 
        !            26:        /* put x in range 1 <= x < 10 */
        !            27:        if(x > 0.0) {
        !            28:                if(x<1e-32){
        !            29:                        x=x* 1e33; 
        !            30:                        ie -= 33;
        !            31:                }
        !            32:                if(x<1e-16){
        !            33:                        x=x* 1e17; 
        !            34:                        ie -= 17;
        !            35:                }
        !            36:                if(x<1e-8){
        !            37:                        x=x* 1e9; 
        !            38:                        ie -= 9;
        !            39:                }
        !            40:                if(x<1e-4){
        !            41:                        x=x* 1e5; 
        !            42:                        ie -= 5;
        !            43:                }
        !            44:                if(x<1e-2){
        !            45:                        x=x* 1e3; 
        !            46:                        ie -= 3;
        !            47:                }
        !            48:                if(x<1e-1){
        !            49:                        x=x* 1e2; 
        !            50:                        ie -= 2;
        !            51:                }
        !            52:                while(x<1.0){ 
        !            53:                        x=x* 10.; 
        !            54:                        ie -= 1;
        !            55:                }
        !            56:        }
        !            57:        if(x>=1e32){
        !            58:                x=x/1e32; 
        !            59:                ie += 32;
        !            60:        }
        !            61:        if(x>=1e16){
        !            62:                x=x/1e16; 
        !            63:                ie += 16;
        !            64:        }
        !            65:        if(x>=1e8){
        !            66:                x=x/1e8; 
        !            67:                ie += 8;
        !            68:        }
        !            69:        if(x>=1e4){
        !            70:                x=x/1e4; 
        !            71:                ie += 4;
        !            72:        }
        !            73:        if(x>=1e2){
        !            74:                x=x/1e2; 
        !            75:                ie += 2;
        !            76:        }
        !            77:        while(x>=10.){
        !            78:                x=x/10.;
        !            79:                ie++;
        !            80:        }
        !            81: 
        !            82:        /* in f format, number of digits is related to size */
        !            83:        if (fstyle) ndig += ie;
        !            84: 
        !            85:        /* round. x is between 1 and 10 and ndig will be printed to
        !            86:           right of decimal point so rounding is ... */
        !            87:        y = 10.0;
        !            88:        for (i = 0; i < ndig; i++)
        !            89:                y=y/10.0;
        !            90:        x=x+.5*y;
        !            91:        if (x >= 10.0) /* repair rounding disasters */
        !            92:        {
        !            93:                x = 1.0;
        !            94:                ie++;
        !            95:                if (fstyle) ndig++;
        !            96:        }
        !            97:        /* now loop.  put out a digit (obtain by multiplying by
        !            98:          10, truncating, subtracting) until enough digits out */
        !            99:        /* if fstyle, and leading zeros, they go out special */
        !           100:        if (fstyle && ie < 0)
        !           101:        {
        !           102:                *str++ = '0';
        !           103:                i = (ndig > 0) ? -1 : (ndig-1);
        !           104:                i -= ie;
        !           105:                if (i > 0)
        !           106:                {
        !           107:                        *str++ = '.';
        !           108:                        while (i--)
        !           109:                                *str++ = '0';
        !           110:                }
        !           111:        }
        !           112:        for (i=0; i < ndig; i++)
        !           113:        {
        !           114:                if (i == (fstyle ? ie+1 : 1)) /* where is decimal point */
        !           115:                        *str++ = '.';
        !           116:                k = x;
        !           117:                *str++ = k + '0';
        !           118:                x -= (y=k);
        !           119:                x=x*10.;
        !           120:        }
        !           121: 
        !           122:        /* now, in estyle,  put out exponent if not zero */
        !           123:        if (!fstyle && ie != 0)
        !           124:        {
        !           125:                *str++ = 'E';
        !           126:                if (ie < 0)
        !           127:                {
        !           128:                        ie = -ie;
        !           129:                        *str++ = '-';
        !           130:                }
        !           131:                for (k=100; k > ie; k /=10);
        !           132:                for (; k > 0; k /=10)
        !           133:                {
        !           134:                        *str++ = ie/k + '0';
        !           135:                        ie = ie%k;
        !           136:                }
        !           137:        }
        !           138:        *str = '\0';
        !           139:        return;
        !           140: }
        !           141: 
        !           142: /* printf here to force loading of pfloat is printf is called. */
        !           143: 
        !           144: printf(arg0)
        !           145: char *arg0;
        !           146: {
        !           147:        extern int putchar();
        !           148:        return(format(putchar,"%r",&arg0));
        !           149: }
        !           150: 
        !           151: __dummy() {
        !           152:        extern pfloat(); /* force pfloat and pscien to be loaded */
        !           153:        pfloat();
        !           154: }

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