  We describe a numerical algorithm for computing spin glass ground states with
a high level of reliability. The method uses a population based search and
applies optimization on multiple scales. Benchmarks are given leading to
estimates of the performance on large lattices.
