← Back to branch summary

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

« back to all changes in this revision

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

  • Committer: Bazaar Package Importer
  • Author(s): Michael Banck
  • Date: 2006-09-10 14:01:33 UTC
  • Revision ID: james.westby@ubuntu.com-20060910140133-ib2j86trekykfsfv
Tags: upstream-3.2.3
Import upstream version 3.2.3

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/bin/ccdensity/V.c
 
1
#include <stdio.h>
 
2
#include <libdpd/dpd.h>
 
3
#define EXTERN
 
4
#include "globals.h"
 
5
 
 
6
void V_build(void)
 
7
{
 
8
  dpdbuf4 V, L, T;
 
9
 
 
10
  if(params.ref == 0 || params.ref == 1) { /** RHF/ROHF **/
 
11
 
 
12
    dpd_buf4_init(&V, CC_MISC, 0, 2, 2, 2, 2, 0, "VMNIJ");
 
13
    dpd_buf4_init(&T, CC_TAMPS, 0, 2, 7, 2, 7, 0, "tauIJAB");
 
14
    dpd_buf4_init(&L, CC_LAMPS, 0, 2, 7, 2, 7, 0, "LIJAB");
 
15
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 0.0);
 
16
    dpd_buf4_close(&L);
 
17
    dpd_buf4_close(&T);
 
18
    dpd_buf4_close(&V);
 
19
 
 
20
    dpd_buf4_init(&V, CC_MISC, 0, 2, 2, 2, 2, 0, "Vmnij");
 
21
    dpd_buf4_init(&T, CC_TAMPS, 0, 2, 7, 2, 7, 0, "tauijab");
 
22
    dpd_buf4_init(&L, CC_LAMPS, 0, 2, 7, 2, 7, 0, "Lijab");
 
23
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 0.0);
 
24
    dpd_buf4_close(&L);
 
25
    dpd_buf4_close(&T);
 
26
    dpd_buf4_close(&V);
 
27
 
 
28
    dpd_buf4_init(&V, CC_MISC, 0, 0, 0, 0, 0, 0, "VMnIj");
 
29
    dpd_buf4_init(&T, CC_TAMPS, 0, 0, 5, 0, 5, 0, "tauIjAb");
 
30
    dpd_buf4_init(&L, CC_LAMPS, 0, 0, 5, 0, 5, 0, "LIjAb");
 
31
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 0.0);
 
32
    dpd_buf4_close(&L);
 
33
    dpd_buf4_close(&T);
 
34
    dpd_buf4_close(&V);
 
35
 
 
36
    dpd_buf4_init(&V, CC_MISC, 0, 10, 10, 10, 10, 0, "VIAJB");
 
37
    dpd_buf4_init(&T, CC_TAMPS, 0, 10, 10, 10, 10, 0, "tIAJB");
 
38
    dpd_buf4_init(&L, CC_LAMPS, 0, 10, 10, 10, 10, 0, "LIAJB");
 
39
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 0.0);
 
40
    dpd_buf4_close(&L);
 
41
    dpd_buf4_close(&T);
 
42
    dpd_buf4_init(&T, CC_TAMPS, 0, 10, 10, 10, 10, 0, "tIAjb");
 
43
    dpd_buf4_init(&L, CC_LAMPS, 0, 10, 10, 10, 10, 0, "LIAjb");
 
44
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 1.0);
 
45
    dpd_buf4_close(&L);
 
46
    dpd_buf4_close(&T);
 
47
    dpd_buf4_close(&V);
 
48
 
 
49
    dpd_buf4_init(&V, CC_MISC, 0, 10, 10, 10, 10, 0, "Viajb");
 
50
    dpd_buf4_init(&T, CC_TAMPS, 0, 10, 10, 10, 10, 0, "tiajb");
 
51
    dpd_buf4_init(&L, CC_LAMPS, 0, 10, 10, 10, 10, 0, "Liajb");
 
52
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 0.0);
 
53
    dpd_buf4_close(&L);
 
54
    dpd_buf4_close(&T);
 
55
    dpd_buf4_init(&T, CC_TAMPS, 0, 10, 10, 10, 10, 0, "tIAjb");
 
56
    dpd_buf4_init(&L, CC_LAMPS, 0, 10, 10, 10, 10, 0, "LIAjb");
 
57
    dpd_contract444(&T, &L, &V, 1, 1, 1.0, 1.0);
 
58
    dpd_buf4_close(&L);
 
59
    dpd_buf4_close(&T);
 
60
    dpd_buf4_close(&V);
 
61
 
 
62
    dpd_buf4_init(&V, CC_MISC, 0, 10, 10, 10, 10, 0, "VIAjb");
 
63
    dpd_buf4_init(&T, CC_TAMPS, 0, 10, 10, 10, 10, 0, "tIAjb");
 
64
    dpd_buf4_init(&L, CC_LAMPS, 0, 10, 10, 10, 10, 0, "Liajb");
 
65
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 0.0);
 
66
    dpd_buf4_close(&L);
 
