search

trilinos / Trilinos

Primary repository for the Trilinos Project
https://trilinos.org/
Other
1.19k stars 565 forks source link

Adding support for UMFPACK in Amesos2 #88

Closed dridzal closed 7 years ago

dridzal commented 8 years ago

@trilinos/amesos2: The @trilinos/rol team would like to use Amesos2 as the interface to a direct solver that is well suited for solving PDEs in two spatial dimensions. In our experience, UMFPACK works well for these problems. This capability would support the PDE-OPT application development kit and test suite.

srajama1 commented 8 years ago

The target date for this is mid february.

egboman commented 8 years ago

SuperLU, MUMPS and MKL-Pardiso are already supported. Do these not work well on PDEs?

dridzal commented 8 years ago

I've used SuperLU a few years ago, and found it (1) slightly slower than UMFPACK and (2) significantly less robust. The latter issue is why I prefer UMFPACK. Four or five years ago, I managed to compile Amesos with MUMPS, and it seemed to segfault on 50% of problems that I gave it, especially on KKT-type systems. It is possible that this has improved now. I have not tried Pardiso. I guess a question to @trilinos/amesos2 is if my observations are consistent with their knowledge of the solvers. I am open to exploring other options if adding UMFPACK is not worth it. Having said that, it seems a little strange that a well-recognized solver like UMFPACK is left out of the Amesos2 interface.

srajama1 commented 8 years ago

It is a fair request to have an interface for UMFPACK. We will add it.

jdbooth commented 8 years ago

I agree. It is a fair request. I have also experienced problems with robustness in SuperLU and even MUMPS (~4.5, I have not ran a large study with 5.0). Also MKL-Pardiso is hard to build if not using Intel compilers. UMFPACK would also be nice to have for comparisons.

srajama1 commented 8 years ago

Denis, I'll work on this next week. I am sorry for the delay. There were some unexpected work on other projects that blocked the work on this.

srajama1 commented 8 years ago

This has become higher priority, Moving to in progress.

srajama1 commented 7 years ago

PR #1069 fixed this. This is still not tested nightly though. See #1255 for details of nightly testing