This is the one of the requirements of advanced pole placement procedures that are being implemented. Also, there seems to be no implementation of such routine anywhere other than the excessive rank checks of
[B, A@B, A^2@B...]
which are still error prone because of the powers of A.
Computation of controllability indices based on the QR decomposition of B matrix and then Arnoldi iteration over the A^k @ B vectors.
[x] Basic functionality and good-weather tests
[x] Write a test suite sufficiently ill-conditioned
This is the one of the requirements of advanced pole placement procedures that are being implemented. Also, there seems to be no implementation of such routine anywhere other than the excessive rank checks of
which are still error prone because of the powers of
A.Computation of controllability indices based on the QR decomposition of B matrix and then Arnoldi iteration over the
A^k @ Bvectors.