sbf8.c 488 B

123456789101112131415161718192021222324252627
  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. /* $Id$ */
  6. /*
  7. SUBTRACT TWO FLOATS - DOUBLE Precision (SBF 8)
  8. */
  9. #include "FP_types.h"
  10. void
  11. sbf8(s2,s1)
  12. DOUBLE s1,s2;
  13. {
  14. EXTEND e1, e2;
  15. if (s2.d[0] == 0 && s2.d[1] == 0) {
  16. return;
  17. }
  18. extend(&s1.d[0],&e1,sizeof(DOUBLE));
  19. extend(&s2.d[0],&e2,sizeof(DOUBLE));
  20. sub_ext(&e1,&e2);
  21. compact(&e1,&s1.d[0],sizeof(DOUBLE));
  22. }