hto_cmplx_rootz Module Reference

Data Types
interface	operator(.crz.)

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

Function/Subroutine Documentation

crz()

real*8 function, dimension(2) hto_cmplx_rootz::crz	(	real*8, intent(in)	x,
		real*8, intent(in)	y
	)

Definition at line 686 of file CALLING_cpHTO.f.

       USE hto_units
       IMPLICIT NONE
       real*8, dimension(2) :: crz
       real*8, INTENT(IN) :: x,y
       real*8 rs,r,u,v,sax
       IF(abs(y/x).lt.1.d-8) THEN
        IF(x.gt.0.d0) THEN
         crz(1)= sqrt(x)*(1.d0+0.5d0*y*y/(x*x))
         crz(2)= 0.5d0*y/sqrt(x)
        ELSEIF(x.lt.0.d0) THEN
         sax= sqrt(abs(x))
         crz(1)= 0.5d0*abs(y)/sax
         crz(2)= sign(one,y)*sax*(1.d0+0.5d0*y*y/(x*x))
        ENDIF
       ELSE
        rs= x*x+y*y
        r= sqrt(rs)
        u= 0.5d0*(r+x)
        crz(1)= sqrt(u)
        IF(x > 0.d0.and.y.eq.0.d0) THEN
         crz(2)= 0.d0
        ELSE
         v= 0.5d0*(r-x)
         crz(2)= sign(one,y)*sqrt(v)
        ENDIF
       ENDIF