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

« back to all changes in this revision

Viewing changes to src/bin/cclambda/cache.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:
 
1
#include <stdio.h>
 
2
#include <stdlib.h>
 
3
#include <libciomr/libciomr.h>
 
4
#include <ccfiles.h>
 
5
#include <psifiles.h>
 
6
 
 
7
void cache_abcd_rhf(int **cachelist);
 
8
void cache_iabc_rhf(int **cachelist);
 
9
void cache_ijab_rhf(int **cachelist);
 
10
void cache_iajb_rhf(int **cachelist);
 
11
void cache_ijka_rhf(int **cachelist);
 
12
void cache_ijkl_rhf(int **cachelist);
 
13
 
 
14
void cache_abcd_uhf(int **cachelist);
 
15
void cache_iabc_uhf(int **cachelist);
 
16
void cache_ijab_uhf(int **cachelist);
 
17
void cache_iajb_uhf(int **cachelist);
 
18
void cache_ijka_uhf(int **cachelist);
 
19
void cache_ijkl_uhf(int **cachelist);
 
20
 
 
21
int **cacheprep_uhf(int level, int *cachefiles)
 
22
{
 
23
  int **cachelist;
 
24
 
 
25
  /* The listing of CC files whose entries may be cached */
 
26
  cachefiles[CC_AINTS] = 1;
 
27
  cachefiles[CC_CINTS] = 1;
 
28
  cachefiles[CC_DINTS] = 1;
 
29
  cachefiles[CC_EINTS] = 1;
 
30
  cachefiles[CC_DENOM] = 1;
 
31
  cachefiles[CC_TAMPS] = 1;
 
32
  cachefiles[CC_LAMBDA] = 1;
 
33
  cachefiles[CC_HBAR] = 1;
 
34
 
 
35
  /* The listing of DPD patterns which may be cached */
 
36
  cachelist = init_int_matrix(32,32);
 
37
 
 
38
  if(level == 0) return cachelist;
 
39
  else if(level == 1) {
 
40
 
 
41
      /*** Cache oooo and ooov ***/
 
42
      cache_ijkl_uhf(cachelist);
 
43
      cache_ijka_uhf(cachelist);
 
44
 
 
45
      return cachelist;
 
46
    }
 
47
  else if(level == 2) {
 
48
 
 
49
      /*** Cache oooo, ooov, oovv, and ovov ***/
 
50
      cache_ijkl_uhf(cachelist);
 
51
      cache_ijka_uhf(cachelist);
 
52
      cache_ijab_uhf(cachelist);
 
53
      cache_iajb_uhf(cachelist);
 
54
 
 
55
      return cachelist;
 
56
    }
 
57
  else if(level == 3) {
 
58
 
 
59
      /*** Cache, oooo, oov, oovv, ovov, and ovvv ***/
 
60
 
 
61
      cache_ijkl_uhf(cachelist);
 
62
      cache_ijka_uhf(cachelist);
 
63
      cache_ijab_uhf(cachelist);
 
64
      cache_iajb_uhf(cachelist);
 
65
      cache_iabc_uhf(cachelist);
 
66
 
 
67
      return cachelist;
 
68
    }
 
69
  else if(level == 4) {
 
70
 
 
71
      /*** Cache everything ***/
 
72
      cache_ijkl_uhf(cachelist);
 
73
      cache_ijka_uhf(cachelist);
 
74
      cache_ijab_uhf(cachelist);
 
75
      cache_iajb_uhf(cachelist);
 
76
      cache_iabc_uhf(cachelist);
 
77
      cache_abcd_uhf(cachelist);
 
78
 
 
79
      return cachelist;
 
80
    }
 
81
  else { printf("Error: invalid cache level!\n"); exit(PSI_RETURN_FAILURE); }
 
82
}
 
83
 
 
84
int **cacheprep_rhf(int level, int *cachefiles)
 
