← Back to branch summary

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

« back to all changes in this revision

Viewing changes to src/bin/ccdensity/sortone_ROHF.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/ccdensity/sortone_ROHF.c
21
21
** TDC, 1/03
22
22
*/
23
23
 
24
 
void sortone_ROHF(void)
 
24
void sortone_ROHF(struct RHO_Params rho_params)
25
25
{
26
26
  int h, nirreps, nmo, nfzv, nfzc, nclsd, nopen;
27
27
  int row, col, i, j, I, J, a, b, A, B, p, q;
44
44
  openpi = moinfo.openpi;
45
45
  qt_occ = moinfo.qt_occ; qt_vir = moinfo.qt_vir;
46
46
 
47
 
  O = block_matrix(nmo-nfzc,nmo-nfzc);
 
47
  /* O = block_matrix(nmo-nfzc,nmo-nfzc); */
 
48
  O = block_matrix(nmo-nfzv, nmo-nfzv);
48
49
 
49
50
  /* Sort A components first */
50
 
  dpd_file2_init(&D, CC_OEI, 0, 0, 0, "DIJ");
 
51
  dpd_file2_init(&D, CC_OEI, 0, 0, 0, rho_params.DIJ_lbl);
51
52
  dpd_file2_mat_init(&D);
52
53
  dpd_file2_mat_rd(&D);
53
54
  for(h=0; h < nirreps; h++) {
62
63
  dpd_file2_mat_close(&D);
63
64
  dpd_file2_close(&D);
64
65
 
65
 
  dpd_file2_init(&D, CC_OEI, 0, 1, 1, "DAB");
 
66
  dpd_file2_init(&D, CC_OEI, 0, 1, 1, rho_params.DAB_lbl);
66
67
  dpd_file2_mat_init(&D);
67
68
  dpd_file2_mat_rd(&D);
68
69
  for(h=0; h < nirreps; h++) {
79
80
  dpd_file2_close(&D);
80
81
 
81
82
  /* Note that this component of the density is stored occ-vir */
82
 
  dpd_file2_init(&D, CC_OEI, 0, 0, 1, "DAI");
 
83
  dpd_file2_init(&D, CC_OEI, 0, 0, 1, rho_params.DAI_lbl);
83
84
  dpd_file2_mat_init(&D);
84
85
  dpd_file2_mat_rd(&D);
85
86
  for(h=0; h < nirreps; h++) {
95
96
  dpd_file2_mat_close(&D);
96
97
  dpd_file2_close(&D);
97
98
 
98
 
  dpd_file2_init(&D, CC_OEI, 0, 0, 1, "DIA");
 
99
  dpd_file2_init(&D, CC_OEI, 0, 0, 1, rho_params.DIA_lbl);
99
100
  dpd_file2_mat_init(&D);
100
101
  dpd_file2_mat_rd(&D);
101
102
  for(h=0; h < nirreps; h++) {
112
113
  dpd_file2_close(&D);
113
114
 
114
115
  /* Sort B components */
115
 
  dpd_file2_init(&D, CC_OEI, 0, 0, 0, "Dij");
 
116
  dpd_file2_init(&D, CC_OEI, 0, 0, 0, rho_params.Dij_lbl);
116
117
  dpd_file2_mat_init(&D); 
117
118
  dpd_file2_mat_rd(&D);
118
119
  for(h=0; h < nirreps; h++) {
127
128
  dpd_file2_mat_close(&D);
128
129
  dpd_file2_close(&D);
129
130
 
130
 
  dpd_file2_init(&D, CC_OEI, 0, 1, 1, "Dab");
 
131
  dpd_file2_init(&D, CC_OEI, 0, 1, 1, rho_params.Dab_lbl);
131
132
  dpd_file2_mat_init(&D);
132
133
  dpd_file2_mat_rd(&D);
133
134
  for(h=0; h < nirreps; h++) {
144
145
  dpd_file2_close(&D);
145
146
 
146
147
  /* Note that this component of the density is stored occ-vir */
147
 
  dpd_file2_init(&D, CC_OEI, 0, 0, 1, "Dai");
 
148
  dpd_file2_init(&D, CC_OEI, 0, 0, 1, rho_params.Dai_lbl);
148
149
  dpd_file2_mat_init(&D);
149
150
  dpd_file2_mat_rd(&D);
150
151
  for(h=0; h < nirreps; h++) {
160
161
  dpd_file2_mat_close(&D);
161
162
  dpd_file2_close(&D);
162
163
 
163
 
  dpd_file2_init(&D, CC_OEI, 0, 0, 1, "Dia");
 
164
  dpd_file2_init(&D, CC_OEI, 0, 0, 1, rho_params.Dia_lbl);
164
165
  dpd_file2_mat_init(&D);
165
166
  dpd_file2_mat_rd(&D);
166
167
  for(h=0; h < nirreps; h++) {
185
186
        }
186
187
    }
187
188
 
 
189
  /*
 
190
  for (i=0;i<(nmo-nfzv);++i) {
 
191
    for (j=0;j<(nmo-nfzv); ++j) {
 
192
      fprintf(outfile,"%15.10lf\n", O[i][j]);
 
193
    }
 
194
    fprintf(outfile,"\n");
 
195
  }
 
196
  fprintf(outfile,"End of sortone\n");
 
197
        */
 
198
  /*
 
199
  for (i=0;i<2;++i) {
 
200
    for (j=0;j<5;++j)
 
201
      fprintf(outfile,"%15.10lf",O[i][j]);
 
202
    fprintf(outfile,"\n");
 
203
  }
 
204
  print_mat(O,nmo-nfzv,nmo-nfzv,outfile);
 
205
  */
 
206
 
188
207
  moinfo.opdm = O;
189
208
 
190
209
}