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

« back to all changes in this revision

Viewing changes to src/bin/ccdensity/dump_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
* 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:
6
6
#define EXTERN
7
7
#include "globals.h"
8
8
 
 
9
EXTERN dpd_gbl dpd_main;
 
10
 
 
11
 
9
12
/* DUMP_ROHF(): Mulliken-order the ROHF-CCSD two-electron density and
10
13
** dump it to a file for subsequent backtransformation.  Basically all
11
14
** we have to do is swap indices two and three, e.g.
28
31
** I really need to give an example of this problem using specific
29
32
** elements of GIJKA so that the code below will be clearer.*/
30
33
 
31
 
void dump_ROHF(struct iwlbuf *OutBuf)
 
34
void dump_ROHF(struct iwlbuf *OutBuf, struct RHO_Params rho_params)
32
35
{
33
36
  int nirreps, nmo, nfzv;
34
37
  int *qt_occ, *qt_vir;
42
45
  nfzv = moinfo.nfzv;
43
46
 
44
47
  psio_open(PSIF_MO_OPDM, PSIO_OPEN_OLD);
45
 
  psio_write_entry(PSIF_MO_OPDM, "MO-basis OPDM", (char *) moinfo.opdm[0],
46
 
                   sizeof(double)*nmo*nmo);
 
48
 /*  psio_write_entry(PSIF_MO_OPDM, "MO-basis OPDM", (char *) moinfo.opdm[0], */
 
49
  psio_write_entry(PSIF_MO_OPDM, rho_params.opdm_lbl, (char *) moinfo.opdm[0],
 
50
                   sizeof(double)*(nmo-nfzv)*(nmo-nfzv));
47
51
  psio_close(PSIF_MO_OPDM, 1);
48
52
 
 
53
if (!params.onepdm) {
49
54
  psio_open(PSIF_MO_LAG, PSIO_OPEN_OLD);
50
55
  psio_write_entry(PSIF_MO_LAG, "MO-basis Lagrangian", (char *) moinfo.I[0],
51
56
                   sizeof(double)*nmo*nmo);
150
155
 
151
156
  dpd_buf4_dump(&G, OutBuf, qt_vir, qt_vir, qt_occ, qt_vir, 0, 0);
152
157
  dpd_buf4_close(&G);
153
 
 
154
158
  dpd_buf4_init(&G, CC_GAMMA, 0, 5, 5, 5, 5, 0, "GAbCd");
155
 
  dpd_buf4_sort_ooc(&G, CC_TMP0, prqs, 5, 5, "G(AC,BD)");
 
159
  dpd_buf4_sort(&G, CC_TMP0, prqs, 5, 5, "G(AC,BD)");
156
160
  dpd_buf4_close(&G);
157
161
  dpd_buf4_init(&G, CC_TMP0, 0, 5, 5, 5, 5, 0, "G(AC,BD)");
158
162
  dpd_buf4_dump(&G, OutBuf, qt_vir, qt_vir, qt_vir, qt_vir, 1, 0);
159
163
  dpd_buf4_close(&G);
 
164
 
 
165
  }
160
166
}