public / ack.git / commitdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | patch | inline | side by side (parent: b948e4b)
avoid use of modf; not portable
authorceriel <none@none>
Fri, 22 Apr 1988 14:11:55 +0000 (14:11 +0000)
committerceriel <none@none>
Fri, 22 Apr 1988 14:11:55 +0000 (14:11 +0000)
mach/i86/cg/mach.c patch | blob | history
mach/m68020/ncg/mach.c patch | blob | history
mach/m68k2/cg/mach.c patch | blob | history
mach/m68k2/ncg/mach.c patch | blob | history
mach/m68k4/ncg/mach.c patch | blob | history
mach/moon3/ncg/mach.c patch | blob | history
mach/pdp/cg/mach.c patch | blob | history
mach/pdp/ncg/mach.c patch | blob | history

diff --git a/mach/i86/cg/mach.c b/mach/i86/cg/mach.c
index 431a870..71e5a4f 100644 (file)
--- a/mach/i86/cg/mach.c
+++ b/mach/i86/cg/mach.c
@@ -64,8 +64,7 @@ con_float() {
        double atof();
        int i;
 #ifndef OWNFLOAT
-       double f1;
-       double frexp(), modf();
+       double frexp();
        int j;
        int sign = 0;
        int fraction[4] ;
@@ -107,15 +106,18 @@ con_float() {
                f += f;
                i --;
        }
-       f = modf(2 * f, &f1); /* hidden bit */
+       f = 2*f - 1.0;          /* hidden bit */
 #ifdef IEEEFLOAT
        if (argval == 4) {
 #endif IEEEFLOAT
                i = (i + 128) & 0377;
                fraction[0] = (sign << 15) | (i << 7);
                for (j = 6; j>= 0; j--) {
-                       if (f >= 0.5) fraction[0] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               f -= 1.0;
+                               fraction[0] |= (1 << j);
+                       }
                }
 #ifdef IEEEFLOAT
        }
@@ -123,16 +125,22 @@ con_float() {
                i = (i + 1024) & 03777;
                fraction[0] = (sign << 15) | (i << 4);
                for (j = 3; j>= 0; j--) {
-                       if (f >= 0.5) fraction[0] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               fraction[0] |= (1 << j);
+                               f -= 1.0;
+                       }
                }
        }
 #endif IEEEFLOAT
        for (i = 1; i < argval / 2; i++) {
                fraction[i] = 0;
                for (j = 15; j>= 0; j--) {
-                       if (f >= 0.5) fraction[i] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               fraction[i] |= (1 << j);
+                               f -= 1.0;
+                       }
                }
        }
        if (f >= 0.5) {
@@ -159,6 +167,7 @@ con_float() {
        putc('\n', codefile);
 }
 #endif
+
 /*
 
 string holstr(n) word n; {
diff --git a/mach/m68020/ncg/mach.c b/mach/m68020/ncg/mach.c
index cf1056e..1d5aa9c 100644 (file)
--- a/mach/m68020/ncg/mach.c
+++ b/mach/m68020/ncg/mach.c
@@ -63,8 +63,7 @@ con_float() {
        double atof();
        int i;
 #ifndef OWNFLOAT
-       double f1;
-       double frexp(), modf();
+       double frexp();
        int j;
        int sign = 0;
        int fraction[4] ;
@@ -106,15 +105,18 @@ con_float() {
                f += f;
                i --;
        }
-       f = modf(2 * f, &f1); /* hidden bit */
+       f = 2*f - 1.0;          /* hidden bit */
 #ifdef IEEEFLOAT
        if (argval == 4) {
 #endif IEEEFLOAT
                i = (i + 128) & 0377;
                fraction[0] = (sign << 15) | (i << 7);
                for (j = 6; j>= 0; j--) {
-                       if (f >= 0.5) fraction[0] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               f -= 1.0;
+                               fraction[0] |= (1 << j);
+                       }
                }
 #ifdef IEEEFLOAT
        }
