  We show that the common local Monte Carlo rules used to simulate the motion
of driven flux lines in disordered media cannot capture the interplay between
elasticity and disorder which lies at the heart of these systems. We therefore
discuss a class of generalized Monte Carlo algorithms where an arbitrary number
of line elements may move at the same time. We prove that all these dynamical
rules have the same value of the critical force and possess phase spaces made
up of a single ergodic component. A variant Monte Carlo algorithm allows to
compute the critical force of a sample in a single pass through the system. We
establish dynamical scaling properties and obtain precise values for the
critical force, which is finite even for an unbounded distribution of the
disorder. Extensions to higher dimensions are outlined.
