Improved convergence in the volume-integral method (VIM) of calculating electromagnetic scattering properties of a structure is achieved by numerically solving a volume integral equation for a vector field, F, rather than the electric field, E. The electric field, E, is determined from the vector field, F, after solving of the volume integral equation. The vector field, F, may be related to the electric field, E, by a change of basis, and may be continuous at material boundaries where the electric field, E, has discontinuities. Convolutions of the vector field, F, are performed using convolution operators according to the finite Laurent rule, which allows for efficient matrix-vector products using Fast Fourier Transforms. An invertible convolution-and-change-of-basis operator, C, is configured to transform the vector field, F, to the electric field, E, by performing a change of basis according to material and geometric properties of the periodic structure.