compute_pipjm.f90
View Options

			subroutine compute_pipjm(p,j,m,XOUT)
			! compute \Pi^{pjm}(b_a)
			use prec_const
			use basic
			use array
			use basis_transformation
			use FMZM
			use coeff
			!
			!
			implicit none
			!
			integer,intent(in) :: p,j,m
			TYPE(FM), intent(inout) :: XOUT

			! local vars.
			integer :: h,s,s1
			CHARACTER(len=T5len) :: T5str_
			TYPE(FM),SAVE :: T5,dmshs1_,kerns_,bapow_,kerns1_
			!
			CALL FM_ENTER_USER_ROUTINE
			!
			XOUT = TO_FM('0')
			!
			bapow_ = TO_FM(bafm_/2)**m
			!
			DO h =0, j + floor(real(p + m )/2)
			T5STR_ = GET_T5(p,j,m,0,h)
			T5 = TO_FM(T5STR_)
			IF(T5 .ne. TO_FM('0')) THEN
			DO s=0,NFLR
			kerns_ = kerneln(bafm_,s)
			DO s1 =0,s + m + h
			kerns1_ = kerneln(bafm_,s1)
			dmshs1_ = ALLX2L(s,h,s1,m)
			XOUT = XOUT + T5dmshs1_FACTORIAL(TO_FM(s))/FACTORIAL(TO_FM(s + m))kerns_kerns1_*bapow_
			ENDDO
			ENDDO
			ENDIF
			ENDDO
			!
			CALL FM_EXIT_USER_ROUTINE
			!
			end subroutine compute_pipjm

cosolver/gkcoulomb_suite/Pipjm/src/compute_pipjm.f90138972b69f9dnew_coulomb