Added parabolic_stop into .fit() function.
Allows to use huge n_components and later do not think about their number: BIC search stage will stop automatically as soon as it will found parabolic-like minimum in BIC.
Default value of parabolic_stop is -1 (i.e. do not use this mode)
Recommended value when using huge n_components number : 0.7