<div dir="ltr"><div><div><div>Hi,<br><br></div>I run the code src/snes/examples/tutorials/ex5.c with options: <br><br><i> mpirun -n 8 ./ex5 -pc_type mg \<br> -ksp_monitor\<br> -pc_mg_levels 3 \<br> -pc_mg_galerkin \<br> -da_grid_x 17 \<br>
-da_grid_y 17 \<br> -mg_levels_ksp_norm_type unpreconditioned \<br> -snes_monitor \<br> -mg_levels_ksp_chebyshev_estimate_eigenvalues 0.5,1.1 \<br> -mg_levels_pc_type sor \<br> -pc_mg_type \</i><br><br><br></div>I run this script several times, but at each time, the residual norms have some tiny differences ( these should not happen). For example:<br>
<br></div>Case 1:<br><br><i> 0 SNES Function norm 1.188788066192e+00 <br> 0 KSP Residual norm 1.573384253521e+00 <br> 1 KSP Residual norm 3.616321708396e-02 <br> 2 KSP Residual norm 2.780221563755e-04 <br> 3 KSP Residual norm 2.194662354037e-06 <br>
1 SNES Function norm 5.125240595190e-03 <br> 0 KSP Residual norm 1.217284338266e-01 <br> 1 KSP Residual norm 1.774247017346e-04 <br> 2 KSP Residual norm 2.557118591292e-06 <br> 3 KSP Residual norm 1.622173269367e-08 <br>
2 SNES Function norm 3.922335995111e-05 <br> 0 KSP Residual norm 9.960745924876e-04 <br> 1 KSP Residual norm 1.273916336665e-06 <br> 2 KSP Residual norm 1.571259383270e-08 <br> 3 KSP Residual norm 1.250266145356e-10 <br>
3 SNES Function norm 2.662898279023e-09 </i><br><div><br><br></div><div>Case 2: <br><br><br><i> 0 SNES Function norm 1.188788066192e+00 <br> 0 KSP Residual norm 1.573384253521e+00 <br> 1 KSP Residual norm 3.616321708396e-02 <br>
2 KSP Residual norm 2.780221563755e-04 <br> 3 KSP Residual norm 2.194662354037e-06 <br> 1 SNES Function norm 5.125240595190e-03 <br> 0 KSP Residual norm 1.217284338266e-01 <br> 1 KSP Residual norm 1.774247017347e-04 <br>
2 KSP Residual norm 2.557118591292e-06 <br> 3 KSP Residual norm 1.622173269367e-08 <br> 2 SNES Function norm 3.922335995<span style="background-color:rgb(255,0,0)"><span style></span></span><span style="background-color:rgb(255,0,0)">108</span>e-05 <br>
0 KSP Residual norm 9.960745924<span style="background-color:rgb(255,0,0)">862</span>e-04 <br> 1 KSP Residual norm 1.273916336<span style="background-color:rgb(255,0,0)">6<span style="background-color:rgb(255,0,0)"><span style></span></span>54</span>e-06 <br>
2 KSP Residual norm 1.571259383<span style="background-color:rgb(255,0,0)">257</span>e-08 <br> 3 KSP Residual norm 1.250266145<span style="background-color:rgb(255,0,0)">346</span>e-10 <br> 3 SNES Function norm 2.6628982<span style="background-color:rgb(255,0,0)">85038</span>e-09 </i><br>
<br></div><div>These differences are marked by red color. So here, I was wondering there are any explanations why these happen?<br><br></div><div>Thanks,<br><br></div><div>Fande,<br></div><div><br></div><div><br></div></div>