We solve the full vectorial Maxwell equations with the medium's nonlinear response evolving alongside the field.
Standard time-domain tools model nonlinearity as a perturbative polarization expansion on a linear host. That treatment holds at low conversion and slow envelopes.
It stops holding when the drive is strong, the pump depletes, the spectrum is broad, or the pulse is few-cycle. We solve the equations those models reduce from.