Number of unstable roots in DR1

For modification, we need number of unstable roots in DR1 after calls to qzdiv.

where/how can we get that information?

After looking around a bit, I guess this may be achieved using

nnz(abs(eigenvalues_) > options_.qz_criterium)

would that be the best way?

Yes, that is the idea, but don’t use eigenvalues_, it is deprecated and returned only by the “check” command.

Use dr.eigval instead

Kind regards

Michel

Thanks. Is then nba variable in DR1 what we should use with higher resilience since nba may have different value if mjdgges.dll is called instead the qzdiv?

Normally, it shouldn’t matter

Best

michel