hto_a_cmplx Module Reference

Functions/Subroutines
real *8 function, dimension(2)	hto_b021_dm_cp (scal, psi, ps0i, xmsi, xm0i)

Function/Subroutine Documentation

hto_b021_dm_cp()

real*8 function, dimension(2) hto_a_cmplx::hto_b021_dm_cp	(	real*8	scal,
		real*8, dimension(2)	psi,
		real*8	ps0i,
		real*8, dimension(2,2), intent(in)	xmsi,
		real*8, dimension(2), intent(in)	xm0i
	)

Definition at line 3360 of file CALLING_cpHTO.f.

      USE hto_cmplx_root
      USE hto_cmplx_rootz
      USE hto_cmplx_srs_root
      USE hto_ln_2_riemann
      USE hto_acmplx_pro
      USE hto_acmplx_rat
      USE hto_full_ln
      USE hto_units
*
      IMPLICIT NONE
*
      INTEGER i,nci,nc
      real*8 scal,scals,ps0i,ps0,xm1,xm2
      real*8, intent(in), dimension(2) :: xm0i
      real*8, intent(in), dimension(2,2) :: xmsi
      real*8, dimension(2) :: value,xm1c,xm2c
      real*8, dimension(2,2) :: xms
      real*8, dimension(2) :: xm0
      real*8, dimension(2) :: psi,ps,aroot,root,lambdasc,lambdas,
     #        lambdac,lambda,argc,arg,llam,l1,l2 
*
      ps= psi
      ps0= ps0i
      xms= xmsi
      xm0= xm0i
*
      xm1c(1:2)= xms(1,1:2)
      xm2c(1:2)= xms(2,1:2)
      xm1c= xm1c.cq.ps
      xm2c= xm2c.cq.ps
      xm1= xm0(1)*xm0(1)/ps0
      xm2= xm0(2)*xm0(2)/ps0
*
      aroot= xm1c.cp.xm2c
      root= (aroot(1).crz.aroot(2))
*
      lambdasc= co+(xm1c.cp.xm1c)+(xm2c.cp.xm2c)-2.d0*(
     #          xm1c+xm2c+(xm1c.cp.xm2c))
      lambdas(1)= 1.d0+xm1*xm1+xm2*xm2-2.d0*(
     #            xm1+xm2+xm1*xm2)
      lambdas(2)= -eps
      lambdac= (lambdasc(1)).crz.(lambdasc(2))
      lambda= (lambdas(1)).cr.(lambdas(2))
      IF(lambda(2).eq.0.d0) lambda(2)= -eps
*
      argc= 0.5d0*((-co+xm1c+xm2c-lambdac).cq.root)
*
      arg(1)= 0.5d0*(-1.d0+xm1+xm2-lambda(1))/sqrt(xm1*xm2)
      arg(2)= eps
*
      llam= argc.lnsrs.arg
*
      l1= xm1c(1).fln.xm1c(2)
      l2= xm2c(1).fln.xm2c(2)
*
      value= (((xm1c-xm2c-co).cq.lambdac).cp.llam)+0.5d0*(l1-l2) 
*
      RETURN