1
#include <libdpd/dpd.h>
10
dpd_file2_init(&D1, CC_OEI, 0, 0, 1, "DAI");
11
dpd_file2_init(&D2, CC_OEI, 0, 1, 0, "D(orb)(A,I)");
12
dpd_file2_axpy(&D2, &D1, 1.0, 1);
16
dpd_file2_init(&D1, CC_OEI, 0, 0, 1, "DIA");
17
dpd_file2_init(&D2, CC_OEI, 0, 1, 0, "D(orb)(A,I)");
18
dpd_file2_axpy(&D2, &D1, 1.0, 1);
22
dpd_file2_init(&D1, CC_OEI, 0, 0, 1, "Dai");
23
dpd_file2_init(&D2, CC_OEI, 0, 1, 0, "D(orb)(a,i)");
24
dpd_file2_axpy(&D2, &D1, 1.0, 1);
28
dpd_file2_init(&D1, CC_OEI, 0, 0, 1, "Dia");
29
dpd_file2_init(&D2, CC_OEI, 0, 1, 0, "D(orb)(a,i)");
30
dpd_file2_axpy(&D2, &D1, 1.0, 1);
34
else if(params.ref == 2) {
36
dpd_file2_init(&D1, CC_OEI, 0, 0, 1, "DAI");
37
dpd_file2_init(&D2, CC_OEI, 0, 1, 0, "D(orb)(A,I)");
38
dpd_file2_axpy(&D2, &D1, 1.0, 1);
42
dpd_file2_init(&D1, CC_OEI, 0, 0, 1, "DIA");
43
dpd_file2_init(&D2, CC_OEI, 0, 1, 0, "D(orb)(A,I)");
44
dpd_file2_axpy(&D2, &D1, 1.0, 1);
48
dpd_file2_init(&D1, CC_OEI, 0, 2, 3, "Dai");
49
dpd_file2_init(&D2, CC_OEI, 0, 3, 2, "D(orb)(a,i)");
50
dpd_file2_axpy(&D2, &D1, 1.0, 1);
54
dpd_file2_init(&D1, CC_OEI, 0, 2, 3, "Dia");
55
dpd_file2_init(&D2, CC_OEI, 0, 3, 2, "D(orb)(a,i)");
56
dpd_file2_axpy(&D2, &D1, 1.0, 1);