I've aplied Pierre's patch with minor modifications. See http://projects.scipy.org/scipy/numpy/changeset/2267 . For the time being trace is defined simply as diagonal().sum(). Thanks Pierre!