lib/float/sbf8.c

   1 /*
   2   (c) copyright 1988 by the Vrije Universiteit, Amsterdam, The Netherlands.
   3   See the copyright notice in the ACK home directory, in the file "Copyright".
   4 */
   5
   6 /* $Header$ */
   7
   8 /*
   9         SUBTRACT TWO FLOATS - DOUBLE Precision (SBF 8)
  10 */
  11
  12 #include        "FP_types.h"
  13
  14 void
  15 sbf8(s2,s1)
  16 DOUBLE  s1,s2;
  17 {
  18         EXTEND e1, e2;
  19
  20         if (s2.d[0] == 0 && s2.d[1] == 0) {
  21                 return;
  22         }
  23         extend(&s1.d[0],&e1,sizeof(DOUBLE));
  24         extend(&s2.d[0],&e2,sizeof(DOUBLE));
  25         sub_ext(&e1,&e2);
  26         compact(&e1,&s1.d[0],sizeof(DOUBLE));
  27 }