adds initial value optimization

NFFT · Jun 16, 2021 · a547a45 · a547a45 · mischmi96 · Jun 16, 2021
1 parent 36ef144
commit a547a45
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 2 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "ANOVAapprox"
 uuid = "5e027bd6-ab01-4733-8320-e0223e929ebb"
 authors = ["Michael Schmischke <michael.schmischke@math.tu-chemnitz.de> and Felix Bartel <felix.bartel@math.tu-chemnitz.de>"]
-version = "1.0.23"
+version = "1.0.24"
 
 [deps]
 Distributed = "8ba89e20-285c-5b6f-9357-94700520ee1b"

diff --git a/src/nperiodic_approx_scat_lsqr.jl b/src/nperiodic_approx_scat_lsqr.jl
@@ -48,7 +48,8 @@ function approximate( approx::nperiodic_approx_scat_lsqr{d,ds}; max_iter::Int64=
         W = ones( M )
     end
 
-    for L in lambda 
+    for i = 1:length(lambda) 
+        L = lambda[i] 
         if verbose
             println( "lambda = ", L )
         end
@@ -59,6 +60,11 @@ function approximate( approx::nperiodic_approx_scat_lsqr{d,ds}; max_iter::Int64=
             M, n )
 
         tmp = zeros( ComplexF64, n )
+
+        if i > 1 
+            tmp = approx.fc[lambda[i-1]].data
+        end
+
         lsqr!( tmp, F, W .* approx.y, maxiter = max_iter, verbose=verbose, damp=sqrt(L) )
         approx.fc[L] = GroupedCoeff(approx.trafo.setting, tmp)
     end