      SUBROUTINE WMUL (AR, AI, BR, BI, CR, CI)
      IMPLICIT NONE
C
C C = A*B
C
      DOUBLE PRECISION AR, AI, BR, BI, CR, CI
C
      DOUBLE PRECISION T
C
      DOUBLE PRECISION FLOP
C
C
      T = AR*BI+AI*BR
      IF (T.NE.0.0D0) T = FLOP (T)
      CR = FLOP (AR*BR-AI*BI)
      CI = T
C
      RETURN
      END
