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".
|
|
|
|
*/
|
|
|
|
|
|
|
|
/* $Header$ */
|
|
|
|
|
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-01-05 12:06:58 +00:00
|
|
|
DOUBLE
|
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-01-05 12:06:58 +00:00
|
|
|
DOUBLE *result = &s1; /* s1 may not be in a register! */
|
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) {
|
1988-07-22 20:54:49 +00:00
|
|
|
return s1;
|
|
|
|
}
|
1989-10-25 17:15:37 +00:00
|
|
|
s2 = ngf8(s2);
|
1988-08-02 11:57:37 +00:00
|
|
|
*result = adf8(s2,s1); /* add and return result */
|
1988-04-07 10:57:49 +00:00
|
|
|
return(s1);
|
|
|
|
}
|