rho, rho_bar, sigma, u, u_bar, v, v_bar, w, w_bar, z, z_bar, sigma_bar, 
-2*u^2*v^2*z_bar^2-2*u^2*v_bar^2*w_bar^2+2*u_bar^2*v^2*w^2+2*u_bar^2*v_bar^2*z^2-rho*u_bar^2+rho_bar*u^2
-2*u^2*v^2*z_bar^2+2*u^2*v_bar^2*w_bar^2-2*u_bar^2*v^2*w^2+2*u_bar^2*v_bar^2*z^2-sigma*v_bar^2+v^2*sigma_bar
2*u^2*v^2*z_bar^2+2*u^2*v_bar^2*w_bar^2-2*u_bar^2*v^2*w^2-2*u_bar^2*v_bar^2*z^2+rho*u_bar^2-rho_bar*u^2
2*u^2*v^2*z_bar^2-2*u^2*v_bar^2*w_bar^2+2*u_bar^2*v^2*w^2-2*u_bar^2*v_bar^2*z^2+sigma*v_bar^2-v^2*sigma_bar
u*u_bar-1
v*v_bar-1
