<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">Georgios Samaras :<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><div>We (my professor (cc) and I) would like to solve a linear system involving a dense matrix, but we have not yet think of the number of right hand sides, but we are working on it.<br></div></div></div></blockquote><div> </div><div>As Barry mentioned, you can either use elemental directly, or use petsc-elemental interface. An example can be found at </div><div>~petsc/src/mat/examples/tests/ex145.c </div><div><br></div><div>You may use petsc KSP interface instead. I just modified ~petsc/src/ksp/ksp/examples/tutorials/ex2.c </div><div><br></div><div><a href="https://bitbucket.org/petsc/petsc/commits/1cc583424d1965a1a245571240baab92605aba58">https://bitbucket.org/petsc/petsc/commits/1cc583424d1965a1a245571240baab92605aba58</a><br></div><div><br></div><div>so this example can be run with elemental with runtime options</div><div><div>$ mpiexec -n 3 ./ex2 -pc_type lu -pc_factor_mat_solver_package elemental -mat_type elemental</div><div>Norm of error 2.81086e-15 iterations 1</div></div><div><br></div><div>Suggest using petsc-dev (master branch) for petsc-elemental interface.</div><div><br></div><div>Hong</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><div><br></div>Thanks,<br></div>Georgios Samaras<br></div><div class=""><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Apr 9, 2015 at 4:31 PM, Patrick Sanan <span dir="ltr"><<a href="mailto:patrick.sanan@gmail.com" target="_blank">patrick.sanan@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Do you want an explicit inverse, or just to solve linear systems involving a dense matrix? If so, how many right hand sides?<br>
<br>
PETSc itself focuses on sparse linear algebra, but it includes an interface to the dense linear algebra library Elemental.<br>
<div><div><br>
Il giorno 09/apr/2015, alle ore 03:06 PM, Georgios Samaras <<a href="mailto:georgesamarasdit@gmail.com" target="_blank">georgesamarasdit@gmail.com</a>> ha scritto:<br>
<br>
> Dear all,<br>
><br>
>   I am looking for a C/C++ linear algebra package for performing inversion of a large dense matrix in a distributed environment. Is your package suitable for this goal?<br>
><br>
> Thanks in advance,<br>
> Georgios Samaras<br>
<br>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div></div>