X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;f=src%2Fmath%2Flinreg.c;h=8028697c6b7f3410b12927a0fe7ae1d82c5a6b4b;hb=07db00919d6f067fd5dd6c6c1c2c2fba4f42cf21;hp=f5ae33e5823f4f79cfaccd73637316baa2f96d0f;hpb=4aa2ac7ebf21370681326bb03182dac973f23a29;p=pspp-builds.git diff --git a/src/math/linreg.c b/src/math/linreg.c index f5ae33e5..8028697c 100644 --- a/src/math/linreg.c +++ b/src/math/linreg.c @@ -669,6 +669,7 @@ rearrange_covariance_matrix (const struct covariance_matrix *cm, pspp_linreg_cac size_t i; size_t j; size_t k; + size_t n_coeffs = 0; assert (cm != NULL); cov = covariance_to_design (cm); @@ -676,11 +677,11 @@ rearrange_covariance_matrix (const struct covariance_matrix *cm, pspp_linreg_cac assert (c != NULL); assert (cov->m->size1 > 0); assert (cov->m->size2 == cov->m->size1); - permutation = xnmalloc (1 + c->n_indeps, sizeof (*permutation)); model_vars = xnmalloc (1 + c->n_indeps, sizeof (*model_vars)); /* Put the model variables in the right order in MODEL_VARS. + Count the number of coefficients. */ for (i = 0; i < c->n_indeps; i++) { @@ -688,6 +689,8 @@ rearrange_covariance_matrix (const struct covariance_matrix *cm, pspp_linreg_cac } model_vars[i] = c->depvar; result = covariance_matrix_create (1 + c->n_indeps, model_vars); + permutation = xnmalloc (design_matrix_get_n_cols (result), sizeof (*permutation)); + for (j = 0; j < cov->m->size2; j++) { k = 0; @@ -742,4 +745,3 @@ double pspp_linreg_mse (const pspp_linreg_cache *c) assert (c != NULL); return (c->sse / c->dfe); } -