← Back to branch summary

~ubuntu-branches/ubuntu/precise/psicode/precise

« back to all changes in this revision

Viewing changes to src/bin/ccdensity/ltdensity_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/ltdensity_intermediates.cc
 
1
/*! \file
 
2
    \ingroup CCDENSITY
 
3
    \brief Enter brief description of file here 
 
4
*/
 
5
#include "MOInfo.h"
 
6
#include "Params.h"
 
7
#include "Frozen.h"
 
8
#define EXTERN
 
9
#include "globals.h"
 
10
 
 
11
namespace psi { namespace ccdensity {
 
12
 
 
13
/* build left transition moment for oscillator strengths 
 
14
   and rotational strengths */
 
15
 
 
16
/* LR1_OO  =  LIE * RJE */
 
17
/* LR1_oo  = Lie * Rje */
 
18
/* LR1_VV = LMA * RMB */
 
19
/* LR1_vv = Lma * Rmb */
 
20
/* LT1_OO  =  LIE * TJE */
 
21
/* LT1_oo  = Lie * Tje */
 
22
/* LT1_VV = LMA * TMB */
 
23
/* LT1_vv = Lma * Tmb */
 
24
/* L2R1_OV = RME * LIMAE + Rme + LImAe */
 
25
/* L2R1_ov = Rme * Limae + RME + LiMaE */
 
26
/* L1R2_OV = LME * RIMAE + Lme * RImAe */
 
27
/* L1R2_ov = Lme * Rimae + LME * RiMaE */
 
28
/* L1T2_OV = LME * TIMAE + Lme * TImAe */
 
29
/* L1T2_ov = Lme * Timae + LME * TiMaE */
 
30
/* LR2_OO  = 0.5 * LIMEF * RJMEF + LImEf * RJmEf */
 
31
/* LR2_oo  = 0.5 * Limef * Rjmef + LiMeF * RjMeF */
 
32
/* LR2_VV = 0.5 * LMNEA * RMNEB + LmNeA * RmNeB */
 
33
/* LR2_vv = 0.5 * Lmnea * Rmneb + LMnEa * RMnEb */
 
34
/* LT2_OO = 0.5 * LIMEF * TJMEF + LImEf * TJmEf */
 
35
/* LT2_oo = 0.5 * Limef * Tjmef + LiMeF * TjMeF */
 
36
/* LT2_VV = 0.5 * LMNEA * TMNEB + LmNeA * TmNeB */
 
37
/* LT2_vv = 0.5 * Lmnea * Tmneb + LMnEa * TMnEb */
 
38
/* LR_OO = LR1_OO + LR2_OO */
 
39
/* LR_oo = LR1_oo + LR2_oo */
 
40
/* LR_VV = LR1_VV + LR2_VV */
 
41
/* LR_vv = LR1_vv + LR2_vv */
 
42
/* LT_OO = LT1_OO + LT2_OO */
 
43
/* LT_oo = LT1_oo + LT2_oo */
 
44
 
 
45
void ltdensity_intermediates(struct TD_Params S)
 
46
{
 
47
  dpdfile2 L1, R1, T1, I, LR1, LR2, LT1, LT2;
 
48
  dpdbuf4 L2, T2, R2;
 
49
  int rohf = 0;
 
50
 
 
51
  if ( (params.ref == 0) || (params.ref == 1) ) rohf = 1;
 
52
 
 
53
  /* LR1_OO(I,J)  =  LIE * RJE */
 
54
  dpd_file2_init(&I, EOM_TMP, S.irrep, 0, 0, "LR1_OO");
 
55
  dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "LIA");
 
56
  dpd_file2_init(&R1, CC_GR, S.irrep, 0, 1, "RIA");
 
57
  dpd_contract222(&L1, &R1, &I, 0, 0, 1.0, 0.0);
 
58
  dpd_file2_close(&R1);
 
59
  dpd_file2_close(&L1);
 
60
  dpd_file2_close(&I);
 
61
 
 
62
  /* LR1_oo(i,j)  = Lia * Rje */
 
63
 
 
64
  if (rohf) {
 
65
    dpd_file2_init(&I, EOM_TMP, S.irrep, 0, 0, "LR1_oo");
 
66
    dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "Lia");
 
67
    dpd_file2_init(&R1, CC_GR, S.irrep, 0, 1, "Ria");
 
68
    dpd_contract222(&L1, &R1, &I, 0, 0, 1.0, 0.0);
 
69
    dpd_file2_close(&R1);
 
70
    dpd_file2_close(&L1);
 
71
    dpd_file2_close(&I);
 
72
  }
 
73
  else {
 
74
    dpd_file2_init(&I, EOM_TMP, S.irrep, 2, 2, "LR1_oo");
 
75
    dpd_file2_init(&L1, CC_GLG, 0, 2, 3, "Lia");
 
76
    dpd_file2_init(&R1, CC_GR, S.irrep, 2, 3, "Ria");
 
77
    dpd_contract222(&L1, &R1, &I, 0, 0, 1.0, 0.0);
 
78
    dpd_file2_close(&R1);
 
79
    dpd_file2_close(&L1);
 
80
    dpd_file2_close(&I);
 
81
  }
 
82
 
 
83
  /* LR1_VV(A,B) = LMA * RMB */
 
84
 
 
85
  dpd_file2_init(&I, EOM_TMP, S.irrep, 1, 1, "LR1_VV");
 
86
  dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "LIA");
 
87
  dpd_file2_init(&R1, CC_GR, S.irrep, 0, 1, "RIA");
 
88
  dpd_contract222(&L1, &R1, &I, 1, 1, 1.0, 0.0);
 
89
  dpd_file2_close(&R1);
 
90
  dpd_file2_close(&L1);
 
91
  dpd_file2_close(&I);
 
92
 
 
93
  /* LR1_vv(a,b) = Lma * Rmb */
 
94
 
 
95
  if (rohf) {
 
96
    dpd_file2_init(&I, EOM_TMP, S.irrep, 1, 1, "LR1_vv");
 
97
    dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "Lia");
 
98
    dpd_file2_init(&R1, CC_GR, S.irrep, 0, 1, "Ria");
 
99
    dpd_contract222(&L1, &R1, &I, 1, 1, 1.0, 0.0);
 
