Raise error when helmholtz_solver is called with a real field

[astanziola]

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.