3
\brief Enter brief description of file here
10
/* dpd_buf4_dot_self(): Evaluates the sum of the squares of the elements of a
14
** dpdbuf4 *BufX: A pointer to the dpdbuf4.
17
double dpd_buf4_dot_self(dpdbuf4 *BufX)
19
int h, nirreps, my_irrep;
23
nirreps = BufX->params->nirreps;
24
my_irrep = BufX->file.my_irrep;
26
for(h=0; h < nirreps; h++) {
28
dpd_buf4_mat_irrep_init(BufX, h);
29
dpd_buf4_mat_irrep_rd(BufX, h);
31
for(row=0; row < BufX->params->rowtot[h]; row++)
32
for(col=0; col < BufX->params->coltot[h^my_irrep]; col++)
33
alpha += BufX->matrix[h][row][col] * BufX->matrix[h][row][col];
35
dpd_buf4_mat_irrep_close(BufX, h);