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
|
|
|
/*
|
1993-01-05 12:06:58 +00:00
|
|
|
* Multiply Double Precision Float (MLF 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
|
|
|
mlf8(s2,s1)
|
1993-01-05 12:06:58 +00:00
|
|
|
DOUBLE s1,s2;
|
1988-04-07 10:57:49 +00:00
|
|
|
{
|
|
|
|
|
EXTEND e1,e2;
|
|
|
|
|
|
1993-01-05 12:06:58 +00:00
|
|
|
extend(&s1.d[0],&e1,sizeof(DOUBLE));
|
|
|
|
|
extend(&s2.d[0],&e2,sizeof(DOUBLE));
|
1988-04-07 10:57:49 +00:00
|
|
|
/* do a multiply */
|
|
|
|
|
mul_ext(&e1,&e2);
|
1993-01-05 12:06:58 +00:00
|
|
|
compact(&e1,&s1.d[0],sizeof(DOUBLE));
|
1988-04-07 10:57:49 +00:00
|
|
|
}
|
