hto_ln_2_riemann Module Reference

Data Types
interface	operator(.lnsrs.)

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

Function/Subroutine Documentation

lnsrs()

real*8 function, dimension(2) hto_ln_2_riemann::lnsrs	(	real*8, dimension(2), intent(in)	x,
		real*8, dimension(2), intent(in)	y
	)

Definition at line 523 of file CALLING_cpHTO.f.

       USE hto_riemann
       USE hto_full_ln
       USE hto_units
       IMPLICIT NONE
       real*8, dimension(2) :: lnsrs,x,y,olnsrs
       real*8 ax,ay,teta,tnteta,ilnx,sr,si,sx,sy,xms
       INTENT(IN) x,y
*
       IF(abs(y(2)).ne.1.d0) THEN
*
        xms= x(1)*x(1)+x(2)*x(2)
        IF(abs(1.d0-sqrt(xms)).lt.1.d-12) THEN
         teta= atan(abs(x(2)/x(1)))
         sr= x(1)/abs(x(1))
         si= x(2)/abs(x(2))
         lnsrs(1)= 0.d0 
         IF(sr > 0.d0) THEN
          lnsrs(2)= si*teta
         ELSE
          lnsrs(2)= si*(pi-teta)
         ENDIF
         RETURN
        ELSE 
         lnsrs= x(1).fln.x(2)
         olnsrs= y(1).fln.y(2)
         RETURN
        ENDIF
       ENDIF 
*
       lnsrs(1)= 0.5d0*log(x(1)*x(1)+x(2)*x(2))
       ax= abs(x(1))
       ay= abs(x(2))
       IF(x(1).eq.0.d0) THEN
        IF(x(2) > 0.d0) THEN
           teta= pi/2.d0
        ELSE
           teta= -pi/2.d0
        ENDIF
        ilnx= teta
       ELSE IF(x(2).eq.0.d0) THEN 
        IF(x(1) > 0.d0) THEN
         teta= 0.d0
        ELSE
         teta= pi
        ENDIF
        ilnx= teta
       ELSE
        tnteta= ay/ax
        teta= atan(tnteta)
        sr= x(1)/ax
        si= x(2)/ay
        IF(sr > 0.d0) THEN
         ilnx= si*teta
        ELSE
         ilnx= si*(pi-teta)
        ENDIF
       ENDIF
*
       IF(x(1) > 0.d0) THEN
        IF(y(1) < 0.d0) THEN
         IF((x(2) > 0.d0.and.y(2) > 0).or.
     #      (x(2) < 0.d0.and.y(2) < 0)) THEN
          lnsrs(2)= ilnx
         ELSEIF(x(2) > 0.d0.and.y(2) < 0) THEN
          lnsrs(2)= ilnx-pi
         ELSE
          print*,'+++++++++++++++++++++++++++++++++++'
          print*,' anomaly ln '
          print*,x
          print*,y
          print*,'+++++++++++++++++++++++++++++++++++'
         ENDIF
        ELSE
         lnsrs(2)= ilnx
        ENDIF
       ELSE
        IF(abs(y(2)).ne.1.d0) THEN
         print*,' Wrong argument for LNSRS '
         stop
        ELSE
*
         IF((x(2)*y(2)) < 0.d0) THEN
          lnsrs(2)= ilnx+2.d0*y(2)*pi
         ELSE
          lnsrs(2)= ilnx
         ENDIF
*
        ENDIF
       ENDIF       
       RETURN