We develop a computational framework for identifying bounds to light-matter interactions, originating from polarization-current-based formulations of local conservation laws embedded in Maxwell's equations. We propose an iterative method for imposing only the maximally violated constraints, enabling rapid convergence to global bounds. Our framework can identify bounds to the minimum size of any scatterer that encodes a specific linear operator, given only its material properties, as we demonstrate for the optical computation of a discrete Fourier transform. It further resolves bounds on far-field scattering properties over any arbitrary bandwidth, where previous bounds diverge.