  We present results of the application of the anisotropic hydrodynamics
(aHydro) framework to (2+1)-dimensional boost invariant systems. The necessary
aHydro dynamical equations are derived by taking moments of the Boltzmann
equation using a momentum-space anisotropic one-particle distribution function.
We present a derivation of the necessary equations and then proceed to
numerical solutions of the resulting partial differential equations using both
realistic smooth Glauber initial conditions and fluctuating Monte-Carlo Glauber
initial conditions. For this purpose we have developed two numerical
implementations: one which is based on straightforward integration of the
resulting partial differential equations supplemented by a two-dimensional
weighted Lax-Friedrichs smoothing in the case of fluctuating initial
conditions; and another that is based on the application of the Kurganov-Tadmor
central scheme. For our final results we compute the collective flow of the
matter via the lab-frame energy-momentum tensor eccentricity as a function of
the assumed shear viscosity to entropy ratio, proper time, and impact
parameter.
