  To date the most precise estimations of the critical exponent for the
Anderson transition have been made using the transfer matrix method. This
method involves the simulation of extremely long quasi one-dimensional systems.
The method is inherently serial and is not well suited to modern massively
parallel supercomputers. The obvious alternative is to simulate a large
ensemble of hypercubic systems and average. While this permits taking full
advantage of both OpenMP and MPI on massively parallel supercomputers, a
straight forward implementation results in data that does not scale. We show
that this problem can be avoided by generating random sets of orthogonal
starting vectors with an appropriate stationary probability distribution. We
have applied this method to the Anderson transition in the three-dimensional
orthogonal universality class and been able to increase the largest $L\times L$
cross section simulated from $L=24$ (New J. Physics, 16, 015012 (2014)) to
$L=64$ here. This permits an estimation of the critical exponent with improved
precision and without the necessity of introducing an irrelevant scaling
variable. In addition, this approach is better suited to simulations with
correlated random potentials such as is needed in quantum Hall or cold atom
systems.
