3
\brief Enter brief description of file here
7
#include <libdpd/dpd.h>
13
namespace psi { namespace ccenergy {
15
void local_filter_T1(dpdfile2 *T1);
18
/* apply denominators to t1 and t2 */
22
dpdfile2 newtIA, dIA, tIA, newtia, dia, tia;
24
if (params.ref == 0) {
25
dpd_file2_init(&newtIA, CC_OEI, 0, 0, 1, "New tIA");
26
dpd_file2_copy(&newtIA, CC_OEI, "New tIA Increment");
27
dpd_file2_close(&newtIA);
29
dpd_file2_init(&newtIA, CC_OEI, 0, 0, 1, "New tIA Increment");
30
if(params.local && local.filter_singles) {
31
local_filter_T1(&newtIA);
34
dpd_file2_init(&dIA, CC_OEI, 0, 0, 1, "dIA");
35
dpd_file2_dirprd(&dIA, &newtIA);
36
dpd_file2_close(&dIA);
38
dpd_file2_close(&newtIA);
40
/* Add the new increment to the old tIA to get the New tIA */
41
dpd_file2_init(&tIA, CC_OEI, 0, 0, 1, "tIA");
42
dpd_file2_copy(&tIA, CC_OEI, "New tIA");
43
dpd_file2_close(&tIA);
44
dpd_file2_init(&newtIA, CC_OEI, 0, 0, 1, "New tIA");
45
dpd_file2_init(&tIA, CC_OEI, 0, 0, 1, "New tIA Increment");
46
dpd_file2_axpy(&tIA, &newtIA, 1, 0);
47
dpd_file2_close(&tIA);
48
dpd_file2_close(&newtIA);
50
else if (params.ref == 1) {
51
dpd_file2_init(&newtIA, CC_OEI, 0, 0, 1, "New tIA");
52
dpd_file2_init(&dIA, CC_OEI, 0, 0, 1, "dIA");
53
dpd_file2_dirprd(&dIA, &newtIA);
54
dpd_file2_close(&dIA);
55
dpd_file2_close(&newtIA);
57
dpd_file2_init(&newtia, CC_OEI, 0, 0, 1, "New tia");
58
dpd_file2_init(&dia, CC_OEI, 0, 0, 1, "dia");
59
dpd_file2_dirprd(&dia, &newtia);
60
dpd_file2_close(&dia);
61
dpd_file2_close(&newtia);
63
else if (params.ref == 2) {
64
dpd_file2_init(&newtIA, CC_OEI, 0, 0, 1, "New tIA");
65
dpd_file2_init(&dIA, CC_OEI, 0, 0, 1, "dIA");
66
dpd_file2_dirprd(&dIA, &newtIA);
67
dpd_file2_close(&dIA);
68
dpd_file2_close(&newtIA);
70
dpd_file2_init(&newtia, CC_OEI, 0, 2, 3, "New tia");
71
dpd_file2_init(&dia, CC_OEI, 0, 2, 3, "dia");
72
dpd_file2_dirprd(&dia, &newtia);
73
dpd_file2_close(&dia);
74
dpd_file2_close(&newtia);
81
}} // namespace psi::ccenergy