← Back to branch summary

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

« back to all changes in this revision

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