cov->n_cm = (n_vars * (n_vars - 1) ) / 2;
- if (cov->n_cm > 0)
- cov->cm = xcalloc (sizeof *cov->cm, cov->n_cm);
+
+ cov->cm = xcalloc (cov->n_cm, sizeof *cov->cm);
cov->categoricals = NULL;
return cov;
static gsl_matrix *
covariance_calculate_double_pass_unnormalized (struct covariance *cov)
{
- size_t i, j;
- for (i = 0 ; i < cov->dim; ++i)
- {
- for (j = 0 ; j < cov->dim; ++j)
- {
- int idx;
- double *x = gsl_matrix_ptr (cov->moments[MOMENT_VARIANCE], i, j);
-
- idx = cm_idx (cov, i, j);
- if ( idx >= 0)
- {
- x = &cov->cm [idx];
- }
- }
- }
-
return cm_to_gsl (cov);
}