ENH: added optional R^2 output map
ENH: adding the capability to re-do the fit after discarding outliers

Benefits of doing this are not clear at the moment - need to look more.
ENH: use different initialization for the optimizer

following the suggestion from Stephan Maier, BWH