lib/ack/libfp/mlf8.fc

   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  * Multiply Double Precision Float (MLF 8)
  10  */
  11
  12 #include        "FP_types.h"
  13
  14 void
  15 mlf8(s2,s1)
  16 DOUBLE  s1,s2;
  17 {
  18         EXTEND  e1,e2;
  19
  20         extend(&s1.d[0],&e1,sizeof(DOUBLE));
  21         extend(&s2.d[0],&e2,sizeof(DOUBLE));
  22                 /* do a multiply */
  23         mul_ext(&e1,&e2);
  24         compact(&e1,&s1.d[0],sizeof(DOUBLE));
  25 }