Now that we have complex Cholesky decomposition, it should be possible to
use cholesky_solve() in lu_solve() even for complex overdetermined systems. 
However it does not seem to work if you remove the special case logic for
such systems.