Highly optimized simulations on single- and multi-GPU systems of 3D Ising spin glass