even simpler handling of delta distribution

Param/Distrib/RangedDistributions.cpp

@@ -13,6 +13,7 @@
 //  ************************************************************************************************
 
 #include "Param/Distrib/RangedDistributions.h"
+#include "Base/Util/Assert.h"
 #include "Param/Distrib/Distributions.h"
 #include "Param/Distrib/ParameterSample.h"
 #include <cmath>
@@ -64,22 +65,16 @@ IRangedDistribution::~IRangedDistribution() = default;
 
 std::vector<ParameterSample> IRangedDistribution::generateSamples(double mean, double stddev) const
 {
-    auto generator = distribution(mean, stddev);
-    if (!generator->isDelta())
-        return generator->equidistantSamples(m_n_samples, m_sigma_factor, m_limits);
+    if (stddev == 0) // it's a delta distribution
+        return {{mean, 1.0}};
+    ASSERT(stddev > 0);
 
-    // handling the case of delta distributions
-    auto samples = generator->equidistantSamples(1, m_sigma_factor, m_limits);
-    ParameterSample& sample = samples[0]; // there is only one element in the vector
-    sample.weight = 1.0;
-    return std::vector<ParameterSample>(1, sample);
+    auto generator = distribution(mean, stddev);
+    return generator->equidistantSamples(m_n_samples, m_sigma_factor, m_limits);
 }
 
 IDistribution1D* IRangedDistribution::distribution(double mean, double stddev) const
 {
-    if (stddev < 0.0)
-        throw std::runtime_error(
-            "Error in IRangedDistribution::distribution: standard deviation is less than zero");
     return distribution_impl(mean, stddev);
 }