| Re: [eigen] 64 bit support on umfpack |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
On 27.10.2014 18:52, Zhexuan Zhang wrote:
I am trying to solve large size sparse linear system using umfpack. In
As Sameer already responded, you can directly call the umfpack routines
(this also gives you more control on the settings).
UmfPackSupport.h, it only has implementations on "di" and "zi" mode, but
not "dl" and "zl" mode. This causes failure when I try to solve large size
problem. I am not a C++ expert. I tried to edit the .h file but no success.
Does anyone know how to force the "dl" mode? If anyone can send me a .h
file that has "dl" mode enabled, it is greatly appreciated.
Basically, what is required is to overload the umfpack_* methods in
UmfPackSupport.h, with SuiteSparse_long instead of int and l instead i
in the names of the called methods. To avoid code duplication, I
recommend using some macros similar to EIGEN_MKL_LU_PARTPIV in
src/LU/PartialPivLU_MKL.h
Then inside UmfPackLU all int which refer to indexes must be replaced by
Index.
Cheers,
Christoph
--
----------------------------------------------
Dipl.-Inf., Dipl.-Math. Christoph Hertzberg
Cartesium 0.049
Universität Bremen
Enrique-Schmidt-Straße 5
28359 Bremen
Tel: +49 (421) 218-64252
----------------------------------------------