1
subroutine dft_roks_hessv( acc, g_x, g_ax )
7
c $Id: dft_roks_hessv.F 21176 2011-10-10 06:35:49Z d3y133 $
12
integer gtype,grow,gcol,growp,gcolp
16
if (crohf_init_flag.ne.1)
17
$ call errquit('dft_roks_hessv: ROKS internal block invalid',0,
19
call ga_inquire(g_x,gtype,grow,gcol)
20
if (grow.ne.crohf_vlen)
21
$ call errquit('dft_roks_hessv: invalid vector length',0,
23
call ga_inquire(g_ax,gtype,growp,gcolp)
24
if (growp.ne.crohf_vlen)
25
$ call errquit('dft_roks_hessv: invalid vector length',0,
28
$ call errquit('dft_roks_hessv: invalid no. of vectors',0,
31
c Call internal routine
33
call dft_roks_hessv_xx( basis, geom, nbf, nmo,
35
$ pflg, g_movecs, oskel,
36
$ crohf_g_fcv, crohf_g_fpv, crohf_g_fcp,
37
$ acc, lshift, g_x, g_ax )
39
c Zap numbers much smaller than acc to ensure hard zeroes
40
c remain unpolluted ... cannot use a threshold larger than the
41
c integral accuracy since can break symmetry in non-abelian groups
42
c Also must ensure that the threshold tends to zero to permit
45
call ga_screen(g_ax, min(acc*acc,acc*0.01d0,1d-12))