New to the CUDA forums and unfortunately, never really had a chance to play with it. However, someone from our team in Shanghai has worked with CUDA for a while and come up with something interesting. I don’t want to spam the forums with a product announcement but I would like to do a quick survey of what is currently available to make sure we are doing something really new.
Our main application is semiconductor physics and we do a lot of FEM modeling of very non-linear equations. Because of that, iterative solvers have never been sufficiently powerful to converge reliably with our badly conditioned asymmetric matrices and we rely on direct solvers. In order to parallelize as much as possible, we had our own multi-frontal solver and recently implemented the well-known MUMPS (http://en.wikipedia.org/wiki/MUMPS) solver.
When we started considering using a CUDA version of a multi-frontal solver a few months ago, it seemed like none was available. To the best of my knowledge, only direct solvers for full/banded matrices (CULA LAPACK) and iterative sparse solvers (http://www.acceleware.com/) exist. I know that ANSYS has also ported its own mechanical FEM solver but I do not know if they are using a direct solver or an iterative one.
So what can the experienced CUDA developers out there tell me ? Has anyone else ported a multi-frontal direct solver to CUDA ? I’d like to believe our developer’s claims that he is the first but a little due diligence never hurt anyone.
If anyone is interested, I can release some benchmark comparisons to MUMPS.