← Back to branch summary

~ubuntu-branches/ubuntu/quantal/psicode/quantal

« back to all changes in this revision

Viewing changes to src/bin/ccdensity/x_te_intermediates.cc

  • Committer: Bazaar Package Importer
  • Author(s): Michael Banck, Michael Banck, Daniel Leidert
  • Date: 2009-02-23 00:12:02 UTC
  • mfrom: (1.1.2 upstream)
  • Revision ID: james.westby@ubuntu.com-20090223001202-rutldoy3dimfpesc
Tags: 3.4.0-1
* New upstream release.

[ Michael Banck ]
* debian/patches/01_DESTDIR.dpatch: Refreshed.
* debian/patches/02_FHS.dpatch: Removed, applied upstream.
* debian/patches/03_debian_docdir: Likewise.
* debian/patches/04_man.dpatch: Likewise.
* debian/patches/06_466828_fix_gcc_43_ftbfs.dpatch: Likewise.
* debian/patches/07_464867_move_executables: Fixed and refreshed.
* debian/patches/00list: Adjusted.
* debian/control: Improved description.
* debian/patches-held: Removed.
* debian/rules (install/psi3): Do not ship the ruby bindings for now.

[ Daniel Leidert ]
* debian/rules: Fix txtdir via DEB_MAKE_INSTALL_TARGET.
* debian/patches/01_DESTDIR.dpatch: Refreshed.

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/bin/ccdensity/x_te_intermediates.cc
 
1
/*! \file
 
2
    \ingroup CCDENSITY
 
3
    \brief Enter brief description of file here 
 
4
*/
 
5
#include <cstdio>
 
6
#include <libdpd/dpd.h>
 
7
#include "MOInfo.h"
 
8
#include "Params.h"
 
9
#include "Frozen.h"
 
10
#define EXTERN
 
11
#include "globals.h"
 
12
 
 
13
namespace psi { namespace ccdensity {
 
14
/* replace R2L2_oOoO and tau2l2_oOoO with sort later */
 
15
/* just removed cleaning  - shouldn't need this */
 
16
 
 
17
extern void c_clean_CIJAB(dpdbuf4 *CMNEF);
 
18
 
 
19
void x_te_intermediates(void)
 
20
{
 
21
  dpdfile2 R1, L1;
 
22
  dpdbuf4 V, L, R, T2;
 
23
  int G_irr, L_irr, R_irr;
 
24
  G_irr = params.G_irr; L_irr = params.L_irr; R_irr = params.R_irr;
 
25
 
 
26
  if(params.ref == 0 || params.ref == 1) { /** RHF/ROHF **/
 
27
 
 
28
    /* R2L2_OOOO = 0.5 Rijef Lklef */
 
29
    dpd_buf4_init(&V, EOM_TMP, G_irr, 2, 2, 2, 2, 0, "R2L2_OOOO");
 
30
    dpd_buf4_init(&R, CC_GR, R_irr, 2, 7, 2, 7, 0, "RIJAB");
 
31
    dpd_buf4_init(&L, CC_GL, L_irr, 2, 7, 2, 7, 0, "LIJAB");
 
32
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 0.0);
 
33
    dpd_buf4_close(&L);
 
34
    dpd_buf4_close(&R);
 
35
    dpd_buf4_close(&V);
 
36
    dpd_buf4_init(&V, EOM_TMP, G_irr, 2, 2, 2, 2, 0, "R2L2_oooo");
 
37
    dpd_buf4_init(&R, CC_GR, R_irr, 2, 7, 2, 7, 0, "Rijab");
 
38
    dpd_buf4_init(&L, CC_GL, L_irr, 2, 7, 2, 7, 0, "Lijab");
 
39
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 0.0);
 
40
    dpd_buf4_close(&L);
 
41
    dpd_buf4_close(&R);
 
42
    dpd_buf4_close(&V);
 
43
    dpd_buf4_init(&V, EOM_TMP, 0, 0, 0, 0, 0, 0, "R2L2_OoOo");
 
44
    dpd_buf4_init(&R, CC_GR, R_irr, 0, 5, 0, 5, 0, "RIjAb");
 
45
    dpd_buf4_init(&L, CC_GL, L_irr, 0, 5, 0, 5, 0, "LIjAb");
 
46
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 0.0);
 
47
    dpd_buf4_close(&L);
 
48
    dpd_buf4_close(&R);
 
49
    dpd_buf4_sort(&V, EOM_TMP, qpsr, 0, 0, "R2L2_oOoO");
 
50
    dpd_buf4_close(&V);
 
51
 
 
52
    /* Tau2L2_OOOO = 0.5 Tau_ijef Lklef */
 
53
    dpd_buf4_init(&V, EOM_TMP, L_irr, 2, 2, 2, 2, 0, "Tau2L2_OOOO");
 
54
    dpd_buf4_init(&T2, CC_TAMPS, 0, 2, 7, 2, 7, 0, "tauIJAB");
 
55
    dpd_buf4_init(&L, CC_GL, L_irr, 2, 7, 2, 7, 0, "LIJAB");
 
56
    dpd_contract444(&T2, &L, &V, 0, 0, 1.0, 0.0);
 
57
    dpd_buf4_close(&L);
 
58
    dpd_buf4_close(&T2);
 
59
    dpd_buf4_close(&V);
 
60
    dpd_buf4_init(&V, EOM_TMP, L_irr, 2, 2, 2, 2, 0, "Tau2L2_oooo");
 
61
    dpd_buf4_init(&T2, CC_TAMPS, 0, 2, 7, 2, 7, 0, "tauijab");
 
62
    dpd_buf4_init(&L, CC_GL, L_irr, 2, 7, 2, 7, 0, "Lijab");
 
63
    dpd_contract444(&T2, &L, &V, 0, 0, 1.0, 0.0);
 
64
    dpd_buf4_close(&L);
 
65
    dpd_buf4_close(&T2);
 
66
    dpd_buf4_close(&V);
 
67
    dpd_buf4_init(&V, EOM_TMP, L_irr, 0, 0, 0, 0, 0, "Tau2L2_OoOo");
 
68
    dpd_buf4_init(&T2, CC_TAMPS, 0, 0, 5, 0, 5, 0, "tauIjAb");
 
69
    dpd_buf4_init(&L, CC_GL, L_irr, 0, 5, 0, 5, 0, "LIjAb");
 
70
    dpd_contract444(&T2, &L, &V, 0, 0, 1.0, 0.0);
 
71
    dpd_buf4_close(&L);
 
72
    dpd_buf4_close(&T2);
 
73
    dpd_buf4_sort(&V, EOM_TMP, qpsr, 0, 0, "Tau2L2_oOoO");
 
74
    dpd_buf4_close(&V);
 
75
 
 
76
    /* R2L2_OVOV = Rimae Ljmbe */
 
77
    dpd_buf4_init(&V, EOM_TMP, G_irr, 10, 10, 10, 10, 0, "R2L2_OVOV");
 
78
    dpd_buf4_init(&R, CC_GR, R_irr, 10, 10, 10, 10, 0, "RIAJB");
 
79
    dpd_buf4_init(&L, CC_GL, L_irr, 10, 10, 10, 10, 0, "LIAJB");
 
80
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 0.0);
 
