DEHam/src/nt2.irp.f

14 lines
379 B
FortranFixed
Raw Permalink Normal View History

2018-01-29 00:10:47 +01:00
use iso_c_binding
BEGIN_PROVIDER [integer(C_SIZE_T), nt2]
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(1.0*(natom-ntrou+1))/((gamma(1.0*(nalpha+1))*gamma(1.0*(nbeta+1)))),selected_int_kind(16))
if(mpiid==0)then
print *,"nt2=",nt2
endif
2016-12-19 23:26:16 +01:00
END_PROVIDER