<div class="gmail_quote">On Sat, Jul 9, 2011 at 17:43, TAY wee-beng <span dir="ltr">&lt;<a href="mailto:zonexo@gmail.com">zonexo@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div bgcolor="#ffffff" text="#000000">Sorry about that. I just run interactively to get the output.<br>
    <br>
    Now I issued :<br>
    <br>
    mpiexec -np 1 ./a.out  -ksp_type gmres -ksp_gmres_restart 1000
    -ksp_gmres_modifiedgramschmidt -ksp_monitor_true_residual -ksp_view
    -pc_type asm -sub_pc_type lu &gt; log1</div></blockquote></div><br><div>Your matrix is singular which implies a bug in your assembly code, probably in the way you impose boundary conditions. Switch to a very small problem size (e.g. 4 cells) and debug it until the matrix is correct and non-singular.</div>
<div><br></div><div>You should also check the error codes so the code doesn&#39;t try to keep going after a fatal error.</div>