85
{
 
86
  int **cachelist;
 
87
 
 
88
  /* The listing of CC files whose entries may be cached */
 
89
  cachefiles[CC_AINTS] = 1;
 
90
  cachefiles[CC_CINTS] = 1;
 
91
  cachefiles[CC_DINTS] = 1;
 
92
  cachefiles[CC_EINTS] = 1;
 
93
  cachefiles[CC_DENOM] = 1;
 
94
  cachefiles[CC_TAMPS] = 1;
 
95
  cachefiles[CC_LAMBDA] = 1;
 
96
  cachefiles[CC_HBAR] = 1;
 
97
 
 
98
  /* The listing of DPD patterns which may be cached */
 
99
  cachelist = init_int_matrix(12,12);
 
100
 
 
101
  if(level == 0) return cachelist;
 
102
  else if(level == 1) {
 
103
 
 
104
      /*** Cache oooo and ooov ***/
 
105
      cache_ijkl_rhf(cachelist);
 
106
      cache_ijka_rhf(cachelist);
 
107
 
 
108
      return cachelist;
 
109
    }
 
110
  else if(level == 2) {
 
111
 
 
112
      /*** Cache oooo, ooov, oovv, and ovov ***/
 
113
      cache_ijkl_rhf(cachelist);
 
114
      cache_ijka_rhf(cachelist);
 
115
      cache_ijab_rhf(cachelist);
 
116
      cache_iajb_rhf(cachelist);
 
117
 
 
118
      return cachelist;
 
119
    }
 
120
  else if(level == 3) {
 
121
 
 
122
      /*** Cache, oooo, oov, oovv, ovov, and ovvv ***/
 
123
 
 
124
      cache_ijkl_rhf(cachelist);
 
125
      cache_ijka_rhf(cachelist);
 
126
      cache_ijab_rhf(cachelist);
 
127
      cache_iajb_rhf(cachelist);
 
128
      cache_iabc_rhf(cachelist);
 
129
 
 
130
      return cachelist;
 
131
    }
 
132
  else if(level == 4) {
 
133
 
 
134
      /*** Cache everything ***/
 
135
      cache_ijkl_rhf(cachelist);
 
136
      cache_ijka_rhf(cachelist);
 
137
      cache_ijab_rhf(cachelist);
 
138
      cache_iajb_rhf(cachelist);
 
139
      cache_iabc_rhf(cachelist);
 
140
      cache_abcd_rhf(cachelist);
 
141
 
 
142
      return cachelist;
 
143
    }
 
144
  else { printf("Error: invalid cache level!\n"); exit(PSI_RETURN_FAILURE); }
 
145
}
 
146
 
 
147
void cache_abcd_uhf(int **cachelist)
 
148
{
 
149
  /* <ab|cd> */
 
150
  cachelist[5][5] = 1;
 
151
  cachelist[5][6] = 1;
 
152
  cachelist[5][7] = 1;
 
153
  cachelist[5][8] = 1;
 
154
  cachelist[5][9] = 1;
 
155
  cachelist[6][5] = 1;
 
156
  cachelist[6][6] = 1;
 
157
  cachelist[6][7] = 1;
 
158
  cachelist[6][8] = 1;
 
159
  cachelist[6][9] = 1;
 
160
  cachelist[7][5] = 1;
 
161
  cachelist[7][6] = 1;
 
162
  cachelist[7][7] = 1;
 
163
  cachelist[7][8] = 1;
 
164
  cachelist[7][9] = 1;
 
165
  cachelist[8][5] = 1;
 
166
  cachelist[8][6] = 1;
 
167
  cachelist[8][7] = 1;
 
168
  cachelist[8][8] = 1;
 
169
  cachelist[8][9] = 1;
 
170
  cachelist[9][5] = 1;
 
171
  cachelist[9][6] = 1;
 
172
  cachelist[9][7] = 1;
 
173
  cachelist[9][8] = 1;
 
174
  cachelist[9][9] = 1;
 
175
  /* <AB|CD> */
 
176
  cachelist[15][15] = 1;
 
177
  cachelist[15][16] = 1;
 
178
  cachelist[15][17] = 1;
 
179
  cachelist[15][18] = 1;
 
180
  cachelist[15][19] = 1;
 
181
  cachelist[16][15] = 1;
 
182
  cachelist[16][16] = 1;
 
183
  cachelist[16][17] = 1;
 
184
  cachelist[16][18] = 1;
 
185
  cachelist[16][19] = 1;
 
186
  cachelist[17][15] = 1;
 
187
  cachelist[17][16] = 1;
 
188
  cachelist[17][17] = 1;
 
189
  cachelist[17][18] = 1;
 
190
  cachelist[17][19] = 1;
 
191
  cachelist[18][15] = 1;
 
192
  cachelist[18][16] = 1;
 
193
  cachelist[18][17] = 1;
 
194
  cachelist[18][18] = 1;
 
195
  cachelist[18][19] = 1;
 
196
  cachelist[19][15] = 1;
 
197
  cachelist[19][16] = 1;
 
198
  cachelist[19][17] = 1;
 
199
  cachelist[19][18] = 1;
 
200
  cachelist[19][19] = 1;
 
201
  /* <Ab|Cd> */
 
202
  cachelist[28][28] = 1;
 
203
  cachelist[29][29] = 1;
 
204
  cachelist[28][29] = 1;
 
205
  cachelist[29][28] = 1;
 
206
}
 
