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

Reply via email to