cleaned up a bit, and made add_ext and sub_ext proper entry points,
without preconditions on the arguments
This commit is contained in:
@@ -27,27 +27,7 @@ _float s1,s2;
|
||||
}
|
||||
extend((_double *)&s1,&e1,sizeof(SINGLE));
|
||||
extend((_double *)&s2,&e2,sizeof(SINGLE));
|
||||
/* if signs differ do subtraction */
|
||||
/* result in e1 */
|
||||
if (e1.sign ^ e2.sign) {
|
||||
/* set sign of e1 to sign of largest */
|
||||
swap = (e2.exp > e1.exp) ? 1 : (e2.exp < e1.exp) ? 0 :
|
||||
(e2.m1 > e1.m1 ) ? 1 : 0;
|
||||
/* adjust mantissas to equal powers */
|
||||
sft_ext(&e1,&e2);
|
||||
/* subtract the extended formats */
|
||||
if (swap) {
|
||||
sub_ext(&e2,&e1);
|
||||
e1 = e2;
|
||||
}
|
||||
else
|
||||
sub_ext(&e1,&e2);
|
||||
}
|
||||
else {
|
||||
/* adjust mantissas to equal powers */
|
||||
sft_ext(&e1,&e2);
|
||||
add_ext(&e1,&e2);
|
||||
}
|
||||
add_ext(&e1,&e2);
|
||||
compact(&e1,(_double *)&s1,sizeof(SINGLE));
|
||||
return(s1);
|
||||
}
|
||||
|
||||
Reference in New Issue
Block a user