Some eigenpairs from a large, sparse, nonsymmetric matrix: Julia vs Matlab

Specific Domains Numerics
, ,
21

With SuiteSparse merged into SparseArrays, UMFPACK.umf_ctrl doesn’t appeared defined? How would one do this now?

22

Refinement is now (Julia v1.9) omitted by default in sparse LU; see this issue for discussion.

It can be enabled by setting element JL_UMFPACK_IRSTEP in the default table returned by get_umfpack_control() (both are in SparseArrays.UMFPACK). The old value was 2.0.

2 Likes
23

This is great information, but we need to be very careful about looking at MATLAB source code. It is not free (libre) to copy from.

I would advise anyone trying to implement Krylov-Schur methods not to look at eigs.m.

1 Like