To speed up solution to quasistatic approximation to Maxwell's equations implement as a staggered grid finite difference solution in the frequency domain. This will likely require implementing a domain decomposition solution.