Nanophotonic structures are capable of generating field hotspots, which can enhance quantum light-matter interactions by many orders of magnitude. However, numerical simulations for applications such as radiative heat transfer, electron energy loss spectroscopy, van der Waals forces, Purcell factor throughout a volume, and many others are challenging and often computationally prohibitive. Common to these simulations is that the Green’s function or local photonic density of states must be known at each point across a volume of space, necessitating the solution of Maxwell’s equations perhaps many thousands of times.