view libcruft/amos/zmlt.f @ 12019:e1695d65a5ca release-3-2-x

fix result class calculation in pt-mat.cc
author Jaroslav Hajek <highegg@gmail.com>
date Thu, 25 Jun 2009 06:36:24 +0200
parents 8b0cb8f79fdc
children
line wrap: on
line source

      SUBROUTINE ZMLT(AR, AI, BR, BI, CR, CI)
C***BEGIN PROLOGUE  ZMLT
C***REFER TO  ZBESH,ZBESI,ZBESJ,ZBESK,ZBESY,ZAIRY,ZBIRY
C
C     DOUBLE PRECISION COMPLEX MULTIPLY, C=A*B.
C
C***ROUTINES CALLED  (NONE)
C***END PROLOGUE  ZMLT
      DOUBLE PRECISION AR, AI, BR, BI, CR, CI, CA, CB
      CA = AR*BR - AI*BI
      CB = AR*BI + AI*BR
      CR = CA
      CI = CB
      RETURN
      END