hto_imag_lnz Module Reference

Data Types
interface	operator(.ilnz.)

Functions/Subroutines
real *8 function	ilnz (x, y)

Function/Subroutine Documentation

ilnz()

real*8 function hto_imag_lnz::ilnz	(	real*8, intent(in)	x,
		real*8, intent(in)	y
	)

Definition at line 385 of file CALLING_cpHTO.f.

      USE hto_riemann
      IMPLICIT NONE
      real*8, INTENT(IN) :: x,y
      real*8 teta,tnteta,sr,si,ax,ay,ilnz
      ax= abs(x)
      ay= abs(y)
      IF(x.eq.0.d0) THEN
       IF(y > 0.d0) THEN
          teta= pi/2.d0
       ELSE
          teta= -pi/2.d0
       ENDIF
       ilnz= teta
       RETURN
      ELSE IF(y.eq.0.d0) THEN 
       IF(x > 0.d0) THEN
          teta= 0.d0
       ELSE
          teta= pi
       ENDIF
       ilnz= teta
       RETURN
      ELSE
       tnteta= ay/ax
       teta= atan(tnteta)
       sr= x/ax
       si= y/ay
       IF(sr > 0.d0) THEN
        ilnz= si*teta
       ELSE
        ilnz= si*(pi-teta)
       ENDIF
       RETURN
      ENDIF