Closed wrs28 closed 2 years ago
Many, many thanks for this. I need a little time to review and test but will do asap!
Do take your time. I'm sure there is something that will break that I missed, and it's important to get it right.
@dpo have you had a chance to look at this? There’s probably some updating to do to make compatible with Julia 1, but it can probably wait till the initial port is done
@wrs28 Sorry for the long delay. I've been overloaded. I really do want to merge this but it'll have to wait until January.
Merged in #64.
Per Issue #46, here is a first pass at upgrading
MUMPS.jl
to using the MUMS C-API. See commit comments for details, but essentially this PR keeps both the original code and portsMUMPS3.jl
. For now there are twoinclude
lines inMUMPS.jl
which govern which set of codes is used. It is easy for a developer to toggle between the two for testing.The exported methods of this PR have the same signature as the original codes, and no new methods are exported. All tests pass with both codes, with one test being slightly modified as noted in the commit comment, as does the example in
examples/mumps_mpi.jl
, again, with suitable updating. This is so that anyone usingMUMPS.jl
will experience minimal changes.As it stands, sparse rhs is now supported. Other features, such as computing and accessing the Schur complement are implemented inside the code, but not exported.