summaryrefslogtreecommitdiff
path: root/test/math/etodec.c
diff options
context:
space:
mode:
authorEric Andersen <andersen@codepoet.org>2002-04-03 10:26:12 +0000
committerEric Andersen <andersen@codepoet.org>2002-04-03 10:26:12 +0000
commit915950ede281243b2c7a5400980ef16681cf3ab4 (patch)
treefb5ee2cd0ee875b4251cc441fb5f3c4ccae3bcd5 /test/math/etodec.c
parente53310b756c8e0e02ed41737dc3573cad33bc083 (diff)
run dos2unix on these files
Diffstat (limited to 'test/math/etodec.c')
-rw-r--r--test/math/etodec.c362
1 files changed, 181 insertions, 181 deletions
diff --git a/test/math/etodec.c b/test/math/etodec.c
index 29fee34cd..22545d6fb 100644
--- a/test/math/etodec.c
+++ b/test/math/etodec.c
@@ -1,181 +1,181 @@
-#include "ehead.h"
-void emovi(), emovo(), ecleaz(), eshdn8(), emdnorm();
-void todec();
-/*
-; convert DEC double precision to e type
-; double d;
-; short e[NE];
-; dectoe( &d, e );
-*/
-void dectoe( d, e )
-unsigned short *d;
-unsigned short *e;
-{
-unsigned short y[NI];
-register unsigned short r, *p;
-
-ecleaz(y); /* start with a zero */
-p = y; /* point to our number */
-r = *d; /* get DEC exponent word */
-if( *d & (unsigned int )0x8000 )
- *p = 0xffff; /* fill in our sign */
-++p; /* bump pointer to our exponent word */
-r &= 0x7fff; /* strip the sign bit */
-if( r == 0 ) /* answer = 0 if high order DEC word = 0 */
- goto done;
-
-
-r >>= 7; /* shift exponent word down 7 bits */
-r += EXONE - 0201; /* subtract DEC exponent offset */
- /* add our e type exponent offset */
-*p++ = r; /* to form our exponent */
-
-r = *d++; /* now do the high order mantissa */
-r &= 0177; /* strip off the DEC exponent and sign bits */
-r |= 0200; /* the DEC understood high order mantissa bit */
-*p++ = r; /* put result in our high guard word */
-
-*p++ = *d++; /* fill in the rest of our mantissa */
-*p++ = *d++;
-*p = *d;
-
-eshdn8(y); /* shift our mantissa down 8 bits */
-done:
-emovo( y, e );
-}
-
-
-
-/*
-; convert e type to DEC double precision
-; double d;
-; short e[NE];
-; etodec( e, &d );
-*/
-#if 0
-static unsigned short decbit[NI] = {0,0,0,0,0,0,0200,0};
-void etodec( x, d )
-unsigned short *x, *d;
-{
-unsigned short xi[NI];
-register unsigned short r;
-int i, j;
-
-emovi( x, xi );
-*d = 0;
-if( xi[0] != 0 )
- *d = 0100000;
-r = xi[E];
-if( r < (EXONE - 128) )
- goto zout;
-i = xi[M+4];
-if( (i & 0200) != 0 )
- {
- if( (i & 0377) == 0200 )
- {
- if( (i & 0400) != 0 )
- {
- /* check all less significant bits */
- for( j=M+5; j<NI; j++ )
- {
- if( xi[j] != 0 )
- goto yesrnd;
- }
- }
- goto nornd;
- }
-yesrnd:
- eaddm( decbit, xi );
- r -= enormlz(xi);
- }
-
-nornd:
-
-r -= EXONE;
-r += 0201;
-if( r < 0 )
- {
-zout:
- *d++ = 0;
- *d++ = 0;
- *d++ = 0;
- *d++ = 0;
- return;
- }
-if( r >= 0377 )
- {
- *d++ = 077777;
- *d++ = -1;
- *d++ = -1;
- *d++ = -1;
- return;
- }
-r &= 0377;
-r <<= 7;
-eshup8( xi );
-xi[M] &= 0177;
-r |= xi[M];
-*d++ |= r;
-*d++ = xi[M+1];
-*d++ = xi[M+2];
-*d++ = xi[M+3];
-}
-#else
-
-extern int rndprc;
-
-void etodec( x, d )
-unsigned short *x, *d;
-{
-unsigned short xi[NI];
-long exp;
-int rndsav;
-
-emovi( x, xi );
-exp = (long )xi[E] - (EXONE - 0201); /* adjust exponent for offsets */
-/* round off to nearest or even */
-rndsav = rndprc;
-rndprc = 56;
-emdnorm( xi, 0, 0, exp, 64 );
-rndprc = rndsav;
-todec( xi, d );
-}
-
-void todec( x, y )
-unsigned short *x, *y;
-{
-unsigned short i;
-unsigned short *p;
-
-p = x;
-*y = 0;
-if( *p++ )
- *y = 0100000;
-i = *p++;
-if( i == 0 )
- {
- *y++ = 0;
- *y++ = 0;
- *y++ = 0;
- *y++ = 0;
- return;
- }
-if( i > 0377 )
- {
- *y++ |= 077777;
- *y++ = 0xffff;
- *y++ = 0xffff;
- *y++ = 0xffff;
- return;
- }
-i &= 0377;
-i <<= 7;
-eshup8( x );
-x[M] &= 0177;
-i |= x[M];
-*y++ |= i;
-*y++ = x[M+1];
-*y++ = x[M+2];
-*y++ = x[M+3];
-}
-#endif
+#include "ehead.h"
+void emovi(), emovo(), ecleaz(), eshdn8(), emdnorm();
+void todec();
+/*
+; convert DEC double precision to e type
+; double d;
+; short e[NE];
+; dectoe( &d, e );
+*/
+void dectoe( d, e )
+unsigned short *d;
+unsigned short *e;
+{
+unsigned short y[NI];
+register unsigned short r, *p;
+
+ecleaz(y); /* start with a zero */
+p = y; /* point to our number */
+r = *d; /* get DEC exponent word */
+if( *d & (unsigned int )0x8000 )
+ *p = 0xffff; /* fill in our sign */
+++p; /* bump pointer to our exponent word */
+r &= 0x7fff; /* strip the sign bit */
+if( r == 0 ) /* answer = 0 if high order DEC word = 0 */
+ goto done;
+
+
+r >>= 7; /* shift exponent word down 7 bits */
+r += EXONE - 0201; /* subtract DEC exponent offset */
+ /* add our e type exponent offset */
+*p++ = r; /* to form our exponent */
+
+r = *d++; /* now do the high order mantissa */
+r &= 0177; /* strip off the DEC exponent and sign bits */
+r |= 0200; /* the DEC understood high order mantissa bit */
+*p++ = r; /* put result in our high guard word */
+
+*p++ = *d++; /* fill in the rest of our mantissa */
+*p++ = *d++;
+*p = *d;
+
+eshdn8(y); /* shift our mantissa down 8 bits */
+done:
+emovo( y, e );
+}
+
+
+
+/*
+; convert e type to DEC double precision
+; double d;
+; short e[NE];
+; etodec( e, &d );
+*/
+#if 0
+static unsigned short decbit[NI] = {0,0,0,0,0,0,0200,0};
+void etodec( x, d )
+unsigned short *x, *d;
+{
+unsigned short xi[NI];
+register unsigned short r;
+int i, j;
+
+emovi( x, xi );
+*d = 0;
+if( xi[0] != 0 )
+ *d = 0100000;
+r = xi[E];
+if( r < (EXONE - 128) )
+ goto zout;
+i = xi[M+4];
+if( (i & 0200) != 0 )
+ {
+ if( (i & 0377) == 0200 )
+ {
+ if( (i & 0400) != 0 )
+ {
+ /* check all less significant bits */
+ for( j=M+5; j<NI; j++ )
+ {
+ if( xi[j] != 0 )
+ goto yesrnd;
+ }
+ }
+ goto nornd;
+ }
+yesrnd:
+ eaddm( decbit, xi );
+ r -= enormlz(xi);
+ }
+
+nornd:
+
+r -= EXONE;
+r += 0201;
+if( r < 0 )
+ {
+zout:
+ *d++ = 0;
+ *d++ = 0;
+ *d++ = 0;
+ *d++ = 0;
+ return;
+ }
+if( r >= 0377 )
+ {
+ *d++ = 077777;
+ *d++ = -1;
+ *d++ = -1;
+ *d++ = -1;
+ return;
+ }
+r &= 0377;
+r <<= 7;
+eshup8( xi );
+xi[M] &= 0177;
+r |= xi[M];
+*d++ |= r;
+*d++ = xi[M+1];
+*d++ = xi[M+2];
+*d++ = xi[M+3];
+}
+#else
+
+extern int rndprc;
+
+void etodec( x, d )
+unsigned short *x, *d;
+{
+unsigned short xi[NI];
+long exp;
+int rndsav;
+
+emovi( x, xi );
+exp = (long )xi[E] - (EXONE - 0201); /* adjust exponent for offsets */
+/* round off to nearest or even */
+rndsav = rndprc;
+rndprc = 56;
+emdnorm( xi, 0, 0, exp, 64 );
+rndprc = rndsav;
+todec( xi, d );
+}
+
+void todec( x, y )
+unsigned short *x, *y;
+{
+unsigned short i;
+unsigned short *p;
+
+p = x;
+*y = 0;
+if( *p++ )
+ *y = 0100000;
+i = *p++;
+if( i == 0 )
+ {
+ *y++ = 0;
+ *y++ = 0;
+ *y++ = 0;
+ *y++ = 0;
+ return;
+ }
+if( i > 0377 )
+ {
+ *y++ |= 077777;
+ *y++ = 0xffff;
+ *y++ = 0xffff;
+ *y++ = 0xffff;
+ return;
+ }
+i &= 0377;
+i <<= 7;
+eshup8( x );
+x[M] &= 0177;
+i |= x[M];
+*y++ |= i;
+*y++ = x[M+1];
+*y++ = x[M+2];
+*y++ = x[M+3];
+}
+#endif