Return to atof.c CVS log

Up to [Power 6/32 Unix Tahoe 4.2BSD] / cci / usr / src / usr.bin / sf77

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

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