
<html dir="ltr">

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" id="owaParaStyle"></style>

</head>

<body fpstyle="1" ocsi="0">

<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">I'm using KSP for 2D Navier Stokes, finite difference. I have DM object created with DMDA_STENTIL_STAR. I'm only using one single process to run the code. It works fine, the linear

 system converges to the correct solution with correct order, as expected. However, if I change it to DMDA_STENTIL_BOX, it does not converge, the residue \|Ax-b\|_2 remains large. 

<div><br>

</div>

<div>I have checked that the rhs for both cases are exactly the same, and the matrices also seem to be the same because I have checked the output of A*y, where y is the analytic solution, and both cases give exactly the same answer. If the linear system is

 exactly the same, then it seems very strange to me that <span style="font-size: 10pt;">DMDA_STENTIL_STAR works while </span><span style="font-size: 10pt;">DMDA_STENTIL_BOX does not work, for I use the same solver, there shouldn't be any difference. </span></div>

<div><span style="font-size: 10pt;"><br>

</span></div>

<div>If you happen to encounter this kind of bug before and have some idea what's going on, please let me know. Thank you very much! If you want to look at the code I've written, feel free to ask. </div>

</div>

</body>

</html>