<div dir="ltr">I think the first thing that I need to do is get ASM working on the CPU in my code as a direct solver.<div><br></div><div>FieldSplit was mentioned the way to do this. How would I do that?</div><div><br></div><div>Thanks,</div><div>Mark</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Dec 30, 2020 at 11:11 PM Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</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">Barry Smith <<a href="mailto:bsmith@petsc.dev" target="_blank">bsmith@petsc.dev</a>> writes:<br>
<br>
>   I was assuming the triangular systems were all solvable and setup in the correct format so the triangular solver would not be the place to first detect problems with the factor. Though not completely generally that seems to cover most cases.<br>
<br>
You could write a solver that way, but that doesn't appear to be what cusparse has done. Perhaps it's because the function takes the raw arrays as inputs (rather than an opaque structure that can only be created by an analysis phase that confirms validity). For all we know, they might do 90% of the work synchronously and just leave the last kernel in the stream. <br>
</blockquote></div>