81
    dpd_buf4_close(&L);
 
82
    dpd_buf4_close(&R);
 
83
    dpd_buf4_init(&R, CC_GR, R_irr, 10, 10, 10, 10, 0, "RIAjb");
 
84
    dpd_buf4_init(&L, CC_GL, L_irr, 10, 10, 10, 10, 0, "LIAjb");
 
85
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 1.0);
 
86
    dpd_buf4_close(&L);
 
87
    dpd_buf4_close(&R);
 
88
    dpd_buf4_close(&V);
 
89
 
 
90
    dpd_buf4_init(&V, EOM_TMP, G_irr, 10, 10, 10, 10, 0, "R2L2_ovov");
 
91
    dpd_buf4_init(&R, CC_GR, R_irr, 10, 10, 10, 10, 0, "Riajb");
 
92
    dpd_buf4_init(&L, CC_GL, L_irr, 10, 10, 10, 10, 0, "Liajb");
 
93
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 0.0);
 
94
    dpd_buf4_close(&L);
 
95
    dpd_buf4_close(&R);
 
96
    dpd_buf4_init(&R, CC_GR, R_irr, 10, 10, 10, 10, 0, "RIAjb");
 
97
    dpd_buf4_init(&L, CC_GL, L_irr, 10, 10, 10, 10, 0, "LIAjb");
 
98
    dpd_contract444(&R, &L, &V, 1, 1, 1.0, 1.0);
 
99
    dpd_buf4_close(&L);
 
100
    dpd_buf4_close(&R);
 
101
    dpd_buf4_close(&V);
 
102
 
 
103
    dpd_buf4_init(&V, EOM_TMP, G_irr, 10, 10, 10, 10, 0, "R2L2_OVov");
 
104
    dpd_buf4_init(&R, CC_GR, R_irr, 10, 10, 10, 10, 0, "RIAjb");
 
105
    dpd_buf4_init(&L, CC_GL, L_irr, 10, 10, 10, 10, 0, "Liajb");
 
106
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 0.0);
 
107
    dpd_buf4_close(&L);
 
108
    dpd_buf4_close(&R);
 
109
    dpd_buf4_init(&R, CC_GR, R_irr, 10, 10, 10, 10, 0, "RIAJB");
 
110
    dpd_buf4_init(&L, CC_GL, L_irr, 10, 10, 10, 10, 0, "LIAjb");
 
111
    dpd_contract444(&R, &L, &V, 0, 1, 1.0, 1.0);
 
112
    dpd_buf4_close(&L);
 
113
    dpd_buf4_close(&R);
 
114
    dpd_buf4_close(&V);
 
115
 
 
116
    dpd_buf4_init(&V, EOM_TMP, G_irr, 10, 10, 10, 10, 0, "R2L2_ovOV");
 
117
    dpd_buf4_init(&R, CC_GR, R_irr, 10, 10, 10, 10, 0, "RiaJB");
 
118
    dpd_buf4_init(&L, CC_GL, L_irr, 10, 10, 10, 10, 0, "LIAJB");
 
119
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 0.0);
 
120
    dpd_buf4_close(&L);
 
121
    dpd_buf4_close(&R);
 
122
    dpd_buf4_init(&R, CC_GR, R_irr, 10, 10, 10, 10, 0, "Riajb");
 
123
    dpd_buf4_init(&L, CC_GL, L_irr, 10, 10, 10, 10, 0, "LIAjb");
 
124
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 1.0);
 
125
    dpd_buf4_close(&L);
 
126
    dpd_buf4_close(&R);
 
127
    dpd_buf4_close(&V);
 
128
 
 
129
    dpd_buf4_init(&V, EOM_TMP, G_irr, 10, 10, 10, 10, 0, "R2L2_oVoV");
 
130
    dpd_buf4_init(&R, CC_GR, R_irr, 10, 10, 10, 10, 0, "RjAIb");
 
131
    dpd_buf4_init(&L, CC_GL, L_irr, 10, 10, 10, 10, 0, "LIbjA");
 
132
    dpd_contract444(&R, &L, &V, 0, 1, 1.0, 0.0);
 
133
    dpd_buf4_close(&L);
 
134
    dpd_buf4_close(&R);
 
135
    dpd_buf4_close(&V);
 
136
 
 
137
    dpd_buf4_init(&V, EOM_TMP, G_irr, 10, 10, 10, 10, 0, "R2L2_OvOv");
 
138
    dpd_buf4_init(&R, CC_GR, R_irr, 10, 10, 10, 10, 0, "RIbjA");
 
139
    dpd_buf4_init(&L, CC_GL, L_irr, 10, 10, 10, 10, 0, "LjAIb");
 
140
    dpd_contract444(&R, &L, &V, 0, 1, 1.0, 0.0);
 
141
    dpd_buf4_close(&L);
 
142
    dpd_buf4_close(&R);
 
143
    dpd_buf4_close(&V);
 
144
 
 
145
    /* L2R1_OOVO = Lijae Rke */
 
146
    dpd_buf4_init(&V, EOM_TMP, G_irr, 2, 11, 2, 11, 0, "L2R1_OOVO");
 
147
    dpd_buf4_init(&L, CC_GL, L_irr, 2, 5, 2, 7, 0, "LIJAB");
 
148
    dpd_file2_init(&R1, CC_GR, R_irr, 0, 1, "RIA");
 
149
    dpd_contract424(&L, &R1, &V, 3, 1, 0, 1.0, 0.0);
 
150
    dpd_file2_close(&R1);
 
