<div dir="ltr">On Tue, Jul 16, 2013 at 1:13 PM, Harshad Sahasrabudhe <span dir="ltr"><<a href="mailto:hsahasra@purdue.edu" target="_blank">hsahasra@purdue.edu</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Jed,<br>
<br>
Thanks for your reply.<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
You're on your own for storage of factors.  Alternatively, you could add<br>
library support so that you could use PCLU and<br>
'-pc_factor_mat_solver_package magma' (or PCFactorSetMatSolverPackage).<br>
Doing this is not a priority for us, but we can provide guidance if you<br>
want to tackle it.<br>
</blockquote>
<br>
I would definitely like to start working on adding library support. I think this is the most efficient way to go about it. Can you give me certain details such as:<br>
<br>
1) How should I start going about it?<br></blockquote><div><br></div><div>Read the UMFPACK implementation</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

2) How will I check-in the changes to Petsc?<br></blockquote><div><br></div><div>Using Git</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
3) What version of Petsc will the changes be reflected in if I started working on it right now?<br></blockquote><div><br></div><div>A branch of 'master'</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

4) How many hours does it generally take to get this done?<br></blockquote><div><br></div><div>How many licks does it take to get to the center of a Tootsie Roll Pop?</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

5) How is the peer review done?<br></blockquote><div><br></div><div>Through a pull request on BitBucket.</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">

Thanks,<br>
Harshad<br>
<br>
On 07/13/2013 12:43 PM, Jed Brown wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
"<a href="mailto:hsahasra@purdue.edu" target="_blank">hsahasra@purdue.edu</a>" <<a href="mailto:hsahasra@purdue.edu" target="_blank">hsahasra@purdue.edu</a>> writes:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi,<br>
<br>
I am working on solving a system of linear equations with square<br>
matrix. I'm first factoring the matrix using LU decomposition.<br>
</blockquote>
I assume you're solving a dense problem because that is all MAGMA does.<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I want to do the LU decomposition step using MAGMA on GPUs. MAGMA<br>
library implements LAPACK functions on a CPU+GPU based system.<br>
<br>
So my question is, how do I extract the data from a Petsc Mat so that<br>
it can be sent to the dgetrf routine in MAGMA.<br>
</blockquote>
MatDenseGetArray<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Is there any need for duplicating the data for this step?<br>
</blockquote>
You're on your own for storage of factors.  Alternatively, you could add<br>
library support so that you could use PCLU and<br>
'-pc_factor_mat_solver_package magma' (or PCFactorSetMatSolverPackage).<br>
Doing this is not a priority for us, but we can provide guidance if you<br>
want to tackle it.<br>
</blockquote>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>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>