207
 
 
208
void cache_abcd_rhf(int **cachelist)
 
209
{
 
210
  /* <ab|cd> */
 
211
  cachelist[5][5] = 1;
 
212
  cachelist[5][6] = 1;
 
213
  cachelist[5][7] = 1;
 
214
  cachelist[5][8] = 1;
 
215
  cachelist[5][9] = 1;
 
216
  cachelist[6][5] = 1;
 
217
  cachelist[6][6] = 1;
 
218
  cachelist[6][7] = 1;
 
219
  cachelist[6][8] = 1;
 
220
  cachelist[6][9] = 1;
 
221
  cachelist[7][5] = 1;
 
222
  cachelist[7][6] = 1;
 
223
  cachelist[7][7] = 1;
 
224
  cachelist[7][8] = 1;
 
225
  cachelist[7][9] = 1;
 
226
  cachelist[8][5] = 1;
 
227
  cachelist[8][6] = 1;
 
228
  cachelist[8][7] = 1;
 
229
  cachelist[8][8] = 1;
 
230
  cachelist[8][9] = 1;
 
231
  cachelist[9][5] = 1;
 
232
  cachelist[9][6] = 1;
 
233
  cachelist[9][7] = 1;
 
234
  cachelist[9][8] = 1;
 
235
  cachelist[9][9] = 1;
 
236
}
 
237
 
 
238
void cache_iabc_rhf(int **cachelist)
 
239
{
 
240
  /* <ia|bc> */
 
241
  cachelist[10][5] = 1;
 
242
  cachelist[10][6] = 1;
 
243
  cachelist[10][7] = 1;
 
244
  cachelist[10][8] = 1;
 
245
  cachelist[10][9] = 1;
 
246
  cachelist[11][5] = 1;
 
247
  cachelist[11][6] = 1;
 
248
  cachelist[11][7] = 1;
 
249
  cachelist[11][8] = 1;
 
250
  cachelist[11][9] = 1;
 
251
  /* <ab|ci> */
 
252
  cachelist[5][10] = 1;
 
253
  cachelist[5][11] = 1;
 
254
  cachelist[6][10] = 1;
 
255
  cachelist[6][11] = 1;
 
256
  cachelist[7][10] = 1;
 
257
  cachelist[7][11] = 1;
 
258
  cachelist[8][10] = 1;
 
259
  cachelist[8][11] = 1;
 
260
  cachelist[9][10] = 1;
 
261
  cachelist[9][11] = 1;
 
262
}
 
263
 
 
264
void cache_iabc_uhf(int **cachelist)
 
