
<html>


<head>


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


</head>


<body>


<p><font size="4"><font face="monospace">Dear all,</font></font></p>


<p><font size="4"><font face="monospace">I don't know if this is a bug, but I observed that when using GMRES with AGG-PCGAMG as preconditioner Dirichlet boundary conditions don't seem to be exactly fulfilled.</font></font></p>


<p><font size="4"><font face="monospace">My Matrix has zero rows and cols with 1 on the diagonal where I have dirichlet-bcs in my FE-mesh and I would expect that the eqs in this rows can be exactly fulfilled (as u_i = g_i) there.</font></font></p>


<p><font size="4"><font face="monospace">However, when I solve A*x = b with the above solver I only get u_i = g_i + error in that part of the vector. Switching from pc_gamg_type agg to pc_gamg_type classical cures this problem, but the classical is not advertised


 in the user manual.</font></font></p>


<p>These are the options I'm currently using:</p>


<p>-ksp_type gmres<br>


-ksp_pc_side right<br>


-pc_type gamg<br>


-pc_gamg_type agg [or classical]<br>


-pc_gamg_sym_graph 1<br>


-pc_gamg_square_graph 1<br>


-pc_gamg_agg_nsmooths 1<br>


-pc_gamg_threshold 0.01<br>


-pc_mg_cycles v<br>


</p>


<p>Iteration counts are basically the same.</p>


<p>Best regards</p>


<p>Elias<br>


</p>


<pre class="moz-signature" cols="72">-- 


Dr. Elias Karabelas


Research Associate


University of Graz


Institute of Mathematics and Scientific Computing


Heinrichstraße 36


A-8010 Graz


Austria


Phone: +43 316 380 8546


Email: <a class="moz-txt-link-abbreviated" href="mailto:elias.karabelas@uni-graz.at">elias.karabelas@uni-graz.at</a>


Web:  <a class="moz-txt-link-freetext" href="https://ccl.medunigraz.at/">https://ccl.medunigraz.at/</a></pre>


</body>


</html>