Use prob.b instead of prob.A to promote abstol and reltol in case A is an operator

Author: ma-sadeghi

src/common.jl

@@ -115,8 +115,8 @@ function SciMLBase.init(prob::LinearProblem, alg::SciMLLinearSolveAlgorithm,
     args...;
     alias_A = default_alias_A(alg, prob.A, prob.b),
     alias_b = default_alias_b(alg, prob.A, prob.b),
-    abstol = default_tol(eltype(prob.A)),
-    reltol = default_tol(eltype(prob.A)),
+    abstol = default_tol(eltype(prob.b)),
+    reltol = default_tol(eltype(prob.b)),
     maxiters::Int = length(prob.b),
     verbose::Bool = false,
     Pl = IdentityOperator(size(prob.A)[1]),