6
#include<libipv1/ip_lib.h>
7
#include<libciomr/libciomr.h>
8
#include<libpsio/psio.h>
9
#include<libint/libint.h>
17
#include"read_scf_opdm.h"
18
#include"read_scf_occ_evec.h"
19
#include"shell_block_matrix.h"
23
#include"xc_grad_fock.h"
25
pthread_mutex_t fock_mutex; /* Lock on the global AO matrix */
36
double *Gtri, *Gtri_o; /* Total and open-shell G matrices
37
and lower triagonal form in SO basis */
39
/*----------------------------------------
40
Read in the difference HF/DFT densities
41
----------------------------------------*/
45
/*-------------------------------------------
46
Compute HF contribution to the Fock matrix
47
-------------------------------------------*/
51
/*-----------------------------------
52
Do numerical interation for KS DFT
53
-----------------------------------*/
54
if(UserOptions.make_dft){
55
/*--- Read in the SCF eigenvector density ---*/
57
/*-- Compute exch+corr contribution to the Fock matrix ---*/
58
if(UserOptions.reftype == rhf)
60
else if(UserOptions.reftype == uhf)
63
punt("\nUnrecognized Kohn-Sham DFT reference");