hto_olas Module Reference

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

Function/Subroutine Documentation

hto_b0af_em()

real*8 function, dimension(2) hto_olas::hto_b0af_em	(	real*8	scal,
		real*8, dimension(2)	psi,
		real*8	ps0i,
		real*8, dimension(2)	xmsi,
		real*8	xm0i
	)

Definition at line 2800 of file CALLING_cpHTO.f.

      USE hto_acmplx_pro
      USE hto_acmplx_rat
      USE hto_cmplx_root
      USE hto_cmplx_rootz
      USE hto_ln_2_riemann
      USE hto_full_ln
      USE hto_units
*
      IMPLICIT NONE
*
      real*8 scal,ps0i,xm0i,scals,ps0,xm0
      real*8, dimension(2) :: value,psi,ps,xmsi,xms,betasc,betas,
     #        betac,beta,argc,arg,lbet
*
      scals= scal*scal
      ps= psi/scals
      ps0= ps0i/scals
      xms= xmsi/scals
      xm0= xm0i/scal
*
      betasc= co-4.d0*(xms.cq.ps)
      IF(betasc(2).eq.0.d0) THEN
       betasc(2)= -eps
       betac= (betasc(1)).cr.(betasc(2))
      ELSE
       betac= (betasc(1)).crz.(betasc(2))
      ENDIF 
      argc= (betac+co).cq.(betac-co)
      IF(argc(2).eq.0.d0) THEN
       argc(2)= eps
       lbet= argc(1).fln.argc(2)
      ELSE
       betas(1)= 1.d0-4.d0*(xm0*xm0)/ps0
       betas(2)= -eps
       beta= (betas(1)).cr.(betas(2))
       arg= (beta+co).cq.(beta-co)
       IF(arg(2).eq.0.d0) arg(2)= eps
       lbet= argc.lnsrs.arg
      ENDIF
*
      value= 2.d0*co-(betac.cp.lbet)
*
      RETURN
*