<div dir="ltr">Hi PETSc,<br><div><br></div><div>I am trying to implement the Time stepping routines in my model, i have a working runge-kutta time scheme that goes to the following steps:</div><div><br></div><div><ul><li>Interpolate u,v,w to the time advancing variable position.</li><li>Calculate nonlinear coefficients and advect velocities with a forward-backward shock capturing scheme.</li><li>Calculate the variable laplacian </li><li>Sum terms to create RHS (nonlinear advected velocities + laplacian)</li><li>Finally, the runge kutta integration is done in a typical way that looks like:</li></ul><div>                   newtemp(t+1) = prevtemp(t) + dt*RHS</div></div><div><br></div><div><br></div><div>So my questions are:</div><div><ul><li>I think my problem is nonlinear, but is being made linearized by the advecting scheme, is this correct? this is to know if i should use the linear or nonlinear routines for TS.</li><li>How do i know what are the appropriate routines i should be using here? so far i think i should use the following:</li></ul><div>call TSCreate(PETSC_COMM_WORLD,ts,ierr)<br>call TSSetProblemType(ts,TS_LINEAR,ierr)<br>call TSSetTimeStep(ts,dt,ierr)<br><br>call TSSetFromOptions(ts,ierr)<br></div><div><br></div><div>call TSSetRHSFunction(ts,NULL,TSComputeRHSFunctionLinear,NULL,ierr)<br>call TSSolve(ts,loctemperature,ierr)<br><br></div></div><div>Should i use call TSSetRHSJacobian for the temperature jacobian in this case? </div><div><br></div><div>I am using <a href="https://www.mcs.anl.gov/petsc/petsc-current/src/ts/examples/tutorials/ex4.c.html">https://www.mcs.anl.gov/petsc/petsc-current/src/ts/examples/tutorials/ex4.c.html</a> as a  general guide, is there a more appropriate example?</div><div><br></div><div>The dt value and total timesteps are controlled by the model,</div><div><br></div><div>Thanks for your help,</div><div><br></div><div><br></div></div>