update A and b in solve

Author: pjaap

src/solvers.jl

@@ -381,14 +381,22 @@ function solve_linear_system!(A, b, sol, soltemp, residual, unknowns, freedofs,
             end
         end
 
-        if !SC.parameters[:initialized]
-            ## init solver if not done before
+
+        if SC.parameters[:initialized]
+            # set/update the linear system
+            SC.linsolver.b = linsolve_b
+            if linsolve_needs_matrix
+                SC.linsolver.A = linsolve_A
+            end
+        else
+            # init solver if not done before
             @timeit timer "linear solver" begin
                 init_linear_solver!(SC, linsolve_A, linsolve_b, timer)
             end
             SC.parameters[:initialized] = true
         end
 
+        # now the linear solver is definitely ready
         linsolve = SC.linsolver
 
         # Solve linear system