Added smoothing of probability distributions to improve robustness of results

Added Cauchy smoothing to 2D factors and Gussian smoothing to 1D ones to
reduce significant overfitting in the high-dimensional scenario that we
approached by setting the discetization resolution to 50 for each
variable. This greatly improves the robustness of results especially in
respect with the D-statistics null distribution calculation.