<div dir="ltr"><div>Hi, <br><br></div>I have a few questions about ex13 from TS explaining how to solve the following PDE in 2D :<pre width="80"><font color="#B22222">u_t = uxx + uyy</font></pre>

<div><div>1) Can I get a reference/explanation for how the RHSJacobian function is derived ?<br><br></div><div>2) If the problem is linear, why is<a name="line147"> TSSetProblemType never set to TS_LINEAR?<br><br></a></div><div>3) Since this is just a 2D version of ex3, can't this be set up by only providing a function for the evaluation for TSSetRHSJacobian with the TSSetRHSFunction set to TSComputeRHSFunctionLinear ? <br></div><div><br>If this is the case and a matrix could be written that equivalently states the problem as u_t = A(t)*u, then how would one make such a matrix A(t) ? (In essence, what I'm asking is how does petsc define 2D vectors. Would A(t) be the familiar block tridiagonal matrix ?)<br>

</div><div><br></div><div>4) In lines 153/154, why does the outer loop iterate over columns and the inner loop iterate over rows? Does this again have something to do with how petsc stores 2d vectors?<br></div><div><br></div><div>Thank You,<br></div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div style="font-size:12.8px">Sajid Ali<br></div><div style="font-size:12.8px">Applied Physics<br></div><div style="font-size:12.8px">Northwestern University</div></div></div></div></div></div>