be identical. If missing values are involved, then element (i,j)
is the moment of the i th variable, when paired with the j th variable.
*/
-const gsl_matrix *
+gsl_matrix *
covariance_moments (const struct covariance *cov, int m)
{
return cov->moments[m];
bool centered)
{
size_t i;
- struct covariance *cov = xzalloc (sizeof *cov);
+ struct covariance *cov = XZALLOC (struct covariance);
cov->centered = centered;
cov->passes = 1;