Return to atof.c CVS log

Up to [Power 6/32 Unix Tahoe 4.2BSD] / cci / usr / src / lib / libc / tahoe / gen

Annotation of cci/usr/src/lib/libc/tahoe/gen/atof.c, revision 1.1.1.1

1.1       root        1: /*     atof.s  4.2     83/06/30        */
                      2: 
                      3: /*
                      4:        C library - ascii to floating
                      5: */
                      6: 
                      7: #include <ctype.h>
                      8: #define LOGHUGE        39
                      9: 
                     10: double
                     11: atof(p)
                     12: register char *p;
                     13: {
                     14:        register c;
                     15:        double fl, flexp, exp5;
                     16:        static double big = 72057594037927936.;  /*2^56 (Dorit)*/
                     17:        double ldexp();
                     18:        int nd;
                     19:        register eexp, exp, neg, negexp, bexp;
                     20: 
                     21:        neg = 1;
                     22:        while((c = *p++) == ' ')
                     23:                ;
                     24:        if (c == '-')
                     25:                neg = -1;
                     26:        else if (c=='+')
                     27:                ;
                     28:        else
                     29:                --p;
                     30: 
                     31:        exp = 0;
                     32:        fl = 0;
                     33:        nd = 0;
                     34:        while ((c = *p++), isdigit(c)) {
                     35:                if (fl<big)
                     36:                        fl = 10*fl + (c-'0');
                     37:                else
                     38:                        exp++;
                     39:                nd++;
                     40:        }
                     41: 
                     42:        if (c == '.') {
                     43:                while ((c = *p++), isdigit(c)) {
                     44:                        if (fl<big) {
                     45:                                fl = 10*fl + (c-'0');
                     46:                                exp--;
                     47:                        }
                     48:                nd++;
                     49:                }
                     50:        }
                     51: 
                     52:        negexp = 1;
                     53:        eexp = 0;
                     54:        if ((c == 'E') || (c == 'e')) {
                     55:                if ((c= *p++) == '+')
                     56:                        ;
                     57:                else if (c=='-')
                     58:                        negexp = -1;
                     59:                else
                     60:                        --p;
                     61: 
                     62:                while ((c = *p++), isdigit(c)) {
                     63:                        eexp = 10*eexp+(c-'0');
                     64:                }
                     65:                if (negexp<0)
                     66:                        eexp = -eexp;
                     67:                exp = exp + eexp;
                     68:        }
                     69: 
                     70:        negexp = 1;
                     71:        if (exp<0) {
                     72:                negexp = -1;
                     73:                exp = -exp;
                     74:        }
                     75: 
                     76: 
                     77:        if((nd+exp*negexp) < -LOGHUGE){
                     78:                fl = 0;
                     79:                exp = 0;
                     80:        }
                     81:        flexp = 1;
                     82:        exp5 = 5;
                     83:        bexp = exp;
                     84:        for (;;) {
                     85:                if (exp&01)
                     86:                        flexp *= exp5;
                     87:                exp >>= 1;
                     88:                if (exp==0)
                     89:                        break;
                     90:                exp5 *= exp5;
                     91:        }
                     92:        if (negexp<0)
                     93:                fl /= flexp;
                     94:        else
                     95:                fl *= flexp;
                     96:        fl = ldexp(fl, negexp*bexp);
                     97:        if (neg<0)
                     98:                fl = -fl;
                     99:        return(fl);
                    100: }

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