← Back to branch summary

~ubuntu-branches/ubuntu/karmic/psicode/karmic

« back to all changes in this revision

Viewing changes to src/bin/cchbar/norm_HET1.c

  • 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/cchbar/norm_HET1.c
1
 
#include <stdio.h>
2
 
#include <math.h>
3
 
#define EXTERN
4
 
#include "globals.h"
5
 
 
6
 
void norm_HET1(void) {
7
 
  int i;
8
 
  double dot;
9
 
  dpdfile2 F;
10
 
  dpdbuf4 W;
11
 
 
12
 
  fprintf(outfile,"Dots of (HeT1)c in names \"CC3 Wxxx\" in CC3_HET1 \n");
13
 
 
14
 
  if (params.ref == 0) { /* RHF */
15
 
    /** Wamef **/
16
 
    dpd_buf4_init(&W, CC3_HET1, 0, 11, 5, 11, 5, 0, "CC3 WAmEf (Am,Ef)");
17
 
    dot = dpd_buf4_dot_self(&W);
18
 
    dpd_buf4_close(&W);
19
 
    fprintf(outfile,"<WAmEf (Am,Ef) | WAmEf (Am,Ef)> = %15.10lf\n", dot);
20
 
 
21
 
    /** Wmnie **/
22
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 10, 0, 10, 0, "CC3 WMnIe (Mn,Ie)");
23
 
    dot = dpd_buf4_dot_self(&W);
24
 
    dpd_buf4_close(&W);
25
 
    fprintf(outfile,"<WMnIe (Mn,Ie) | WMnIe (Mn,Ie)> = %15.10lf\n", dot);
26
 
 
27
 
    /** Wmnij **/
28
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 0, 0, 0, 0, "CC3 WMnIj (Mn,Ij)");
29
 
    dot = dpd_buf4_dot_self(&W);
30
 
    dpd_buf4_close(&W);
31
 
    fprintf(outfile,"<WMnIj (Mn,Ij) | WMnIj (Mn,Ij)> = %15.10lf\n", dot);
32
 
 
33
 
    /** Wmbij **/
34
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 10, 0, 10, 0, "CC3 WMbIj (Ij,Mb)");
35
 
    dot = dpd_buf4_dot_self(&W);
36
 
    dpd_buf4_close(&W);
37
 
    fprintf(outfile,"<WMbIj (Ij,Mb) | WMbIj (Ij,Mb)> = %15.10lf\n", dot);
38
 
 
39
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 0, 10, 0, 0, "CC3 WMbIj (Mb,Ij)");
40
 
    dot = dpd_buf4_dot_self(&W);
41
 
    dpd_buf4_close(&W);
42
 
    fprintf(outfile,"<WMbIj (Mb,Ij) | WMbIj (Mb,Ij)> = %15.10lf\n", dot);
43
 
 
44
 
    /** Wmbej **/
45
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 10, 10, 10, 0, "CC3 WMbEj (ME,jb)");
46
 
    dot = dpd_buf4_dot_self(&W);
47
 
    dpd_buf4_close(&W);
48
 
    fprintf(outfile,"<WMbEj (ME,jb) | WMbEj (ME,jb)> = %15.10lf\n", dot);
49
 
 
50
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 10, 10, 10, 0, "CC3 WMbeJ (Me,Jb)");
51
 
    dot = dpd_buf4_dot_self(&W);
52
 
    dpd_buf4_close(&W);
53
 
    fprintf(outfile,"<WMbeJ (Me,Jb) | WMbeJ (Me,Jb)> = %15.10lf\n", dot);
54
 
 
55
 
    /** WABEI **/
56
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 5, 10, 5, 0, "CC3 WAbEi (iE,bA)");
57
 
    dot = dpd_buf4_dot_self(&W);
58
 
    dpd_buf4_close(&W);
59
 
    fprintf(outfile,"<WAbEi (Ie,bA) | WAbEi (Ie,bA)> = %15.10lf\n", dot);
60
 
 
61
 
  }
