Calculate internal fields inside scatterers for HoloPy

For some applications (topological derivatives, visualizations, pedagogy) it may be useful to calculate the internal field inside a scatterer. Currently, trying to calculate a field inside a scatterer gives something unphysical.