3
#include <libdpd/dpd.h>
8
/* spinad_amps(): For RHF references, build the T2 AA and BB amplitudes from
9
** the existing T2 AB amplitudes and copy the existing T1 A amplitudes
12
** T2(IJ,AB) = T2(ij,ab) = T2(Ij,Ab) - T2(Ij,Ba)
18
void spinad_amps(void)
23
if(params.ref == 0) { /** RHF **/
25
dpd_file2_init(&T1, CC_LAMBDA, 0, 0, 1, "LIA");
26
dpd_file2_copy(&T1, CC_LAMBDA, "Lia");
29
dpd_buf4_init(&T2, CC_LAMBDA, 0, 2, 7, 0, 5, 1, "LIjAb");
30
dpd_buf4_copy(&T2, CC_LAMBDA, "LIJAB");
31
dpd_buf4_copy(&T2, CC_LAMBDA, "Lijab");