Use projected gradient solver in transform to support sparse matrices
Use same parameters when solving the transform
Use scipy.nnls.optimize for dense data