Use complex-valued linear solvers even when operator is real-valued when the RHS is complex

[sebastiangrimberg]

In the limit which both systems are solved exactly, this has no effect. But it does affect the solution for some non-zero convergence tolerance.