{
ordered_cols[i] = i;
}
- for (i = last_col + 1; i < A->size1 - 1; i++)
+ for (i = last_col + 1; i < A->size1; i++)
{
ordered_cols[i - 1] = i;
}
}
}
gsl_matrix_free (B);
- return GSL_SUCCESS;
free (ordered_cols);
+ return GSL_SUCCESS;
}
return GSL_ENOTSQR;
}