Raise error when `helmholtz_solver` is called with a real field

**Describe the bug**
It does not make sense to call the Helmholtz operator with a real source field, as it is defined on complex fields. 

One option would be to convert the field to complex, but this could be confusing, so it is probably better to just raise an error when one attempts to do so.