
<html>


<head>


<meta http-equiv="Content-Type" content="text/html; charset=utf-8">


</head>


<body>


<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">


<div class="PlainText">Hi,<br>


<br>


I am getting close to finish porting a code from petsc 3.3 / sieve to main / dmplex, but am now encountering difficulties


<br>


I am reasonably sure that the Jacobian and residual are correct. The codes handle boundary conditions differently (MatZeroRowCols vs dmplex constraints) so it is not trivial to compare them. Running with snes_type ksponly pc_type Jacobi or hyper gives me the


 same results in roughly the same number of iterations.<br>


<br>


In my old code, gamg would work out of the box. When using petsc-main, -pc_type gamg -pc_gamg_type agg works for _some_ problems using P1-Lagrange elements, but never for P2-Lagrange. The typical error message is in gamg_agg.txt<br>


<br>


When using -pc_type classical, a problem where the KSP would converge in 47 iteration in 3.3 now takes 1400.  ksp_view_3.3.txt and ksp_view_main.txt show the output of -ksp_view for both versions. I don’t notice anything obvious.<br>


<br>


Strangely, removing the call to PCSetCoordinates does not have any impact on the convergence.<br>


<br>


I am sure that I am missing something, or not passing the right options. What’s a good starting point for 3D elasticity?<br>


Regards,<br>


Blaise<br>


<br>


<br>


</div>


</span></font></div>


<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">


<div class="PlainText"><br>


— <br>


Canada Research Chair in Mathematical and Computational Aspects of Solid Mechanics (Tier 1)<br>


Professor, Department of Mathematics & Statistics<br>


Hamilton Hall room 409A, McMaster University<br>


1280 Main Street West, Hamilton, Ontario L8S 4K1, Canada <br>


<a href="https://www.math.mcmaster.ca/bourdin">https://www.math.mcmaster.ca/bourdin</a> | +1 (905) 525 9140 ext. 27243<br>


<br>


</div>


</span></font></div>


</body>


</html>