Exact equal time statistics of Orszag–McLaughlin dynamics investigated using the Hopf characteristic functional approach