1 files changed, 52 insertions, 39 deletions

diff --git a/qfunctionconfigurator/functionconfig.cpp b/qfunctionconfigurator/functionconfig.cpp
index 38a9e703..dc50a4d1 100644
--- a/qfunctionconfigurator/functionconfig.cpp
+++ b/qfunctionconfigurator/functionconfig.cpp
@@ -19,20 +19,26 @@ Map::Map() :
     _mutex(QMutex::Recursive),
     activep(false),
     max_x(0),
-    max_y(0)
+    max_y(0),
+    lazy_reload(true)
 {
 }
 
 Map::num Map::getValue(Map::num x) {
     QMutexLocker foo(&_mutex);
+    if (lazy_reload)
+    {
+        lazy_reload = false;
+        reload();
+    }
     num q  = x * precision();
     int xi = (int)q;
     num yi = getValueInternal(xi);
-    num yiplus1 = getValueInternal(xi+1);
+    num yiplus1 = getValueInternal(xi + (x < 0 ? -1 : 1));
     num f = (q-xi);
-    num ret = yiplus1 * f + yi * (1.0f - f); // at least do a linear interpolation.
-    last_input_value.setX(x);
-    last_input_value.setY(ret);
+    num ret = yiplus1 * f + yi * (1 - f); // at least do a linear interpolation.
+    last_input_value.setX(std::abs(x));
+    last_input_value.setY(std::abs(ret));
     return ret;
 }
 
