Tolerande for the solver normalized by input as in sfcf

parent 073a92b8
......@@ -42,6 +42,8 @@ function CG!(si, U, A, dpar::DiracParam, lp::SpaceParm, dws::DiracWorkspace{T},
fill!(si,zero(eltype(si)))
err = 0.0
tol = tol * norm
iterations = 0
sumf = scalar_field(Complex{T}, lp)
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment