Fix Kaczmarz tikhonov test

JuliaImageRecon · Nov 26, 2024 · 178a581 · 178a581
1 parent 6e5c1db
commit 178a581
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/test/testKaczmarz.jl b/test/testKaczmarz.jl
@@ -50,11 +50,11 @@ Random.seed!(12345)
 
         # @show A, x, regMatrix
         # use regularization matrix
-        S = createLinearSolver(solver, arrayType(A), iterations=100, reg=[L2Regularization(arrayType(regMatrix))])
+        S = createLinearSolver(Kaczmarz, arrayType(A), iterations=100, reg=[L2Regularization(arrayType(regMatrix))])
         x_matrix = Array(solve!(S, arrayType(b)))
 
         # use standard reconstruction
-        S = createLinearSolver(solver, arrayType(A * Diagonal(1 ./ sqrt.(regMatrix))), reg = [L2Regularization(1.0)], iterations=100)
+        S = createLinearSolver(Kaczmarz, arrayType(A * Diagonal(1 ./ sqrt.(regMatrix))), reg = [L2Regularization(1.0)], iterations=100)
         x_approx = Array(solve!(S, arrayType(b))) ./ sqrt.(regMatrix)
 
         # test
@@ -63,10 +63,10 @@ Random.seed!(12345)
 
         # Compare reg. matrix of equal elements to standard reco
         λ = rand()
-        S = createLinearSolver(solver, arrayType(A), iterations=100, reg=[L2Regularization(λ)])
+        S = createLinearSolver(Kaczmarz, arrayType(A), iterations=100, reg=[L2Regularization(λ)])
         x_standard = Array(solve!(S, arrayType(b)))
 
-        S = createLinearSolver(solver, arrayType(A), iterations=100, reg=[L2Regularization(fill(λ, N))])
+        S = createLinearSolver(Kaczmarz, arrayType(A), iterations=100, reg=[L2Regularization(fill(λ, N))])
         x_matrix = Array(solve!(S, arrayType(b)))
         @test isapprox(x_standard, x_matrix)
       end