Multiadaptive run seems to work quite well for the bimodal case
Fixed improper usage of numeric_limits in EMPeriodicGaussian and samplePoints in EMGaussian