hto_full_ln Module Reference

Data Types
interface	operator(.fln.)

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

Function/Subroutine Documentation

fln()

real*8 function, dimension(2) hto_full_ln::fln	(	real*8, intent(in)	x,
		real*8, intent(in)	y
	)

Definition at line 470 of file CALLING_cpHTO.f.

       USE hto_riemann
       IMPLICIT NONE
       real*8 ax,ay,teta,tnteta,sr,si
       real*8, INTENT(IN) :: x,y
       real*8, dimension(2) :: fln
       IF(abs(y).eq.1.d0) THEN
        fln(1)= log(abs(x))
        IF(x > 0.d0) THEN
         fln(2)= 0.d0
        ELSE
         fln(2)= y*pi
        ENDIF 
       ELSE
        fln(1)= 0.5d0*log(x*x+y*y)
        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
         fln(2)= teta
        ELSE IF(y.eq.0.d0) THEN 
         IF(x > 0.d0) THEN
          teta= 0.d0
         ELSE
          teta= pi
         ENDIF
         fln(2)= teta
        ELSE
         tnteta= ay/ax
         teta= atan(tnteta)
         sr= x/ax
         si= y/ay
         IF(sr > 0.d0) THEN
          fln(2)= si*teta
         ELSE
          fln(2)= si*(pi-teta)
         ENDIF
        ENDIF
       ENDIF