Return to frexp.c CVS log

Up to [CSRG BSD Unix] / 3BSD / libc / gen

Annotation of 3BSD/libc/gen/frexp.c, revision 1.1.1.1

1.1       root        1: /*
                      2:        the call
                      3:                x = frexp(arg,&exp);
                      4:        must return a double fp quantity x which is <1.0
                      5:        and the corresponding binary exponent "exp".
                      6:        such that
                      7:                arg = x*2^exp
                      8: */
                      9: 
                     10: double
                     11: frexp(x,i)
                     12: double x;
                     13: int *i;
                     14: {
                     15:        int neg;
                     16:        int j;
                     17:        j = 0;
                     18:        neg = 0;
                     19:        if(x<0){
                     20:                x = -x;
                     21:                neg = 1;
                     22:                }
                     23:        if(x>1.0)
                     24:                while(x>1){
                     25:                        j = j+1;
                     26:                        x = x/2;
                     27:                        }
                     28:        else if(x<0.5)
                     29:                while(x<0.5){
                     30:                        j = j-1;
                     31:                        x = 2*x;
                     32:                        }
                     33:        *i = j;
                     34:        if(neg) x = -x;
                     35:        return(x);
                     36:        }

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