<div dir="ltr"><div><br></div><div>So Baky and I have been at the Brookhaven GPU Hackathon now for three days, talking to everyone. We have also been emailing with people who will respond to us from the hypre team and the PETSc team, as well as reading every blog post and mail archive and message board and from what we can tell, a distributed AMG preconditioner will not be available for us on a Summit platform for the foreseeable future.</div><div><br></div><div>There is a hypre build for CUDA, but it has a problem with it's use of CUSP, and nobody seems to be working on it.</div><div><br></div><div>PETSc has some .cu cuda files for the SpMV and Vector operations but the preconditioners are limited to point Jacobi and similar simple operations and a version of ILU. Neither works for our stiff projection in the embedded boundary algorithms. We built it and ran it and PETSc takes several hundred iterations to get the residual down by a factor of 6. We need to get down to more like 10e-11 for this solve.</div><div><br></div><div>The AMG being worked on by the NVIDIA team is not targeted for multi-node solving, and I haven't heard back from them in months.<br></div><div><br></div><div> We are left with two options as I see it to meet our ECP Milestones:</div><div><br></div><div>1. Build yet another interface, this time to see if there is a distributed GPU AMG preconditioner in Trilinos<br></div><div><br></div><div> 2. Implement our own special-purpose EB-GMG solver written in Chombo.<br></div><div><br></div><div>I would love to be wrong about all this.</div><div><br></div><div>Brian</div><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><span style="font-family:monospace,monospace">Brian Van Straalen Lawrence Berkeley Lab<br><a href="mailto:BVStraalen@lbl.gov" target="_blank">BVStraalen@lbl.gov</a> Computational Research<br>(510) 486-4976 Division (<a href="http://crd.lbl.gov" target="_blank">crd.lbl.gov</a>)</span><br></div></div>
</div>