Function nnls 

pub fn nnls(a: &[Vec<f64>], b: &[f64]) -> Vec<f64>

a is column-major: a[j] is the j-th column (length m, matching b’s length). Returns the n = a.len()-length non-negative solution.