@brief Gauss-Hermite quadrature used in one-electron integral code @author Vladimir Mironov @date Sep, 2018 Initial release
@brief Gauss-Hermite quadrature using minimum point formula @details Compute: xint = sum( w(1:npts,npts) * (h(1:npts,npts)t+dxi)(ni-1) * (h(1:npts,npts)t+dxj)(nj-1) ) yint = sum( w(1:npts,npts) * (h(1:npts,npts)*t+dyi)(ni-1) * (h(1:npts,npts)t+dyj)(nj-1) ) zint = sum( w(1:npts,npts) * (h(1:npts,npts)t+dzi)(ni-1) * (h(1:npts,npts)*t+dzj)(nj-1) )
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=dp), | intent(out) | :: | tint(3) | |||
real(kind=dp), | intent(in) | :: | t | |||
real(kind=dp), | intent(in) | :: | rij(3) | |||
real(kind=dp), | intent(in) | :: | ri(3) | |||
real(kind=dp), | intent(in) | :: | rj(3) | |||
integer, | intent(in) | :: | ni | |||
integer, | intent(in) | :: | nj |
@brief Gauss-Hermite quadrature using minimum point formula @details Compute: xint = sum( w(1:npts,npts) * (h(1:npts,npts)t+dxi)(ni-1) * (h(1:npts,npts)t+dxj)(nj-1) ) yint = sum( w(1:npts,npts) * (h(1:npts,npts)*t+dyi)(ni-1) * (h(1:npts,npts)t+dyj)(nj-1) ) zint = sum( w(1:npts,npts) * (h(1:npts,npts)t+dzi)(ni-1) * (h(1:npts,npts)*t+dzj)(nj-1) )
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=dp), | intent(out) | :: | tint(3,0:*) | |||
real(kind=dp), | intent(in) | :: | t | |||
real(kind=dp), | intent(in) | :: | rij(3) | |||
real(kind=dp), | intent(in) | :: | ri(3) | |||
real(kind=dp), | intent(in) | :: | rj(3) | |||
real(kind=dp), | intent(in) | :: | r(3) | |||
integer, | intent(in) | :: | ni | |||
integer, | intent(in) | :: | nj | |||
integer, | intent(in) | :: | m |