irpjast/jastrow_provider.irp.f

16 lines
373 B
Fortran

BEGIN_PROVIDER [ double precision, jastrow_full ]
implicit none
BEGIN_DOC
! Complete jastrow factor
END_DOC
integer :: i, j
print *, "J_ee = ", factor_ee
print *, "J_en = ", factor_en
print *, "J_een = ", factor_een
print *, "J = J_ee + J_en + J_een = ", factor_ee + factor_en + factor_een
jastrow_full = dexp(factor_ee + factor_en + factor_een)
END_PROVIDER