100
    dpd_file2_close(&R1);
 
101
    dpd_file2_close(&L1);
 
102
    dpd_file2_close(&I);
 
103
  }
 
104
  else {
 
105
    dpd_file2_init(&I, EOM_TMP, S.irrep, 3, 3, "LR1_vv");
 
106
    dpd_file2_init(&L1, CC_GLG, 0, 2, 3, "Lia");
 
107
    dpd_file2_init(&R1, CC_GR, S.irrep, 2, 3, "Ria");
 
108
    dpd_contract222(&L1, &R1, &I, 1, 1, 1.0, 0.0);
 
109
    dpd_file2_close(&R1);
 
110
    dpd_file2_close(&L1);
 
111
    dpd_file2_close(&I);
 
112
  }
 
113
  
 
114
  /* LT1_OO(I,J)  =  LIE * TJE */
 
115
 
 
116
  dpd_file2_init(&I, EOM_TMP, 0, 0, 0, "LT1_OO");
 
117
  dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "LIA");
 
118
  dpd_file2_init(&T1, CC_OEI, 0, 0, 1, "tIA");
 
119
  dpd_contract222(&L1, &T1, &I, 0, 0, 1.0, 0.0);
 
120
  dpd_file2_close(&T1);
 
121
  dpd_file2_close(&L1);
 
122
  dpd_file2_close(&I);
 
123
 
 
124
  /* LT1_oo(i,j)  = Lia * Rje */
 
125
 
 
126
  if (rohf) {
 
127
    dpd_file2_init(&I, EOM_TMP, 0, 0, 0, "LT1_oo");
 
128
    dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "Lia");
 
129
    dpd_file2_init(&T1, CC_OEI, 0, 0, 1, "tia");
 
130
    dpd_contract222(&L1, &T1, &I, 0, 0, 1.0, 0.0);
 
131
    dpd_file2_close(&T1);
 
132
    dpd_file2_close(&L1);
 
133
    dpd_file2_close(&I);
 
134
  }
 
135
  else {
 
136
    dpd_file2_init(&I, EOM_TMP, 0, 2, 2, "LT1_oo");
 
137
    dpd_file2_init(&L1, CC_GLG, 0, 2, 3, "Lia");
 
138
    dpd_file2_init(&T1, CC_OEI, 0, 2, 3, "tia");
 
139
    dpd_contract222(&L1, &T1, &I, 0, 0, 1.0, 0.0);
 
140
    dpd_file2_close(&T1);
 
141
    dpd_file2_close(&L1);
 
142
    dpd_file2_close(&I);
 
143
  }
 
144
 
 
145
  /* LT1_VV(A,B) = LMA * TMB */
 
146
 
 
147
  dpd_file2_init(&I, EOM_TMP, 0, 1, 1, "LT1_VV");
 
148
  dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "LIA");
 
149
  dpd_file2_init(&T1, CC_OEI, 0, 0, 1, "tIA");
 
150
  dpd_contract222(&L1, &T1, &I, 1, 1, 1.0, 0.0);
 
151
  dpd_file2_close(&T1);
 
152
  dpd_file2_close(&L1);
 
153
  dpd_file2_close(&I);
 
154
 
 
155
  /* LT1_vv(a,b) = Lma * Tmb */
 
156
 
 
157
  if (rohf) {
 
158
    dpd_file2_init(&I, EOM_TMP, 0, 1, 1, "LT1_vv");
 
159
    dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "Lia");
 
160
    dpd_file2_init(&T1, CC_OEI, 0, 0, 1, "tia");
 
161
    dpd_contract222(&L1, &T1, &I, 1, 1, 1.0, 0.0);
 
162
    dpd_file2_close(&T1);
 
163
    dpd_file2_close(&L1);
 
164
    dpd_file2_close(&I);
 
165
  }
 
166
  else {
 
167
    dpd_file2_init(&I, EOM_TMP, 0, 3, 3, "LT1_vv");
 
168
    dpd_file2_init(&L1, CC_GLG, 0, 2, 3, "Lia");
 
169
    dpd_file2_init(&T1, CC_OEI, 0, 2, 3, "tia");
 
170
    dpd_contract222(&L1, &T1, &I, 1, 1, 1.0, 0.0);
 
171
    dpd_file2_close(&T1);
 
172
    dpd_file2_close(&L1);
 
173
    dpd_file2_close(&I);
 
174
  }
 
175
 
 
176
  /* L2R1_OV(I,A) = RME * LIMAE + Rme + LImAe */
 
177
 
 
178
  if (rohf) {
 
179
    dpd_file2_init(&I, EOM_TMP, S.irrep, 0, 1, "L2R1_OV");
 
180
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 5, 2, 7, 0, "LIJAB");
 
181
    dpd_file2_init(&R1, CC_GR, S.irrep, 0, 1, "RIA");
 
182
    dpd_dot24(&R1, &L2, &I, 0, 0, 1.0, 0.0);
 
183
    dpd_file2_close(&R1);
 
184
    dpd_buf4_close(&L2);
 
185
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 5, 0, 5, 0, "LIjAb");
 
186
    dpd_file2_init(&R1, CC_GR, S.irrep, 0, 1, "Ria");
 
187
    dpd_dot24(&R1, &L2, &I, 0, 0, 1.0, 1.0);
 
188
    dpd_file2_close(&R1);
 
189
    dpd_buf4_close(&L2);
 
190
    dpd_file2_close(&I);
 
191
  }
 
192
  else {
 
193
    dpd_file2_init(&I, EOM_TMP, S.irrep, 0, 1, "L2R1_OV");
 
194
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 5, 2, 7, 0, "LIJAB");
 
195
    dpd_file2_init(&R1, CC_GR, S.irrep, 0, 1, "RIA");
 
196
    dpd_dot24(&R1, &L2, &I, 0, 0, 1.0, 0.0);
 
197
    dpd_file2_close(&R1);
 
198
    dpd_buf4_close(&L2);
 
199
    dpd_buf4_init(&L2, CC_GLG, 0, 22, 28, 22, 28, 0, "LIjAb");
 
200
    dpd_file2_init(&R1, CC_GR, S.irrep, 2, 3, "Ria");
 
201
    dpd_dot24(&R1, &L2, &I, 0, 0, 1.0, 1.0);
 
