[petsc-users] Simple question regarding dense matrix algebra using elemental

Preyas Shah shah.preyas at gmail.com
Sun Apr 12 18:57:10 CDT 2015


Hi,

I have been recently investigating the use of Petsc for solving a PDE
related to my research and web search suggests that I should use petsc with
elemental.

So far, I was required to solve a matrix equation Ax=b where A was dense
(with number of non zeros =0) but the size of the matrix was order
5000x5000. I employed the standard serial LU solver from Gnu Scientific
Library and obtained a decent runtime that served my needs.

Now I am investigating the same problem in a particularly singular limit of
one parameters in my PDE. As a result, to obtain grid convergence on the
physical domain, I am forced to go to sizes of A beyond 30000x30000. I am
trying to find a good library that can solve such dense systems in
**parallel**. Petsc says its capable of doing dense linear algebra but my
web search hasn't shown me any examples where dense equations are solved in
**parallel**. A webpage showing a minimum working example would be enough.
Or any other advice :)

Thanks for your time!
​~Preyas
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150412/1c8ce005/attachment.html>


More information about the petsc-users mailing list