151
    dpd_buf4_close(&L);
 
152
    dpd_buf4_close(&V);
 
153
 
 
154
    dpd_buf4_init(&V, EOM_TMP, G_irr, 2, 11, 2, 11, 0, "L2R1_oovo");
 
155
    dpd_buf4_init(&L, CC_GL, L_irr, 2, 5, 2, 7, 0, "Lijab");
 
156
    dpd_file2_init(&R1, CC_GR, R_irr, 0, 1, "Ria");
 
157
    dpd_contract424(&L, &R1, &V, 3, 1, 0, 1.0, 0.0);
 
158
    dpd_file2_close(&R1);
 
159
    dpd_buf4_close(&L);
 
160
    dpd_buf4_close(&V);
 
161
 
 
162
    dpd_buf4_init(&V, EOM_TMP, G_irr, 0, 11, 0, 11, 0, "L2R1_OoVo");
 
163
    dpd_buf4_init(&L, CC_GL, L_irr, 0, 5, 0, 5, 0, "LIjAb");
 
164
    dpd_file2_init(&R1, CC_GR, R_irr, 0, 1, "Ria");
 
165
    dpd_contract424(&L, &R1, &V, 3, 1, 0, 1.0, 0.0);
 
166
    dpd_file2_close(&R1);
 
167
    dpd_buf4_close(&L);
 
168
    dpd_buf4_close(&V);
 
169
 
 
170
    dpd_buf4_init(&V, EOM_TMP, G_irr, 0, 11, 0, 11, 0, "L2R1_OovO");
 
171
    dpd_buf4_init(&L, CC_GL, L_irr, 0, 5, 0, 5, 0, "LIjaB");
 
172
    dpd_file2_init(&R1, CC_GR, R_irr, 0, 1, "RIA");
 
173
    dpd_contract424(&L, &R1, &V, 3, 1, 0, 1.0, 0.0);
 
174
    dpd_file2_close(&R1);
 
175
    dpd_buf4_close(&L);
 
176
    dpd_buf4_close(&V);
 
177
 
 
178
    dpd_buf4_init(&V, EOM_TMP, G_irr, 2, 11, 2, 11, 0, "L2R1_OOVO");
 
179
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 2, 10, "L2R1_OOVO(pqsr)");
 
180
    dpd_buf4_close(&V);
 
181
    dpd_buf4_init(&V, EOM_TMP, G_irr, 2, 11, 2, 11, 0, "L2R1_oovo");
 
182
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 2, 10, "L2R1_oovo(pqsr)");
 
183
    dpd_buf4_close(&V);
 
184
 
 
185
    dpd_buf4_init(&V, EOM_TMP, G_irr, 0, 11, 0, 11, 0, "L2R1_OoVo");
 
186
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 0, 10, "L2R1_OoVo(pqsr)");
 
187
    dpd_buf4_sort(&V, EOM_TMP, qprs, 0, 11, "L2R1_OoVo(qprs)");
 
188
    dpd_buf4_close(&V);
 
189
    dpd_buf4_init(&V, EOM_TMP, G_irr, 0, 10, 0, 10, 0, "L2R1_OoVo(pqsr)");
 
190
    dpd_buf4_sort(&V, EOM_TMP, qprs, 0, 10, "L2R1_OoVo(qpsr)");
 
191
    dpd_buf4_close(&V);
 
192
 
 
193
    dpd_buf4_init(&V, EOM_TMP, G_irr, 0, 11, 0, 11, 0, "L2R1_OovO");
 
194
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 0, 10, "L2R1_OovO(pqsr)");
 
195
    dpd_buf4_sort(&V, EOM_TMP, qprs, 0, 11, "L2R1_OovO(qprs)");
 
196
    dpd_buf4_close(&V);
 
197
 
 
198
    /* L1R2_OOVO = Rijae Lke */
 
199
    dpd_buf4_init(&V, EOM_TMP, G_irr, 2, 11, 2, 11, 0, "L1R2_OOVO");
 
200
    dpd_buf4_init(&R, CC_GR, R_irr, 2, 5, 2, 7, 0, "RIJAB");
 
201
    dpd_file2_init(&L1, CC_GL, L_irr, 0, 1, "LIA");
 
202
    dpd_contract424(&R, &L1, &V, 3, 1, 0, 1.0, 0.0);
 
203
    dpd_file2_close(&L1);
 
204
    dpd_buf4_close(&R);
 
205
    dpd_buf4_close(&V);
 
206
 
 
207
    dpd_buf4_init(&V, EOM_TMP, G_irr, 2, 11, 2, 11, 0, "L1R2_oovo");
 
208
    dpd_buf4_init(&R, CC_GR, R_irr, 2, 5, 2, 7, 0, "Rijab");
 
209
    dpd_file2_init(&L1, CC_GL, L_irr, 0, 1, "Lia");
 
210
    dpd_contract424(&R, &L1, &V, 3, 1, 0, 1.0, 0.0);
 
211
    dpd_file2_close(&L1);
 
212
    dpd_buf4_close(&R);
 
213
    dpd_buf4_close(&V);
 
214
 
 
215
    dpd_buf4_init(&V, EOM_TMP, G_irr, 0, 11, 0, 11, 0, "L1R2_OoVo");
 
216
    dpd_buf4_init(&R, CC_GR, R_irr, 0, 5, 0, 5, 0, "RIjAb");
 
217
    dpd_file2_init(&L1, CC_GL, L_irr, 0, 1, "Lia");
 
218
    dpd_contract424(&R, &L1, &V, 3, 1, 0, 1.0, 0.0);
 
219
    dpd_file2_close(&L1);
 
220
    dpd_buf4_close(&R);
 
221
    dpd_buf4_close(&V);
 
222
 
 
223
    dpd_buf4_init(&V, EOM_TMP, G_irr, 0, 11, 0, 11, 0, "L1R2_OovO");
 
224
    dpd_buf4_init(&R, CC_GR, R_irr, 0, 5, 0, 5, 0, "RIjaB");
 
225
    dpd_file2_init(&L1, CC_GL, L_irr, 0, 1, "LIA");
 
226
    dpd_contract424(&R, &L1, &V, 3, 1, 0, 1.0, 0.0);
 
227
    dpd_file2_close(&L1);
 
228
    dpd_buf4_close(&R);
 
229
    dpd_buf4_close(&V);
 
