2
c wrapper utility function to print ao integrals
3
subroutine print_integrals(basis, oskel)
5
C$Id: print_integrals.F 23464 2013-01-18 20:22:20Z niri $
10
#include "mafdecls.fh"
22
if (util_print('kinetic',print_never)) then
23
if (ga_nodeid().eq.0) then
24
write(luout,*) "===================================="
25
write(luout,*) "Begin kinetic 1-e integrals"
26
write(luout,*) "===================================="
28
call print_1eints(basis, basis, 'kinetic', oskel)
29
if (ga_nodeid().eq.0) then
30
write(luout,*) "===================================="
31
write(luout,*) "End kinetic 1-e integrals"
32
write(luout,*) "===================================="
36
if (util_print('potential',print_never)) then
37
if (ga_nodeid().eq.0) then
38
write(luout,*) "===================================="
39
write(luout,*) "Begin potential 1-e integrals"
40
write(luout,*) "===================================="
42
call print_1eints(basis, basis, 'potential', oskel)
43
if (ga_nodeid().eq.0) then
44
write(luout,*) "===================================="
45
write(luout,*) "End potential 1-e integrals"
46
write(luout,*) "===================================="
50
if (util_print('overlap',print_never)) then
51
if (ga_nodeid().eq.0) then
52
write(luout,*) "===================================="
53
write(luout,*) "Begin overlap 1-e integrals"
54
write(luout,*) "===================================="
56
call print_1eints(basis, basis, 'overlap', oskel)
57
if (ga_nodeid().eq.0) then
58
write(luout,*) "===================================="
59
write(luout,*) "End overlap 1-e integrals"
60
write(luout,*) "===================================="
64
if (util_print('ao2eints',print_never)) then
65
if (ga_nodeid().eq.0) then
66
write(luout,*) "===================================="
67
write(luout,*) "Begin all 2-e integrals"
68
write(luout,*) "===================================="
70
call print_aoints(basis)
71
if (ga_nodeid().eq.0) then
72
write(luout,*) "===================================="
73
write(luout,*) "End all 2-e integrals"
74
write(luout,*) "===================================="
78
if (util_print('allaoints',print_never)) ! Debug
79
$ call print_aoints(basis)