Author: Hans Terlouw <gipsy@astro.rug.nl>
Function to search for gaussian components in a profile.
Parameters: 


Returns:  a list with up to ncomp tuples of which each tuple contains the amplitude, the centre and the dispersion of the gaussian, in that order. If the argument flat is True, a ‘flat’ list with ncomp * 3 numbers is returned which may directly be used as initial estimates for kmpfit.Fitter. 
In this function the second derivative of the profile in the signal region is calculated by fitting a second degree polynomal. The smoothing parameter q determines the number of points used for this (\(2q+1\)). The gaussians are then estimated as described by [Schwarz1968].
See also
Tutorial chapter “Automatic initial estimates for profiles with multi component Gaussians.”
[Schwarz1968]  Schwarz, U.J., 1968. Analysis of an Observed Function into Components, using its Second Derivative, Communication from the Netherlands Foundation for Radio Astronomy and the Kapteyn Astronomical Laboratory at Groningen, 19 405. (local copy) 