62
 
 
63
 
  else if (params.ref == 1) { /* ROHF */
64
 
 
65
 
    fprintf(outfile,"Wamef terms\n");
66
 
    /** WAMEF **/
67
 
    dpd_buf4_init(&W, CC3_HET1, 0, 11, 5, 11, 7, 0, "CC3 WAMEF (AM,E>F)");
68
 
    dot = dpd_buf4_dot_self(&W);
69
 
    dpd_buf4_close(&W);
70
 
    fprintf(outfile,"\t<WAMEF|WAMEF> = %15.10lf\n", dot); 
71
 
 
72
 
    /** Wamef **/
73
 
    dpd_buf4_init(&W, CC3_HET1, 0, 11, 5, 11, 7, 0, "CC3 Wamef (am,e>f)");
74
 
    dot = dpd_buf4_dot_self(&W);
75
 
    dpd_buf4_close(&W);
76
 
    fprintf(outfile,"\t<Wamef|Wamef> = %15.10lf\n", dot); 
77
 
 
78
 
    /** WAmEf **/
79
 
    dpd_buf4_init(&W, CC3_HET1, 0, 11, 5, 11, 5, 0, "CC3 WAmEf (Am,Ef)");
80
 
    dot = dpd_buf4_dot_self(&W);
81
 
    dpd_buf4_close(&W);
82
 
    fprintf(outfile,"\t<WAmEf|WAmEf> = %15.10lf\n", dot); 
83
 
 
84
 
    /** WaMeF **/
85
 
    dpd_buf4_init(&W, CC3_HET1, 0, 11, 5, 11, 5, 0, "CC3 WaMeF (aM,eF)");
86
 
    dot = dpd_buf4_dot_self(&W);
87
 
    dpd_buf4_close(&W);
88
 
    fprintf(outfile,"\t<WaMeF|WaMeF> = %15.10lf\n", dot); 
89
 
 
90
 
    /** WMNIE **/
91
 
    fprintf(outfile,"Wmnie terms\n");
92
 
 
93
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 10, 2, 10, 0, "CC3 WMNIE (M>N,IE)");
94
 
    dot = dpd_buf4_dot_self(&W);
95
 
    dpd_buf4_close(&W);
96
 
    fprintf(outfile,"\t<WMNIE|WMNIE> = %15.10lf\n", dot); 
97
 
 
98
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 10, 2, 10, 0, "CC3 Wmnie (m>n,ie)");
99
 
    dot = dpd_buf4_dot_self(&W);
100
 
    dpd_buf4_close(&W);
101
 
    fprintf(outfile,"\t<Wmnie|Wmnie> = %15.10lf\n", dot); 
102
 
 
103
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 10, 0, 10, 0, "CC3 WMnIe (Mn,Ie)");
104
 
    dot = dpd_buf4_dot_self(&W);
105
 
    dpd_buf4_close(&W);
106
 
    fprintf(outfile,"\t<WMnIe|WMnIe> = %15.10lf\n", dot); 
107
 
 
108
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 10, 0, 10, 0, "CC3 WmNiE (mN,iE)");
109
 
    dot = dpd_buf4_dot_self(&W);
110
 
    dpd_buf4_close(&W);
111
 
    fprintf(outfile,"\t<WmNiE|WmNiE> = %15.10lf\n", dot); 
112
 
 
113
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 11, 2, 11, 0, "CC3 WMNIE (M>N,EI)");
114
 
    dot = dpd_buf4_dot_self(&W);
115
 
    dpd_buf4_close(&W);
116
 
    fprintf(outfile,"\t<WMNIE(M>N,EI)|WMNIE(M>N,EI)> = %15.10lf\n", dot); 
117
 
 
118
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 11, 2, 11, 0, "CC3 Wmnie (m>n,ei)");
119
 
    dot = dpd_buf4_dot_self(&W);
120
 
    dpd_buf4_close(&W);
121
 
    fprintf(outfile,"\t<Wmnie(m>n,ei)|Wmnie(m>n,ei)> = %15.10lf\n", dot); 
122
 
 
123
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 11, 0, 11, 0, "CC3 WMnIe (Mn,eI)");
124
 
    dot = dpd_buf4_dot_self(&W);
125
 
    dpd_buf4_close(&W);
126
 
    fprintf(outfile,"\t<WMnIe(Mn,eI)|WMnIe(Mn,eI)> = %15.10lf\n", dot); 
127
 
 
128
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 11, 0, 11, 0, "CC3 WmNiE (mN,Ei)");
129
 
    dot = dpd_buf4_dot_self(&W);
130
 
    dpd_buf4_close(&W);
131
 
    fprintf(outfile,"\t<WmNiE(mN,Ei)|WmNiE(mN,Ei)> = %15.10lf\n", dot); 
132
 
 
133
 
    /** WMNIJ **/
