First, solve a linear system without a loop in time. Is the solution the same? Second, you problem are not the same (different values for mu1
, mu2
, and mu3
). Third, do not use fancy preconditioner until you sort this out and stick to an exact factorization instead of GAMG.
2 Likes