<div dir="ltr"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br><br>
<br>
If you increase resolution at differentiating a discontinuity, you'd<br>
expect to become more and more confident that the error is large.<br></blockquote><div><br></div><div>I am trying to use smooth functions but my first function might have had a discontinuity. I'm now just using a sin function that is perfectly smooth.</div><div><br></div><div>In looking at this some more, the dynamics/stiffness of the system, when these (cold) sources are added, is just stiff and this behavior seems reasonable. I've found that with a dt-min of 0.1 and ts_tol 1.e-1 it works OK and the results are the same in the eyeball norm as dt-min 1e-3 and ts_tol 1.e-3. So I'm gonna call that converged.</div><div><br></div><div>Thanks,<br></div><div>Mark</div><div><br></div><div> </div></div></div>