inline SGVector< float64_t > get_mean(int32_t c) const
{
    return SGVector< float64_t >(m_means.get_column_vector(c), m_dim);
}

T* get_column_vector(index_t col) const
{
    return &matrix[col*num_rows];
}