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

[Antonio Stanziola]

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.