daeeb5aca3
It seems that someone wanted to build flt_arith with a compiler that had long but not unsigned long. This required extra code to accomplish unsigned right shift, unsigned division, and unsigned comparison using the signed operations. Now that we use uint32_t, we can simply use the unsigned operations and remove the ucmp() function. We have similar code in mach/proto/fp/ and in lang/cem/libcc.ansi/stdlib/ext_comp.c where we use the unsigned operations. Some long variables become uint32_t, and some masks with 0xFFFFFFFF disappear because uint32_t has only 32 bits. Update flt_arith.3 to show that mantissa uses uint32_t. Provide a target to install modules/src/flt_arith/test.c as flt_test so I can run the tests.
84 lines
1.8 KiB
C
84 lines
1.8 KiB
C
/*
|
|
(c) copyright 1988 by the Vrije Universiteit, Amsterdam, The Netherlands.
|
|
See the copyright notice in the ACK home directory, in the file "Copyright".
|
|
*/
|
|
|
|
/* $Id$ */
|
|
|
|
#include "flt_misc.h"
|
|
|
|
void
|
|
flt_add(e1,e2,e3)
|
|
register flt_arith *e1,*e2,*e3;
|
|
{
|
|
/* Add two extended numbers e1 and e2, and put the result
|
|
in e3
|
|
*/
|
|
flt_arith ce1, ce2;
|
|
int diff;
|
|
|
|
flt_status = 0;
|
|
if ((e2->m1 | e2->m2) == 0L) {
|
|
*e3 = *e1;
|
|
return;
|
|
}
|
|
if ((e1->m1 | e1->m2) == 0L) {
|
|
*e3 = *e2;
|
|
return;
|
|
}
|
|
ce2 = *e2;
|
|
ce1 = *e1;
|
|
e1 = &ce1;
|
|
e2 = &ce2;
|
|
|
|
/* adjust mantissas to equal power */
|
|
diff = e2->flt_exp - e1->flt_exp;
|
|
if (diff < 0) {
|
|
diff = -diff;
|
|
e2->flt_exp += diff;
|
|
flt_b64_sft(&(e2->flt_mantissa), diff);
|
|
}
|
|
else if (diff > 0) {
|
|
e1->flt_exp += diff;
|
|
flt_b64_sft(&(e1->flt_mantissa), diff);
|
|
}
|
|
if (e1->flt_sign != e2->flt_sign) {
|
|
/* e2 + e1 = e2 - (-e1) */
|
|
if (e1->m1 > e2->m1 ||
|
|
(e1->m1 == e2->m1 && e1->m2 > e2->m2)) {
|
|
/* abs(e1) > abs(e2) */
|
|
if (e1->m2 < e2->m2) {
|
|
e1->m1 -= 1; /* carry in */
|
|
}
|
|
e1->m1 -= e2->m1;
|
|
e1->m2 -= e2->m2;
|
|
*e3 = *e1;
|
|
}
|
|
else {
|
|
if (e1->m2 > e2->m2)
|
|
e2->m1 -= 1; /* carry in */
|
|
e2->m1 -= e1->m1;
|
|
e2->m2 -= e1->m2;
|
|
*e3 = *e2;
|
|
}
|
|
}
|
|
else {
|
|
*e3 = *e2;
|
|
if (flt_b64_add(&e3->flt_mantissa,&e1->flt_mantissa)) {/* addition carry */
|
|
flt_b64_sft(&e3->flt_mantissa, 1);
|
|
e3->m1 |= 0x80000000L; /* set max bit */
|
|
e3->flt_exp++; /* increase the exponent */
|
|
}
|
|
}
|
|
flt_nrm(e3);
|
|
flt_chk(e3);
|
|
}
|
|
|
|
void
|
|
flt_sub(e1,e2,e3)
|
|
flt_arith *e1,*e2,*e3;
|
|
{
|
|
e2->flt_sign = ! e2->flt_sign;
|
|
flt_add(e1,e2,e3);
|
|
e2->flt_sign = ! e2->flt_sign;
|
|
}
|