134
 
    fprintf(outfile,"Doing Wmnij terms.\n");
135
 
 
136
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 0, 2, 2, 0, "CC3 WMNIJ (M>N,I>J)");
137
 
    dot = dpd_buf4_dot_self(&W);
138
 
    dpd_buf4_close(&W);
139
 
    fprintf(outfile,"\t<WMNIJ (M>N,IJ)|WMNIJ> = %15.10lf\n", dot); 
140
 
 
141
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 0, 2, 2, 0, "CC3 Wmnij (m>n,i>j)");
142
 
    dot = dpd_buf4_dot_self(&W);
143
 
    dpd_buf4_close(&W);
144
 
    fprintf(outfile,"\t<Wmnij (m>n,ij)|Wmnij> = %15.10lf\n", dot); 
145
 
 
146
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 0, 0, 0, 0, "CC3 WMnIj (Mn,Ij)");
147
 
    dot = dpd_buf4_dot_self(&W);
148
 
    dpd_buf4_close(&W);
149
 
    fprintf(outfile,"\t<WMnIj (Mn,Ij)|WMnIj> = %15.10lf\n", dot); 
150
 
 
151
 
    fprintf(outfile,"Doing Wmbij terms.\n");
152
 
    /** WMBIJ **/
153
 
 
154
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 0, 10, 2, 0, "CC3 WMBIJ (MB,I>J)");
155
 
    dot = dpd_buf4_dot_self(&W);
156
 
    dpd_buf4_close(&W);
157
 
    fprintf(outfile,"\t<WMBIJ (MB,I>J)|WMBIJ> = %15.10lf\n", dot);
158
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 10, 2, 10, 0, "CC3 WMBIJ (I>J,MB)");
159
 
    dot = dpd_buf4_dot_self(&W);
160
 
    dpd_buf4_close(&W);
161
 
    fprintf(outfile,"\t<WMBIJ (I>J,MB)|WMBIJ> = %15.10lf\n", dot);
162
 
 
163
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 0, 10, 2, 0, "CC3 Wmbij (mb,i>j)");
164
 
    dot = dpd_buf4_dot_self(&W);
165
 
    dpd_buf4_close(&W);
166
 
    fprintf(outfile,"\t<Wmbij (mb,i>j)|Wmbij> = %15.10lf\n", dot);
167
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 10, 2, 10, 0, "CC3 Wmbij (i>j,mb)");
168
 
    dot = dpd_buf4_dot_self(&W);
169
 
    dpd_buf4_close(&W);
170
 
    fprintf(outfile,"\t<Wmbij (i>j,mb)|Wmbij> = %15.10lf\n", dot);
171
 
 
172
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 0, 10, 0, 0, "CC3 WMbIj (Mb,Ij)");
173
 
    dot = dpd_buf4_dot_self(&W);
174
 
    dpd_buf4_close(&W);
175
 
    fprintf(outfile,"\t<WMbIj (Mb,Ij)|WMbIj> = %15.10lf\n", dot);
176
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 0, 10, 0, 0, "CC3 WMbIj (Ij,Mb)");
177
 
    dot = dpd_buf4_dot_self(&W);
178
 
    dpd_buf4_close(&W);
179
 
    fprintf(outfile,"\t<WMbIj (Ij,Mb)|WMbIj> = %15.10lf\n", dot);
180
 
 
181
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 0, 10, 0, 0, "CC3 WmBiJ (mB,iJ)");
182
 
    dot = dpd_buf4_dot_self(&W);
183
 
    dpd_buf4_close(&W);
184
 
    fprintf(outfile,"\t<WmBiJ (mB,iJ)|WmBiJ> = %15.10lf\n", dot);
185
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 0, 10, 0, 0, "CC3 WmBiJ (iJ,mB)");
186
 
    dot = dpd_buf4_dot_self(&W);
187
 
    dpd_buf4_close(&W);
188
 
    fprintf(outfile,"\t<WmBiJ (iJ,mB)|WmBiJ> = %15.10lf\n", dot);
189
 
 
190
 
    /** WMBEJ **/
191
 
    fprintf(outfile,"Doing Wmbej terms.\n");
192
 
 
193
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 10, 10, 10, 0, "CC3 WMBEJ (ME,JB)");
194
 
    dot = dpd_buf4_dot_self(&W);
195
 
    dpd_buf4_close(&W);
