start of iterative implementation
Optimized working version of AENMF
Optimized NMF solver using ALS with dual coordinate descent and initialization using convex hulls