number of unstable roots in DR1

[perendia]

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

where/how can we get that information?

[perendia]

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

nnz(abs(eigenvalues_) > options_.qz_criterium)

would that be the best way?

[MichelJuillard]

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

[perendia]

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?

[MichelJuillard]

Normally, it shouldn't matter

Best

michel