@@ -44,13 +50,13 @@ bool Map::getLastPoint(QPointF& point ) {
 
 Map::num Map::getValueInternal(int x) {
     num sign = x < 0 ? -1 : 1;
-    x = abs(x);
+    x = std::abs(x);
     num ret;
     int sz = cur.data.size();
     if (sz == 0)
         ret = 0;
     else
-        ret = cur.data[std::min<unsigned>(x, sz-1)];
+        ret = cur.data[std::min<unsigned>(x, sz-1)] * max_y / integral_max;
     return ret * sign;
 }
 
@@ -75,29 +81,29 @@ void Map::reload() {
         QList<QPointF> input = cur.input;
         auto& data = cur.data;
 
-        data = std::vector<num>(value_count);
+        data = std::vector<integral>(value_count);
         const int mult = precision();
         
         const int sz = data.size();
         
         for (int i = 0; i < sz; i++)
-            data[i] = -1;
+            data[i] = integral_max;
         
         if (input.size() == 1)
         {
             for (int k = 0; k < input[0].x() * mult; k++) {
                 if (k < sz)
-                    data[k] = input[0].y() * k / (input[0].x() * mult);
+                    data[k] = input[0].y() * k * integral_max / (input[0].x() * mult) / max_y ;
             }
         }
         else if (input[0].x() > 1e-2)
             input.prepend(QPointF(0, 0));
         
         for (int i = 0; i < sz; i++) {
-            QPointF p0 = ensureInBounds(input, i - 1);
-            QPointF p1 = ensureInBounds(input, i);
-            QPointF p2 = ensureInBounds(input, i + 1);
-            QPointF p3 = ensureInBounds(input, i + 2);
+            const QPointF p0 = ensureInBounds(input, i - 1);
+            const QPointF p1 = ensureInBounds(input, i);
+            const QPointF p2 = ensureInBounds(input, i + 1);
+            const QPointF p3 = ensureInBounds(input, i + 2);
 
             using n = double;
             const n p0_x = p0.x(), p1_x = p1.x(), p2_x = p2.x(), p3_x = p3.x();
@@ -105,35 +111,39 @@ void Map::reload() {
             
             // multiplier helps fill in all the x's needed
             const int mult_ = mult * 20;
-            int end = std::min<int>(sz, p2.x() * mult_);
-            int start = p1.x() * mult;
+            const int end = std::min<int>(sz, p2.x() * mult_);
+            const int start = p1.x() * mult;
             const n max = end - start;
             
             for (int j = start; j < end; j++) {
-                n t = (j - start) / max;
-                n t2 = t*t;
-                n t3 = t*t*t;
-
-                int x = .5 * ((2. * p1_x) +
-                              (-p0_x + p2_x) * t +
-                              (2. * p0_x - 5. * p1_x + 4. * p2_x - p3_x) * t2 +
-                              (-p0_x + 3. * p1_x - 3. * p2_x + p3_x) * t3)
-                        * mult;
+                const n t = (j - start) / max;
+                const n t2 = t*t;
+                const n t3 = t*t*t;
+
+                const int x = .5 * ((2. * p1_x) +
+                                    (-p0_x + p2_x) * t +
+                                    (2. * p0_x - 5. * p1_x + 4. * p2_x - p3_x) * t2 +
+                                    (-p0_x + 3. * p1_x - 3. * p2_x + p3_x) * t3)
+                                 * mult;
+
+                if (x < 0 || x >= sz || data[x] != integral_max)
+                    continue;
                 
-                num y = .5 * ((2. * p1_y) +
-                              (-p0_y + p2_y) * t +
-                              (2. * p0_y - 5. * p1_y + 4. * p2_y - p3_y) * t2 +
-                              (-p0_y + 3. * p1_y - 3. * p2_y + p3_y) * t3);
-                
-                if (x >= 0 && x < sz)
-                    data[x] = y;
+                const n y = .5 * ((2. * p1_y) +
+                                  (-p0_y + p2_y) * t +
+                                  (2. * p0_y - 5. * p1_y + 4. * p2_y - p3_y) * t2 +
+                                  (-p0_y + 3. * p1_y - 3. * p2_y + p3_y) * t3);
+
+                const n y_ = std::min<n>(max_y, std::max<n>(y, 0));
+
+                data[x] = y_ / max_y * integral_max;
             }
         }
         
         num last = 0;
         for (int i = 0; i < sz; i++)
         {
-            if (data[i] < 0)
+            if (data[i] == integral_max)
                 data[i] = last;
             last = data[i];
         }
@@ -147,14 +157,15 @@ void Map::removePoint(int i) {
     if (i >= 0 && i < cur.input.size())
     {
         cur.input.removeAt(i);
-        reload();
+        lazy_reload = true;
     }
 }
 
 void Map::addPoint(QPointF pt) {
     QMutexLocker foo(&_mutex);
     cur.input.append(pt);
-    reload();
+    lazy_reload = true;
+    qStableSort(cur.input.begin(), cur.input.end(), sortFn);
 }
 
 void Map::movePoint(int idx, QPointF pt) {
@@ -162,7 +173,8 @@ void Map::movePoint(int idx, QPointF pt) {
     if (idx >= 0 && idx < cur.input.size())
     {
         cur.input[idx] = pt;
-        reload();
+        lazy_reload = true;
+        // we don't allow points to be reodered, so no sort here
     }
 }
 
@@ -173,8 +185,9 @@ const QList<QPointF> Map::getPoints() {
 
 void Map::invalidate_unsaved_settings()
 {
+    QMutexLocker foo(&_mutex);
     cur = saved;
-    reload();
+    lazy_reload = true;
 }
 
 void Map::loadSettings(QSettings& settings, const QString& title) {
@@ -203,7 +216,7 @@ void Map::loadSettings(QSettings& settings, const QString& title) {
         points.append(QPointF(maxInput(), maxOutput()));
     
     cur.input = points;
-    reload();
+    lazy_reload = true;
     saved = cur;
 }
 
@@ -235,6 +248,6 @@ void Map::saveSettings(QSettings& settings, const QString& title) {
 
 int Map::precision() const {
     if (cur.input.size())
-        return value_count / std::max<num>(1.f, (cur.input[cur.input.size() - 1].x()));
+        return value_count / std::max<num>(1, (cur.input[cur.input.size() - 1].x()));
     return 1;
 }