  Observations, mainly of outbursts in dwarf novae, imply that the anomalous
viscosity in highly ionized accretion discs is magnetic in origin, and requires
that the plasma $\beta \sim 1$. Until now most simulations of the magnetic
dynamo in accretion discs have used a local approximation (known as the
shearing box). While these simulations demonstrate the possibility of a
self-sustaining dynamo, the magnetic activity generated in these models
saturates at $\beta \gg 1$. This long-standing discrepancy has previously been
attributed to the local approximation itself. There have been recent attempts
at simulating magnetic activity in global accretion discs with parameters
relevant to the dwarf novae. These too find values of $\beta \gg 1$. We
speculate that the tension between these models and the observations may be
caused by numerical magnetic diffusivity. As a pedagogical example, we present
exact time-dependent solutions for the evolution of weak magnetic fields in an
incompressible fluid subject to linear shear and magnetic diffusivity. We find
that the maximum factor by which the initial magnetic energy can be increased
depends on the magnetic Reynolds number as ${\mathcal R}_{\rm m}^{2/3}$. We
estimate that current global numerical simulations of dwarf nova discs have
numerical magnetic Reynolds numbers around 6 orders of magnitude less than the
physical value found in dwarf nova discs of ${\mathcal R}_{\rm m} \sim
10^{10}$. We suggest that, given the current limitations on computing power,
expecting to be able to compute realistic dynamo action in observable accretion
discs using numerical MHD is, for the time being, a step too far.
