1
c $Id: int_vstat1e.F 19696 2010-10-29 16:53:42Z d3y133 $
4
C> \brief Compute the nuclear attaction integrals for pseudo nucleii
6
C> This subroutine is modified from int_1epe. It evaluates the nuclear
7
C> attraction integrals for each pseudo nucleus at positions specified by the
8
C> array qxyz, which are just the quardrature points. The effective charges
9
C> of these pseudo nuclei are determined by the transition densities and/or
10
C> the quardrature weights.
12
c $Id: int_vstat1e.F 23398 2013-01-08 08:18:46Z d3y133 $
2
13
subroutine int_vstat1e(i_basis,ish,j_basis,jsh,lscr,scr,nqpts,
3
14
& vint,qxyz,zq,dens)
35
46
external cando_nw_1e
38
integer i_basis ! [input] basis set handle for ish
39
integer ish ! [input] i shell/contraction
40
integer j_basis ! [input] basis set handle for jsh
41
integer jsh ! [input] j shell/contraction
42
integer lscr ! [input] length of scratch array
43
double precision scr(lscr) ! [scratch] scratch array
44
integer nqpts ! [input] length of potential buffer
45
double precision vint(nqpts) ! [output] potential integrals
46
double precision qxyz(3,nqpts) ! [input] coordinates of q points
47
double precision zq(nqpts) ! [input] effective charges on q points
48
double precision dens(*) ! [input] elements of density matrix
49
integer i_basis !< [Input] basis set handle for ish
50
integer ish !< [Input] i shell/contraction
51
integer j_basis !< [Input] basis set handle for jsh
52
integer jsh !< [Input] j shell/contraction
53
integer lscr !< [Input] length of scratch array
54
double precision scr(lscr) !< [Scratch] scratch array
55
integer nqpts !< [Input] length of potential buffer
56
double precision vint(nqpts) !< [Output] potential integrals
57
double precision qxyz(3,nqpts) !< [Input] coordinates of q points
58
double precision zq(nqpts) !< [Input] effective charges on q points
59
double precision dens(*) !< [Input] elements of density matrix
50
61
logical shells_ok, orel, oirel, ojrel, oNR
51
62
integer i_geom, j_geom, ibas, jbas, ucont