202
    dpd_file2_close(&R1);
 
203
    dpd_buf4_close(&L2);
 
204
    dpd_file2_close(&I);
 
205
  }
 
206
 
 
207
  /* L2R1_OV(i,a) = Rme * Limae + RME + LiMaE */
 
208
 
 
209
  if (rohf) {
 
210
    dpd_file2_init(&I, EOM_TMP, S.irrep, 0, 1, "L2R1_ov");
 
211
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 5, 2, 7, 0, "Lijab");
 
212
    dpd_file2_init(&R1, CC_GR, S.irrep, 0, 1, "Ria");
 
213
    dpd_dot24(&R1, &L2, &I, 0, 0, 1.0, 0.0);
 
214
    dpd_file2_close(&R1);
 
215
    dpd_buf4_close(&L2);
 
216
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 5, 0, 5, 0, "LiJaB");
 
217
    dpd_file2_init(&R1, CC_GR, S.irrep, 0, 1, "RIA");
 
218
    dpd_dot24(&R1, &L2, &I, 0, 0, 1.0, 1.0);
 
219
    dpd_file2_close(&R1);
 
220
    dpd_buf4_close(&L2);
 
221
    dpd_file2_close(&I);
 
222
  }
 
223
  else {
 
224
    dpd_file2_init(&I, EOM_TMP, S.irrep, 2, 3, "L2R1_ov");
 
225
    dpd_buf4_init(&L2, CC_GLG, 0, 10, 15, 12, 17, 0, "Lijab");
 
226
    dpd_file2_init(&R1, CC_GR, S.irrep, 2, 3, "Ria");
 
227
    dpd_dot24(&R1, &L2, &I, 0, 0, 1.0, 0.0);
 
228
    dpd_file2_close(&R1);
 
229
    dpd_buf4_close(&L2);
 
230
    dpd_buf4_init(&L2, CC_GLG, 0, 23, 29, 23, 29, 0, "LiJaB");
 
231
    dpd_file2_init(&R1, CC_GR, S.irrep, 0, 1, "RIA");
 
232
    dpd_dot24(&R1, &L2, &I, 0, 0, 1.0, 1.0);
 
233
    dpd_file2_close(&R1);
 
234
    dpd_buf4_close(&L2);
 
235
    dpd_file2_close(&I);
 
236
  }
 
237
 
 
238
  /* L1R2_OV(I,A) = LME * RIMAE + Lme * RImAe */
 
239
 
 
240
  if (rohf) {
 
241
    dpd_file2_init(&I, EOM_TMP, S.irrep, 0, 1, "L1R2_OV");
 
242
    dpd_buf4_init(&R2, CC_GR, S.irrep, 0, 5, 2, 7, 0, "RIJAB");
 
243
    dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "LIA");
 
244
    dpd_dot24(&L1, &R2, &I, 0, 0, 1.0, 0.0);
 
245
    dpd_file2_close(&L1);
 
246
    dpd_buf4_close(&R2);
 
247
    dpd_buf4_init(&R2, CC_GR, S.irrep, 0, 5, 0, 5, 0, "RIjAb");
 
248
    dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "Lia");
 
249
    dpd_dot24(&L1, &R2, &I, 0, 0, 1.0, 1.0);
 
250
    dpd_file2_close(&L1);
 
251
    dpd_buf4_close(&R2);
 
252
    dpd_file2_close(&I);
 
253
  }
 
254
  else {
 
255
    dpd_file2_init(&I, EOM_TMP, S.irrep, 0, 1, "L1R2_OV");
 
256
    dpd_buf4_init(&R2, CC_GR, S.irrep, 0, 5, 2, 7, 0, "RIJAB");
 
257
    dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "LIA");
 
258
    dpd_dot24(&L1, &R2, &I, 0, 0, 1.0, 0.0);
 
259
    dpd_file2_close(&L1);
 
260
    dpd_buf4_close(&R2);
 
261
    dpd_buf4_init(&R2, CC_GR, S.irrep, 22, 28, 22, 28, 0, "RIjAb");
 
262
    dpd_file2_init(&L1, CC_GLG, 0, 2, 3, "Lia");
 
263
    dpd_dot24(&L1, &R2, &I, 0, 0, 1.0, 1.0);
 
264
    dpd_file2_close(&L1);
 
265
    dpd_buf4_close(&R2);
 
266
    dpd_file2_close(&I);
 
267
  }
 
268
 
 
269
  /* L1R2_ov(i,a) = Lme * Rimae + LME * RiMaE */
 
270
 
 
271
  if (rohf) {
 
272
    dpd_file2_init(&I, EOM_TMP, S.irrep, 0, 1, "L1R2_ov");
 
273
    dpd_buf4_init(&R2, CC_GR, S.irrep, 0, 5, 2, 7, 0, "Rijab");
 
274
    dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "Lia");
 
275
    dpd_dot24(&L1, &R2, &I, 0, 0, 1.0, 0.0);
 
276
    dpd_file2_close(&L1);
 
277
    dpd_buf4_close(&R2);
 
278
    dpd_buf4_init(&R2, CC_GR, S.irrep, 0, 5, 0, 5, 0, "RiJaB");
 
279
    dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "LIA");
 
280
    dpd_dot24(&L1, &R2, &I, 0, 0, 1.0, 1.0);
 
281
    dpd_file2_close(&L1);
 
282
    dpd_buf4_close(&R2);
 
283
    dpd_file2_close(&I);
 
284
  }
 
285
  else {
 
286
    dpd_file2_init(&I, EOM_TMP, S.irrep, 2, 3, "L1R2_ov");
 
287
    dpd_buf4_init(&R2, CC_GR, S.irrep, 10, 15, 12, 17, 0, "Rijab");
 
288
    dpd_file2_init(&L1, CC_GLG, 0, 2, 3, "Lia");
 
289
    dpd_dot24(&L1, &R2, &I, 0, 0, 1.0, 0.0);
 
290
    dpd_file2_close(&L1);
 
291
    dpd_buf4_close(&R2);
 
292
    dpd_buf4_init(&R2, CC_GR, S.irrep, 23, 29, 23, 29, 0, "RiJaB");
 
293
    dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "LIA");
 
