			program main
			! test T4T5 library
			use fmzm
			use t4t5

			implicit none
			integer :: lmax, kmax, pmax, jmax ,mmax
			integer :: l, k, p, j, m
			type(fm) :: NX0, NX1
			character(100) :: str1, str2
			CALL FM_ENTER_USER_ROUTINE

			call init_T4T5(.true., 10 )
			lmax = 4
			kmax = 3
			pmax = 4
			jmax = 3
			mmax = 2
			print*, 'test T4/T5 coefficients '

			l = lmax
			k = kmax
			p = pmax
			m = mmax
			j= jmax
			! do l =0, lmax
			! do k=0, kmax
			! do p = 0, pmax
			! do m=0, p
			! do j =0, jmax

			! NX0 = get_T5(l, k, m , p ,j ) ! T_{lkm}^{pj}
			! CALL FMFORM('F60.40',NX0 ,str1)
			! print*,l,k, m , p, j
			! print*, 'T5 =', trim(adjustl(str1))

			! NX1 = get_T5inv(l, k, p, j, m ) ! (T^{-1})_{lk}^{pjm}
			NX1 = get_T5inv_prec(l, k, p, j, m ) ! (T^{-1})_{lk}^{pjm}
			CALL FMFORM('F60.40',NX1 ,str2)
			print*, l,k, p,j,m
			print*, 'T5inv =', trim(adjustl(str2))


			! enddo
			! enddo
			! enddo
			! enddo
			! enddo
			CALL FM_EXIT_USER_ROUTINE
			end program main

cosolver/basis_transformation/T4T5/main.f90
aa929621a7cbmaster

main.f90
View Options