<div dir="ltr"><div dir="ltr">On Mon, Dec 26, 2022 at 10:29 AM Edoardo Centofanti <<a href="mailto:edoardo.centofanti01@universitadipavia.it">edoardo.centofanti01@universitadipavia.it</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Thank you for your answer. Can you provide me the full path of the example you have in mind? The one I found does not seem to exploit the algebraic multigrid, but just the geometric one.</div></blockquote><div><br></div><div>cd $PETSC_DIR/src/snes/tutorials/ex5</div><div>./ex5 -da_grid_x 64 -da_grid_y 64 -mms 3 -pc_type gang</div><div><br></div><div>and for GPUs I think you need the options to move things over</div><div><br></div><div>  -dm_vec_type cuda -dm_mat_type aijcusparse</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</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 dir="ltr"><div>Thanks,</div><div>Edoardo</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno lun 26 dic 2022 alle ore 15:39 Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> ha scritto:<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 dir="ltr"><div dir="ltr">On Mon, Dec 26, 2022 at 4:41 AM Edoardo Centofanti <<a href="mailto:edoardo.centofanti01@universitadipavia.it" target="_blank">edoardo.centofanti01@universitadipavia.it</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi PETSc Users,<div><br></div><div>I am experimenting some issues with the GAMG precondtioner when used with GPU.</div><div>In particular, it seems to go out of memory very easily (around 5000 dofs are enough to make it throw the "[0]PETSC ERROR: cuda error 2 (cudaErrorMemoryAllocation) : out of memory" error). </div><div>I have these issues both with single and multiple GPUs (on the same or on different nodes). The exact same problems work like a charm with HYPRE BoomerAMG on GPUs.</div><div>With both preconditioners I exploit the device acceleration by giving the usual command line options "-dm_vec_type cuda" and "-dm_mat_type aijcusparse" (I am working with structured meshes). My PETSc version is 3.17.</div><div><br></div><div>Is this a known issue of the GAMG preconditioner? </div></div></blockquote><div><br></div><div>No. Can you get it to do this with a PETSc example? Say SNES ex5?</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</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 dir="ltr"><div>Thank you in advance,</div><div>Edoardo </div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>
</blockquote></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>