← Back to branch summary

~ubuntu-branches/ubuntu/precise/psicode/precise

« back to all changes in this revision

Viewing changes to src/bin/ccenergy/diis_RHF.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/ccenergy/diis_RHF.c
27
27
{
28
28
  int nvector=8;  /* Number of error vectors to keep */
29
29
  int h, nirreps;
30
 
  int row, col, word, p, q;
 
30
  int row, col, word, t1_word, p, q;
31
31
  int diis_cycle;
32
32
  int vector_length=0;
33
33
  int errcod, *ipiv;
72
72
  dpd_file2_mat_close(&T1b);
73
73
  dpd_file2_close(&T1b);
74
74
 
 
75
  t1_word = word;
75
76
  dpd_buf4_init(&T2a, CC_TAMPS, 0, 0, 5, 0, 5, 0, "New tIjAb");
 
77
  for(h=0; h < nirreps; h++) {
 
78
    dpd_buf4_mat_irrep_init(&T2a, h);
 
79
    dpd_buf4_mat_irrep_rd(&T2a, h);
 
80
    for(row=0; row < T2a.params->rowtot[h]; row++)
 
81
      for(col=0; col < T2a.params->coltot[h]; col++)
 
82
        error[0][word++] = T2a.matrix[h][row][col];
 
83
    dpd_buf4_mat_irrep_close(&T2a, h);
 
84
  }
 
85
  dpd_buf4_close(&T2a);
 
86
 
 
87
  word = t1_word;
76
88
  dpd_buf4_init(&T2b, CC_TAMPS, 0, 0, 5, 0, 5, 0, "tIjAb");
77
89
  for(h=0; h < nirreps; h++) {
78
 
    dpd_buf4_mat_irrep_init(&T2a, h);
79
 
    dpd_buf4_mat_irrep_rd(&T2a, h);
80
90
    dpd_buf4_mat_irrep_init(&T2b, h);
81
91
    dpd_buf4_mat_irrep_rd(&T2b, h);
82
 
    for(row=0; row < T2a.params->rowtot[h]; row++)
83
 
      for(col=0; col < T2a.params->coltot[h]; col++)
84
 
        error[0][word++] = T2a.matrix[h][row][col] - T2b.matrix[h][row][col];
85
 
    dpd_buf4_mat_irrep_close(&T2a, h);
 
92
    for(row=0; row < T2b.params->rowtot[h]; row++)
 
93
      for(col=0; col < T2b.params->coltot[h]; col++)
 
94
        error[0][word++] -= T2b.matrix[h][row][col];
86
95
    dpd_buf4_mat_irrep_close(&T2b, h);
87
96
  }
88
 
  dpd_buf4_close(&T2a);
89
97
  dpd_buf4_close(&T2b);
90
98
 
91
99
  start = psio_get_address(PSIO_ZERO, diis_cycle*vector_length*sizeof(double));