#2564 Regression in solve?

[kcrisman]

In Maxima 5.26:

(%i6) solve([x^2*y*z=18,x*y^3*z=24,x*y*z^4=6],[x,y,z]);
(%o6) [[x = 3, y = 2, z = 1], [x = 1.3372150673296066 - 2.6854898740652002 %i, y = 1.0528643257547123 %i - 1.7004342714592284, z = .9324722294043555 - .3612416661871516 %i],     
<snip>
[x = - .5512485534497119 %i - 2.9489192990517044, y = .18453671892660384 - 1.9914683525900694 %i, z = .8951632913550623 %i + .4457383557765383]]

In Maxima 5.29.1:

(%i5) solve([x^2*y*z=18,x*y^3*z=24,x*y*z^4=6],[x,y,z]);

Maxima encountered a Lisp error:

 In function <, the value of the first argument is
  ((MPLUS SIMP) -7.105427357601002e-15
   ((MTIMES SIMP) -1.0658141036401503e-14 $%I))
which is not of the expected type REAL

Automatically continuing.
To enable the Lisp debugger set *debugger-hook* to nil.

[Dan Gildea]

was able to reproduce in 5.29.1,
but works for me in current git