I have a question. Why is the default value of `reltol` often picked as `sqrt(eps)` and not e.g. `eps`? Is there a particular reason? https://github.com/JuliaLinearAlgebra/IterativeSolvers.jl/blob/40d9e9db7e0c6ee2047d121125005b217023baf9/src/cg.jl#L122 cc @Elise-palethorpe