230
 
 
231
    dpd_buf4_init(&V, EOM_TMP, G_irr, 2, 11, 2, 11, 0, "L1R2_OOVO");
 
232
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 2, 10, "L1R2_OOVO(pqsr)");
 
233
    dpd_buf4_close(&V);
 
234
    dpd_buf4_init(&V, EOM_TMP, G_irr, 2, 11, 2, 11, 0, "L1R2_oovo");
 
235
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 2, 10, "L1R2_oovo(pqsr)");
 
236
    dpd_buf4_close(&V);
 
237
    dpd_buf4_init(&V, EOM_TMP, G_irr, 0, 11, 0, 11, 0, "L1R2_OoVo");
 
238
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 0, 10, "L1R2_OoVo(pqsr)");
 
239
    dpd_buf4_close(&V);
 
240
    dpd_buf4_init(&V, EOM_TMP, G_irr, 0, 10, 0, 10, 0, "L1R2_OoVo(pqsr)");
 
241
    dpd_buf4_sort(&V, EOM_TMP, qprs, 0, 10, "L1R2_OoVo(qpsr)");
 
242
    dpd_buf4_close(&V);
 
243
    dpd_buf4_init(&V, EOM_TMP, G_irr, 0, 11, 0, 11, 0, "L1R2_OovO");
 
244
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 0, 10, "L1R2_OovO(pqsr)");
 
245
    dpd_buf4_close(&V);
 
246
 
 
247
    /* L2R1_VVOV = Limab Rmc */
 
248
    dpd_buf4_init(&V, EOM_TMP, G_irr, 7, 10, 7, 10, 0, "L2R1_VVOV");
 
249
    dpd_buf4_init(&L, CC_GL, L_irr, 0, 7, 2, 7, 0, "LIJAB");
 
250
    dpd_file2_init(&R1, CC_GR, R_irr, 0, 1, "RIA");
 
251
    dpd_contract424(&L, &R1, &V, 1, 0, 0, 1.0, 0.0);
 
252
    dpd_file2_close(&R1);
 
253
    dpd_buf4_close(&L);
 
254
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 7, 11, "L2R1_VVOV(pqsr)");
 
255
    dpd_buf4_close(&V);
 
256
 
 
257
    dpd_buf4_init(&V, EOM_TMP, G_irr, 7, 10, 7, 10, 0, "L2R1_vvov");
 
258
    dpd_buf4_init(&L, CC_GL, L_irr, 0, 7, 2, 7, 0, "Lijab");
 
259
    dpd_file2_init(&R1, CC_GR, R_irr, 0, 1, "Ria");
 
260
    dpd_contract424(&L, &R1, &V, 1, 0, 0, 1.0, 0.0);
 
261
    dpd_file2_close(&R1);
 
262
    dpd_buf4_close(&L);
 
263
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 7, 11, "L2R1_vvov(pqsr)");
 
264
    dpd_buf4_close(&V);
 
265
 
 
266
    dpd_buf4_init(&V, EOM_TMP, G_irr, 5, 10, 5, 10, 0, "L2R1_VvOv");
 
267
    dpd_buf4_init(&L, CC_GL, L_irr, 0, 5, 0, 5, 0, "LIjAb");
 
268
    dpd_file2_init(&R1, CC_GR, R_irr, 0, 1, "Ria");
 
269
    dpd_contract424(&L, &R1, &V, 1, 0, 0, 1.0, 0.0);
 
270
    dpd_file2_close(&R1);
 
271
    dpd_buf4_close(&L);
 
272
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 5, 11, "L2R1_VvOv(pqsr)");
 
273
    dpd_buf4_close(&V);
 
274
 
 
275
    dpd_buf4_init(&V, EOM_TMP, G_irr, 5, 10, 5, 10, 0, "L2R1_VvoV");
 
276
    dpd_buf4_init(&L, CC_GL, L_irr, 0, 5, 0, 5, 0, "LiJAb");
 
277
    dpd_file2_init(&R1, CC_GR, R_irr, 0, 1, "RIA");
 
278
    dpd_contract424(&L, &R1, &V, 1, 0, 0, 1.0, 0.0);
 
279
    dpd_file2_close(&R1);
 
280
    dpd_buf4_close(&L);
 
281
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 5, 11, "L2R1_VvoV(pqsr)");
 
282
    dpd_buf4_close(&V);
 
283
 
 
284
    /* R2L1_VVOV = Rimab Lmc */
 
285
    dpd_buf4_init(&V, EOM_TMP, G_irr, 7, 10, 7, 10, 0, "R2L1_VVOV");
 
286
    dpd_buf4_init(&R, CC_GR, R_irr, 0, 7, 2, 7, 0, "RIJAB");
 
287
    dpd_file2_init(&L1, CC_GL, L_irr, 0, 1, "LIA");
 
288
    dpd_contract424(&R, &L1, &V, 1, 0, 0, 1.0, 0.0);
 
289
    dpd_file2_close(&L1);
 
290
    dpd_buf4_close(&R);
 
291
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 7, 11, "R2L1_VVOV(pqsr)");
 
292
    dpd_buf4_close(&V);
 
293
 
 
294
    dpd_buf4_init(&V, EOM_TMP, G_irr, 7, 10, 7, 10, 0, "R2L1_vvov");
 
295
    dpd_buf4_init(&R, CC_GR, R_irr, 0, 7, 2, 7, 0, "Rijab");
 
296
    dpd_file2_init(&L1, CC_GL, L_irr, 0, 1, "Lia");
 
297
    dpd_contract424(&R, &L1, &V, 1, 0, 0, 1.0, 0.0);
 
298
    dpd_file2_close(&L1);
 
299
    dpd_buf4_close(&R);
 
300
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 7, 11, "R2L1_vvov(pqsr)");
 
301
    dpd_buf4_close(&V);
 
302
 
 
303
    dpd_buf4_init(&V, EOM_TMP, G_irr, 5, 10, 5, 10, 0, "R2L1_VvOv");
 
304
    dpd_buf4_init(&R, CC_GR, R_irr, 0, 5, 0, 5, 0, "RIjAb");
 
305
    dpd_file2_init(&L1, CC_GL, L_irr, 0, 1, "Lia");
 
306
    dpd_contract424(&R, &L1, &V, 1, 0, 0, 1.0, 0.0);
 
307
    dpd_file2_close(&L1);
 
308
    dpd_buf4_close(&R);
 
309
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 5, 11, "R2L1_VvOv(pqsr)");
 
310
    dpd_buf4_close(&V);
 
311
 
 
312
    dpd_buf4_init(&V, EOM_TMP, G_irr, 5, 10, 5, 10, 0, "R2L1_VvoV");
 
313
    dpd_buf4_init(&R, CC_GR, R_irr, 0, 5, 0, 5, 0, "RiJAb");
 
314
    dpd_file2_init(&L1, CC_GL, L_irr, 0, 1, "LIA");
 
315
    dpd_contract424(&R, &L1, &V, 1, 0, 0, 1.0, 0.0);
 
316
    dpd_file2_close(&L1);
 
317
    dpd_buf4_close(&R);
 
318
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 5, 11, "R2L1_VvoV(pqsr)");
 
319
    dpd_buf4_close(&V);
 
320
 
 
321
    dpd_buf4_init(&T2, CC_TAMPS, 0, 0, 5, 0, 5, 0, "tauIjAb");
 
322
    dpd_buf4_sort(&T2, CC_TAMPS, qprs, 0, 5, "taujIAb");
 
323
    dpd_buf4_close(&T2);
 
324
 
 
325
  }
 
326
  else if(params.ref == 2) { /** UHF **/
 
327
    dpd_buf4_init(&V, EOM_TMP, G_irr, 2, 2, 2, 2, 0, "R2L2_OOOO");
 
328
    dpd_buf4_init(&R, CC_GR, R_irr, 2, 7, 2, 7, 0, "RIJAB");
 
329
    dpd_buf4_init(&L, CC_GL, L_irr, 2, 7, 2, 7, 0, "LIJAB");
 
330
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 0.0);
 
331
    dpd_buf4_close(&L);
 
332
    dpd_buf4_close(&R);
 
333
    dpd_buf4_close(&V);
 
334
    dpd_buf4_init(&V, EOM_TMP, G_irr, 12, 12, 12, 12, 0, "R2L2_oooo");
 
335
    dpd_buf4_init(&R, CC_GR, R_irr, 12, 17, 12, 17, 0, "Rijab");
 
336
    dpd_buf4_init(&L, CC_GL, L_irr, 12, 17, 12, 17, 0, "Lijab");
 
337
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 0.0);
 
338
    dpd_buf4_close(&L);
 
339
    dpd_buf4_close(&R);
 
340
    dpd_buf4_close(&V);
 
341
    dpd_buf4_init(&V, EOM_TMP, G_irr, 22, 22, 22, 22, 0, "R2L2_OoOo");
 
342
    dpd_buf4_init(&R, CC_GR, R_irr, 22, 28, 22, 28, 0, "RIjAb");
 
343
    dpd_buf4_init(&L, CC_GL, L_irr, 22, 28, 22, 28, 0, "LIjAb");
 
344
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 0.0);
 
345
    dpd_buf4_close(&L);
 
346
    dpd_buf4_close(&R);
 
347
    dpd_buf4_sort(&V, EOM_TMP, qpsr, 23, 23, "R2L2_oOoO");
 
348
    dpd_buf4_close(&V);
 
349
 
 
350
    /* Tau2L2_OOOO = 0.5 Tau_ijef Lklef */
 
351
    dpd_buf4_init(&V, EOM_TMP, L_irr, 2, 2, 2, 2, 0, "Tau2L2_OOOO");
 
352
    dpd_buf4_init(&T2, CC_TAMPS, 0, 2, 7, 2, 7, 0, "tauIJAB");
 
353
    dpd_buf4_init(&L, CC_GL, L_irr, 2, 7, 2, 7, 0, "LIJAB");
 
354
    dpd_contract444(&T2, &L, &V, 0, 0, 1.0, 0.0);
 
355
    dpd_buf4_close(&L);
 
356
    dpd_buf4_close(&T2);
 
357
    dpd_buf4_close(&V);
 
358
    dpd_buf4_init(&V, EOM_TMP, L_irr, 12, 12, 12, 12, 0, "Tau2L2_oooo");
 
359
    dpd_buf4_init(&T2, CC_TAMPS, 0, 12, 17, 12, 17, 0, "tauijab");
 
360
    dpd_buf4_init(&L, CC_GL, L_irr, 12, 17, 12, 17, 0, "Lijab");
 
361
    dpd_contract444(&T2, &L, &V, 0, 0, 1.0, 0.0);
 
362
    dpd_buf4_close(&L);
 
363
    dpd_buf4_close(&T2);
 
364
    dpd_buf4_close(&V);
 
365
    dpd_buf4_init(&V, EOM_TMP, L_irr, 22, 22, 22, 22, 0, "Tau2L2_OoOo");
 
366
    dpd_buf4_init(&T2, CC_TAMPS, 0, 22, 28, 22, 28, 0, "tauIjAb");
 
367
    dpd_buf4_init(&L, CC_GL, L_irr, 22, 28, 22, 28, 0, "LIjAb");
 
368
    dpd_contract444(&T2, &L, &V, 0, 0, 1.0, 0.0);
 
369
    dpd_buf4_close(&L);
 
370
    dpd_buf4_close(&T2);
 
371
    dpd_buf4_sort(&V, EOM_TMP, qpsr, 23, 23, "Tau2L2_oOoO");
 
372
    dpd_buf4_close(&V);
 
373
 
 
374
    /* R2L2_OVOV = Rimae Ljmbe */
 
375
    dpd_buf4_init(&V, EOM_TMP, G_irr, 20, 20, 20, 20, 0, "R2L2_OVOV");
 
376
    dpd_buf4_init(&R, CC_GR, R_irr, 20, 20, 20, 20, 0, "RIAJB");
 
377
    dpd_buf4_init(&L, CC_GL, L_irr, 20, 20, 20, 20, 0, "LIAJB");
 
378
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 0.0);
 
379
    dpd_buf4_close(&L);
 
380
    dpd_buf4_close(&R);
 
381
    dpd_buf4_init(&R, CC_GR, R_irr, 20, 30, 20, 30, 0, "RIAjb");
 
382
    dpd_buf4_init(&L, CC_GL, L_irr, 20, 30, 20, 30, 0, "LIAjb");
 
383
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 1.0);
 
384
    dpd_buf4_close(&L);
 
385
    dpd_buf4_close(&R);
 
386
    dpd_buf4_close(&V);
 
387
 
 
388
    dpd_buf4_init(&V, EOM_TMP, G_irr, 30, 30, 30, 30, 0, "R2L2_ovov");
 
389
    dpd_buf4_init(&R, CC_GR, R_irr, 30, 30, 30, 30, 0, "Riajb");
 
390
    dpd_buf4_init(&L, CC_GL, L_irr, 30, 30, 30, 30, 0, "Liajb");
 
391
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 0.0);
 
392
    dpd_buf4_close(&L);
 
393
    dpd_buf4_close(&R);
 
394
    dpd_buf4_init(&R, CC_GR, R_irr, 20, 30, 20, 30, 0, "RIAjb");
 
395
    dpd_buf4_init(&L, CC_GL, L_irr, 20, 30, 20, 30, 0, "LIAjb");
 
396
    dpd_contract444(&R, &L, &V, 1, 1, 1.0, 1.0);
 
397
    dpd_buf4_close(&L);
 
398
    dpd_buf4_close(&R);
 
399
    dpd_buf4_close(&V);
 
400
 
 
401
    dpd_buf4_init(&V, EOM_TMP, G_irr, 20, 30, 20, 30, 0, "R2L2_OVov");
 
402
    dpd_buf4_init(&R, CC_GR, R_irr, 20, 30, 20, 30, 0, "RIAjb");
 
403
    dpd_buf4_init(&L, CC_GL, L_irr, 30, 30, 30, 30, 0, "Liajb");
 
404
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 0.0);
 
405
    dpd_buf4_close(&L);
 
406
    dpd_buf4_close(&R);
 
407
    dpd_buf4_init(&R, CC_GR, R_irr, 20, 20, 20, 20, 0, "RIAJB");
 
408
    dpd_buf4_init(&L, CC_GL, L_irr, 20, 30, 20, 30, 0, "LIAjb");
 
409
    dpd_contract444(&R, &L, &V, 0, 1, 1.0, 1.0);
 
410
    dpd_buf4_close(&L);
 
411
    dpd_buf4_close(&R);
 
412
    dpd_buf4_close(&V);
 
413
 
 
414
    dpd_buf4_init(&V, EOM_TMP, G_irr, 30, 20, 30, 20, 0, "R2L2_ovOV");
 
415
    dpd_buf4_init(&R, CC_GR, R_irr, 30, 20, 30, 20, 0, "RiaJB");
 
416
    dpd_buf4_init(&L, CC_GL, L_irr, 20, 20, 20, 20, 0, "LIAJB");
 
417
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 0.0);
 
418
    dpd_buf4_close(&L);
 
419
    dpd_buf4_close(&R);
 
420
    dpd_buf4_init(&R, CC_GR, R_irr, 30, 30, 30, 30, 0, "Riajb");
 
421
    dpd_buf4_init(&L, CC_GL, L_irr, 20, 30, 20, 30, 0, "LIAjb");
 
422
    dpd_contract444(&R, &L, &V, 0, 0, 1.0, 1.0);
 
423
    dpd_buf4_close(&L);
 
424
    dpd_buf4_close(&R);
 
425
    dpd_buf4_close(&V);
 
426
 
 
427
    dpd_buf4_init(&V, EOM_TMP, G_irr, 27, 27, 27, 27, 0, "R2L2_oVoV");
 
428
    dpd_buf4_init(&L, CC_GL, L_irr, 24, 27, 24, 27, 0, "LIbjA");
 
429
    dpd_buf4_init(&R, CC_GR, R_irr, 27, 24, 27, 24, 0, "RjAIb");
 
430
    dpd_contract444(&R, &L, &V, 0, 1, 1.0, 0.0);
 
431
    dpd_buf4_close(&R);
 
432
    dpd_buf4_close(&L);
 
433
    dpd_buf4_close(&V);
 
434
 
 
435
    dpd_buf4_init(&V, EOM_TMP, G_irr, 24, 24, 24, 24, 0, "R2L2_OvOv");
 
436
    dpd_buf4_init(&L, CC_GL, L_irr, 27, 24, 27, 24, 0, "LjAIb");
 
437
    dpd_buf4_init(&R, CC_GR, R_irr, 24, 27, 24, 27, 0, "RIbjA");
 
438
    dpd_contract444(&R, &L, &V, 0, 1, 1.0, 0.0);
 
439
    dpd_buf4_close(&R);
 
440
    dpd_buf4_close(&L);
 
441
    dpd_buf4_close(&V);
 
442
 
 
443
    /* L2R1_OOVO = Lijae Rke */
 
444
    dpd_buf4_init(&V, EOM_TMP, G_irr, 2, 21, 2, 21, 0, "L2R1_OOVO");
 
445
    dpd_buf4_init(&L, CC_GL, L_irr, 2, 5, 2, 7, 0, "LIJAB");
 
446
    dpd_file2_init(&R1, CC_GR, R_irr, 0, 1, "RIA");
 
447
    dpd_contract424(&L, &R1, &V, 3, 1, 0, 1.0, 0.0);
 
448
    dpd_file2_close(&R1);
 
449
    dpd_buf4_close(&L);
 
450
    dpd_buf4_close(&V);
 
451
 
 
452
    dpd_buf4_init(&V, EOM_TMP, G_irr, 12, 31, 12, 31, 0, "L2R1_oovo");
 
453
    dpd_buf4_init(&L, CC_GL, L_irr, 12, 15, 12, 17, 0, "Lijab");
 
454
    dpd_file2_init(&R1, CC_GR, R_irr, 2, 3, "Ria");
 
455
    dpd_contract424(&L, &R1, &V, 3, 1, 0, 1.0, 0.0);
 
456
    dpd_file2_close(&R1);
 
457
    dpd_buf4_close(&L);
 
458
    dpd_buf4_close(&V);
 
459
 
 
460
    dpd_buf4_init(&V, EOM_TMP, G_irr, 22, 26, 22, 26, 0, "L2R1_OoVo");
 
461
    dpd_buf4_init(&L, CC_GL, L_irr, 22, 28, 22, 28, 0, "LIjAb");
 
462
    dpd_file2_init(&R1, CC_GR, R_irr, 2, 3, "Ria");
 
463
    dpd_contract424(&L, &R1, &V, 3, 1, 0, 1.0, 0.0);
 
464
    dpd_file2_close(&R1);
 
465
    dpd_buf4_close(&L);
 
