pycdiv.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pycdiv.f

 
C*********************************************************************
 
C...PYCDIV
C...Auxiliary to PYCMQR
C
C     COMPLEX DIVISION, (CR,CI) = (AR,AI)/(BR,BI)
C
 
      SUBROUTINE pycdiv(AR,AI,BR,BI,CR,CI)
 
      DOUBLE PRECISION ar,ai,br,bi,cr,ci
      DOUBLE PRECISION s,ars,ais,brs,bis
 
      s = dabs(br) + dabs(bi)
      ars = ar/s
      ais = ai/s
      brs = br/s
      bis = bi/s
      s = brs**2 + bis**2
      cr = (ars*brs + ais*bis)/s
      ci = (ais*brs - ars*bis)/s
      RETURN
      END