safer convolution

lib/LatAnalyze/Core/Plot.cpp

@@ -511,6 +511,27 @@ PlotImpulses::PlotImpulses(const DVec &x, const DVec &y)
     setCommand("'" + tmpFileName + "' u 1:2 w impulses");
 }
 
+// PlotSteps constructor ////////////////////////////////////////////////////
+PlotSteps::PlotSteps(const DVec &x, const DVec &y)
+{
+    if (x.rows() != y.rows())
+    {
+        LATAN_ERROR(Size, "x and y vector does not have the same size");
+    }
+
+    DMat   d(x.rows(), 2);
+    string tmpFileName;
+
+    for (Index i = 0; i < x.rows(); ++i)
+    {
+        d(i, 0) = x(i);
+        d(i, 1) = y(i);
+    }
+    tmpFileName = dumpToTmpFile(d);
+    pushTmpFile(tmpFileName);
+    setCommand("'" + tmpFileName + "' u 1:2 w steps");
+}
+
 // PlotGrid constructor ////////////////////////////////////////////////////////
 PlotGrid::PlotGrid(const DVec &x, const DVec &y, const DMat &value)
 {

lib/LatAnalyze/Core/Plot.hpp

@@ -191,6 +191,15 @@ public:
     virtual ~PlotHistogram(void) = default;
 };
 
+class PlotSteps: public PlotObject
+{
+public:
+    // constructor
+    PlotSteps(const DVec &x, const DVec &y);
+    // destructor
+    virtual ~PlotSteps(void) = default;
+};
+
 class PlotImpulses: public PlotObject
 {
 public:

lib/LatAnalyze/Numerical/DWT.cpp

@@ -38,6 +38,10 @@ void filterConvolution(MatType &out, const MatType &data,
 {
     Index n = data.rows(), nf = n*filter.size();
 
+    if (&out == &data)
+    {
+        LATAN_ERROR(Argument, "filter convolution does not support in-place operation");
+    }
     out.resizeLike(data);
     out.fill(0.);
     for (unsigned int i = 0; i < filter.size(); ++i)