Extend polarization benchmark

[redeboer]

Some ideas:

• Full expression, with all angles substituted (only arguments are $\sigma1, \sigma2$)
• Full expression, but now with simplfication of the trigonomatric functions
  1. SymPy's simplification of WIgner-D functions (already there)
  2. Wigner-D of a acos (simplfication of trigonometric functions)
  3. ~Rewrite numerator of cos(zeta) expressions as Kibble function~
• Lambdification with and without use_cse=True
• JIT versus not-JITted?

• 4