294
    dpd_dot24(&L1, &R2, &I, 0, 0, 1.0, 1.0);
 
295
    dpd_file2_close(&L1);
 
296
    dpd_buf4_close(&R2);
 
297
    dpd_file2_close(&I);
 
298
  }
 
299
 
 
300
  /* L1T2_OV = LME * TIMAE + Lme * TImAe */
 
301
 
 
302
  if (rohf) {
 
303
    dpd_file2_init(&I, EOM_TMP, 0, 0, 1, "L1T2_OV");
 
304
    dpd_buf4_init(&T2, CC_TAMPS, 0, 0, 5, 2, 7, 0, "tIJAB");
 
305
    dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "LIA");
 
306
    dpd_dot24(&L1, &T2, &I, 0, 0, 1.0, 0.0);
 
307
    dpd_file2_close(&L1);
 
308
    dpd_buf4_close(&T2);
 
309
    dpd_buf4_init(&T2, CC_TAMPS, 0, 0, 5, 0, 5, 0, "tIjAb");
 
310
    dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "Lia");
 
311
    dpd_dot24(&L1, &T2, &I, 0, 0, 1.0, 1.0);
 
312
    dpd_file2_close(&L1);
 
313
    dpd_buf4_close(&T2);
 
314
    dpd_file2_close(&I);
 
315
  }
 
316
  else {
 
317
    dpd_file2_init(&I, EOM_TMP, 0, 0, 1, "L1T2_OV");
 
318
    dpd_buf4_init(&T2, CC_TAMPS, 0, 0, 5, 2, 7, 0, "tIJAB");
 
319
    dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "LIA");
 
320
    dpd_dot24(&L1, &T2, &I, 0, 0, 1.0, 0.0);
 
321
    dpd_file2_close(&L1);
 
322
    dpd_buf4_close(&T2);
 
323
    dpd_buf4_init(&T2, CC_TAMPS, 0, 22, 28, 22, 28, 0, "tIjAb");
 
324
    dpd_file2_init(&L1, CC_GLG, 0, 2, 3, "Lia");
 
325
    dpd_dot24(&L1, &T2, &I, 0, 0, 1.0, 1.0);
 
326
    dpd_file2_close(&L1);
 
327
    dpd_buf4_close(&T2);
 
328
    dpd_file2_close(&I);
 
329
  }
 
330
 
 
331
  /* L1T2_ov = Lme * Timae + LME * TiMaE */
 
332
 
 
333
  if (rohf) {
 
334
    dpd_file2_init(&I, EOM_TMP, 0, 0, 1, "L1T2_ov");
 
335
    dpd_buf4_init(&T2, CC_TAMPS, 0, 0, 5, 2, 7, 0, "tijab");
 
336
    dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "Lia");
 
337
    dpd_dot24(&L1, &T2, &I, 0, 0, 1.0, 0.0);
 
338
    dpd_file2_close(&L1);
 
339
    dpd_buf4_close(&T2);
 
340
    dpd_buf4_init(&T2, CC_TAMPS, 0, 0, 5, 0, 5, 0, "tiJaB");
 
341
    dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "LIA");
 
342
    dpd_dot24(&L1, &T2, &I, 0, 0, 1.0, 1.0);
 
343
    dpd_file2_close(&L1);
 
344
    dpd_buf4_close(&T2);
 
345
    dpd_file2_close(&I);
 
346
  }
 
347
  else {
 
348
    dpd_file2_init(&I, EOM_TMP, 0, 2, 3, "L1T2_ov");
 
349
    dpd_buf4_init(&T2, CC_TAMPS, 0, 10, 15, 12, 17, 0, "tijab");
 
350
    dpd_file2_init(&L1, CC_GLG, 0, 2, 3, "Lia");
 
351
    dpd_dot24(&L1, &T2, &I, 0, 0, 1.0, 0.0);
 
352
    dpd_file2_close(&L1);
 
353
    dpd_buf4_close(&T2);
 
354
    dpd_buf4_init(&T2, CC_TAMPS, 0, 23, 29, 23, 29, 0, "tiJaB");
 
355
    dpd_file2_init(&L1, CC_GLG, 0, 0, 1, "LIA");
 
356
    dpd_dot24(&L1, &T2, &I, 0, 0, 1.0, 1.0);
 
357
    dpd_file2_close(&L1);
 
358
    dpd_buf4_close(&T2);
 
359
    dpd_file2_close(&I);
 
360
  }
 
361
 
 
362
  /* LR2_OO(I,J)  = 0.5 * LIMEF * RJMEF + LImEf * RJmEf */
 
363
 
 
364
  if (rohf) {
 
365
    dpd_file2_init(&I, EOM_TMP, S.irrep, 0, 0, "LR2_OO");
 
366
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 7, 2, 7, 0, "LIJAB");
 
367
    dpd_buf4_init(&R2, CC_GR, S.irrep, 0, 7, 2, 7, 0, "RIJAB");
 
368
    dpd_contract442(&L2, &R2, &I, 0, 0, 1.0, 0.0);
 
369
    dpd_buf4_close(&R2);
 
370
    dpd_buf4_close(&L2);
 
371
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 5, 0, 5, 0, "LIjAb");
 
372
    dpd_buf4_init(&R2, CC_GR, S.irrep, 0, 5, 0, 5, 0, "RIjAb");
 
373
    dpd_contract442(&L2, &R2, &I, 0, 0, 1.0, 1.0);
 
374
    dpd_buf4_close(&R2);
 
375
    dpd_buf4_close(&L2);
 
376
    dpd_file2_close(&I);
 
377
  }
 
378
  else {
 
379
    dpd_file2_init(&I, EOM_TMP, S.irrep, 0, 0, "LR2_OO");
 
380
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 7, 2, 7, 0, "LIJAB");
 
381
    dpd_buf4_init(&R2, CC_GR, S.irrep, 0, 7, 2, 7, 0, "RIJAB");
 
382
    dpd_contract442(&L2, &R2, &I, 0, 0, 1.0, 0.0);
 
383
    dpd_buf4_close(&R2);
 
384
    dpd_buf4_close(&L2);
 
385
    dpd_buf4_init(&L2, CC_GLG, 0, 22, 28, 22, 28, 0, "LIjAb");
 
