Added a script to create a mock galaxy spectrum and fit using Imfit.
Changed fitting code to Imfit.

Cosmetic changes to plots.
Added more options for PSF.

Changed Imfit calls to be compatible to latest dev version.
Better flags and mask handling for single-lambda spectra slices.

Added support to fat spectra slices.
Fixed bug when not setting PSF beta parameter.
Fixed bug where using single-lambda slices causes a crash because of bad values for noise.
Refactored fitting to support different initial models for every slice, allowing multiple passes.
Added optional 2-pass fitting with smoothing of center, PA and b/a.
Save plots to a directory.
Fixed bug with unflagged/unmasked pixels with noise==0.
Changed plot names.
Set grid name for easier debugging of starlight runs.
Skip fitting if there are less than 1000 pixels for fitting.
Added Sérsic index to bulge fitting.

Changed initial guess (must find a better way).

Better plotting limits and mask values.
Plot radial profiles.
Faster v_d convolution from pystarlight, no need to create cache files any more.
Starting to separate the "documentation code" from the library.
Small changes to the simple mock galaxy.
Fixed bug creating planes image.
Initial fit using Differential Evolution.

Added option to retry a fit using Nelder-Mead simplex if least squares fails.