  We present Cahn-Hilliard and Allen-Cahn numerical integration algorithms that
are unconditionally stable and so provide significantly faster
accuracy-controlled simulation. Our stability analysis is based on Eyre's
theorem and unconditional von Neumann stability analysis, both of which we
present. Numerical tests confirm the accuracy of the von Neumann approach,
which is straightforward and should be widely applicable in phase-field
modeling. We show that accuracy can be controlled with an unbounded time step
Delta-t that grows with time t as Delta-t ~ t^alpha. We develop a
classification scheme for the step exponent alpha and demonstrate that a class
of simple linear algorithms gives alpha=1/3. For this class the speed up
relative to a fixed time step grows with the linear size of the system as N/log
N, and we estimate conservatively that an 8192^2 lattice can be integrated 300
times faster than with the Euler method.