386
    dpd_buf4_init(&R2, CC_GR, S.irrep, 22, 28, 22, 28, 0, "RIjAb");
 
387
    dpd_contract442(&L2, &R2, &I, 0, 0, 1.0, 1.0);
 
388
    dpd_buf4_close(&R2);
 
389
    dpd_buf4_close(&L2);
 
390
    dpd_file2_close(&I);
 
391
  }
 
392
 
 
393
  /* LR2_oo(i,j)  = 0.5 * Limef * Rjmef + LiMeF * RjMeF */
 
394
 
 
395
  if (rohf) {
 
396
    dpd_file2_init(&I, EOM_TMP, S.irrep, 0, 0, "LR2_oo");
 
397
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 7, 2, 7, 0, "Lijab");
 
398
    dpd_buf4_init(&R2, CC_GR, S.irrep, 0, 7, 2, 7, 0, "Rijab");
 
399
    dpd_contract442(&L2, &R2, &I, 0, 0, 1.0, 0.0);
 
400
    dpd_buf4_close(&R2);
 
401
    dpd_buf4_close(&L2);
 
402
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 5, 0, 5, 0, "LiJaB");
 
403
    dpd_buf4_init(&R2, CC_GR, S.irrep, 0, 5, 0, 5, 0, "RiJaB");
 
404
    dpd_contract442(&L2, &R2, &I, 0, 0, 1.0, 1.0);
 
405
    dpd_buf4_close(&R2);
 
406
    dpd_buf4_close(&L2);
 
407
    dpd_file2_close(&I);
 
408
  }
 
409
  else {
 
410
    dpd_file2_init(&I, EOM_TMP, S.irrep, 2, 2, "LR2_oo");
 
411
    dpd_buf4_init(&L2, CC_GLG, 0, 10, 17, 12, 17, 0, "Lijab");
 
412
    dpd_buf4_init(&R2, CC_GR, S.irrep, 10, 17, 12, 17, 0, "Rijab");
 
413
    dpd_contract442(&L2, &R2, &I, 0, 0, 1.0, 0.0);
 
414
    dpd_buf4_close(&R2);
 
415
    dpd_buf4_close(&L2);
 
416
    dpd_buf4_init(&L2, CC_GLG, 0, 23, 29, 23, 29, 0, "LiJaB");
 
417
    dpd_buf4_init(&R2, CC_GR, S.irrep, 23, 29, 23, 29, 0, "RiJaB");
 
418
    dpd_contract442(&L2, &R2, &I, 0, 0, 1.0, 1.0);
 
419
    dpd_buf4_close(&R2);
 
420
    dpd_buf4_close(&L2);
 
421
    dpd_file2_close(&I);
 
422
  }
 
423
 
 
424
  /* LR2_VV(A,B) = 0.5 * LMNEA * RMNEB + LmNeA * RmNeB */
 
425
 
 
426
  if (rohf) {
 
427
    dpd_file2_init(&I, EOM_TMP, S.irrep, 1, 1, "LR2_VV");
 
428
    dpd_buf4_init(&L2, CC_GLG, 0, 2, 5, 2, 7, 0, "LIJAB");
 
429
    dpd_buf4_init(&R2, CC_GR, S.irrep, 2, 5, 2, 7, 0, "RIJAB");
 
430
    dpd_contract442(&L2, &R2, &I, 3, 3, 1.0, 0.0);
 
431
    dpd_buf4_close(&R2);
 
432
    dpd_buf4_close(&L2);
 
433
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 5, 0, 5, 0, "LiJaB");
 
434
    dpd_buf4_init(&R2, CC_GR, S.irrep, 0, 5, 0, 5, 0, "RiJaB");
 
435
    dpd_contract442(&L2, &R2, &I, 3, 3, 1.0, 1.0);
 
436
    dpd_buf4_close(&R2);
 
437
    dpd_buf4_close(&L2);
 
438
    dpd_file2_close(&I);
 
439
  }
 
440
  else {
 
441
    dpd_file2_init(&I, EOM_TMP, S.irrep, 1, 1, "LR2_VV");
 
442
    dpd_buf4_init(&L2, CC_GLG, 0, 2, 5, 2, 7, 0, "LIJAB");
 
443
    dpd_buf4_init(&R2, CC_GR, S.irrep, 2, 5, 2, 7, 0, "RIJAB");
 
444
    dpd_contract442(&L2, &R2, &I, 3, 3, 1.0, 0.0);
 
445
    dpd_buf4_close(&R2);
 
446
    dpd_buf4_close(&L2);
 
447
    dpd_buf4_init(&L2, CC_GLG, 0, 23, 29, 23, 29, 0, "LiJaB");
 
448
    dpd_buf4_init(&R2, CC_GR, S.irrep, 23, 29, 23, 29, 0, "RiJaB");
 
449
    dpd_contract442(&L2, &R2, &I, 3, 3, 1.0, 1.0);
 
450
    dpd_buf4_close(&R2);
 
451
    dpd_buf4_close(&L2);
 
452
    dpd_file2_close(&I);
 
453
  }
 
454
 
 
455
  /* LR2_vv(a,b) = 0.5 * Lmnea * Rmneb + LMnEa * RMnEb */
 
456
 
 
457
  if (rohf) {
 
458
    dpd_file2_init(&I, EOM_TMP, S.irrep, 1, 1, "LR2_vv");
 
459
    dpd_buf4_init(&L2, CC_GLG, 0, 2, 5, 2, 7, 0, "Lijab");
 
460
    dpd_buf4_init(&R2, CC_GR, S.irrep, 2, 5, 2, 7, 0, "Rijab");
 
461
    dpd_contract442(&L2, &R2, &I, 3, 3, 1.0, 0.0);
 
462
    dpd_buf4_close(&R2);
 
463
    dpd_buf4_close(&L2);
 
464
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 5, 0, 5, 0, "LIjAb");
 
465
    dpd_buf4_init(&R2, CC_GR, S.irrep, 0, 5, 0, 5, 0, "RIjAb");
 
466
    dpd_contract442(&L2, &R2, &I, 3, 3, 1.0, 1.0);
 
467
    dpd_buf4_close(&R2);
 
468
    dpd_buf4_close(&L2);
 
469
    dpd_file2_close(&I);
 
470
  }
 
