1988-04-07 11:40:46 +00:00
|
|
|
/*
|
|
|
|
(c) copyright 1988 by the Vrije Universiteit, Amsterdam, The Netherlands.
|
|
|
|
See the copyright notice in the ACK home directory, in the file "Copyright".
|
|
|
|
*/
|
|
|
|
|
1994-06-24 14:02:31 +00:00
|
|
|
/* $Id$ */
|
1988-04-07 11:40:46 +00:00
|
|
|
|
1988-04-07 10:57:49 +00:00
|
|
|
/*
|
1988-07-25 10:46:15 +00:00
|
|
|
SUBTRACT TWO FLOATS - DOUBLE Precision (SBF 8)
|
1988-04-07 10:57:49 +00:00
|
|
|
*/
|
|
|
|
|
|
|
|
#include "FP_types.h"
|
|
|
|
|
1993-11-24 15:37:27 +00:00
|
|
|
void
|
1988-04-07 10:57:49 +00:00
|
|
|
sbf8(s2,s1)
|
1993-01-05 12:06:58 +00:00
|
|
|
DOUBLE s1,s2;
|
1988-04-07 10:57:49 +00:00
|
|
|
{
|
1993-11-24 15:37:27 +00:00
|
|
|
EXTEND e1, e2;
|
1988-04-07 10:57:49 +00:00
|
|
|
|
1993-01-05 12:06:58 +00:00
|
|
|
if (s2.d[0] == 0 && s2.d[1] == 0) {
|
1993-11-24 15:37:27 +00:00
|
|
|
return;
|
1988-07-22 20:54:49 +00:00
|
|
|
}
|
1993-11-24 15:37:27 +00:00
|
|
|
extend(&s1.d[0],&e1,sizeof(DOUBLE));
|
|
|
|
extend(&s2.d[0],&e2,sizeof(DOUBLE));
|
|
|
|
sub_ext(&e1,&e2);
|
|
|
|
compact(&e1,&s1.d[0],sizeof(DOUBLE));
|
1988-04-07 10:57:49 +00:00
|
|
|
}
|