Rewrite complex coefficients that cannot be represented using single precision

Loopy stops with an exception in this case, asking the user to explicitly cast to np.complex64/128

This is needed for 3D Helmholtz expansions of order >= 10