471
  else {
 
472
    dpd_file2_init(&I, EOM_TMP, S.irrep, 3, 3, "LR2_vv");
 
473
    dpd_buf4_init(&L2, CC_GLG, 0, 12, 15, 12, 17, 0, "Lijab");
 
474
    dpd_buf4_init(&R2, CC_GR, S.irrep, 12, 15, 12, 17, 0, "Rijab");
 
475
    dpd_contract442(&L2, &R2, &I, 3, 3, 1.0, 0.0);
 
476
    dpd_buf4_close(&R2);
 
477
    dpd_buf4_close(&L2);
 
478
    dpd_buf4_init(&L2, CC_GLG, 0, 22, 28, 22, 28, 0, "LIjAb");
 
479
    dpd_buf4_init(&R2, CC_GR, S.irrep, 22, 28, 22, 28, 0, "RIjAb");
 
480
    dpd_contract442(&L2, &R2, &I, 3, 3, 1.0, 1.0);
 
481
    dpd_buf4_close(&R2);
 
482
    dpd_buf4_close(&L2);
 
483
    dpd_file2_close(&I);
 
484
  }
 
485
 
 
486
  /* LT2_OO(I,J) = 0.5 * LIMEF * TJMEF + LImEf * TJmEf */
 
487
 
 
488
  if (rohf) {
 
489
    dpd_file2_init(&I, EOM_TMP, 0, 0, 0, "LT2_OO");
 
490
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 7, 2, 7, 0, "LIJAB");
 
491
    dpd_buf4_init(&T2, CC_TAMPS, 0, 0, 7, 2, 7, 0, "tIJAB");
 
492
    dpd_contract442(&L2, &T2, &I, 0, 0, 1.0, 0.0);
 
493
    dpd_buf4_close(&T2);
 
494
    dpd_buf4_close(&L2);
 
495
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 5, 0, 5, 0, "LIjAb");
 
496
    dpd_buf4_init(&T2, CC_TAMPS, 0, 0, 5, 0, 5, 0, "tIjAb");
 
497
    dpd_contract442(&L2, &T2, &I, 0, 0, 1.0, 1.0);
 
498
    dpd_buf4_close(&T2);
 
499
    dpd_buf4_close(&L2);
 
500
    dpd_file2_close(&I);
 
501
  }
 
502
  else {
 
503
    dpd_file2_init(&I, EOM_TMP, 0, 0, 0, "LT2_OO");
 
504
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 7, 2, 7, 0, "LIJAB");
 
505
    dpd_buf4_init(&T2, CC_TAMPS, 0, 0, 7, 2, 7, 0, "tIJAB");
 
506
    dpd_contract442(&L2, &T2, &I, 0, 0, 1.0, 0.0);
 
507
    dpd_buf4_close(&T2);
 
508
    dpd_buf4_close(&L2);
 
509
    dpd_buf4_init(&L2, CC_GLG, 0, 22, 28, 22, 28, 0, "LIjAb");
 
510
    dpd_buf4_init(&T2, CC_TAMPS, 0, 22, 28, 22, 28, 0, "tIjAb");
 
511
    dpd_contract442(&L2, &T2, &I, 0, 0, 1.0, 1.0);
 
512
    dpd_buf4_close(&T2);
 
513
    dpd_buf4_close(&L2);
 
514
    dpd_file2_close(&I);
 
515
  }
 
516
  
 
517
  /* LT2_oo(i,j) = 0.5 * Limef * Tjmef + LiMeF * TjMeF */
 
518
 
 
519
  if (rohf) {
 
520
    dpd_file2_init(&I, EOM_TMP, 0, 0, 0, "LT2_oo");
 
521
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 7, 2, 7, 0, "Lijab");
 
522
    dpd_buf4_init(&T2, CC_TAMPS, 0, 0, 7, 2, 7, 0, "tijab");
 
523
    dpd_contract442(&L2, &T2, &I, 0, 0, 1.0, 0.0);
 
524
    dpd_buf4_close(&T2);
 
525
    dpd_buf4_close(&L2);
 
526
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 5, 0, 5, 0, "LiJaB");
 
527
    dpd_buf4_init(&T2, CC_TAMPS, 0, 0, 5, 0, 5, 0, "tiJaB");
 
528
    dpd_contract442(&L2, &T2, &I, 0, 0, 1.0, 1.0);
 
529
    dpd_buf4_close(&T2);
 
530
    dpd_buf4_close(&L2);
 
531
    dpd_file2_close(&I);
 
532
  }
 
533
  else {
 
534
    dpd_file2_init(&I, EOM_TMP, 0, 2, 2, "LT2_oo");
 
535
    dpd_buf4_init(&L2, CC_GLG, 0, 10, 17, 12, 17, 0, "Lijab");
 
536
    dpd_buf4_init(&T2, CC_TAMPS, 0, 10, 17, 12, 17, 0, "tijab");
 
537
    dpd_contract442(&L2, &T2, &I, 0, 0, 1.0, 0.0);
 
538
    dpd_buf4_close(&T2);
 
539
    dpd_buf4_close(&L2);
 
540
    dpd_buf4_init(&L2, CC_GLG, 0, 23, 29, 23, 29, 0, "LiJaB");
 
541
    dpd_buf4_init(&T2, CC_TAMPS, 0, 23, 29, 23, 29, 0, "tiJaB");
 
542
    dpd_contract442(&L2, &T2, &I, 0, 0, 1.0, 1.0);
 
543
    dpd_buf4_close(&T2);
 
544
    dpd_buf4_close(&L2);
 
545
    dpd_file2_close(&I);
 
546
  }
 
547
 
 
548
  /* LT2_VV(A,B) = 0.5 * LMNEA * TMNEB + LmNeA * TmNeB */
 
549
 
 
550
  if (rohf) {
 
551
    dpd_file2_init(&I, EOM_TMP, 0, 1, 1, "LT2_VV");
 
552
    dpd_buf4_init(&L2, CC_GLG, 0, 2, 5, 2, 7, 0, "LIJAB");
 
553
    dpd_buf4_init(&T2, CC_TAMPS, 0, 2, 5, 2, 7, 0, "tIJAB");
 
554
    dpd_contract442(&L2, &T2, &I, 3, 3, 1.0, 0.0);
 
555
    dpd_buf4_close(&T2);
 
556
    dpd_buf4_close(&L2);
 
557
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 5, 0, 5, 0, "LiJaB");
 
