[petsc-users] PETSc for distributed computing?

[Hong]

Georgios Samaras :

> 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.
>

As Barry mentioned, you can either use elemental directly, or use
petsc-elemental interface. An example can be found at
~petsc/src/mat/examples/tests/ex145.c

You may use petsc KSP interface instead. I just modified
~petsc/src/ksp/ksp/examples/tutorials/ex2.c

https://bitbucket.org/petsc/petsc/commits/1cc583424d1965a1a245571240baab92605aba58

so this example can be run with elemental with runtime options
$ mpiexec -n 3 ./ex2 -pc_type lu -pc_factor_mat_solver_package elemental
-mat_type elemental
Norm of error 2.81086e-15 iterations 1

Suggest using petsc-dev (master branch) for petsc-elemental interface.

Hong


> Thanks,
> Georgios Samaras
>
> On Thu, Apr 9, 2015 at 4:31 PM, Patrick Sanan <patrick.sanan at gmail.com>
> wrote:
>
>> Do you want an explicit inverse, or just to solve linear systems
>> involving a dense matrix? If so, how many right hand sides?
>>
>> PETSc itself focuses on sparse linear algebra, but it includes an
>> interface to the dense linear algebra library Elemental.
>>
>> Il giorno 09/apr/2015, alle ore 03:06 PM, Georgios Samaras <
>> georgesamarasdit at gmail.com> ha scritto:
>>
>> > Dear all,
>> >
>> >   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?
>> >
>> > Thanks in advance,
>> > Georgios Samaras
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150410/1f42e43e/attachment.html>