From 21c2135588a6a76b6dcc4defc6e64a3d9032db59 Mon Sep 17 00:00:00 2001 From: Stanislaw Halik Date: Sun, 26 Jul 2015 10:50:36 +0200 Subject: qfc: interpolate the right way for negative values --- qfunctionconfigurator/functionconfig.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'qfunctionconfigurator/functionconfig.cpp') diff --git a/qfunctionconfigurator/functionconfig.cpp b/qfunctionconfigurator/functionconfig.cpp index 38a9e703..825aa3c1 100644 --- a/qfunctionconfigurator/functionconfig.cpp +++ b/qfunctionconfigurator/functionconfig.cpp @@ -28,7 +28,7 @@ Map::num Map::getValue(Map::num x) { 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); -- cgit v1.2.3 From b34bb2dc35e2d11bd6473a75356a49303bbeb556 Mon Sep 17 00:00:00 2001 From: Stanislaw Halik Date: Sun, 26 Jul 2015 10:51:08 +0200 Subject: qfc: we need positive retval here --- qfunctionconfigurator/functionconfig.cpp | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) (limited to 'qfunctionconfigurator/functionconfig.cpp') diff --git a/qfunctionconfigurator/functionconfig.cpp b/qfunctionconfigurator/functionconfig.cpp index 825aa3c1..b4806e19 100644 --- a/qfunctionconfigurator/functionconfig.cpp +++ b/qfunctionconfigurator/functionconfig.cpp @@ -30,9 +30,9 @@ Map::num Map::getValue(Map::num x) { num yi = getValueInternal(xi); 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; } -- cgit v1.2.3 From ff290993484e6708ec57fbc4c4d4eeaaae379c4f Mon Sep 17 00:00:00 2001 From: Stanislaw Halik Date: Sun, 26 Jul 2015 10:53:33 +0200 Subject: qfc: use integers as storage, not floats Float mantissa has 23 bits. 9 bits get wasted per value. Instead, use uint16_t that doesn't lose much precision, but saves half the space. --- qfunctionconfigurator/functionconfig.cpp | 62 +++++++++++++++++--------------- qfunctionconfigurator/functionconfig.h | 6 +++- 2 files changed, 38 insertions(+), 30 deletions(-) (limited to 'qfunctionconfigurator/functionconfig.cpp') diff --git a/qfunctionconfigurator/functionconfig.cpp b/qfunctionconfigurator/functionconfig.cpp index b4806e19..31e3641d 100644 --- a/qfunctionconfigurator/functionconfig.cpp +++ b/qfunctionconfigurator/functionconfig.cpp @@ -44,13 +44,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(x, sz-1)]; + ret = cur.data[std::min(x, sz-1)] * max_y / integral_max; return ret * sign; } @@ -75,29 +75,29 @@ void Map::reload() { QList input = cur.input; auto& data = cur.data; - data = std::vector(value_count); + data = std::vector(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 / (input[0].x() * mult) / max_y * integral_max; } } 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 +105,39 @@ void Map::reload() { // multiplier helps fill in all the x's needed const int mult_ = mult * 20; - int end = std::min(sz, p2.x() * mult_); - int start = p1.x() * mult; + const int end = std::min(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(max_y, std::max(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]; } @@ -235,6 +239,6 @@ void Map::saveSettings(QSettings& settings, const QString& title) { int Map::precision() const { if (cur.input.size()) - return value_count / std::max(1.f, (cur.input[cur.input.size() - 1].x())); + return value_count / std::max(1, (cur.input[cur.input.size() - 1].x())); return 1; } diff --git a/qfunctionconfigurator/functionconfig.h b/qfunctionconfigurator/functionconfig.h index 446c73ac..68ec2869 100644 --- a/qfunctionconfigurator/functionconfig.h +++ b/qfunctionconfigurator/functionconfig.h @@ -14,16 +14,19 @@ #include #include #include +#include #include "opentrack-compat/qcopyable-mutex.hpp" class Map { private: static constexpr int value_count = 5000; using num = float; + using integral = std::uint16_t; + static constexpr integral integral_max = std::numeric_limits::max(); struct State { QList input; - std::vector data; + std::vector data; }; int precision() const; @@ -45,6 +48,7 @@ public: { setMaxInput(maxx); setMaxOutput(maxy); + reload(); } num getValue(num x); -- cgit v1.2.3 From f54f96a67572d770d3eceb33d5df18f0b64016ce Mon Sep 17 00:00:00 2001 From: Stanislaw Halik Date: Sun, 26 Jul 2015 11:24:52 +0200 Subject: qfc: add missing mutex --- qfunctionconfigurator/functionconfig.cpp | 1 + 1 file changed, 1 insertion(+) (limited to 'qfunctionconfigurator/functionconfig.cpp') diff --git a/qfunctionconfigurator/functionconfig.cpp b/qfunctionconfigurator/functionconfig.cpp index 31e3641d..cf2516ca 100644 --- a/qfunctionconfigurator/functionconfig.cpp +++ b/qfunctionconfigurator/functionconfig.cpp @@ -177,6 +177,7 @@ const QList Map::getPoints() { void Map::invalidate_unsaved_settings() { + QMutexLocker foo(&_mutex); cur = saved; reload(); } -- cgit v1.2.3 From de900bb1159af9f85d6f874925648844566363e3 Mon Sep 17 00:00:00 2001 From: Stanislaw Halik Date: Sun, 26 Jul 2015 11:26:13 +0200 Subject: qfc: reload spline data lazily --- qfunctionconfigurator/functionconfig.cpp | 22 +++++++++++++++------- qfunctionconfigurator/functionconfig.h | 5 +++-- 2 files changed, 18 insertions(+), 9 deletions(-) (limited to 'qfunctionconfigurator/functionconfig.cpp') diff --git a/qfunctionconfigurator/functionconfig.cpp b/qfunctionconfigurator/functionconfig.cpp index cf2516ca..dc50a4d1 100644 --- a/qfunctionconfigurator/functionconfig.cpp +++ b/qfunctionconfigurator/functionconfig.cpp @@ -19,12 +19,18 @@ 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); @@ -87,7 +93,7 @@ void Map::reload() { { for (int k = 0; k < input[0].x() * mult; k++) { if (k < sz) - data[k] = input[0].y() * k / (input[0].x() * mult) / max_y * integral_max; + data[k] = input[0].y() * k * integral_max / (input[0].x() * mult) / max_y ; } } else if (input[0].x() > 1e-2) @@ -151,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) { @@ -166,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 } } @@ -179,7 +187,7 @@ void Map::invalidate_unsaved_settings() { QMutexLocker foo(&_mutex); cur = saved; - reload(); + lazy_reload = true; } void Map::loadSettings(QSettings& settings, const QString& title) { @@ -208,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; } diff --git a/qfunctionconfigurator/functionconfig.h b/qfunctionconfigurator/functionconfig.h index 68ec2869..74e20cf0 100644 --- a/qfunctionconfigurator/functionconfig.h +++ b/qfunctionconfigurator/functionconfig.h @@ -40,6 +40,7 @@ private: double max_y; State cur, saved; + bool lazy_reload; public: double maxInput() const { return max_x; } double maxOutput() const { return max_y; } @@ -48,7 +49,7 @@ public: { setMaxInput(maxx); setMaxOutput(maxy); - reload(); + lazy_reload = true; } num getValue(num x); @@ -57,7 +58,7 @@ public: void removeAllPoints() { QMutexLocker foo(&_mutex); cur.input.clear(); - reload(); + lazy_reload = true; } void addPoint(QPointF pt); -- cgit v1.2.3