196
 
    fprintf(outfile,"\t<WMBEJ (ME,JB)|WMBEJ> = %15.10lf\n", dot);
197
 
 
198
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 10, 10, 10, 0, "CC3 Wmbej (me,jb)");
199
 
    dot = dpd_buf4_dot_self(&W);
200
 
    dpd_buf4_close(&W);
201
 
    fprintf(outfile,"\t<Wmbej (me,jb)|Wmbej> = %15.10lf\n", dot);
202
 
 
203
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 10, 10, 10, 0, "CC3 WMbEj (ME,jb)");
204
 
    dot = dpd_buf4_dot_self(&W);
205
 
    dpd_buf4_close(&W);
206
 
    fprintf(outfile,"\t<WMbEj (ME,jb)|WMbEj> = %15.10lf\n", dot);
207
 
 
208
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 10, 10, 10, 0, "CC3 WmBeJ (me,JB)");
209
 
    dot = dpd_buf4_dot_self(&W);
210
 
    dpd_buf4_close(&W);
211
 
    fprintf(outfile,"\t<WmBeJ (me,JB)|WmBeJ> = %15.10lf\n", dot);
212
 
 
213
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 10, 10, 10, 0, "CC3 WMbeJ (Me,Jb)");
214
 
    dot = dpd_buf4_dot_self(&W);
215
 
    dpd_buf4_close(&W);
216
 
    fprintf(outfile,"\t<WMbeJ (Me,Jb)|WMbeJ> = %15.10lf\n", dot);
217
 
 
218
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 10, 10, 10, 0, "CC3 WmBEj (mE,jB)");
219
 
    dot = dpd_buf4_dot_self(&W);
220
 
    dpd_buf4_close(&W);
221
 
    fprintf(outfile,"\t<WmBEj (mE,jB)|WmBEj> = %15.10lf\n", dot);
222
 
 
223
 
    /** WABEI **/
224
 
 
225
 
    fprintf(outfile,"Doing Wabei terms.\n");
226
 
 
227
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 5, 10, 7, 0, "CC3 WABEI (IE,B>A)");
228
 
    dot = dpd_buf4_dot_self(&W);
229
 
    dpd_buf4_close(&W);
230
 
    fprintf(outfile,"\t<WABEI (IE,B>A) |WABEI> = %15.10lf\n", dot);
231
 
 
232
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 5, 10, 7, 0, "CC3 Wabei (ie,b>a)");
233
 
    dot = dpd_buf4_dot_self(&W);
234
 
    dpd_buf4_close(&W);
235
 
    fprintf(outfile,"\t<Wabei (ie,b>a) |Wabei> = %15.10lf\n", dot);
236
 
 
237
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 5, 10, 5, 0, "CC3 WAbEi (iE,bA)");
238
 
    dot = dpd_buf4_dot_self(&W);
239
 
    dpd_buf4_close(&W);
240
 
    fprintf(outfile,"\t<WAbEi (iE,bA) |WAbEi> = %15.10lf\n", dot);
241
 
 
242
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 5, 10, 5, 0, "CC3 WaBeI (Ie,Ba)");
243
 
    dot = dpd_buf4_dot_self(&W);
244
 
    dpd_buf4_close(&W);
245
 
    fprintf(outfile,"\t<WaBeI (Ie,Ba) |WaBeI> = %15.10lf\n", dot);
246
 
  }
