A bit different results for almost the same conditions but recomposed input

I have some expression to solve. But if try to rewrite or recompose the expression, then the result a more or less different.

Here is the source, where solution is used as a prove for some statement:
https://github.com/andry81/bittools/blob/trunk/src/bitsync/correlation.cpp

##Input 1:

(a + c) * c + (b + d) * d <= max((a + c) ^ 2 + d ^ 2, (b + d) ^ 2 + c ^ 2), a > 0, b > 0, c >= 1, d >= 1

Solve

##Solution:

{a c + b d + c^2 + d^2<=1/2 a^2 sgn(-c^2 + (a + c)^2 + d^2 - (b + d)^2) + a^2/2 + 1/2 b^2 sgn(c^2 - (a + c)^2 - d^2 + (b + d)^2) + a c sgn(-c^2 + (a + c)^2 + d^2 - (b + d)^2) + 1/2 c^2 sgn(-c^2 + (a + c)^2 + d^2 - (b + d)^2) + 1/2 d^2 sgn(-c^2 + (a + c)^2 + d^2 - (b + d)^2) + 1/2 c^2 sgn(c^2 - (a + c)^2 - d^2 + (b + d)^2) + 1/2 d^2 sgn(c^2 - (a + c)^2 - d^2 + (b + d)^2) + b d sgn(c^2 - (a + c)^2 - d^2 + (b + d)^2) + a c + b^2/2 + b d + c^2 + d^2, a>0, b>0, c>=1, d>=1}

Input 2:

(x + u) * u + (y + v) * v <= max((x + u) ^ 2 + v ^ 2, (y + v) ^ 2 + u ^ 2), x > 0, y > 0, u >= 1, v >= 1

Solve

Solution:

{u^2 + u x + v^2 + v y<=1/2 x^2 sgn(-u^2 + v^2 + (u + x)^2 - (v + y)^2) + 1/2 y^2 sgn(u^2 - v^2 - (u + x)^2 + (v + y)^2) + 1/2 u^2 sgn(-u^2 + v^2 + (u + x)^2 - (v + y)^2) + 1/2 u^2 sgn(u^2 - v^2 - (u + x)^2 + (v + y)^2) + u x sgn(-u^2 + v^2 + (u + x)^2 - (v + y)^2) + 1/2 v^2 sgn(-u^2 + v^2 + (u + x)^2 - (v + y)^2) + 1/2 v^2 sgn(u^2 - v^2 - (u + x)^2 + (v + y)^2) + v y sgn(u^2 - v^2 - (u + x)^2 + (v + y)^2) + u^2 + u x + v^2 + v y + x^2/2 + y^2/2, x>0, y>0, u>=1, v>=1}

Input 3:

Solve [(x + z) * z + (y + w) * w <= max((x + z) ^ 2 + w ^ 2, (y + w) ^ 2 + z ^ 2), x > 0, y > 0, z >= 1, w >= 1]

Solve

Result (not solution):

solve (x + z) z + (y + w) w<=max((x + z)^2 + w^2, (y + w)^2 + z^2)
x>0
y>0
z>=1
w>=1

If try to change a bit the 3d input:

Input 3 (changed):

Solve [(x + z) * z + (y + w) * y <= max((x + z) ^ 2 + w ^ 2, (y + w) ^ 2 + z ^ 2), x > 0, y > 0, z >= 1, w >= 1]

Solve

Solution:

w>=1 and x>0 and 0<y<=sqrt(2 w^2 + x^2) and z>=1
w>=1 and x>0 and y>sqrt(2 w^2 + x^2) and 1<=z<=(w^2 + w y)/x
w>=1 and x>0 and y>sqrt(2 w^2 + x^2) and z>=(-w^2 + w y - x^2 + y^2)/x
  1. Is that supposed to work like that?
  2. What the input and solution is correct?
  3. Why the 3d solution is not complete?