265
{
 
266
  /* <IA|BC> */
 
267
  cachelist[20][5] = 1;
 
268
  cachelist[20][6] = 1;
 
269
  cachelist[20][7] = 1;
 
270
  cachelist[20][8] = 1;
 
271
  cachelist[20][9] = 1;
 
272
  cachelist[21][5] = 1;
 
273
  cachelist[21][6] = 1;
 
274
  cachelist[21][7] = 1;
 
275
  cachelist[21][8] = 1;
 
276
  cachelist[21][9] = 1;
 
277
  /* <AB|CI> */
 
278
  cachelist[5][20] = 1;
 
279
  cachelist[5][21] = 1;
 
280
  cachelist[6][20] = 1;
 
281
  cachelist[6][21] = 1;
 
282
  cachelist[7][20] = 1;
 
283
  cachelist[7][21] = 1;
 
284
  cachelist[8][20] = 1;
 
285
  cachelist[8][21] = 1;
 
286
  cachelist[9][20] = 1;
 
287
  cachelist[9][21] = 1;
 
288
 
 
289
  /* <ia|bc> */
 
290
  cachelist[30][15] = 1;
 
291
  cachelist[30][16] = 1;
 
292
  cachelist[30][17] = 1;
 
293
  cachelist[30][18] = 1;
 
294
  cachelist[30][19] = 1;
 
295
  cachelist[31][15] = 1;
 
296
  cachelist[31][16] = 1;
 
297
  cachelist[31][17] = 1;
 
298
  cachelist[31][18] = 1;
 
299
  cachelist[31][19] = 1;
 
300
  /* <ab|ci> */
 
301
  cachelist[15][30] = 1;
 
302
  cachelist[15][31] = 1;
 
303
  cachelist[16][30] = 1;
 
304
  cachelist[16][31] = 1;
 
305
  cachelist[17][30] = 1;
 
306
  cachelist[17][31] = 1;
 
307
  cachelist[18][30] = 1;
 
308
  cachelist[18][31] = 1;
 
309
  cachelist[19][30] = 1;
 
310
  cachelist[19][31] = 1;
 
311
 
 
312
  /* <Ia|Bc> */
 
313
  cachelist[24][28] = 1;
 
314
  cachelist[24][29] = 1;
 
315
  cachelist[25][28] = 1;
 
316
  cachelist[25][29] = 1;
 
317
 
 
318
  /* <Ab|Ci> */
 
319
  cachelist[28][24] = 1;
 
320
  cachelist[28][25] = 1;
 
321
  cachelist[29][24] = 1;
 
322
  cachelist[29][25] = 1;
 
323
}
 
324
 
 
325
void cache_ijab_rhf(int **cachelist)
 
326
{
 
327
  /* <ij|ab> */
 
328
  cachelist[0][5] = 1;
 
329
  cachelist[0][6] = 1;
 
330
  cachelist[0][7] = 1;
 
331
  cachelist[0][8] = 1;
 
332
  cachelist[0][9] = 1;
 
333
  cachelist[1][5] = 1;
 
334
  cachelist[1][6] = 1;
 
335
  cachelist[1][7] = 1;
 
336
  cachelist[1][8] = 1;
 
337
  cachelist[1][9] = 1;
 
338
  cachelist[2][5] = 1;
 
339
  cachelist[2][6] = 1;
 
340
  cachelist[2][7] = 1;
 
341
  cachelist[2][8] = 1;
 
342
  cachelist[2][9] = 1;
 
343
  cachelist[3][5] = 1;
 
344
  cachelist[3][6] = 1;
 
345
  cachelist[3][7] = 1;
 
346
  cachelist[3][8] = 1;
 
347
  cachelist[3][9] = 1;
 
348
  cachelist[4][5] = 1;
 
349
  cachelist[4][6] = 1;
 
350
  cachelist[4][7] = 1;
 
351
  cachelist[4][8] = 1;
 
352
  cachelist[4][9] = 1;
 
353
  /* <ab|ij> */
 
354
  cachelist[5][0] = 1;
 
355
  cachelist[5][1] = 1;
 
356
  cachelist[5][2] = 1;
 
357
  cachelist[5][3] = 1;
 
358
  cachelist[5][4] = 1;
 
359
  cachelist[6][0] = 1;
 
360
  cachelist[6][1] = 1;
 
361
  cachelist[6][2] = 1;
 
362
  cachelist[6][3] = 1;
 
363
  cachelist[6][4] = 1;
 
364
  cachelist[7][0] = 1;
 
365
  cachelist[7][1] = 1;
 
366
  cachelist[7][2] = 1;
 
367
  cachelist[7][3] = 1;
 
368
  cachelist[7][4] = 1;
 
369
  cachelist[8][0] = 1;
 
370
  cachelist[8][1] = 1;
 
371
  cachelist[8][2] = 1;
 
372
  cachelist[8][3] = 1;
 
373
  cachelist[8][4] = 1;
 
374
  cachelist[9][0] = 1;
 
375
  cachelist[9][1] = 1;
 
376
  cachelist[9][2] = 1;
 
377
  cachelist[9][3] = 1;
 
378
  cachelist[9][4] = 1;
 
379
}
 
