4 ms·

i'd be curious what matlab would do with finding the eigs for op's matrix. i'd bet a coffee that it (one of the eigs functions) prints a warning. entirely possi

by a-dub 5y ago

i'd be curious what matlab would do with finding the eigs for op's matrix. i'd bet a coffee that it (one of the eigs functions) prints a warning. entirely possible it may even switch search methods if precision loss is detected.

mattkrause 5y ago
I get the same results as the author if I pass the matrix to eig (R2020a, Intel Mac) >> lesp = gallery('lesp', 100) >> any(imag(eig(lesp))) ans = logical 0 >> any(imag(eig(lesp'))) ans = logical 1
- a-dub 5y ago
  what about eigs() with the 'Display' option turned on?
  - mattkrause 5y ago
    No luck there either, I'm afraid. It sets 'flag' to zero (i.e., converged), but A and A^T still differ.