Back to home page

sPhenix code displayed by LXR

 
 

    


File indexing completed on 2025-08-05 08:21:09

0001  
0002 C*********************************************************************
0003  
0004 C...PYCDIV
0005 C...Auxiliary to PYCMQR
0006 C
0007 C     COMPLEX DIVISION, (CR,CI) = (AR,AI)/(BR,BI)
0008 C
0009  
0010       SUBROUTINE PYCDIV(AR,AI,BR,BI,CR,CI)
0011  
0012       DOUBLE PRECISION AR,AI,BR,BI,CR,CI
0013       DOUBLE PRECISION S,ARS,AIS,BRS,BIS
0014  
0015       S = DABS(BR) + DABS(BI)
0016       ARS = AR/S
0017       AIS = AI/S
0018       BRS = BR/S
0019       BIS = BI/S
0020       S = BRS**2 + BIS**2
0021       CR = (ARS*BRS + AIS*BIS)/S
0022       CI = (AIS*BRS - ARS*BIS)/S
0023       RETURN
0024       END