380
 
 
381
void cache_ijab_uhf(int **cachelist)
 
382
{
 
383
  /* <IJ|AB> */
 
384
  cachelist[0][5] = 1;
 
385
  cachelist[0][6] = 1;
 
386
  cachelist[0][7] = 1;
 
387
  cachelist[0][8] = 1;
 
388
  cachelist[0][9] = 1;
 
389
  cachelist[1][5] = 1;
 
390
  cachelist[1][6] = 1;
 
391
  cachelist[1][7] = 1;
 
392
  cachelist[1][8] = 1;
 
393
  cachelist[1][9] = 1;
 
394
  cachelist[2][5] = 1;
 
395
  cachelist[2][6] = 1;
 
396
  cachelist[2][7] = 1;
 
397
  cachelist[2][8] = 1;
 
398
  cachelist[2][9] = 1;
 
399
  cachelist[3][5] = 1;
 
400
  cachelist[3][6] = 1;
 
401
  cachelist[3][7] = 1;
 
402
  cachelist[3][8] = 1;
 
403
  cachelist[3][9] = 1;
 
404
  cachelist[4][5] = 1;
 
405
  cachelist[4][6] = 1;
 
406
  cachelist[4][7] = 1;
 
407
  cachelist[4][8] = 1;
 
408
  cachelist[4][9] = 1;
 
409
  /* <AB|IJ> */
 
410
  cachelist[5][0] = 1;
 
411
  cachelist[5][1] = 1;
 
412
  cachelist[5][2] = 1;
 
413
  cachelist[5][3] = 1;
 
414
  cachelist[5][4] = 1;
 
415
  cachelist[6][0] = 1;
 
416
  cachelist[6][1] = 1;
 
417
  cachelist[6][2] = 1;
 
418
  cachelist[6][3] = 1;
 
419
  cachelist[6][4] = 1;
 
420
  cachelist[7][0] = 1;
 
421
  cachelist[7][1] = 1;
 
422
  cachelist[7][2] = 1;
 
423
  cachelist[7][3] = 1;
 
424
  cachelist[7][4] = 1;
 
425
  cachelist[8][0] = 1;
 
426
  cachelist[8][1] = 1;
 
427
  cachelist[8][2] = 1;
 
428
  cachelist[8][3] = 1;
 
429
  cachelist[8][4] = 1;
 
430
  cachelist[9][0] = 1;
 
431
  cachelist[9][1] = 1;
 
432
  cachelist[9][2] = 1;
 
433
  cachelist[9][3] = 1;
 
434
  cachelist[9][4] = 1;
 
435
 
 
436
  /* <ij|ab> */
 
437
  cachelist[10][15] = 1;
 
438
  cachelist[10][16] = 1;
 
439
  cachelist[10][17] = 1;
 
440
  cachelist[10][18] = 1;
 
441
  cachelist[10][19] = 1;
 
442
  cachelist[11][15] = 1;
 
443
  cachelist[11][16] = 1;
 
444
  cachelist[11][17] = 1;
 
445
  cachelist[11][18] = 1;
 
446
  cachelist[11][19] = 1;
 
447
  cachelist[12][15] = 1;
 
448
  cachelist[12][16] = 1;
 
449
  cachelist[12][17] = 1;
 
450
  cachelist[12][18] = 1;
 
451
  cachelist[12][19] = 1;
 
452
  cachelist[13][15] = 1;
 
453
  cachelist[13][16] = 1;
 
454
  cachelist[13][17] = 1;
 
455
  cachelist[13][18] = 1;
 
456
  cachelist[13][19] = 1;
 
457
  cachelist[14][15] = 1;
 
458
  cachelist[14][16] = 1;
 
459
  cachelist[14][17] = 1;
 
460
  cachelist[14][18] = 1;
 
461
  cachelist[14][19] = 1;
 
462
  /* <ab|ij> */
 
463
  cachelist[15][10] = 1;
 
464
  cachelist[15][11] = 1;
 
465
  cachelist[15][12] = 1;
 
466
  cachelist[15][13] = 1;
 
467
  cachelist[15][14] = 1;
 
468
  cachelist[16][10] = 1;
 
469
  cachelist[16][11] = 1;
 
470
  cachelist[16][12] = 1;
 
471
  cachelist[16][13] = 1;
 
472
  cachelist[16][14] = 1;
 
473
  cachelist[17][10] = 1;
 
474
  cachelist[17][11] = 1;
 
475
  cachelist[17][12] = 1;
 
476
  cachelist[17][13] = 1;
 
477
  cachelist[17][14] = 1;
 
478
  cachelist[18][10] = 1;
 
479
  cachelist[18][11] = 1;
 
480
  cachelist[18][12] = 1;
 
481
  cachelist[18][13] = 1;
 
482
  cachelist[18][14] = 1;
 
483
  cachelist[19][10] = 1;
 
484
  cachelist[19][11] = 1;
 
485
  cachelist[19][12] = 1;
 
486
  cachelist[19][13] = 1;
 
487
  cachelist[19][14] = 1;
 
488
 
 
489
  /* <Ij|Ab> */
 
490
  cachelist[22][28] = 1;
 
491
  cachelist[23][28] = 1;
 
492
  cachelist[22][29] = 1;
 
493
  cachelist[23][29] = 1;
 
494
  /* <Ab|Ij> */
 
495
  cachelist[28][22] = 1;
 
496
  cachelist[28][23] = 1;
 
497
  cachelist[29][22] = 1;
 
498
  cachelist[29][23] = 1;
 
499
}
 
