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