A numerical solution to full-vector electromagnetic scattering by three-dimensional nonlinear bounded dielectrics