← Back to branch summary

~ubuntu-branches/ubuntu/vivid/psicode/vivid

« back to all changes in this revision

Viewing changes to src/bin/cclambda/L2.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/cclambda/L2.c
3
3
#define EXTERN
4
4
#include "globals.h"
5
5
 
6
 
void DL2(int L_irr, int root_L_irr);
 
6
/* L2_build(): 
 
7
 
 
8
*/
 
9
 
 
10
void DL2(struct L_Params L_params);
7
11
void FaeL2(int L_irr);
8
12
void FmiL2(int L_irr);
9
13
void WijmnL2(int L_irr);
17
21
void dijabL2(int L_irr);
18
22
 
19
23
void BL2_AO(int L_irr);
20
 
 
21
 
void L2_build(int L_irr, int root_L_irr) {
22
 
 
23
 
  DL2(L_irr, root_L_irr);
24
 
#ifdef EOM_DEBUG
25
 
check_sum("DL2", L_irr);
26
 
#endif
27
 
 
28
 
  FaeL2(L_irr);
29
 
#ifdef EOM_DEBUG
30
 
check_sum("FaeL2", L_irr);
31
 
#endif
32
 
 
33
 
  FmiL2(L_irr);
34
 
#ifdef EOM_DEBUG
35
 
check_sum("FmiL2", L_irr);
 
24
void status(char *, FILE *);
 
25
 
 
26
void L2_build(struct L_Params L_params) {
 
27
  dpdbuf4 L2;
 
28
  int L_irr;
 
29
  L_irr = L_params.irrep;
 
30
 
 
31
  DL2(L_params);
 
32
  if(params.print & 2) status("<ij||ab> -> L2", outfile);
 
33
 
 
34
#ifdef EOM_DEBUG
 
35
  check_sum("DL2", L_irr);
36
36
#endif
37
37
 
38
38
  WijmnL2(L_irr);
39
39
#ifdef EOM_DEBUG
40
 
check_sum("WijmnL2", L_irr);
 
40
  check_sum("WijmnL2", L_irr);
41
41
#endif
 
42
  if(params.print & 2) status("Wmnij -> L2", outfile);
42
43
 
43
44
  WefabL2(L_irr);
44
45
#ifdef EOM_DEBUG
45
 
check_sum("WefabL2", L_irr);
 
46
  check_sum("WefabL2", L_irr);
46
47
#endif
 
48
  if(params.print & 2) status("Wabef -> L2", outfile);
47
49
 
48
50
  WejabL2(L_irr);
49
51
#ifdef EOM_DEBUG
50
 
check_sum("WejabL2", L_irr);
 
52
  check_sum("WejabL2", L_irr);
51
53
#endif
 
54
  if(params.print & 2) status("Wamef -> L2", outfile);
52
55
 
53
56
  WijmbL2(L_irr);
54
57
#ifdef EOM_DEBUG
55
 
check_sum("WijmbL2", L_irr);
56
 
#endif
57
 
 
58
 
  WmbejL2(L_irr); 
59
 
#ifdef EOM_DEBUG
60
 
check_sum("WmbejL2", L_irr);
61
 
#endif
62
 
 
63
 
  L1FL2(L_irr);
64
 
#ifdef EOM_DEBUG
65
 
check_sum("L1FL2", L_irr);
66
 
#endif
 
58
  check_sum("WijmbL2", L_irr);
 
59
#endif
 
60
  if(params.print & 2) status("Wmnie -> L2", outfile);
67
61
 
68
62
  GaeL2(L_irr);
69
63
#ifdef EOM_DEBUG
70
 
check_sum("GaeL2", L_irr);
 
64
  check_sum("GaeL2", L_irr);
71
65
#endif
72
66
 
73
67
  GmiL2(L_irr);
74
68
#ifdef EOM_DEBUG
75
 
check_sum("GmiL2", L_irr);
76
 
#endif
 
69
  check_sum("GmiL2", L_irr);
 
70
#endif
 
71
  if(params.print & 2) status("G -> L2", outfile);
 
72
 
 
73
  /* For RHF-CCSD response calculations, save all the above
 
74
     contributions to the L2 residual for use in the ccresponse code
 
75
     (specifically, HX1Y1 and LHX1Y1). */
 
76
  if(params.ref == 0 && params.dertype == 3) {
 
77
    dpd_buf4_init(&L2, CC_LAMBDA, L_irr, 0, 5, 0, 5, 0, "New LIjAb");
 
78
    dpd_buf4_copy(&L2, CC_LAMPS, "LHX1Y1 Residual I");
 
79
    dpd_buf4_close(&L2);
 
80
  }
 
81
 
 
82
  FaeL2(L_irr);
 
83
#ifdef EOM_DEBUG
 
84
  check_sum("FaeL2", L_irr);
 
85
#endif
 
86
 
 
87
  FmiL2(L_irr);
 
88
#ifdef EOM_DEBUG
 
89
  check_sum("FmiL2", L_irr);
 
90
#endif
 
91
  if(params.print & 2) status("F -> L2", outfile);
 
92
 
 
93
  WmbejL2(L_irr); 
 
94
#ifdef EOM_DEBUG
 
95
  check_sum("WmbejL2", L_irr);
 
96
#endif
 
97
  if(params.print & 2) status("Wmbej -> L2", outfile);
 
98
 
 
99
  if(!params.sekino) L1FL2(L_irr); /* should be dropped for Sekino-Bartlett modelIII approach */
 
100
#ifdef EOM_DEBUG
 
101
  check_sum("L1FL2", L_irr);
 
102
#endif
 
103
  if(params.print & 2) status("L1*F -> L2", outfile);
77
104
 
78
105
  dijabL2(L_irr);
79
106
#ifdef EOM_DEBUG
80
107
  check_sum("after D2s", L_irr);
81
108
#endif
 
109
  if(params.print & 2) status("L2 amplitudes", outfile);
82
110
}
83
111