<div dir="ltr"><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Sep 18, 2022 at 6:19 PM Barry Smith <<a href="mailto:bsmith@petsc.dev">bsmith@petsc.dev</a>> wrote:<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 style="overflow-wrap: break-word;"><div><br></div>  Mark,<div><br></div><div>    Some how your changes to GAMG in June slowed down the time to compute the MIS dramatically. I cannot figure out what options to use to </div><div>get the exact same performance as an older branch. -mat_coarsen_type mis -pc_gamg_threshold 0 result in longer times than the older code with its default options.</div></div></blockquote><div><br></div><div>The new MIS-2 folds in the square graph with the MIS. Before the square graph was in a separate method that created an squared graph explicitly.  So don't use aggressive coarsening (you will see a PtAP if you use aggressive coarsening in the new code)</div><div>And -pc_gamg_threshold 0 will filter (zeros only). Use < 0 for no filtering.</div><div>The old code also had this optimization to not create a graph for bs==1 and no filter,</div><div><br></div><div>Mark</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 style="overflow-wrap: break-word;"><div><br></div><div><br><div><br><blockquote type="cite"><div>On Sep 18, 2022, at 4:21 PM, Mark Adams <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>> wrote:</div><br><div><div dir="ltr"><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Sep 18, 2022 at 4:02 PM Barry Smith <<a href="mailto:bsmith@petsc.dev" target="_blank">bsmith@petsc.dev</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
  Mark,<br>
<br>
   Do all MIS algorithms  in PETSc require a symmetric graph structure? And parallel ones can hang if not structurally symmetric?<br></blockquote><div><br></div><div>Yes,</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">
<br>
   When used sequentially I guess it never hangs but it may not produce a "correct" MIS if the matrix structure is not symmetric?</blockquote><div><br></div><div>It is fine in serial and it is not necessarily an MIS of the symmetrized graph.</div><div>If there is a one way edge between two vertices and the order of the greedy MIS process picks the root of the edge it is an MIS of the symmetrized graph, otherwise both vertices could get selected.</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"> But like the MIS is fine for GAMG in this circumstance?<br></blockquote><div><br></div><div>It will be fine for GAMG. The MIS is just a heuristic.</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">
<br>
  Barry<br>
<br>
</blockquote></div></div>
</div></blockquote></div><br></div></div></blockquote></div></div>