@@ -122,16 +124,22 @@ con_float() {
                i = (i + 1024) & 03777;
                fraction[0] = (sign << 15) | (i << 4);
                for (j = 3; j>= 0; j--) {
-                       if (f >= 0.5) fraction[0] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               fraction[0] |= (1 << j);
+                               f -= 1.0;
+                       }
                }
        }
 #endif IEEEFLOAT
        for (i = 1; i < argval / 2; i++) {
                fraction[i] = 0;
                for (j = 15; j>= 0; j--) {
-                       if (f >= 0.5) fraction[i] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               fraction[i] |= (1 << j);
+                               f -= 1.0;
+                       }
                }
        }
        if (f >= 0.5) {
diff --git a/mach/m68k2/cg/mach.c b/mach/m68k2/cg/mach.c
index 2fefd3f..9757e2a 100644 (file)
--- a/mach/m68k2/cg/mach.c
+++ b/mach/m68k2/cg/mach.c
@@ -62,8 +62,7 @@ con_float() {
        double atof();
        int i;
 #ifndef OWNFLOAT
-       double f1;
-       double frexp(), modf();
+       double frexp();
        int j;
        int sign = 0;
        int fraction[4] ;
@@ -105,15 +104,18 @@ con_float() {
                f += f;
                i --;
        }
-       f = modf(2 * f, &f1); /* hidden bit */
+       f = 2*f - 1.0;          /* hidden bit */
 #ifdef IEEEFLOAT
        if (argval == 4) {
 #endif IEEEFLOAT
                i = (i + 128) & 0377;
                fraction[0] = (sign << 15) | (i << 7);
                for (j = 6; j>= 0; j--) {
-                       if (f >= 0.5) fraction[0] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               f -= 1.0;
+                               fraction[0] |= (1 << j);
+                       }
                }
 #ifdef IEEEFLOAT
        }
@@ -121,16 +123,22 @@ con_float() {
                i = (i + 1024) & 03777;
                fraction[0] = (sign << 15) | (i << 4);
                for (j = 3; j>= 0; j--) {
-                       if (f >= 0.5) fraction[0] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               fraction[0] |= (1 << j);
+                               f -= 1.0;
+                       }
                }
        }
 #endif IEEEFLOAT
        for (i = 1; i < argval / 2; i++) {
                fraction[i] = 0;
                for (j = 15; j>= 0; j--) {
-                       if (f >= 0.5) fraction[i] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               fraction[i] |= (1 << j);
+                               f -= 1.0;
+                       }
                }
        }
        if (f >= 0.5) {
diff --git a/mach/m68k2/ncg/mach.c b/mach/m68k2/ncg/mach.c
index cf1056e..1d5aa9c 100644 (file)
--- a/mach/m68k2/ncg/mach.c
+++ b/mach/m68k2/ncg/mach.c
@@ -63,8 +63,7 @@ con_float() {
        double atof();
        int i;
 #ifndef OWNFLOAT
-       double f1;
-       double frexp(), modf();
+       double frexp();
        int j;
        int sign = 0;
        int fraction[4] ;
@@ -106,15 +105,18 @@ con_float() {
                f += f;
                i --;
        }
-       f = modf(2 * f, &f1); /* hidden bit */
+       f = 2*f - 1.0;          /* hidden bit */
 #ifdef IEEEFLOAT
        if (argval == 4) {
 #endif IEEEFLOAT
                i = (i + 128) & 0377;
                fraction[0] = (sign << 15) | (i << 7);
                for (j = 6; j>= 0; j--) {
-                       if (f >= 0.5) fraction[0] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               f -= 1.0;
+                               fraction[0] |= (1 << j);
+                       }
                }
 #ifdef IEEEFLOAT
        }
@@ -122,16 +124,22 @@ con_float() {
                i = (i + 1024) & 03777;
                fraction[0] = (sign << 15) | (i << 4);
                for (j = 3; j>= 0; j--) {
-                       if (f >= 0.5) fraction[0] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               fraction[0] |= (1 << j);
+                               f -= 1.0;
+                       }
                }
        }
 #endif IEEEFLOAT
        for (i = 1; i < argval / 2; i++) {
                fraction[i] = 0;
                for (j = 15; j>= 0; j--) {
-                       if (f >= 0.5) fraction[i] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               fraction[i] |= (1 << j);
+                               f -= 1.0;
+                       }
                }
        }
        if (f >= 0.5) {
diff --git a/mach/m68k4/ncg/mach.c b/mach/m68k4/ncg/mach.c
index cf1056e..1d5aa9c 100644 (file)
--- a/mach/m68k4/ncg/mach.c
+++ b/mach/m68k4/ncg/mach.c
@@ -63,8 +63,7 @@ con_float() {
        double atof();
        int i;
 #ifndef OWNFLOAT
-       double f1;
-       double frexp(), modf();
+       double frexp();
        int j;
        int sign = 0;
        int fraction[4] ;
@@ -106,15 +105,18 @@ con_float() {
                f += f;
                i --;
        }
-       f = modf(2 * f, &f1); /* hidden bit */
+       f = 2*f - 1.0;          /* hidden bit */
 #ifdef IEEEFLOAT
        if (argval == 4) {
 #endif IEEEFLOAT
                i = (i + 128) & 0377;
                fraction[0] = (sign << 15) | (i << 7);
                for (j = 6; j>= 0; j--) {
-                       if (f >= 0.5) fraction[0] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               f -= 1.0;
+                               fraction[0] |= (1 << j);
+                       }
                }
 #ifdef IEEEFLOAT
        }
@@ -122,16 +124,22 @@ con_float() {
                i = (i + 1024) & 03777;
                fraction[0] = (sign << 15) | (i << 4);
                for (j = 3; j>= 0; j--) {
-                       if (f >= 0.5) fraction[0] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               fraction[0] |= (1 << j);
+                               f -= 1.0;
+                       }
                }
        }
 #endif IEEEFLOAT
        for (i = 1; i < argval / 2; i++) {
                fraction[i] = 0;
                for (j = 15; j>= 0; j--) {
-                       if (f >= 0.5) fraction[i] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               fraction[i] |= (1 << j);
+                               f -= 1.0;
+                       }
                }
        }
        if (f >= 0.5) {
diff --git a/mach/moon3/ncg/mach.c b/mach/moon3/ncg/mach.c
index cf1056e..1d5aa9c 100644 (file)
--- a/mach/moon3/ncg/mach.c
+++ b/mach/moon3/ncg/mach.c
@@ -63,8 +63,7 @@ con_float() {
        double atof();
        int i;
 #ifndef OWNFLOAT
-       double f1;
-       double frexp(), modf();
+       double frexp();
        int j;
        int sign = 0;
        int fraction[4] ;
@@ -106,15 +105,18 @@ con_float() {
                f += f;
                i --;
        }
-       f = modf(2 * f, &f1); /* hidden bit */
+       f = 2*f - 1.0;          /* hidden bit */
 #ifdef IEEEFLOAT
        if (argval == 4) {
 #endif IEEEFLOAT
                i = (i + 128) & 0377;
                fraction[0] = (sign << 15) | (i << 7);
                for (j = 6; j>= 0; j--) {
-                       if (f >= 0.5) fraction[0] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               f -= 1.0;
+                               fraction[0] |= (1 << j);
+                       }
                }
 #ifdef IEEEFLOAT
        }
@@ -122,16 +124,22 @@ con_float() {
                i = (i + 1024) & 03777;
                fraction[0] = (sign << 15) | (i << 4);
                for (j = 3; j>= 0; j--) {
-                       if (f >= 0.5) fraction[0] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               fraction[0] |= (1 << j);
+                               f -= 1.0;
+                       }
                }
        }
 #endif IEEEFLOAT
        for (i = 1; i < argval / 2; i++) {
                fraction[i] = 0;
                for (j = 15; j>= 0; j--) {
-                       if (f >= 0.5) fraction[i] |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               fraction[i] |= (1 << j);
+                               f -= 1.0;
+                       }
                }
        }
        if (f >= 0.5) {
diff --git a/mach/pdp/cg/mach.c b/mach/pdp/cg/mach.c
index e53e812..65a4282 100644 (file)
--- a/mach/pdp/cg/mach.c
+++ b/mach/pdp/cg/mach.c
@@ -54,8 +54,8 @@ con_mult(sz) word sz; {
 
 con_float() {
 #ifdef ACK_ASS
-       double f, f1;
-       double atof(), frexp(), modf();
+       double f;
+       double atof(), frexp();
        int i, j;
        int sign = 0;
        int fraction ;
@@ -77,19 +77,25 @@ con_float() {
                f += f;
                i --;
        }
-       f = modf(2 * f, &f1); /* hidden bit */
+       f = 2*f - 1.0;          /* hidden bit */
        i = (i + 128) & 0377;
        fraction = (sign << 15) | (i << 7);
        for (j = 6; j>= 0; j--) {
-               if (f >= 0.5) fraction |= (1 << j);
-               f = modf(2*f, &f1);
+               f *= 2;
+               if (f >= 1.0) {
+                       fraction |= (1 << j);
+                       f -= 1.0;
+               }
        }
        fprintf(codefile, ".data2 0%o", fraction);
        for (i = argval / 2 - 1; i; i--) {
                fraction = 0;
                for (j = 15; j>= 0; j--) {
-                       if (f >= 0.5) fraction |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               fraction |= (1 << j);
+                               f -= 1.0;
+                       }
                }
                fprintf(codefile, ", 0%o", fraction);
        }
diff --git a/mach/pdp/ncg/mach.c b/mach/pdp/ncg/mach.c
index b15dbfc..65a4282 100644 (file)
--- a/mach/pdp/ncg/mach.c
+++ b/mach/pdp/ncg/mach.c
@@ -54,8 +54,8 @@ con_mult(sz) word sz; {
 
 con_float() {
 #ifdef ACK_ASS
-       double f, f1;
-       double atof(), frexp(), modf();
+       double f;
+       double atof(), frexp();
        int i, j;
        int sign = 0;
        int fraction ;
@@ -63,6 +63,11 @@ con_float() {
        if (argval != 4 && argval != 8)
                fatal("bad fcon size");
        f = atof(str);
+       if (f == 0) {
+               if (argval == 8) fprintf(codefile, ".data2 0, 0\n");
+               fprintf(codefile, ".data2 0, 0\n");
+               return;
+       }
        f = frexp(f, &i);
        if (f < 0) {
                f = -f;
@@ -72,19 +77,25 @@ con_float() {
                f += f;
                i --;
        }
-       f = modf(2 * f, &f1); /* hidden bit */
+       f = 2*f - 1.0;          /* hidden bit */
        i = (i + 128) & 0377;
        fraction = (sign << 15) | (i << 7);
        for (j = 6; j>= 0; j--) {
-               if (f >= 0.5) fraction |= (1 << j);
-               f = modf(2*f, &f1);
+               f *= 2;
+               if (f >= 1.0) {
+                       fraction |= (1 << j);
+                       f -= 1.0;
+               }
        }
        fprintf(codefile, ".data2 0%o", fraction);
        for (i = argval / 2 - 1; i; i--) {
                fraction = 0;
                for (j = 15; j>= 0; j--) {
-                       if (f >= 0.5) fraction |= (1 << j);
-                       f = modf(2*f, &f1);
+                       f *= 2;
+                       if (f >= 1.0) {
+                               fraction |= (1 << j);
+                               f -= 1.0;
+                       }
                }
                fprintf(codefile, ", 0%o", fraction);
        }
Unnamed repository; edit this file 'description' to name the repository.