Hello!

I would like to use Optim.jl with the target function defined on a complex Stiefel manifold (unitary matrices) of some large size.

Stiefel manifold is already implemented in Optim.jl. The problem is I can’t use autodiff due to complex numbers, and optimization performs badly without the gradient.

On the other hand unitary matrices can be implemented as arrays of real numbers (of double size), so autodiff will work. But in this case I will miss the unitary constrains (so I’ll have to reimplement them, and I’m not sure how).

Is there any simple workaround on this?