466
    dpd_buf4_close(&V);
 
467
 
 
468
    dpd_buf4_init(&V, EOM_TMP, G_irr, 22, 25, 22, 25, 0, "L2R1_OovO");
 
469
    dpd_buf4_init(&L, CC_GL, L_irr, 22, 29, 22, 29, 0, "LIjaB");
 
470
    dpd_file2_init(&R1, CC_GR, R_irr, 0, 1, "RIA");
 
471
    dpd_contract424(&L, &R1, &V, 3, 1, 0, 1.0, 0.0);
 
472
    dpd_file2_close(&R1);
 
473
    dpd_buf4_close(&L);
 
474
    dpd_buf4_close(&V);
 
475
 
 
476
    dpd_buf4_init(&V, EOM_TMP, G_irr, 2, 21, 2, 21, 0, "L2R1_OOVO");
 
477
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 2, 20, "L2R1_OOVO(pqsr)");
 
478
    dpd_buf4_close(&V);
 
479
    dpd_buf4_init(&V, EOM_TMP, G_irr, 12, 31, 12, 31, 0, "L2R1_oovo");
 
480
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 12, 30, "L2R1_oovo(pqsr)");
 
481
    dpd_buf4_close(&V);
 
482
 
 
483
    dpd_buf4_init(&V, EOM_TMP, G_irr, 22, 26, 22, 26, 0, "L2R1_OoVo");
 
484
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 22, 27, "L2R1_OoVo(pqsr)");
 
485
    dpd_buf4_sort(&V, EOM_TMP, qprs, 23, 26, "L2R1_OoVo(qprs)");
 
486
    dpd_buf4_close(&V);
 
487
    dpd_buf4_init(&V, EOM_TMP, G_irr, 22, 27, 22, 27, 0, "L2R1_OoVo(pqsr)");
 
488
    dpd_buf4_sort(&V, EOM_TMP, qprs, 23, 27, "L2R1_OoVo(qpsr)");
 
489
    dpd_buf4_close(&V);
 
490
 
 
491
    dpd_buf4_init(&V, EOM_TMP, G_irr, 22, 25, 22, 25, 0, "L2R1_OovO");
 
492
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 22, 24, "L2R1_OovO(pqsr)");
 
493
    dpd_buf4_sort(&V, EOM_TMP, qprs, 23, 25, "L2R1_OovO(qprs)");
 
494
    dpd_buf4_close(&V);
 
495
 
 
496
    /* L1R2_OOVO = Rijae Lke */
 
497
    dpd_buf4_init(&V, EOM_TMP, G_irr, 2, 21, 2, 21, 0, "L1R2_OOVO");
 
498
    dpd_buf4_init(&R, CC_GR, R_irr, 2, 5, 2, 7, 0, "RIJAB");
 
499
    dpd_file2_init(&L1, CC_GL, L_irr, 0, 1, "LIA");
 
500
    dpd_contract424(&R, &L1, &V, 3, 1, 0, 1.0, 0.0);
 
501
    dpd_file2_close(&L1);
 
502
    dpd_buf4_close(&R);
 
503
    dpd_buf4_close(&V);
 
504
 
 
505
    dpd_buf4_init(&V, EOM_TMP, G_irr, 12, 31, 12, 31, 0, "L1R2_oovo");
 
506
    dpd_buf4_init(&R, CC_GR, R_irr, 12, 15, 12, 17, 0, "Rijab");
 
507
    dpd_file2_init(&L1, CC_GL, L_irr, 2, 3, "Lia");
 
508
    dpd_contract424(&R, &L1, &V, 3, 1, 0, 1.0, 0.0);
 
509
    dpd_file2_close(&L1);
 
510
    dpd_buf4_close(&R);
 
511
    dpd_buf4_close(&V);
 
512
 
 
513
    dpd_buf4_init(&V, EOM_TMP, G_irr, 22, 26, 22, 26, 0, "L1R2_OoVo");
 
514
    dpd_buf4_init(&R, CC_GR, R_irr, 22, 28, 22, 28, 0, "RIjAb");
 
515
    dpd_file2_init(&L1, CC_GL, L_irr, 2, 3, "Lia");
 
516
    dpd_contract424(&R, &L1, &V, 3, 1, 0, 1.0, 0.0);
 
517
    dpd_file2_close(&L1);
 
518
    dpd_buf4_close(&R);
 
519
    dpd_buf4_close(&V);
 
520
 
 
521
    dpd_buf4_init(&V, EOM_TMP, G_irr, 22, 25, 22, 25, 0, "L1R2_OovO");
 
522
    dpd_buf4_init(&R, CC_GR, R_irr, 22, 29, 22, 29, 0, "RIjaB");
 
523
    dpd_file2_init(&L1, CC_GL, L_irr, 0, 1, "LIA");
 
524
    dpd_contract424(&R, &L1, &V, 3, 1, 0, 1.0, 0.0);
 
525
    dpd_file2_close(&L1);
 
526
    dpd_buf4_close(&R);
 
527
    dpd_buf4_close(&V);
 
528
 
 
529
    dpd_buf4_init(&V, EOM_TMP, G_irr, 2, 21, 2, 21, 0, "L1R2_OOVO");
 
530
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 2, 20, "L1R2_OOVO(pqsr)");
 
531
    dpd_buf4_close(&V);
 
532
    dpd_buf4_init(&V, EOM_TMP, G_irr, 12, 31, 12, 31, 0, "L1R2_oovo");
 
533
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 12, 30, "L1R2_oovo(pqsr)");
 
534
    dpd_buf4_close(&V);
 
535
    dpd_buf4_init(&V, EOM_TMP, G_irr, 22, 26, 22, 26, 0, "L1R2_OoVo");
 
536
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 22, 27, "L1R2_OoVo(pqsr)");
 
537
    dpd_buf4_close(&V);
 
538
    dpd_buf4_init(&V, EOM_TMP, G_irr, 22, 27, 22, 27, 0, "L1R2_OoVo(pqsr)");
 
539
    dpd_buf4_sort(&V, EOM_TMP, qprs, 23, 27, "L1R2_OoVo(qpsr)");
 
540
    dpd_buf4_close(&V);
 
541
    dpd_buf4_init(&V, EOM_TMP, G_irr, 22, 25, 22, 25, 0, "L1R2_OovO");
 
542
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 22, 24, "L1R2_OovO(pqsr)");
 
