pub fn col_scale(a: &CsrMatrix, scales: &[f64]) -> CsrMatrix
Column-scale a CSR matrix: A D_c where D_c = diag(scales).