Add an optimizer option to SciPy to return the whole simplex

[andsor]

As calculating the errors involves the whole final simplex:

pyfssa/fssa/fssa.py

Line 681 in 6d584a8

errors, varco = _neldermead_errors(sim, fsim, func)

[andsor]

PR scipy/scipy#5205