compute_curlyNpjmn.f90
73330a622fd4new_coulomb

compute_curlyNpjmn.f90
View Options

			subroutine compute_curlyNpjmn(p,j,m,n)
			! compute curlyNpjmn
			use prec_const
			use basic
			use array
			use basis_transformation
			use FMZM
			use coeff
			use T5
			!
			!
			implicit none
			!
			integer,intent(in) :: p,j,m,n
			! local vars.
			integer :: g,h,s1,icol
			CHARACTER(len=T5len) :: T5str_
			TYPE(FM),SAVE :: X1,T5_,dmnhs1_
			logical :: is_null = .false.
			!
			CALL FM_ENTER_USER_ROUTINE
			!
			DO g=0,min(Pmax,p + 2*j + m)
			X1 = SQRT(TO_FM(2)*gFACTORIAL(TO_FM(g)))
			DO h =0, j + floor(real(p + m )/2)

			! Read T4 from database
			! T5STR_ = GET_T5(p,j,m,g,h)
			! T5_ = TO_FM(T5STR_)

			is_null = T5zero( p,j,m,g,h)
			T5_ = TO_FM('0.0')
			if(.not. is_null) T5_ =Bruteforce_T5(p,j,m,g,h)

			IF(T5_ .ne. ZEROFM) THEN
			!! IF( n + m + h .le. Jmax) then
			DO s1 =0, min(Jmax,n + m + h)
			! Change order : n + m + h <=Jmax
			dmnhs1_ = ALLX2L(n,h,s1,m)
			icol = numb(g,s1)
			RowNpj(icol) = RowNpj(icol)+ X1T5_dmnhs1_
			ENDDO
			!! ENDIF
			ENDIF
			ENDDO
			ENDDO
			!
			CALL FM_EXIT_USER_ROUTINE
			!
			end subroutine compute_curlyNpjmn