<div class="gmail_quote">On Tue, May 17, 2011 at 01:25, Tian(ICT) <span dir="ltr">&lt;<a href="mailto:rongtian@ncic.ac.cn">rongtian@ncic.ac.cn</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="im"><div><font size="2">pgfem -snes_monitor -snes_mf_operator -mat_mffd_type ds 
-ksp_type fgmres -ksp_monitor_true_residual -ksp_gmres_modifiedgramschmidt 
-ksp_rtol 1e-10 <font color="#ff0000">-pc_type none</font></font></div>
</div><div><font size="2"><font color="#ff0000"><br></font><div class="im">atol=1e-050, rtol=1e-008, 
stol=1e-008, maxit=50, maxf=10000<br>  0 SNES Function norm 
7.071067809722e+001 <br>  0 KSP preconditioned resid norm 
7.071067809722e+001 true resid norm 7.071067809722e+001 ||Ae||/||Ax|| 
1.000000000000e+000<br></div>  1 KSP preconditioned resid norm 
6.628127802131e+001 true resid norm 6.628127802131e+001 ||Ae||/||Ax|| 
9.373588233757e-001<br>  2 KSP preconditioned resid norm 
6.060349268020e+001 true resid norm 8.023787055846e+001 ||Ae||/||Ax|| 
1.134734847941e+000<br>  3 KSP preconditioned resid norm 
5.893436541308e+001 true resid norm 1.661766291099e+002 ||Ae||/||Ax|| 
2.350092427079e+000</font></div></blockquote></div><br><div>This is remarkably early loss of orthogonality. Does it do the same thing if you remove -snes_mf_operator?</div><div><br></div><div>./pgfem -snes_monitor -ksp_type fgmres -ksp_monitor_true_residual -ksp_gmres_modifiedgramschmidt -ksp_rtol 1e-10 -pc_type none</div>
<div><br></div><div>If so, add -ksp_view_binary and send the matrix (probably named &quot;binaryoutput&quot;) to <a href="mailto:petsc-maint@mcs.anl.gov">petsc-maint@mcs.anl.gov</a> .</div><div><br></div><div>./pgfem -snes_monitor -ksp_type fgmres -ksp_monitor_true_residual -ksp_gmres_modifiedgramschmidt -ksp_rtol 1e-10 -pc_type none -ksp_view_binary</div>