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.