Fix broadcasting for vector field source

[mfschubert]

Previously, simulations involving anisotropic eigensolve and batch dimensions (e.g. multiple wavelengths) would have failed. This fixes the issue.

I also modified the sources module, allowing the amplitudes_for_fields function to be used without Brillouin zone integration.