Iterative Krylov solution methods for geophysical electromagnetic simulations on throughput-oriented processing units