DEHam/src/nt2.irp.f

use iso_c_binding
BEGIN_PROVIDER [integer(C_SIZE_T), nt2]
    BEGIN_DOC
    ! calculates the number of det the 1's moving
    END_DOC


!   call combin(idet2(1,nt2+1),natrest,ial0,nt2,32,jrangmax)
    nt2=   nint(gamma(real(natom-ntrou+1,16))/((gamma(real(nalpha+1,16))*gamma(real(nbeta+1,16)))),selected_int_kind(16))
    print *,"nt2=",nt2
END_PROVIDER
looks like now it works for natom > 16 2018-01-29 00:10:47 +01:00			`use iso_c_binding`
			`BEGIN_PROVIDER [integer(C_SIZE_T), nt2]`
initial commit 2016-12-19 23:26:16 +01:00			`BEGIN_DOC`
			`! calculates the number of det the 1's moving`
			`END_DOC`


			`! call combin(idet2(1,nt2+1),natrest,ial0,nt2,32,jrangmax)`
			`nt2= nint(gamma(real(natom-ntrou+1,16))/((gamma(real(nalpha+1,16))*gamma(real(nbeta+1,16)))),selected_int_kind(16))`
			`print *,"nt2=",nt2`
			`END_PROVIDER`