Fast frequency-domain electromagnetic modeling of a 3-D earth using finite differences