500
 
 
501
void cache_iajb_rhf(int **cachelist)
 
502
{
 
503
  /* <ia|jb> */
 
504
  cachelist[10][10] = 1;
 
505
  cachelist[10][11] = 1;
 
506
  cachelist[11][10] = 1;
 
507
  cachelist[11][11] = 1;
 
508
}
 
509
 
 
510
void cache_iajb_uhf(int **cachelist)
 
511
{
 
512
  /* <IA|JB> */
 
513
  cachelist[20][20] = 1;
 
514
  cachelist[20][21] = 1;
 
515
  cachelist[21][20] = 1;
 
516
  cachelist[21][21] = 1;
 
517
  /* <ia|jb> */
 
518
  cachelist[30][30] = 1;
 
519
  cachelist[30][31] = 1;
 
520
  cachelist[31][30] = 1;
 
521
  cachelist[31][31] = 1;
 
522
  /* <Ia|Jb> */
 
523
  cachelist[24][24] = 1;
 
524
  cachelist[24][25] = 1;
 
525
  cachelist[25][24] = 1;
 
526
  cachelist[25][25] = 1;
 
527
}
 
528
 
 
529
void cache_ijka_rhf(int **cachelist)
 
530
{
 
531
  /* <ij|ka> */
 
532
  cachelist[0][10] = 1;
 
533
  cachelist[0][11] = 1;
 
534
  cachelist[1][10] = 1;
 
535
  cachelist[1][11] = 1;
 
536
  cachelist[2][10] = 1;
 
537
  cachelist[2][11] = 1;
 
538
  cachelist[3][10] = 1;
 
539
  cachelist[3][11] = 1;
 
540
  cachelist[4][10] = 1;
 
541
  cachelist[4][11] = 1;
 
542
  /* <ia|jk> */
 
543
  cachelist[10][0] = 1;
 
544
  cachelist[10][1] = 1;
 
545
  cachelist[10][2] = 1;
 
546
  cachelist[10][3] = 1;
 
547
  cachelist[10][4] = 1;
 
548
  cachelist[11][0] = 1;
 
549
  cachelist[11][1] = 1;
 
550
  cachelist[11][2] = 1;
 
551
  cachelist[11][3] = 1;
 
552
  cachelist[11][4] = 1;
 
553
}
 
