<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr; font-family: Tahoma; color: rgb(0, 0, 0); font-size: 13px;">
<div style=""> Hi,<br>
<br>
Last time I was suggested by Jed using the DAE solver to build my nonlinear CFD code based on TS. For incompressible N-S equations, I think it's difficult to derive the analytical Jacobian formulation, and we need a finite difference approximation or matrix-free
method. <br>
<br>
What I previously done was integrating the incompressible N-S equations into tutorial code ts/ex7.c . It simply uses an argument "TSDefaultComputeJacobianColor" to finish the Jacobian issue. When I switch to TSSetIFunction() and TSSetIJacobian(), this argument
doesn't work anymore. And if I remove all Jacobian-related lines and use -snes_mf command option, it prompts out error message "Must Set Jacobian!". Do you know how to resolve this for both methods?<br>
<br>
<br>
Thank you.<br>
<br>
<br>
Regards,<br>
<br>
<br>
<br>
Zhisong Li<br>
</div>
</div>
</body>
</html>