247
 
 
248
 
 
249
 
 
250
 
 
251
 
 
252
 
  else { /******************** UHF */
253
 
 
254
 
    fprintf(outfile,"Wamef terms\n");
255
 
    /** WAMEF **/
256
 
    dpd_buf4_init(&W, CC3_HET1, 0, 21, 5, 21, 7, 0, "CC3 WAMEF (AM,E>F)");
257
 
    dot = dpd_buf4_dot_self(&W);
258
 
    dpd_buf4_close(&W);
259
 
    fprintf(outfile,"\t<WAMEF|WAMEF> = %15.10lf\n", dot); 
260
 
 
261
 
    /** Wamef **/
262
 
    dpd_buf4_init(&W, CC3_HET1, 0, 31, 15, 31, 17, 0, "CC3 Wamef (am,e>f)");
263
 
    dot = dpd_buf4_dot_self(&W);
264
 
    dpd_buf4_close(&W);
265
 
    fprintf(outfile,"\t<Wamef|Wamef> = %15.10lf\n", dot); 
266
 
 
267
 
    /** WAmEf **/
268
 
    dpd_buf4_init(&W, CC3_HET1, 0, 26, 28, 26, 28, 0, "CC3 WAmEf (Am,Ef)");
269
 
    dot = dpd_buf4_dot_self(&W);
270
 
    dpd_buf4_close(&W);
271
 
    fprintf(outfile,"\t<WAmEf|WAmEf> = %15.10lf\n", dot); 
272
 
 
273
 
    /** WaMeF **/
274
 
    dpd_buf4_init(&W, CC3_HET1, 0, 25, 29, 25, 29, 0, "CC3 WaMeF (aM,eF)");
275
 
    dot = dpd_buf4_dot_self(&W);
276
 
    dpd_buf4_close(&W);
277
 
    fprintf(outfile,"\t<WaMeF|WaMeF> = %15.10lf\n", dot); 
278
 
 
279
 
    /** WMNIE **/
280
 
    fprintf(outfile,"Wmnie terms\n");
281
 
 
282
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 20, 2, 20, 0, "CC3 WMNIE (M>N,IE)");
283
 
    dot = dpd_buf4_dot_self(&W);
284
 
    dpd_buf4_close(&W);
285
 
    fprintf(outfile,"\t<WMNIE (M>N,IE)|WMNIE> = %15.10lf\n", dot); 
286
 
 
287
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 30, 12, 30, 0, "CC3 Wmnie (m>n,ie)");
288
 
    dot = dpd_buf4_dot_self(&W);
289
 
    dpd_buf4_close(&W);
290
 
    fprintf(outfile,"\t<Wmnie (m>n,ie)|Wmnie> = %15.10lf\n", dot); 
291
 
 
292
 
    dpd_buf4_init(&W, CC3_HET1, 0, 22, 24, 22, 24, 0, "CC3 WMnIe (Mn,Ie)");
293
 
    dot = dpd_buf4_dot_self(&W);
294
 
    dpd_buf4_close(&W);
295
 
    fprintf(outfile,"\t<WMnIe (Mn,Ie)|WMnIe> = %15.10lf\n", dot); 
296
 
 
297
 
    dpd_buf4_init(&W, CC3_HET1, 0, 23, 27, 23, 27, 0, "CC3 WmNiE (mN,iE)");
298
 
    dot = dpd_buf4_dot_self(&W);
299
 
    dpd_buf4_close(&W);
300
 
    fprintf(outfile,"\t<WmNiE (mN,iE)|WmNiE> = %15.10lf\n", dot); 
301
 
 
302
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 21, 2, 21, 0, "CC3 WMNIE (M>N,EI)");
303
 
    dot = dpd_buf4_dot_self(&W);
304
 
    dpd_buf4_close(&W);
305
 
    fprintf(outfile,"\t<WMNIE(M>N,EI)|WMNIE(M>N,EI)> = %15.10lf\n", dot); 
306
 
 
307
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 31, 12, 31, 0, "CC3 Wmnie (m>n,ei)");
308
 
    dot = dpd_buf4_dot_self(&W);
309
 
    dpd_buf4_close(&W);
310
 
    fprintf(outfile,"\t<Wmnie(m>n,ei)|Wmnie(m>n,ei)> = %15.10lf\n", dot); 
311
 
 
312
 
    dpd_buf4_init(&W, CC3_HET1, 0, 22, 25, 22, 25, 0, "CC3 WMnIe (Mn,eI)");
313
 
    dot = dpd_buf4_dot_self(&W);
314
 
    dpd_buf4_close(&W);
315
 
    fprintf(outfile,"\t<WMnIe(Mn,eI)|WMnIe(Mn,eI)> = %15.10lf\n", dot); 
316
 
 
317
 
    dpd_buf4_init(&W, CC3_HET1, 0, 23, 26, 23, 26, 0, "CC3 WmNiE (mN,Ei)");
318
 
    dot = dpd_buf4_dot_self(&W);
319
 
    dpd_buf4_close(&W);
320
 
    fprintf(outfile,"\t<WmNiE(mN,Ei)|WmNiE(mN,Ei)> = %15.10lf\n", dot); 
321
 
 
322
 
    /** WMNIJ **/
323
 
    fprintf(outfile,"Doing Wmnij terms.\n");
324
 
 
325
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 0, 2, 2, 0, "CC3 WMNIJ (M>N,I>J)");
326
 
    dot = dpd_buf4_dot_self(&W);