67
    dpd_buf4_close(&T);
 
68
    dpd_buf4_init(&T, CC_TAMPS, 0, 10, 10, 10, 10, 0, "tIAJB");
 
69
    dpd_buf4_init(&L, CC_LAMPS, 0, 10, 10, 10, 10, 0, "LIAjb");
 
70
    dpd_contract444(&T, &L, &V, 0, 1, 1.0, 1.0);
 
71
    dpd_buf4_close(&L);
 
72
    dpd_buf4_close(&T);
 
73
    dpd_buf4_close(&V);
 
74
 
 
75
    dpd_buf4_init(&V, CC_MISC, 0, 10, 10, 10, 10, 0, "ViaJB");
 
76
    dpd_buf4_init(&T, CC_TAMPS, 0, 10, 10, 10, 10, 0, "tiaJB");
 
77
    dpd_buf4_init(&L, CC_LAMPS, 0, 10, 10, 10, 10, 0, "LIAJB");
 
78
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 0.0);
 
79
    dpd_buf4_close(&L);
 
80
    dpd_buf4_close(&T);
 
81
    dpd_buf4_init(&T, CC_TAMPS, 0, 10, 10, 10, 10, 0, "tiajb");
 
82
    dpd_buf4_init(&L, CC_LAMPS, 0, 10, 10, 10, 10, 0, "LIAjb");
 
83
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 1.0);
 
84
    dpd_buf4_close(&L);
 
85
    dpd_buf4_close(&T);
 
86
    dpd_buf4_close(&V);
 
87
 
 
88
    dpd_buf4_init(&V, CC_MISC, 0, 10, 10, 10, 10, 0, "ViAjB");
 
89
    dpd_buf4_init(&L, CC_LAMPS, 0, 10, 10, 10, 10, 0, "LIbjA");
 
90
    dpd_buf4_init(&T, CC_TAMPS, 0, 10, 10, 10, 10, 0, "tjAIb");
 
91
    dpd_contract444(&T, &L, &V, 0, 1, 1.0, 0.0);
 
92
    dpd_buf4_close(&T);
 
93
    dpd_buf4_close(&L);
 
94
    dpd_buf4_close(&V);
 
95
 
 
96
    dpd_buf4_init(&V, CC_MISC, 0, 10, 10, 10, 10, 0, "VIaJb");
 
97
    dpd_buf4_init(&L, CC_LAMPS, 0, 10, 10, 10, 10, 0, "LjAIb");
 
98
    dpd_buf4_init(&T, CC_TAMPS, 0, 10, 10, 10, 10, 0, "tIbjA");
 
99
    dpd_contract444(&T, &L, &V, 0, 1, 1.0, 0.0);
 
100
    dpd_buf4_close(&T);
 
101
    dpd_buf4_close(&L);
 
102
    dpd_buf4_close(&V);
 
103
  }
 
104
  else if(params.ref == 2) { /** UHF **/
 
105
    dpd_buf4_init(&V, CC_MISC, 0, 2, 2, 2, 2, 0, "VMNIJ");
 
106
    dpd_buf4_init(&T, CC_TAMPS, 0, 2, 7, 2, 7, 0, "tauIJAB");
 
107
    dpd_buf4_init(&L, CC_LAMPS, 0, 2, 7, 2, 7, 0, "LIJAB");
 
108
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 0.0);
 
109
    dpd_buf4_close(&L);
 
110
    dpd_buf4_close(&T);
 
111
    dpd_buf4_close(&V);
 
112
 
 
113
    dpd_buf4_init(&V, CC_MISC, 0, 12, 12, 12, 12, 0, "Vmnij");
 
114
    dpd_buf4_init(&T, CC_TAMPS, 0, 12, 17, 12, 17, 0, "tauijab");
 
115
    dpd_buf4_init(&L, CC_LAMPS, 0, 12, 17, 12, 17, 0, "Lijab");
 
116
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 0.0);
 
117
    dpd_buf4_close(&L);
 
118
    dpd_buf4_close(&T);
 
119
    dpd_buf4_close(&V);
 
120
 
 
121
    dpd_buf4_init(&V, CC_MISC, 0, 22, 22, 22, 22, 0, "VMnIj");
 
122
    dpd_buf4_init(&T, CC_TAMPS, 0, 22, 28, 22, 28, 0, "tauIjAb");
 
123
    dpd_buf4_init(&L, CC_LAMPS, 0, 22, 28, 22, 28, 0, "LIjAb");
 
124
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 0.0);
 
125
    dpd_buf4_close(&L);
 
126
    dpd_buf4_close(&T);
 
127
    dpd_buf4_close(&V);
 
128
 
 
129
    dpd_buf4_init(&V, CC_MISC, 0, 20, 20, 20, 20, 0, "VIAJB");
 
130
    dpd_buf4_init(&T, CC_TAMPS, 0, 20, 20, 20, 20, 0, "tIAJB");
 
131
    dpd_buf4_init(&L, CC_LAMPS, 0, 20, 20, 20, 20, 0, "LIAJB");
 
