Use the C interface of hsl_ma57

[amontoison]

• It is required if we want to use the functions in hsl_ma57_patch.jl.

[goulart-paul]

Is there any update to this issue? Specifically, I would really like to use ma57_get_factors, but can't get it to work on any branch. It seems to want a ma57lfd_ symbol that doesn't appear in the HSL binaries.

Asking because the package is occasionally producing factorisations for me that appear to have a lot more nonzeros in the factors than I think should be there, but it's quite hard to debug without being able to query for the factors directly.

[amontoison]

ma57lfd_ was available with a patch for the Fortran code of ma57 developed by Dominique. By default we can't use these routines.
With the lastest releases of hsl_ma57, they added the features that we wanted but only in the C interface.

So, we need to rewrite our Julia interface for ma57 bacause we interfaced it before that they added the C interface.

[amontoison]

@goulart-paul
I found a fix for ma57_get_factors, we will be able to use this function with the next release of JuliaHSL / libHSL.
If you don't want to wait, send me an email with your platform (Base.BinaryPlatforms.host_triplet()).

[amontoison]

@goulart-paul
ma57_get_factors works with the HSL_jll.jl provided by libHSL.
I was released this week. 🎉