solve_for didn't throw error in some nonlinear equation

[songjhaha]

Consider this simple equation:

using Symbolics
@variables x y
eqs = [x*y ~ 4, x+y~4]
Symbolics.solve_for(eqs, [x,y])

Which gives a result but solution is represented with variable y:

julia> Symbolics.solve_for(eqs, [x,y])
2-element Vector{SymbolicUtils.BasicSymbolic{Real}}:
 4 / y
 (4//1) + -4 / y

Symbolics.islinear gives false in this situation, while Symbolics.linear_expansion gives true:

julia> Symbolics.islinear(x*y, [x, y])
false

julia> Symbolics.linear_expansion([x*y~0], [x,y])
(Num[-y 0], Num[0], true)

Versions:
Julia: 1.7.3
Symbolics: v5.5.1