327
 
    dpd_buf4_close(&W);
328
 
    fprintf(outfile,"\t<WMNIJ (M>N,I>J)|WMNIJ> = %15.10lf\n", dot); 
329
 
 
330
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 10, 12, 12, 0, "CC3 Wmnij (m>n,i>j)");
331
 
    dot = dpd_buf4_dot_self(&W);
332
 
    dpd_buf4_close(&W);
333
 
    fprintf(outfile,"\t<Wmnij (m>n,i>j)|Wmnij> = %15.10lf\n", dot); 
334
 
 
335
 
    dpd_buf4_init(&W, CC3_HET1, 0, 22, 22, 22, 22, 0, "CC3 WMnIj (Mn,Ij)");
336
 
    dot = dpd_buf4_dot_self(&W);
337
 
    dpd_buf4_close(&W);
338
 
    fprintf(outfile,"\t<WMnIj (Mn,Ij)|WMnIj> = %15.10lf\n", dot); 
339
 
 
340
 
    fprintf(outfile,"Doing Wmbij terms.\n");
341
 
    /** WMBIJ **/
342
 
 
343
 
    dpd_buf4_init(&W, CC3_HET1, 0, 20, 0, 20, 2, 0, "CC3 WMBIJ (MB,I>J)");
344
 
    dot = dpd_buf4_dot_self(&W);
345
 
    dpd_buf4_close(&W);
346
 
    fprintf(outfile,"\t<WMBIJ (MB,I>J)|WMBIJ> = %15.10lf\n", dot);
347
 
    dpd_buf4_init(&W, CC3_HET1, 0, 0, 20, 2, 20, 0, "CC3 WMBIJ (I>J,MB)");
348
 
    dot = dpd_buf4_dot_self(&W);
349
 
    dpd_buf4_close(&W);
350
 
    fprintf(outfile,"\t<WMBIJ (I>J,MB)|WMBIJ> = %15.10lf\n", dot);
351
 
 
352
 
    dpd_buf4_init(&W, CC3_HET1, 0, 30, 10, 30, 12, 0, "CC3 Wmbij (mb,i>j)");
353
 
    dot = dpd_buf4_dot_self(&W);
354
 
    dpd_buf4_close(&W);
355
 
    fprintf(outfile,"\t<Wmbij (mb,i>j)|Wmbij> = %15.10lf\n", dot);
356
 
    dpd_buf4_init(&W, CC3_HET1, 0, 10, 30, 12, 30, 0, "CC3 Wmbij (i>j,mb)");
357
 
    dot = dpd_buf4_dot_self(&W);
358
 
    dpd_buf4_close(&W);
359
 
    fprintf(outfile,"\t<Wmbij (i>j,mb)|Wmbij> = %15.10lf\n", dot);
360
 
 
361
 
    dpd_buf4_init(&W, CC3_HET1, 0, 24, 22, 24, 22, 0, "CC3 WMbIj (Mb,Ij)");
362
 
    dot = dpd_buf4_dot_self(&W);
363
 
    dpd_buf4_close(&W);
364
 
    fprintf(outfile,"\t<WMbIj (Mb,Ij)|WMbIj> = %15.10lf\n", dot);
365
 
    dpd_buf4_init(&W, CC3_HET1, 0, 22, 24, 22, 24, 0, "CC3 WMbIj (Ij,Mb)");
366
 
    dot = dpd_buf4_dot_self(&W);
367
 
    dpd_buf4_close(&W);
368
 
    fprintf(outfile,"\t<WMbIj (Ij,Mb)|WMbIj> = %15.10lf\n", dot);
369
 
 
370
 
    dpd_buf4_init(&W, CC3_HET1, 0, 27, 23, 27, 23, 0, "CC3 WmBiJ (mB,iJ)");
371
 
    dot = dpd_buf4_dot_self(&W);
372
 
    dpd_buf4_close(&W);
373
 
    fprintf(outfile,"\t<WmBiJ (mB,iJ)|WmBiJ> = %15.10lf\n", dot);
374
 
    dpd_buf4_init(&W, CC3_HET1, 0, 23, 27, 23, 27, 0, "CC3 WmBiJ (iJ,mB)");
375
 
    dot = dpd_buf4_dot_self(&W);
376
 
    dpd_buf4_close(&W);
377
 
    fprintf(outfile,"\t<WmBiJ (iJ,mB)|WmBiJ> = %15.10lf\n", dot);
