Added storing of models visited during local optimization

Added the possibility to find best mlik per model in subsampling scenarios
Fixed a bug in the greedy algorithm