Try different Non negative matrix factorization decompositions
- different initialization methods in Julia?
- try 10 random initialization points and take the best
- try matching on the results
- try to make the decomposition sparse
- make sure that the input data is non-negative