Optimization of objective function that's the product of unitary matrices

Click For Summary
The discussion centers on optimizing the objective function ||M-M_target||^2, where M is the product of over 100 2D unitary complex matrices and a vector. The user, Bastien, seeks a faster optimization method due to the structured nature of the problem, emphasizing that even a slight reduction in computation time is significant. A suggestion is made to consider using the Strassen algorithm for faster matrix multiplication. Additionally, it is noted that maximizing ||M-M_target|| directly could simplify the process, as squaring is monotonic for non-negative values. Overall, the conversation highlights the complexity of the optimization challenge while exploring potential algorithmic improvements.
bastpg
Messages
2
Reaction score
0
Hi,
I work in NRM and need for some reason to optimize an objective function of the form ||M-M_target||^2 where M is the product of a large number (>100) 2D unitary complex matrices (Qi) and a vector (A), i.e. M=Q1*Q2*...*QN*A, and M_target is a constant complex vector. I can do it directly, like people have done so far, but the problem has so much structure to it that it seems something smarter could be done. Note that this optimization should be fast, so reducing computation time even a little bit (>10%) could be a big deal. Let me know if the pb seems familiar and you have some ideas about it...
Thanks,
Bastien
 
Physics news on Phys.org
I do not see shortcuts, as your M_target and A have nothing in common. You could at best use a faster algorithm for matrix multiplication, in this case the Strassen algorithm should save some time.
 
What is constant and what is a variable here?
 
Only thing I can think of is, since squaring is monotonic for non-negatives, just seek to maximize ||M-M_target|| without concern for the square if you want to maximize or minimize difference if you wish to minimize the full expression.
 
I am studying the mathematical formalism behind non-commutative geometry approach to quantum gravity. I was reading about Hopf algebras and their Drinfeld twist with a specific example of the Moyal-Weyl twist defined as F=exp(-iλ/2θ^(μν)∂_μ⊗∂_ν) where λ is a constant parametar and θ antisymmetric constant tensor. {∂_μ} is the basis of the tangent vector space over the underlying spacetime Now, from my understanding the enveloping algebra which appears in the definition of the Hopf algebra...

Similar threads

  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 43 ·
2
Replies
43
Views
7K
  • · Replies 1 ·
Replies
1
Views
5K
  • · Replies 9 ·
Replies
9
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 11 ·
Replies
11
Views
3K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 6 ·
Replies
6
Views
2K
Replies
18
Views
5K