zdiv.f

Go to the documentation of this file.

      SUBROUTINE zdiv(AR, AI, BR, BI, CR, CI)
C***BEGIN PROLOGUE  ZDIV
C***REFER TO  ZBESH,ZBESI,ZBESJ,ZBESK,ZBESY,ZAIRY,ZBIRY
C
C     DOUBLE PRECISION COMPLEX DIVIDE C=A/B.
C
C***ROUTINES CALLED  XZABS
C***END PROLOGUE  ZDIV
      DOUBLE PRECISION ar, ai, br, bi, cr, ci, bm, ca, cb, cc, cd
      DOUBLE PRECISION xzabs
      bm = 1.0d0/xzabs(br,bi)
      cc = br*bm
      cd = bi*bm
      ca = (ar*cc+ai*cd)*bm
      cb = (ai*cc-ar*cd)*bm
      cr = ca
      ci = cb
      RETURN
      END