fit_double_gaussian fails in input data

[peendebak]

The fit_double_gaussian can fail with an error message reporting nan values. A minimal example:

import numpy as np
import qtt

xx=np.arange(10.)
yy=0*np.arange(10)
yy[3]=80.
yy[6]=.1

qtt.algorithms.fitting.fit_double_gaussian(xx, yy)

Notes

A solution might be to redefine the gaussian used in the model to fit a exp(-alpha x+b) instead of a exp(x/sigma + b)