<html><body><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; color: #000000"><div>Dear all,<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>I am working on a code that solves a nonlinear system of equations G(x)=0 with Gauss-Seidel method. I managed to parallelize it by using DMDA with very good results. The previous week I changed my Gauss-Seidel solver by SNES. The code using SNES gives the same result as before, but I do not obtain the performance that I expected:<br data-mce-bogus="1"></div><div>1) When using the Gauss-Seidel method (-snes_type ngs) the residual G(x) seems not be scallable to the amplitude of x and I have to add the option -snes_secant_h in order to make SNES converge. However, I varied the step from 1.E-1 to 1.E50 and obtained the same result within the same computation time. Is it normal that <!--StartFragment-->snes_secant_h<!--EndFragment--> can vary so many orders of magnitude?<br data-mce-bogus="1"></div><div>2) Compared to my Gauss-Seidel algorithm, SNES does (approximately) the same number of iterations (with the same convergence criterium) but it is about 100 times slower. What can be the reason(s) of this slow performance of SNES solver? <!--StartFragment-->I do not use preconditioner with my algorithm so I did not add one to SNES.<!--EndFragment--> </div><div><br data-mce-bogus="1"></div><div>The main PETSc subroutines that I have included (in this order) are the following:<br data-mce-bogus="1"></div><div>call DMDACreate3D<br data-mce-bogus="1"></div><div>call DMSetUp</div><div>call DMCreateLocalVector<br data-mce-bogus="1"></div><div>call DMCreateGlobalVector<br data-mce-bogus="1"></div><div>call SNESCreate<br data-mce-bogus="1"></div><div>call SNESSetConvergenceTest<br data-mce-bogus="1"></div><div>call SNESSetDM<br data-mce-bogus="1"></div><div>call DMDASNESSetFunctionLocal<br data-mce-bogus="1"></div><div>call SNESSetFromOptions<br data-mce-bogus="1"></div><div>call SNESSolve<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Thanks in advance for you help.<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Best regards,<br data-mce-bogus="1"></div><div>Pedro<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div> <!--StartFragment--><span class="Object" role="link" id="OBJ_PREFIX_DWT363_ZmEmailObjectHandler"><span class="Object" role="link" id="OBJ_PREFIX_DWT369_ZmEmailObjectHandler"><br></span></span><!--EndFragment--> </div></div></body></html>