From caf884a8e15eda6c9c32020a2fd80e9317beaf24 Mon Sep 17 00:00:00 2001
From: cilix42 <cilix42@users.sourceforge.net>
Date: Thu, 3 Jul 2008 09:39:33 +0000
Subject: [PATCH] Partial 2geom update (anticipating changes that will be part
 of the next major update anyway, so overwriting them on next update is not a
 problem)

---
 src/2geom/numeric/linear_system.h |  4 ++--
 src/2geom/piecewise.h             | 18 +++++++++++++++++-
 2 files changed, 19 insertions(+), 3 deletions(-)

diff --git a/src/2geom/numeric/linear_system.h b/src/2geom/numeric/linear_system.h
index 9cb521eb2..7699c5224 100644
--- a/src/2geom/numeric/linear_system.h
+++ b/src/2geom/numeric/linear_system.h
@@ -6,8 +6,8 @@
 
 #include <gsl/gsl_linalg.h>
 
-#include "numeric/matrix.h"
-#include "numeric/vector.h"
+#include "2geom/numeric/matrix.h"
+#include "2geom/numeric/vector.h"
 
 
 namespace Geom { namespace NL {
diff --git a/src/2geom/piecewise.h b/src/2geom/piecewise.h
index 8992c0097..96e64bcf9 100644
--- a/src/2geom/piecewise.h
+++ b/src/2geom/piecewise.h
@@ -704,7 +704,23 @@ Piecewise<T> derivative(Piecewise<T> const &a) {
 }
 
 std::vector<double> roots(Piecewise<SBasis> const &f);
-Piecewise<SBasis> reverse(Piecewise<SBasis> const &f);
+
+template<typename T>
+Piecewise<T> reverse(Piecewise<T> const &f) {
+    Piecewise<T> ret = Piecewise<T>();
+    ret.cuts.resize(f.cuts.size());
+    ret.segs.resize(f.segs.size());
+    double start = f.cuts[0];
+    double end = f.cuts.back();
+    for (unsigned i = 0; i < f.cuts.size(); i++) {
+        double x = f.cuts[f.cuts.size() - 1 - i];
+        ret.cuts[i] = end - (x - start);
+    }
+    for (unsigned i = 0; i < f.segs.size(); i++)
+        ret.segs[i] = reverse(f[f.segs.size() - i - 1]);
+    return ret;
+}
+
 
 }
 
-- 
2.30.2