Enhance initial guess fit speed

The current implementation is extremely slow, despite executing in parallel. This is because it uses DifferentialEvolution to find a global minimum, but its very computationally expensive.

Options:

• GPU fit (https://github.com/gpufit/Gpufit), possible overkill
• Normal gradient descent + better initial guesses
• Curve peeling