558
    dpd_buf4_init(&T2, CC_TAMPS, 0, 0, 5, 0, 5, 0, "tiJaB");
 
559
    dpd_contract442(&L2, &T2, &I, 3, 3, 1.0, 1.0);
 
560
    dpd_buf4_close(&T2);
 
561
    dpd_buf4_close(&L2);
 
562
    dpd_file2_close(&I);
 
563
  }
 
564
  else {
 
565
    dpd_file2_init(&I, EOM_TMP, 0, 1, 1, "LT2_VV");
 
566
    dpd_buf4_init(&L2, CC_GLG, 0, 2, 5, 2, 7, 0, "LIJAB");
 
567
    dpd_buf4_init(&T2, CC_TAMPS, 0, 2, 5, 2, 7, 0, "tIJAB");
 
568
    dpd_contract442(&L2, &T2, &I, 3, 3, 1.0, 0.0);
 
569
    dpd_buf4_close(&T2);
 
570
    dpd_buf4_close(&L2);
 
571
    dpd_buf4_init(&L2, CC_GLG, 0, 23, 29, 23, 29, 0, "LiJaB");
 
572
    dpd_buf4_init(&T2, CC_TAMPS, 0, 23, 29, 23, 29, 0, "tiJaB");
 
573
    dpd_contract442(&L2, &T2, &I, 3, 3, 1.0, 1.0);
 
574
    dpd_buf4_close(&T2);
 
575
    dpd_buf4_close(&L2);
 
576
    dpd_file2_close(&I);
 
577
  }
 
578
 
 
579
  /* LT2_vv(a,b) = 0.5 * Lmnea * Tmneb + LMnEa * TMnEb */
 
580
 
 
581
  if (rohf) {
 
582
    dpd_file2_init(&I, EOM_TMP, 0, 1, 1, "LT2_vv");
 
583
    dpd_buf4_init(&L2, CC_GLG, 0, 2, 5, 2, 7, 0, "Lijab");
 
584
    dpd_buf4_init(&T2, CC_TAMPS, 0, 2, 5, 2, 7, 0, "tijab");
 
585
    dpd_contract442(&L2, &T2, &I, 3, 3, 1.0, 0.0);
 
586
    dpd_buf4_close(&T2);
 
587
    dpd_buf4_close(&L2);
 
588
    dpd_buf4_init(&L2, CC_GLG, 0, 0, 5, 0, 5, 0, "LIjAb");
 
589
    dpd_buf4_init(&T2, CC_TAMPS, 0, 0, 5, 0, 5, 0, "tIjAb");
 
590
    dpd_contract442(&L2, &T2, &I, 3, 3, 1.0, 1.0);
 
591
    dpd_buf4_close(&T2);
 
592
    dpd_buf4_close(&L2);
 
593
    dpd_file2_close(&I);
 
594
  }
 
595
  else {
 
596
    dpd_file2_init(&I, EOM_TMP, 0, 3, 3, "LT2_vv");
 
597
    dpd_buf4_init(&L2, CC_GLG, 0, 12, 15, 12, 17, 0, "Lijab");
 
598
    dpd_buf4_init(&T2, CC_TAMPS, 0, 12, 15, 12, 17, 0, "tijab");
 
599
    dpd_contract442(&L2, &T2, &I, 3, 3, 1.0, 0.0);
 
600
    dpd_buf4_close(&T2);
 
601
    dpd_buf4_close(&L2);
 
602
    dpd_buf4_init(&L2, CC_GLG, 0, 22, 28, 22, 28, 0, "LIjAb");
 
603
    dpd_buf4_init(&T2, CC_TAMPS, 0, 22, 28, 22, 28, 0, "tIjAb");
 
604
    dpd_contract442(&L2, &T2, &I, 3, 3, 1.0, 1.0);
 
605
    dpd_buf4_close(&T2);
 
606
    dpd_buf4_close(&L2);
 
607
    dpd_file2_close(&I);
 
608
  }
 
609
 
 
610
  /* LR_OO = LR1_OO + LR2_OO */
 
611
  /* LR_oo = LR1_oo + LR2_oo */
 
612
  /* LR_VV = LR1_VV + LR2_VV */
 
613
  /* LR_vv = LR1_vv + LR2_vv */
 
614
 
 
615
  dpd_file2_init(&I, EOM_TMP, S.irrep, 0, 0, "LR_OO");
 
616
  dpd_file2_init(&LR1, EOM_TMP, S.irrep, 0, 0, "LR1_OO");
 
617
  dpd_file2_init(&LR2, EOM_TMP, S.irrep, 0, 0, "LR2_OO");
 
618
  dpd_file2_axpbycz(&LR1, &LR2, &I, 1.0, 1.0, 0.0);
 
619
  dpd_file2_close(&LR2);
 
620
  dpd_file2_close(&LR1);
 
621
  dpd_file2_close(&I);
 
622
 
 
623
  if (rohf) {
 
624
    dpd_file2_init(&I, EOM_TMP, S.irrep, 0, 0, "LR_oo");
 
625
    dpd_file2_init(&LR1, EOM_TMP, S.irrep, 0, 0, "LR1_oo");
 
626
    dpd_file2_init(&LR2, EOM_TMP, S.irrep, 0, 0, "LR2_oo");
 
627
    dpd_file2_axpbycz(&LR1, &LR2, &I, 1.0, 1.0, 0.0);
 
628
    dpd_file2_close(&LR2);
 
629
    dpd_file2_close(&LR1);
 
630
    dpd_file2_close(&I);
 
631
  }
 
632
  else {
 
633
    dpd_file2_init(&I, EOM_TMP, S.irrep, 2, 2, "LR_oo");
 
634
    dpd_file2_init(&LR1, EOM_TMP, S.irrep, 2, 2, "LR1_oo");
 
635
    dpd_file2_init(&LR2, EOM_TMP, S.irrep, 2, 2, "LR2_oo");
 
636
    dpd_file2_axpbycz(&LR1, &LR2, &I, 1.0, 1.0, 0.0);
 
637
    dpd_file2_close(&LR2);
 
638
    dpd_file2_close(&LR1);
 
639
    dpd_file2_close(&I);
 
640
  }
 
