answered 2022-06-12 22:50:30 +0100
You can get much shorter expressions for the components of u_gradient by running
u_gradient
u_gradient.apply_map(factor)
which factorizes the components.