554
 
 
555
void cache_ijka_uhf(int **cachelist)
 
556
{
 
557
  /* <IJ|KA> */
 
558
  cachelist[0][20] = 1;
 
559
  cachelist[0][21] = 1;
 
560
  cachelist[1][20] = 1;
 
561
  cachelist[1][21] = 1;
 
562
  cachelist[2][20] = 1;
 
563
  cachelist[2][21] = 1;
 
564
  cachelist[3][20] = 1;
 
565
  cachelist[3][21] = 1;
 
566
  cachelist[4][20] = 1;
 
567
  cachelist[4][21] = 1;
 
568
  /* <IA|JK> */
 
569
  cachelist[20][0] = 1;
 
570
  cachelist[20][1] = 1;
 
571
  cachelist[20][2] = 1;
 
572
  cachelist[20][3] = 1;
 
573
  cachelist[20][4] = 1;
 
574
  cachelist[21][0] = 1;
 
575
  cachelist[21][1] = 1;
 
576
  cachelist[21][2] = 1;
 
577
  cachelist[21][3] = 1;
 
578
  cachelist[21][4] = 1;
 
579
 
 
580
  /* <ij|ka> */
 
581
  cachelist[10][30] = 1;
 
582
  cachelist[10][31] = 1;
 
583
  cachelist[11][30] = 1;
 
584
  cachelist[11][31] = 1;
 
585
  cachelist[12][30] = 1;
 
586
  cachelist[12][31] = 1;
 
587
  cachelist[13][30] = 1;
 
588
  cachelist[13][31] = 1;
 
589
  cachelist[14][30] = 1;
 
590
  cachelist[14][31] = 1;
 
591
  /* <ia|jk> */
 
592
  cachelist[30][10] = 1;
 
593
  cachelist[30][11] = 1;
 
594
  cachelist[30][12] = 1;
 
595
  cachelist[30][13] = 1;
 
596
  cachelist[30][14] = 1;
 
597
  cachelist[31][10] = 1;
 
598
  cachelist[31][11] = 1;
 
599
  cachelist[31][12] = 1;
 
600
  cachelist[31][13] = 1;
 
601
  cachelist[31][14] = 1;
 
602
 
 
603
  /* <Ij|Ka> */
 
604
  cachelist[22][24] = 1;
 
605
  cachelist[22][25] = 1;
 
606
  cachelist[23][24] = 1;
 
607
  cachelist[23][25] = 1;
 
608
  /* <Ka|Ij> */
 
609
  cachelist[24][22] = 1;
 
610
  cachelist[25][22] = 1;
 
611
  cachelist[24][23] = 1;
 
612
  cachelist[25][23] = 1;
 
613
}
 
614
 
 
615
void cache_ijkl_rhf(int **cachelist)
 
