hto_acmplx_rat Module Reference

Data Types
interface	operator(.cq.)

Functions/Subroutines
real *8 function, dimension(2)	cq (x, y)

Function/Subroutine Documentation

cq()

real*8 function, dimension(2) hto_acmplx_rat::cq	(	real*8, dimension(2), intent(in)	x,
		real*8, dimension(2), intent(in)	y
	)

Definition at line 298 of file CALLING_cpHTO.f.

       IMPLICIT NONE
       real*8, dimension(2) :: cq,z1,z2
       real*8, INTENT(IN), dimension(2) :: x,y
       real*8 theta
       IF(x(2).eq.0.d0.and.y(2).eq.0.d0) THEN
        cq(1)= x(1)/y(1)
        cq(2)= 0.d0
        RETURN
       ENDIF
       IF(abs(x(2)).eq.1.d0.and.abs(y(2)).eq.1.d0) THEN
        cq(1)= x(1)/y(1)
        cq(2)= 0.d0
       ELSE
        IF(abs(x(2)).eq.1.d0) then
         z1(1)= x(1)
         z1(2)= 0.d0
         z2(1)= y(1)
         z2(2)= y(2)         
        ELSEIF(abs(y(2)).eq.1.d0) then
         z1(1)= x(1)
         z1(2)= x(2)
         z2(1)= y(1)
         z2(2)= 0.d0         
        ELSE
         z1= x
         z2= y
        ENDIF
        IF(abs(z2(1)) > abs(z2(2))) THEN
         theta= z2(2)/z2(1)
         cq(1)= (z1(1)+theta*z1(2))/(theta*z2(2)+z2(1))
         cq(2)= (z1(2)-theta*z1(1))/(theta*z2(2)+z2(1))
        ELSE
         theta= z2(1)/z2(2)
         cq(1)= (theta*z1(1)+z1(2))/(theta*z2(1)+z2(2))
         cq(2)= (theta*z1(2)-z1(1))/(theta*z2(1)+z2(2)) 
        ENDIF     
       ENDIF