The dispatch mechanism is totally different from MatLUFactorSymbolic/Numeric, with no ability to choose the factorization package. Should it be changed to take a MatSolverPackage and then to dispatch using a MatGetFactor?