616
{
 
617
  /* <ij|kl> */
 
618
  cachelist[0][0] = 1;
 
619
  cachelist[0][1] = 1;
 
620
  cachelist[0][2] = 1;
 
621
  cachelist[0][3] = 1;
 
622
  cachelist[0][4] = 1;
 
623
  cachelist[1][0] = 1;
 
624
  cachelist[1][1] = 1;
 
625
  cachelist[1][2] = 1;
 
626
  cachelist[1][3] = 1;
 
627
  cachelist[1][4] = 1;
 
628
  cachelist[2][0] = 1;
 
629
  cachelist[2][1] = 1;
 
630
  cachelist[2][2] = 1;
 
631
  cachelist[2][3] = 1;
 
632
  cachelist[2][4] = 1;
 
633
  cachelist[3][0] = 1;
 
634
  cachelist[3][1] = 1;
 
635
  cachelist[3][2] = 1;
 
636
  cachelist[3][3] = 1;
 
637
  cachelist[3][4] = 1;
 
638
  cachelist[4][0] = 1;
 
639
  cachelist[4][1] = 1;
 
640
  cachelist[4][2] = 1;
 
641
  cachelist[4][3] = 1;
 
642
  cachelist[4][4] = 1;
 
643
}
 
644
 
 
645
void cache_ijkl_uhf(int **cachelist)
 
646
{
 
647
  /* <IJ|KL> */
 
648
  cachelist[0][0] = 1;
 
649
  cachelist[0][1] = 1;
 
650
  cachelist[0][2] = 1;
 
651
  cachelist[0][3] = 1;
 
652
  cachelist[0][4] = 1;
 
653
  cachelist[1][0] = 1;
 
654
  cachelist[1][1] = 1;
 
655
  cachelist[1][2] = 1;
 
656
  cachelist[1][3] = 1;
 
657
  cachelist[1][4] = 1;
 
658
  cachelist[2][0] = 1;
 
659
  cachelist[2][1] = 1;
 
660
  cachelist[2][2] = 1;
 
661
  cachelist[2][3] = 1;
 
662
  cachelist[2][4] = 1;
 
663
  cachelist[3][0] = 1;
 
664
  cachelist[3][1] = 1;
 
665
  cachelist[3][2] = 1;
 
666
  cachelist[3][3] = 1;
 
667
  cachelist[3][4] = 1;
 
668
  cachelist[4][0] = 1;
 
669
  cachelist[4][1] = 1;
 
670
  cachelist[4][2] = 1;
 
671
  cachelist[4][3] = 1;
 
672
  cachelist[4][4] = 1;
 
673
  /* <ij|kl> */
 
674
  cachelist[10][10] = 1;
 
675
  cachelist[10][11] = 1;
 
676
  cachelist[10][12] = 1;
 
677
  cachelist[10][13] = 1;
 
678
  cachelist[10][14] = 1;
 
679
  cachelist[11][10] = 1;
 
680
  cachelist[11][11] = 1;
 
681
  cachelist[11][12] = 1;
 
682
  cachelist[11][13] = 1;
 
683
  cachelist[11][14] = 1;
 
684
  cachelist[12][10] = 1;
 
685
  cachelist[12][11] = 1;
 
686
  cachelist[12][12] = 1;
 
687
  cachelist[12][13] = 1;
 
688
  cachelist[12][14] = 1;
 
689
  cachelist[13][10] = 1;
 
690
  cachelist[13][11] = 1;
 
691
  cachelist[13][12] = 1;
 
692
  cachelist[13][13] = 1;
 
693
  cachelist[13][14] = 1;
 
694
  cachelist[14][10] = 1;
 
695
  cachelist[14][11] = 1;
 
696
  cachelist[14][12] = 1;
 
697
  cachelist[14][13] = 1;
 
698
  cachelist[14][14] = 1;
 
699
  /* <Ij|Kl> */
 
700
  cachelist[22][22] = 1;
 
701
  cachelist[22][23] = 1;
 
702
  cachelist[23][22] = 1;
 
703
  cachelist[23][23] = 1;
 
704
}
 
705
 
 
706
void cachedone_uhf(int **cachelist)
 
707
{
 
708
  free_int_matrix(cachelist, 32);
 
709
}
 
710
 
 
711
void cachedone_rhf(int **cachelist)
 
712
{
 
713
  free_int_matrix(cachelist, 12);
 
714
}