Fix bug in CGNR with only AHA, b is now used to init x0, similar to F…

…ISTA and co
JuliaImageRecon · Aug 13, 2024 · 6b17ba9 · 6b17ba9
1 parent ef0c58f
commit 6b17ba9
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/CGNR.jl b/src/CGNR.jl
@@ -131,7 +131,7 @@ end
 
 initCGNR(x₀, A, b) = mul!(x₀, adjoint(A), b)
 #initCGNR(x₀, prod::ProdOp{T, <:WeightingOp, matT}, b) where {T, matT} = mul!(x₀, adjoint(prod.B), b)
-initCGNR(x₀, ::Nothing, b) = x₀ .= one(eltype(x₀))
+initCGNR(x₀, ::Nothing, b) = x₀ .= b
 
 solverconvergence(state::CGNRState) = (; :residual => norm(state.x₀))