cleanup

Author: vpuri3

src/common.jl

@@ -56,9 +56,9 @@ SciMLBase.init(prob::LinearProblem, args...; kwargs...) = SciMLBase.init(prob,no
 
 function SciMLBase.init(prob::LinearProblem, alg::Union{SciMLLinearSolveAlgorithm,Nothing}, args...;
                         alias_A = false, alias_b = false,
-                        abstol=√eps(eltype(prob.A)), # TODO handle when eps()
-                        reltol=√eps(eltype(prob.A)), # not defined on eltype
-                        maxiters=length(prob.b),     # like Int of Complex
+                        abstol=√eps(eltype(prob.A)),
+                        reltol=√eps(eltype(prob.A)),
+                        maxiters=length(prob.b),
                         verbose=false,
                         Pl = nothing,
                         Pr = nothing,

src/pardiso.jl

@@ -76,10 +76,11 @@ function init_cacheval(alg::PardisoJL, cache::LinearCache)
         Pardiso.pardiso(solver, u, A, b)
     end
 
-    #   ipram/dpram
+    """
+    use ipram/dpram to set tolerances???
+    """
 #   abstol = cache.abstol
 #   reltol = cache.reltol
-#   kwargs = (abstol=abstol, reltol=reltol)
 
     return solver
 end

test/runtests.jl

@@ -145,8 +145,7 @@ end
                 MKLPardisoIterate(),
                )
 
-        u = solve(prob, alg; verbose=true)
-
+        u = solve(prob, alg; cache_kwargs...)
         @test A * u ≈ b
     end