Hello everyone,
I'm trying to replicate nearly the following in Sage: mathematica.stackexchange.com/questions/252630/how-to-find-a-unitary-transform-between-two-matrices
The linked post is asking for U such that A=UB¯UT but I am looking for the following:
Given 2 lists of vectors say A={a1,...,an}, B= {b1,...,bn}, is there a way to easily solve for a unitary matrix, U, such that {a1,...,an} = {Ub1,...,Ubn}.
I understand I can do the following:
U = H.solve_left(W)
U /= U.norm()
U.is_unitary()
but it doesn't seem every efficient, as there may be other solutions that I cannot see using this method. I am new to Sage and any advice would be super helpful!
Thank you