641
 
 
642
  dpd_file2_init(&I, EOM_TMP, S.irrep, 1, 1, "LR_VV");
 
643
  dpd_file2_init(&LR1, EOM_TMP, S.irrep, 1, 1, "LR1_VV");
 
644
  dpd_file2_init(&LR2, EOM_TMP, S.irrep, 1, 1, "LR2_VV");
 
645
  dpd_file2_axpbycz(&LR1, &LR2, &I, 1.0, 1.0, 0.0);
 
646
  dpd_file2_close(&LR2);
 
647
  dpd_file2_close(&LR1);
 
648
  dpd_file2_close(&I);
 
649
 
 
650
  if (rohf) {
 
651
    dpd_file2_init(&I, EOM_TMP, S.irrep, 1, 1, "LR_vv");
 
652
    dpd_file2_init(&LR1, EOM_TMP, S.irrep, 1, 1, "LR1_vv");
 
653
    dpd_file2_init(&LR2, EOM_TMP, S.irrep, 1, 1, "LR2_vv");
 
654
    dpd_file2_axpbycz(&LR1, &LR2, &I, 1.0, 1.0, 0.0);
 
655
    dpd_file2_close(&LR2);
 
656
    dpd_file2_close(&LR1);
 
657
    dpd_file2_close(&I);
 
658
  }
 
659
  else {
 
660
    dpd_file2_init(&I, EOM_TMP, S.irrep, 3, 3, "LR_vv");
 
661
    dpd_file2_init(&LR1, EOM_TMP, S.irrep, 3, 3, "LR1_vv");
 
662
    dpd_file2_init(&LR2, EOM_TMP, S.irrep, 3, 3, "LR2_vv");
 
663
    dpd_file2_axpbycz(&LR1, &LR2, &I, 1.0, 1.0, 0.0);
 
664
    dpd_file2_close(&LR2);
 
665
    dpd_file2_close(&LR1);
 
666
    dpd_file2_close(&I);
 
667
  }
 
668
 
 
669
  /* LT_OO = LT1_OO + LT2_OO */
 
670
  /* LT_oo = LT1_oo + LT2_oo */
 
671
 
 
672
  dpd_file2_init(&I, EOM_TMP, 0, 0, 0, "LT_OO");
 
673
  dpd_file2_init(&LT1, EOM_TMP, 0, 0, 0, "LT1_OO");
 
674
  dpd_file2_init(&LT2, EOM_TMP, 0, 0, 0, "LT2_OO");
 
675
  dpd_file2_axpbycz(&LT1, &LT2, &I, 1.0, 1.0, 0.0);
 
676
  dpd_file2_close(&LT2);
 
677
  dpd_file2_close(&LT1);
 
678
  dpd_file2_close(&I);
 
679
 
 
680
  if (rohf) {
 
681
    dpd_file2_init(&I, EOM_TMP, 0, 0, 0, "LT_oo");
 
682
    dpd_file2_init(&LT1, EOM_TMP, 0, 0, 0, "LT1_oo");
 
683
    dpd_file2_init(&LT2, EOM_TMP, 0, 0, 0, "LT2_oo");
 
684
    dpd_file2_axpbycz(&LT1, &LT2, &I, 1.0, 1.0, 0.0);
 
685
    dpd_file2_close(&LT2);
 
686
    dpd_file2_close(&LT1);
 
687
    dpd_file2_close(&I);
 
688
  }
 
689
  else {
 
690
    dpd_file2_init(&I, EOM_TMP, 0, 2, 2, "LT_oo");
 
691
    dpd_file2_init(&LT1, EOM_TMP, 0, 2, 2, "LT1_oo");
 
692
    dpd_file2_init(&LT2, EOM_TMP, 0, 2, 2, "LT2_oo");
 
693
    dpd_file2_axpbycz(&LT1, &LT2, &I, 1.0, 1.0, 0.0);
 
694
    dpd_file2_close(&LT2);
 
695
    dpd_file2_close(&LT1);
 
696
    dpd_file2_close(&I);
 
697
  }
 
698
 
 
699
  /* LT_VV = LT1_VV + LT2_VV */
 
700
  /* LT_vv = LT1_vv + LT2_vv */
 
701
 
 
702
  dpd_file2_init(&I, EOM_TMP, 0, 1, 1, "LT_VV");
 
703
  dpd_file2_init(&LT1, EOM_TMP, 0, 1, 1, "LT1_VV");
 
704
  dpd_file2_init(&LT2, EOM_TMP, 0, 1, 1, "LT2_VV");
 
705
  dpd_file2_axpbycz(&LT1, &LT2, &I, 1.0, 1.0, 0.0);
 
706
  dpd_file2_close(&LT2);
 
707
  dpd_file2_close(&LT1);
 
708
  dpd_file2_close(&I);
 
709
 
 
710
  if (rohf) {
 
711
    dpd_file2_init(&I, EOM_TMP, 0, 1, 1, "LT_vv");
 
712
    dpd_file2_init(&LT1, EOM_TMP, 0, 1, 1, "LT1_vv");
 
713
    dpd_file2_init(&LT2, EOM_TMP, 0, 1, 1, "LT2_vv");
 
714
    dpd_file2_axpbycz(&LT1, &LT2, &I, 1.0, 1.0, 0.0);
 
715
    dpd_file2_close(&LT2);
 
716
    dpd_file2_close(&LT1);
 
717
    dpd_file2_close(&I);
 
718
  }
 
719
  else {
 
720
    dpd_file2_init(&I, EOM_TMP, 0, 3, 3, "LT_vv");
 
721
    dpd_file2_init(&LT1, EOM_TMP, 0, 3, 3, "LT1_vv");
 
722
    dpd_file2_init(&LT2, EOM_TMP, 0, 3, 3, "LT2_vv");
 
723
    dpd_file2_axpbycz(&LT1, &LT2, &I, 1.0, 1.0, 0.0);
 
724
    dpd_file2_close(&LT2);
 
725
    dpd_file2_close(&LT1);
 
726
    dpd_file2_close(&I);
 
727
  }
 
728
 
 
729
  return;
 
730
}
 
731
 
 
732
}} // namespace psi::ccdensity