<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Apr 20, 2017 at 5:03 AM, Damian Kaliszan <span dir="ltr"><<a href="mailto:damian@man.poznan.pl" target="_blank">damian@man.poznan.pl</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
Currently  I'm  doing a simple Ax=b with KSP method and MPI by loading<br>
in parallel a data from file<br>
<br>
viewer = petsc4py.PETSc.Viewer().<wbr>createBinary(Afilename, 'r')<br>
A = petsc4py.PETSc.Mat().load(<wbr>viewer).<br>
<br>
(same for b & x vectors) and then calling a solver<br>
<br>
ksp.solve(b, x)<br>
<br>
I would like to do the same using GPU.<br>
I  changed respectively the above to<br>
viewer = petsc4py.PETSc.Viewer().<wbr>createBinary(Afilename, 'r')<br>
A = PETSc.Mat().create(comm=comm)<br>
A.setType(PETSc.Mat.Type.<wbr>MPIAIJCUSPARSE)<br>
A.load(viewer)<br>
...<br>
What else needs to be changed?<br>
Running  the  above and checking nvidia-smi output confirms the python<br>
script and computations are done on CPU, not GPU as I 'd like to...<br></blockquote><div><br></div><div>1) Make sure the Vec types are also for the GPU</div><div><br></div><div>2) Send the output of</div><div><br></div><div>  -ksp_view -log_view</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Any help would be appreciated.<br>
<br>
Best,<br>
Damian<br>
<br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">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></div>