<div dir="ltr">I came across this thing recently, and I couldn't figure out where the issue could be.<div><br></div><div>The problem I'm solving is a simple DG advection, the ode is M*udot = K*u+b, M is diagonal.</div><div><br></div><div>Attached is a MWE that reproduces the problem.<br></div><div><br></div><div><div>The problem is formed in two different cases depending on the command line option -m_lhs</div><div><br></div><div>a) -m_lhs 1 : F(u,udot,t) = M*udot, G(u,t) = K*u+b</div><div><div>b) -m_lhs 0 : F(u,udot,t) = udot, G(u,t) = M^-1(K*u+b)</div></div><div><br></div><div>Using option b) and RK4, the solution is ok.</div><div>If run with any implicit TS method except arkimex, no matter if I'm choosing option a) or b), the solution is always very close (say, final error < 0.05) to the expected one (computed with BDF).</div><div><br></div><div>When using ARKIMEX, case b) gives a good solution, but not case a). In fact, the solution does not seem to be advected at all in this case.</div><div><br></div><div>I was wondering if I'm doing something wrong or there's a bug in the ARKIMEX implementation.</div><div><br></div><div>I also noticed that, using  -ts_arkimex_fully_implicit does not produce the same output for case a) and b). Shouldn't they produce the same method with this option?</div><div><br></div><div>Thanks,</div>-- <br><div class="gmail_signature">Stefano</div>
</div></div>