<div dir="ltr"><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Feb 9, 2021 at 4:31 AM Chen Gang <<a href="mailto:569615491@qq.com">569615491@qq.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>Dear Professor,</div><div><br></div><div>I have a question about the bjacobi+ilu, I'm in the place  </div><div><br></div><div>src/ksp/ksp/tutorials </div><div><br></div><div>and running the ex2 with</div><div><br></div><div><div>mpirun -np 40 -bind-to core ./ex2 -m 2000 -n 2000  -ksp_max_it 1000000 -ksp_rtol 1e-10 -ksp_type gmres -pc_type bjacobi -sub_pc_type ilu  -sub_pc_factor_levels 10 -ksp_monitor_short</div></div><div><br></div><div>on the Intel machine with 40 cores, the petsc is complied with MKL and ICC/MPIICC with O3.</div><div><br></div><div>I found the numerical results are not the same from time to time, although they all right, the iterations are not exactly the same from time to time. Where have I done wrong?</div></blockquote><div><br></div><div>Humm, iterations should rarely change. You can see floating point scale differences. Try jacobi instead of ilu and/or try not using MKL.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div>I'm looking forward to hear you!</div><div><br></div><div>Best,</div><div><br></div><div>Gang Chen</div><div><br></div><div>Sichuan Univertsity,</div><div>School of Mathematics,</div><div>Chengdu, China</div></blockquote></div></div>