1
#include <libdpd/dpd.h>
7
dpdfile2 tIA, tia, fIA, fia, dIA, dia;
8
dpdbuf4 tIJAB, tijab, tIjAb, D, dIJAB, dijab, dIjAb;
10
if(params.ref == 0) { /** RHF **/
11
dpd_buf4_init(&D, CC_DINTS, 0, 0, 5, 0, 5, 0, "D <ij|ab>");
12
dpd_buf4_copy(&D, CC_TAMPS, "tIjAb");
15
dpd_buf4_init(&tIjAb, CC_TAMPS, 0, 0, 5, 0, 5, 0, "tIjAb");
16
dpd_buf4_init(&dIjAb, CC_DENOM, 0, 0, 5, 0, 5, 0, "dIjAb");
17
dpd_buf4_dirprd(&dIjAb, &tIjAb);
18
dpd_buf4_close(&dIjAb);
19
dpd_buf4_close(&tIjAb);
21
else if(params.ref == 2) { /** UHF **/
22
if(params.semicanonical) {
23
dpd_file2_init(&fIA, CC_OEI, 0, 0, 1, "fIA");
24
dpd_file2_copy(&fIA, CC_OEI, "tIA");
25
dpd_file2_close(&fIA);
27
dpd_file2_init(&fia, CC_OEI, 0, 2, 3, "fia");
28
dpd_file2_copy(&fia, CC_OEI, "tia");
29
dpd_file2_close(&fia);
31
dpd_file2_init(&tIA, CC_OEI, 0, 0, 1, "tIA");
32
dpd_file2_init(&dIA, CC_OEI, 0, 0, 1, "dIA");
33
dpd_file2_dirprd(&dIA, &tIA);
34
dpd_file2_close(&tIA);
35
dpd_file2_close(&dIA);
37
dpd_file2_init(&tia, CC_OEI, 0, 2, 3, "tia");
38
dpd_file2_init(&dia, CC_OEI, 0, 2, 3, "dia");
39
dpd_file2_dirprd(&dia, &tia);
40
dpd_file2_close(&tia);
41
dpd_file2_close(&dia);
44
dpd_buf4_init(&D, CC_DINTS, 0, 2, 7, 2, 7, 0, "D <IJ||AB> (I>J,A>B)");
45
dpd_buf4_copy(&D, CC_TAMPS, "tIJAB");
47
dpd_buf4_init(&dIJAB, CC_DENOM, 0, 1, 6, 1, 6, 0, "dIJAB");
48
dpd_buf4_init(&tIJAB, CC_TAMPS, 0, 2, 7, 2, 7, 0, "tIJAB");
49
dpd_buf4_dirprd(&dIJAB, &tIJAB);
50
dpd_buf4_close(&tIJAB);
51
dpd_buf4_close(&dIJAB);
53
dpd_buf4_init(&D, CC_DINTS, 0, 12, 17, 12, 17, 0, "D <ij||ab> (i>j,a>b)");
54
dpd_buf4_copy(&D, CC_TAMPS, "tijab");
56
dpd_buf4_init(&dIJAB, CC_DENOM, 0, 11, 16, 11, 16, 0, "dijab");
57
dpd_buf4_init(&tIJAB, CC_TAMPS, 0, 12, 17, 12, 17, 0, "tijab");
58
dpd_buf4_dirprd(&dIJAB, &tIJAB);
59
dpd_buf4_close(&tIJAB);
60
dpd_buf4_close(&dIJAB);
62
dpd_buf4_init(&D, CC_DINTS, 0, 22, 28, 22, 28, 0, "D <Ij|Ab>");
63
dpd_buf4_copy(&D, CC_TAMPS, "tIjAb");
65
dpd_buf4_init(&dIJAB, CC_DENOM, 0, 22, 28, 22, 28, 0, "dIjAb");
66
dpd_buf4_init(&tIJAB, CC_TAMPS, 0, 22, 28, 22, 28, 0, "tIjAb");
67
dpd_buf4_dirprd(&dIJAB, &tIJAB);
68
dpd_buf4_close(&tIJAB);
69
dpd_buf4_close(&dIJAB);