132
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 0.0);
 
133
    dpd_buf4_close(&L);
 
134
    dpd_buf4_close(&T);
 
135
    dpd_buf4_init(&T, CC_TAMPS, 0, 20, 30, 20, 30, 0, "tIAjb");
 
136
    dpd_buf4_init(&L, CC_LAMPS, 0, 20, 30, 20, 30, 0, "LIAjb");
 
137
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 1.0);
 
138
    dpd_buf4_close(&L);
 
139
    dpd_buf4_close(&T);
 
140
    dpd_buf4_close(&V);
 
141
 
 
142
    dpd_buf4_init(&V, CC_MISC, 0, 30, 30, 30, 30, 0, "Viajb");
 
143
    dpd_buf4_init(&T, CC_TAMPS, 0, 30, 30, 30, 30, 0, "tiajb");
 
144
    dpd_buf4_init(&L, CC_LAMPS, 0, 30, 30, 30, 30, 0, "Liajb");
 
145
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 0.0);
 
146
    dpd_buf4_close(&L);
 
147
    dpd_buf4_close(&T);
 
148
    dpd_buf4_init(&T, CC_TAMPS, 0, 20, 30, 20, 30, 0, "tIAjb");
 
149
    dpd_buf4_init(&L, CC_LAMPS, 0, 20, 30, 20, 30, 0, "LIAjb");
 
150
    dpd_contract444(&T, &L, &V, 1, 1, 1.0, 1.0);
 
151
    dpd_buf4_close(&L);
 
152
    dpd_buf4_close(&T);
 
153
    dpd_buf4_close(&V);
 
154
 
 
155
    dpd_buf4_init(&V, CC_MISC, 0, 20, 30, 20, 30, 0, "VIAjb");
 
156
    dpd_buf4_init(&T, CC_TAMPS, 0, 20, 30, 20, 30, 0, "tIAjb");
 
157
    dpd_buf4_init(&L, CC_LAMPS, 0, 30, 30, 30, 30, 0, "Liajb");
 
158
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 0.0);
 
159
    dpd_buf4_close(&L);
 
160
    dpd_buf4_close(&T);
 
161
    dpd_buf4_init(&T, CC_TAMPS, 0, 20, 20, 20, 20, 0, "tIAJB");
 
162
    dpd_buf4_init(&L, CC_LAMPS, 0, 20, 30, 20, 30, 0, "LIAjb");
 
163
    dpd_contract444(&T, &L, &V, 0, 1, 1.0, 1.0);
 
164
    dpd_buf4_close(&L);
 
165
    dpd_buf4_close(&T);
 
166
    dpd_buf4_close(&V);
 
167
 
 
168
    dpd_buf4_init(&V, CC_MISC, 0, 30, 20, 30, 20, 0, "ViaJB");
 
169
    dpd_buf4_init(&T, CC_TAMPS, 0, 30, 20, 30, 20, 0, "tiaJB");
 
170
    dpd_buf4_init(&L, CC_LAMPS, 0, 20, 20, 20, 20, 0, "LIAJB");
 
171
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 0.0);
 
172
    dpd_buf4_close(&L);
 
173
    dpd_buf4_close(&T);
 
174
    dpd_buf4_init(&T, CC_TAMPS, 0, 30, 30, 30, 30, 0, "tiajb");
 
175
    dpd_buf4_init(&L, CC_LAMPS, 0, 20, 30, 20, 30, 0, "LIAjb");
 
176
    dpd_contract444(&T, &L, &V, 0, 0, 1.0, 1.0);
 
177
    dpd_buf4_close(&L);
 
178
    dpd_buf4_close(&T);
 
179
    dpd_buf4_close(&V);
 
180
 
 
181
    dpd_buf4_init(&V, CC_MISC, 0, 27, 27, 27, 27, 0, "ViAjB");
 
182
    dpd_buf4_init(&L, CC_LAMPS, 0, 24, 27, 24, 27, 0, "LIbjA");
 
183
    dpd_buf4_init(&T, CC_TAMPS, 0, 27, 24, 27, 24, 0, "tjAIb");
 
184
    dpd_contract444(&T, &L, &V, 0, 1, 1.0, 0.0);
 
185
    dpd_buf4_close(&T);
 
186
    dpd_buf4_close(&L);
 
187
    dpd_buf4_close(&V);
 
188
 
 
189
    dpd_buf4_init(&V, CC_MISC, 0, 24, 24, 24, 24, 0, "VIaJb");
 
190
    dpd_buf4_init(&L, CC_LAMPS, 0, 27, 24, 27, 24, 0, "LjAIb");
 
191
    dpd_buf4_init(&T, CC_TAMPS, 0, 24, 27, 24, 27, 0, "tIbjA");
 
192
    dpd_contract444(&T, &L, &V, 0, 1, 1.0, 0.0);
 
193
    dpd_buf4_close(&T);
 
194
    dpd_buf4_close(&L);
 
195
    dpd_buf4_close(&V);
 
196
  }
 
197
}