  An algoritm for the simulation of the 3--dimensional random field Ising model
with a binary distribution of the random fields is presented. It uses
multi-spin coding and simulates 64 physically different systems simultaneously.
On one processor of a Cray YMP it reaches a speed of 184 Million spin updates
per second. For smaller field strength we present a version of the algorithm
that can perform 242 Million spin updates per second on the same machine.
