<div dir="ltr"><div>I noticed that when I am solving a problem with the ts and I am *not* using a da, if I want to use an implicit time stepping routine:</div><div>1. I have to explicitly provide the Jacobian</div><div>2. When I do provide the Jacobian, if I want to access the elements of x(t) to construct f(t,x), I need to use a const PetscScalar and a VecGetArrayRead to get it to work.  <br></div><div>3.  My code works without declaring const when I'm using an explicit scheme.<br></div><div><br></div><div>In contrast, if I solve a problem using a da, my code works, I can use implicit schemes without having to provide the Jacobian, and I don't have to use const anywhere.</div><div><br></div><div>Can someone clarify what is expected/preferred?<br></div><div><br>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature">gideon</div></div></div>