- for (j = 0;;)
- {
- double fij = xt->mat[j + i * n_cols];
-
- if (proc->statistics & CRS_ST_BTAU)
- {
- const double temp = (2. * sqrt (Dr * Dc) * (Cij - Dij)
- + v[3] * (xt->row_tot[i] * Dc
- + xt->col_tot[j] * Dr));
- btau_cum += fij * temp * temp;
- }
-
- {
- const double temp = Cij - Dij;
- ctau_cum += fij * temp * temp;
- }
-
- if (proc->statistics & CRS_ST_GAMMA)
- {
- const double temp = Q * Cij - P * Dij;
- gamma_cum += fij * temp * temp;
- }
-
- if (proc->statistics & CRS_ST_D)
- {
- d_yx_cum += fij * pow2 (Dr * (Cij - Dij)
- - (P - Q) * (xt->total - xt->row_tot[i]));
- d_xy_cum += fij * pow2 (Dc * (Dij - Cij)
- - (Q - P) * (xt->total - xt->col_tot[j]));
- }
-
- if (++j == n_cols)
- break;
- assert (j < n_cols);
-
- Cij -= xt->col_tot[j] - cum[j + i * n_cols];
- Dij += xt->col_tot[j - 1] - cum[j - 1 + i * n_cols];
-
- if (i > 0)
- {
- Cij += cum[j - 1 + (i - 1) * n_cols];
- Dij -= cum[j + (i - 1) * n_cols];
- }
- }
- }
- }
+ if (proc->statistics & CRS_ST_GAMMA)
+ gamma_cum += fij * pow2 (Q * Cij - P * Dij);
+
+ if (proc->statistics & CRS_ST_D)
+ {
+ d_yx_cum += fij * pow2 (Dr * (Cij - Dij)
+ - (P - Q) * (xt->total - xt->row_tot[i]));
+ d_xy_cum += fij * pow2 (Dc * (Dij - Cij)
+ - (Q - P) * (xt->total - xt->col_tot[j]));
+ }
+
+ if (++j >= n_cols)
+ break;
+
+ Cij -= xt->col_tot[j] - cum[j + i * n_cols];
+ Dij += xt->col_tot[j - 1] - cum[j - 1 + i * n_cols];
+
+ if (i > 0)
+ {
+ Cij += cum[j - 1 + (i - 1) * n_cols];
+ Dij -= cum[j + (i - 1) * n_cols];
+ }
+ }
+ }