
<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>

</head>

<body dir="ltr">

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

<p>Hello,</p>

<p><br>

</p>

<p>Our team at Oak Ridge National Lab requires a distributed and GPU-enabled (ideally) iterative solver as part of a new, high-dimensional PDE solver. We are exploring options for software packages with this capability vs. rolling our own, i.e. having some

 of our team members write one.</p>

<p><br>

</p>

<p>Our code already has a distributed, GPU-enabled matrix-vector multiply that we'd like to use for the core of GMRES or similar technique.  I've looked through the PETSc API and found that matrix-free methods are supported, and this: https://www.mcs.anl.gov/petsc/features/gpus.html

 seems to indicate that GPU acceleration is available for iterative solvers. <br>

</p>

<p><br>

</p>

<p>My question is: does PETSc support all of these things together? E.g. is it possible for me to use a distributed, matrix free iterative solver with a preconditioner shell on the GPU with PETSc?</p>

<p><br>

</p>

<p>Best,</p>

<p>Tyler McDaniel</p>

<br>

</div>

</body>

</html>