>>> a = torch.randn(5, 5, dtype=torch.cdouble, requires_grad=True)
>>> a = (a + a.transpose(-2, -1).conj()).div_(2)
>>> e, v = torch.symeig(a, eigenvectors=True)
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
RuntimeError: symeig does not support automatic differentiation for outputs with complex dtype.
cc @ezyang @anjali411 @dylanbespalko @mruberry @aocsa