Hi all,
I understand 'Data types not conformable for operation' but could it be
that if A was a real matrix then Re(A) = A and Im(A) = zeros(.,.) ?
I do not know how difficult that might be.
At the moment it is not possible.
The issue arose as I was experimenting with 'eigensym' and 'eiggen2'.
The latter seems to be way faster than the eigensym for symmetric
matrices (sorting eigenvalues considered). But this is another issue I
will raise with a code example when it's ready.
Yiannis