forked from nwchemgit/nwchem
-
Notifications
You must be signed in to change notification settings - Fork 2
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Merge pull request nwchemgit#952 from edoapra/flaccid-fraction
updates
- Loading branch information
Showing
11 changed files
with
122 additions
and
72 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -13,6 +13,7 @@ | |
calc_x2c_hfc_scalar.o \ | ||
write_x2c_pc_dip.o \ | ||
read_x2c_pc_dip.o \ | ||
x2c_utils.o \ | ||
x2c_scf.o | ||
|
||
|
||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,57 @@ | ||
subroutine x2c_guess(rtdb,nexc,ipol, | ||
t toll_s, | ||
d do_x2c, | ||
g g_vxc,g_rel,g_dens_at, | ||
g g_x2c_u,g_x2c_u2c, | ||
N nbf_ao) | ||
implicit none | ||
#include "rtdb.fh" | ||
#include "global.fh" | ||
#include "stdio.fh" | ||
#include "errquit.fh" | ||
integer rtdb,nexc,ipol | ||
double precision toll_s | ||
logical do_x2c | ||
integer g_vxc(*) ! [input] | ||
integer g_rel(*) ! [output] | ||
integer g_dens_at(*) | ||
integer g_x2c_u(*),g_x2c_u2c,g_x2c_c2u | ||
integer nbf_ao | ||
c | ||
integer me | ||
c | ||
c ------------------------------------------------- | ||
c set up an exact decoupling Hamiltonian (optional) | ||
c ------------------------------------------------- | ||
me=ga_nodeid() | ||
if (do_x2c) then | ||
if (me.eq.0) write (luout,*) 'allocating g_x2c',nbf_ao | ||
|
||
if(.not.ga_duplicate(g_vxc(1),g_rel(1),'x2c 1')) | ||
& call errquit('dft_scf: ga_duplicate failed',666, GA_ERR) | ||
call ga_zero(g_rel(1)) | ||
|
||
if (me.eq.0) write (luout,*) 'calling calc_x2c_1e_scalar' | ||
|
||
c note: g_vxc is used as temp array in the numerical | ||
c integration routines. | ||
c g_x2c_u2c is created in the x2c routine and survives | ||
c if keepU=.T. Same for the transformation matrices g_xc_u(:). | ||
c optionally, g_x2c_c2u is also allocated and kept. | ||
|
||
call calc_x2c_1e_scalar_ga(rtdb, | ||
& nexc, g_dens_at, g_vxc, g_rel(1), g_x2c_u, | ||
& g_x2c_u2c, g_x2c_c2u, toll_s) | ||
|
||
|
||
if (ipol.gt.1) then | ||
if(.not.ga_duplicate(g_vxc(2),g_rel(2),'x2c 2')) | ||
& call errquit('dft_scf: ga_duplicate failed',1, GA_ERR) | ||
call ga_copy(g_rel(1),g_rel(2)) | ||
end if | ||
|
||
|
||
end if ! do_x2c | ||
|
||
return | ||
end |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters