1
c $Id: int_giaobq10.F 21176 2011-10-10 06:35:49Z d3y133 $
1
c $Id: int_giaobq10.F 23378 2013-01-05 22:46:06Z d3y133 $
6
C> \brief Compute the kinetic and nuclear attraction 1-electron GIAO integrals
7
C> perturbed by the nuclear magnetic moments including embedding point charges
9
C> Compute the kinetic energy and nuclear attraction 1-electron GIAO integrals
10
C> perturbed by the nuclear magnetic moments.
11
C> In the nuclear attraction integrals the contributions of any embedding point
12
C> charges are included. For a definition of the GIAO basis functions see the
13
C> `int_giao_2e` routine (see also [1]).
16
C> "New integral transforms for molecular properties and application
17
C> to a massively parallel GIAO-SCF implementation",
18
C> Comp. Phys. Comm. <b>134</b>, 150-166 (2001), DOI:
19
C> <a href="http://dx.doi.org/10.1016/S0010-4655(00)00195-8">
20
C> 10.1016/S0010-4655(00)00195-8</a>.
3
22
c:tex-% this is part of the API Standard Integral routines.
4
23
c:tex-\subsection{int\_giaobq10}
5
24
c:tex-This routine computes the 1-elec GIAO integrals of the kinetic and nuclear attraction
35
54
external cando_hnd_1e_prp
37
56
c:tex-\begin{verbatim}
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 lbq10 ! [input] length of bq10 buffer
45
double precision bq10(lbq10) ! [output] bq10 integrals
46
integer nbq ! [input] number of bq centers
47
double precision qbq(nbq) ! [input] bq charges
48
double precision cbq(3*nbq) ! [input] bq coordinates
57
integer i_basis !< [Input] basis set handle for ish
58
integer ish !< [Input] i shell/contraction
59
integer j_basis !< [Input] basis set handle for jsh
60
integer jsh !< [Input] j shell/contraction
61
integer lscr !< [Input] length of scratch array
62
double precision scr(lscr) !< [Scratch] scratch array
63
integer lbq10 !< [Input] length of bq10 buffer
64
double precision bq10(lbq10) !< [Output] bq10 integrals
65
integer nbq !< [Input] number of bq centers
66
double precision qbq(nbq) !< [Input] bq charges
67
double precision cbq(3*nbq) !< [Input] bq coordinates
49
68
c:tex-\end{verbatim}
51
70
integer igeom, jgeom, ibas, jbas, ucont