Return to float.c CVS log

Up to [Research Unix] / researchv10no / cmd / lcc / etc

Annotation of researchv10no/cmd/lcc/etc/float.c, revision 1.1

1.1     ! root        1: #include <math.h>
        !             2: double ipow();
        !             3: 
        !             4: main()
        !             5: {
        !             6: #if mips || sun
        !             7:        flt(2,  -125,  128, 24, 0x34000000, 0x7f7fffff, 0x00800000);
        !             8:        dbl(2, -1021, 1024, 53, 0x3cb00000, 0x00000000, 0x7fefffff, 0xffffffff, 0x00100000, 0x00000000);
        !             9: #endif
        !            10: #ifdef vax
        !            11:        flt(2,  -127,  127, 24, 0x00003500, 0xffff7fff, 0x00000080);
        !            12:        dbl(2,  -127,  127, 56, 0x00002500, 0x00000000, 0xffff7fff, 0xffffffff, 0x00000080, 0x00000000);
        !            13: #endif
        !            14: }
        !            15: 
        !            16: dbl(b, emin, emax, p, eps, max, min)
        !            17: double eps, max, min;
        !            18: {
        !            19:        printf("#define DBL_DIG         %g\n", floor((p-1)*log10((double)b)+ispower10(b)));
        !            20:        printf("#define DBL_EPSILON     %.20e\n", eps);
        !            21:        printf("#define DBL_EPSILON     %.20e\n", ipow(b, 1-p));
        !            22:        printf("#define DBL_MANT_DIG    %d\n", p);
        !            23:        printf("#define DBL_MAX         %.20e\n", (1-ipow(b, -p))*ipow(b, emax-1)*b);
        !            24:        printf("#define DBL_MAX         %.20e\n", max);
        !            25:        printf("#define DBL_MAX_10_EXP  %g\n", floor(log10((1-ipow(b, -p))*ipow(b, emax-1)*b)));
        !            26:        printf("#define DBL_MAX_EXP     %d\n", emax);
        !            27:        printf("#define DBL_MIN         %.20e\n", ipow(b, emin-1));
        !            28:        printf("#define DBL_MIN         %.20e\n", min);
        !            29:        printf("#define DBL_MIN_10_EXP  %g\n", ceil(log10(ipow(b, emin-1))));
        !            30:        printf("#define DBL_MIN_EXP     %d\n", emin);
        !            31:        printf("\n");
        !            32: }
        !            33: 
        !            34: flt(b, emin, emax, p, eps, max, min)
        !            35: {
        !            36:        printf("#define FLT_DIG         %g\n", floor((p-1)*log10((double)b)+ispower10(b)));
        !            37:        printf("#define FLT_EPSILON     %.20e\n", *(float*)&eps);
        !            38:        printf("#define FLT_EPSILON     %.20e\n", ipow(b, 1-p));
        !            39:        printf("#define FLT_MANT_DIG    %d\n", p);
        !            40:        printf("#define FLT_MAX         %.20e\n", (1-ipow(b, -p))*ipow(b, emax-1)*b);
        !            41:        printf("#define FLT_MAX         %.20e\n", *(float*)&max);
        !            42:        printf("#define FLT_MAX_10_EXP  %g\n", floor(log10((1-ipow(b, -p))*ipow(b, emax-1)*b)));
        !            43:        printf("#define FLT_MAX_EXP     %d\n", emax);
        !            44:        printf("#define FLT_MIN         %.20e\n", *(float*)&min);
        !            45:        printf("#define FLT_MIN         %.20e\n", ipow(b, emin-1));
        !            46:        printf("#define FLT_MIN_10_EXP  %g\n", ceil(log10(ipow(b, emin-1))));
        !            47:        printf("#define FLT_MIN_EXP     %d\n", emin);
        !            48:        printf("\n");
        !            49: }
        !            50: 
        !            51: double ipow(x, n)
        !            52: {
        !            53:        double ac;
        !            54: 
        !            55:        if (n >= 0)
        !            56:                for (ac = 1; n > 0; n--)
        !            57:                        ac *= x;
        !            58:        else
        !            59:                for (ac = 1; n < 0; n++)
        !            60:                        ac /= x;
        !            61:        return ac;
        !            62: }
        !            63: 
        !            64: ispower10(x)
        !            65: {
        !            66:        double p;
        !            67: 
        !            68:        for (p = 1; p < x; p *= 10)
        !            69:                ;
        !            70:        return p == x;
        !            71: }

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