2
#include<libipv1/ip_lib.h>
4
#include<libciomr/libciomr.h>
5
#include<libint/libint.h>
12
#include"compute_scf_opdm.h"
13
#include"read_gen_opdm.h"
14
#include"enuc_deriv1.h"
16
#include"te_deriv1_scf.h"
17
#include"te_deriv1_corr.h"
18
#include"symmetrize_deriv1.h"
22
pthread_mutex_t deriv1_mutex;
27
/* Either contract integrals with the densities ... */
28
if (UserOptions.symm_ints == 0) {
29
/*--- Gradient in the canonical frame ---*/
30
Grad = block_matrix(Molecule.num_atoms,3);
32
if (Molecule.num_atoms != 0) {
33
if (!strcmp(UserOptions.wfn,"SCF")) {
41
if (!strcmp(UserOptions.wfn,"SCF"))
47
if (!strcmp(UserOptions.wfn,"SCF"))
52
chkpt_wt_grad(Grad[0]);
55
/* ... or compute integrals alone */