378
 
 
379
 
    /** WMBEJ **/
380
 
    fprintf(outfile,"Doing Wmbej terms.\n");
381
 
 
382
 
    dpd_buf4_init(&W, CC3_HET1, 0, 20, 20, 20, 20, 0, "CC3 WMBEJ (ME,JB)");
383
 
    dot = dpd_buf4_dot_self(&W);
384
 
    dpd_buf4_close(&W);
385
 
    fprintf(outfile,"\t<WMBEJ (all ME,JB)|WMBEJ> = %15.10lf\n", dot);
386
 
 
387
 
    dpd_buf4_init(&W, CC3_HET1, 0, 30, 30, 30, 30, 0, "CC3 Wmbej (me,jb)");
388
 
    dot = dpd_buf4_dot_self(&W);
389
 
    dpd_buf4_close(&W);
390
 
    fprintf(outfile,"\t<Wmbej|Wmbej> = %15.10lf\n", dot);
391
 
 
392
 
    dpd_buf4_init(&W, CC3_HET1, 0, 24, 26, 24, 26, 0, "CC3 WMbEj (ME,jb)");
393
 
    dot = dpd_buf4_dot_self(&W);
394
 
    dpd_buf4_close(&W);
395
 
    fprintf(outfile,"\t<WMbEj|WMbEj> = %15.10lf\n", dot);
396
 
 
397
 
    dpd_buf4_init(&W, CC3_HET1, 0, 27, 25, 27, 25, 0, "CC3 WmBeJ (me,JB)");
398
 
    dot = dpd_buf4_dot_self(&W);
399
 
    dpd_buf4_close(&W);
400
 
    fprintf(outfile,"\t<WmBeJ|WmBeJ> = %15.10lf\n", dot);
401
 
 
402
 
    dpd_buf4_init(&W, CC3_HET1, 0, 24, 24, 24, 24, 0, "CC3 WMbeJ (Me,Jb)");
403
 
    dot = dpd_buf4_dot_self(&W);
404
 
    dpd_buf4_close(&W);
405
 
    fprintf(outfile,"\t<WMbeJ|WMbeJ> = %15.10lf\n", dot);
406
 
 
407
 
    dpd_buf4_init(&W, CC3_HET1, 0, 27, 27, 27, 27, 0, "CC3 WmBEj (mE,jB)");
408
 
    dot = dpd_buf4_dot_self(&W);
409
 
    dpd_buf4_close(&W);
410
 
    fprintf(outfile,"\t<WmBEj|WmBEj> = %15.10lf\n", dot);
411
 
 
412
 
    /** WABEI **/
413
 
 
414
 
    fprintf(outfile,"Doing Wabei terms.\n");
415
 
 
416
 
    dpd_buf4_init(&W, CC3_HET1, 0, 20, 5, 20, 7, 0, "CC3 WABEI (IE,B>A)");
417
 
    dot = dpd_buf4_dot_self(&W);
418
 
    dpd_buf4_close(&W);
419
 
    fprintf(outfile,"\t<WABEI (IE,B>A) |WABEI> = %15.10lf\n", dot);
420
 
 
421
 
    dpd_buf4_init(&W, CC3_HET1, 0, 30, 15, 30, 17, 0, "CC3 Wabei (ie,b>a)");
422
 
    dot = dpd_buf4_dot_self(&W);
423
 
    dpd_buf4_close(&W);
424
 
    fprintf(outfile,"\t<Wabei (ie,b>a)|Wabei> = %15.10lf\n", dot);
425
 
 
426
 
    dpd_buf4_init(&W, CC3_HET1, 0, 27, 29, 27, 29, 0, "CC3 WAbEi (iE,bA)");
427
 
    dot = dpd_buf4_dot_self(&W);
428
 
    dpd_buf4_close(&W);
429
 
    fprintf(outfile,"\t<WAbEi (iE,bA)|WAbEi> = %15.10lf\n", dot);
430
 
 
431
 
    dpd_buf4_init(&W, CC3_HET1, 0, 24, 28, 24, 28, 0, "CC3 WaBeI (Ie,Ba)");
432
 
    dot = dpd_buf4_dot_self(&W);
433
 
    dpd_buf4_close(&W);
434
 
    fprintf(outfile,"\t<WAbEi (iE,Ba)|WAbEi> = %15.10lf\n", dot);
435
 
  }
436
 
 
437
 
  return;
438
 
}
439