Ask Your Question

Revision history [back]

click to hide/show revision 1
initial version

Solving for a Unitary Matrix

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\overline{U}^{T}$ but I am looking for the following:

Given 2 lists of vectors say $A = ${$a_1, ..., a_n$}, $B =$ {$b_1, ..., b_n$}, is there a way to easily solve for a unitary matrix, $U$, such that {$a_1, ..., a_n$} = {$Ub_1, ..., Ub_n$}.

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

Solving for a Unitary Matrix

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-matriceshttps://mathematica.stackexchange.com/q/252630

The linked post is asking for $U$ such that $A = UB\overline{U}^{T}$ but I am looking for the following:

Given 2 lists of vectors say $A = ${$a_1, ..., a_n$}, $B =$ {$b_1, ..., b_n$}, is there a way to easily solve for a unitary matrix, $U$, such that {$a_1, ..., a_n$} = {$Ub_1, ..., Ub_n$}.

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