543
    dpd_buf4_close(&V);
 
544
 
 
545
    /* L2R1_VVOV = Limab Rmc */
 
546
    dpd_buf4_init(&V, EOM_TMP, G_irr, 7, 20, 7, 20, 0, "L2R1_VVOV");
 
547
    dpd_buf4_init(&L, CC_GL, L_irr, 0, 7, 2, 7, 0, "LIJAB");
 
548
    dpd_file2_init(&R1, CC_GR, R_irr, 0, 1, "RIA");
 
549
    dpd_contract424(&L, &R1, &V, 1, 0, 0, 1.0, 0.0);
 
550
    dpd_file2_close(&R1);
 
551
    dpd_buf4_close(&L);
 
552
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 7, 21, "L2R1_VVOV(pqsr)");
 
553
    dpd_buf4_close(&V);
 
554
 
 
555
    dpd_buf4_init(&V, EOM_TMP, G_irr, 17, 30, 17, 30, 0, "L2R1_vvov");
 
556
    dpd_buf4_init(&L, CC_GL, L_irr, 10, 17, 12, 17, 0, "Lijab");
 
557
    dpd_file2_init(&R1, CC_GR, R_irr, 2, 3, "Ria");
 
558
    dpd_contract424(&L, &R1, &V, 1, 0, 0, 1.0, 0.0);
 
559
    dpd_file2_close(&R1);
 
560
    dpd_buf4_close(&L);
 
561
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 17, 31, "L2R1_vvov(pqsr)");
 
562
    dpd_buf4_close(&V);
 
563
 
 
564
    dpd_buf4_init(&V, EOM_TMP, G_irr, 28, 24, 28, 24, 0, "L2R1_VvOv");
 
565
    dpd_buf4_init(&L, CC_GL, L_irr, 22, 28, 22, 28, 0, "LIjAb");
 
566
    dpd_file2_init(&R1, CC_GR, R_irr, 2, 3, "Ria");
 
567
    dpd_contract424(&L, &R1, &V, 1, 0, 0, 1.0, 0.0);
 
568
    dpd_file2_close(&R1);
 
569
    dpd_buf4_close(&L);
 
570
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 28, 25, "L2R1_VvOv(pqsr)");
 
571
    dpd_buf4_close(&V);
 
572
 
 
573
    dpd_buf4_init(&V, EOM_TMP, G_irr, 28, 27, 28, 27, 0, "L2R1_VvoV");
 
574
    dpd_buf4_init(&L, CC_GL, L_irr, 23, 28, 23, 28, 0, "LiJAb");
 
575
    dpd_file2_init(&R1, CC_GR, R_irr, 0, 1, "RIA");
 
576
    dpd_contract424(&L, &R1, &V, 1, 0, 0, 1.0, 0.0);
 
577
    dpd_file2_close(&R1);
 
578
    dpd_buf4_close(&L);
 
579
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 28, 26, "L2R1_VvoV(pqsr)");
 
580
    dpd_buf4_close(&V);
 
581
 
 
582
    /* R2L1_VVOV = Rimab Lmc */
 
583
    dpd_buf4_init(&V, EOM_TMP, G_irr, 7, 20, 7, 20, 0, "R2L1_VVOV");
 
584
    dpd_buf4_init(&R, CC_GR, R_irr, 0, 7, 2, 7, 0, "RIJAB");
 
585
    dpd_file2_init(&L1, CC_GL, L_irr, 0, 1, "LIA");
 
586
    dpd_contract424(&R, &L1, &V, 1, 0, 0, 1.0, 0.0);
 
587
    dpd_file2_close(&L1);
 
588
    dpd_buf4_close(&R);
 
589
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 7, 21, "R2L1_VVOV(pqsr)");
 
590
    dpd_buf4_close(&V);
 
591
 
 
592
    dpd_buf4_init(&V, EOM_TMP, G_irr, 17, 30, 17, 30, 0, "R2L1_vvov");
 
593
    dpd_buf4_init(&R, CC_GR, R_irr, 10, 17, 12, 17, 0, "Rijab");
 
594
    dpd_file2_init(&L1, CC_GL, L_irr, 2, 3, "Lia");
 
595
    dpd_contract424(&R, &L1, &V, 1, 0, 0, 1.0, 0.0);
 
596
    dpd_file2_close(&L1);
 
597
    dpd_buf4_close(&R);
 
598
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 17, 31, "R2L1_vvov(pqsr)");
 
599
    dpd_buf4_close(&V);
 
600
 
 
601
    dpd_buf4_init(&V, EOM_TMP, G_irr, 28, 24, 28, 24, 0, "R2L1_VvOv");
 
602
    dpd_buf4_init(&R, CC_GR, R_irr, 22, 28, 22, 28, 0, "RIjAb");
 
603
    dpd_file2_init(&L1, CC_GL, L_irr, 2, 3, "Lia");
 
604
    dpd_contract424(&R, &L1, &V, 1, 0, 0, 1.0, 0.0);
 
605
    dpd_file2_close(&L1);
 
606
    dpd_buf4_close(&R);
 
607
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 28, 25, "R2L1_VvOv(pqsr)");
 
608
    dpd_buf4_close(&V);
 
609
 
 
610
    dpd_buf4_init(&V, EOM_TMP, G_irr, 28, 27, 28, 27, 0, "R2L1_VvoV");
 
611
    dpd_buf4_init(&R, CC_GR, R_irr, 23, 28, 23, 28, 0, "RiJAb");
 
612
    dpd_file2_init(&L1, CC_GL, L_irr, 0, 1, "LIA");
 
613
    dpd_contract424(&R, &L1, &V, 1, 0, 0, 1.0, 0.0);
 
614
    dpd_file2_close(&L1);
 
615
    dpd_buf4_close(&R);
 
616
    dpd_buf4_sort(&V, EOM_TMP, pqsr, 28, 26, "R2L1_VvoV(pqsr)");
 
617
    dpd_buf4_close(&V);
 
618
 
 
619
    dpd_buf4_init(&T2, CC_TAMPS, 0, 22, 28, 22, 28, 0, "tauIjAb");
 
620
    dpd_buf4_sort(&T2, CC_TAMPS, qprs, 23, 28, "taujIAb");
 
621
    dpd_buf4_close(&T2);
 
622
  }
 
623
}
 
624
 
 
625
}} // namespace psi::ccdensity