← Back to branch summary

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

« back to all changes in this revision

Viewing changes to src/bin/ccdensity/x_oe_intermediates_rhf.c

  • Committer: Bazaar Package Importer
  • Author(s): Michael Banck
  • Date: 2008-06-07 16:49:57 UTC
  • mfrom: (2.1.2 hardy)
  • Revision ID: james.westby@ubuntu.com-20080607164957-8pifvb133yjlkagn
Tags: 3.3.0-3
http://bugs.debian.org/464867
* debian/rules (DEB_MAKE_CHECK_TARGET): Do not abort test suite on
  failures.
* debian/rules (DEB_CONFIGURE_EXTRA_FLAGS): Set ${bindir} to /usr/lib/psi.
* debian/rules (install/psi3): Move psi3 file to /usr/bin.
* debian/patches/07_464867_move_executables.dpatch: New patch, add
  /usr/lib/psi to the $PATH, so that the moved executables are found.
  (closes: #464867)
* debian/patches/00list: Adjusted.

Show diffs side-by-side

added added

removed removed

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