Fix the bug of overwriting original signal

felixcremer · May 7, 2021 · 2d25f77 · 2d25f77
1 parent 908586a
commit 2d25f77
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/sift.jl b/src/sift.jl
@@ -65,14 +65,14 @@ function iterate(it::SiftIterable, iteration::Int=start(it))
     it.y .-= subs
 
     # Return the residual at item and iteration number as state
-    it.tol, iteration + 1
+    it.residual, iteration + 1
 end
 
 # Utility functions
 
 function sift_iterator(y, x;
                       reltol::Real = 0.1,
-                      maxiter::Int = 10,
+                      maxiter::Int = 6,
                       stop_steps::Int = 4)
     maxs, mins = Int[], Int[]
     crosses = Int[]
@@ -121,7 +121,7 @@ function sift!(imf, y, x;
              verbose::Bool = false,
              kwargs...)
 
-    imf = y
+    imf .= y
     # Actually perform sifting
     iterable = sift_iterator(imf, x; reltol, maxiter, stop_steps)