NAME
        IEEEDPMul -- multiply one double precision IEEE number by another

   SYNOPSIS
          x   = IEEEDPMul(  y  ,  z  );
        d0/d1             d0/d1 d2/d3

        double  x,y,z;

   FUNCTION
        Compute x = y * z in IEEE double precision.

   INPUTS
        y -- IEEE double precision floating point value
        z -- IEEE double precision floating point value

   RESULT
        x -- IEEE double precision floating point value

   BUGS

   SEE ALSO
        ieeedpdiv()