What has Maximize been thinking about for so long?

The command

Maximize[{Sqrt[2 x + 13] + (3 y + 5)^(1/3) + (8 z + 12)^(1/4), x + y + z == 3&&{x,y,z} >= 0}, {x, y, z}] 

is running without any response on my comp for hours (as well as many other Mathematica commands). This optimization problem is quite standard and easily solved by Lagrange multipliers

L=Sqrt[2x+13]+(3y+5)^(1/3)+(8z+12)^(1/4)-t*(x+y+z-3); 

For simplicity the conditions {x,y,z}>=0 are not taken into account.

Reduce[D[L, x] == 0 && D[L, y] == 0 && D[L, z] == 0 && D[L, t] == 0, {x, y, z, t}, Reals] 

x == 3/2 && y == 1 && z == 1/2 && t == 1/4

The Hessian of L at this point is easily found by

M = {{D[L, x, x], D[L, x, y], D[L, x, z], D[L, x, t]}, {D[L, y, x],  D[L, y, y], D[L, y, z], D[L, y, t]}, {D[L, z, x], D[L, z, y],  D[L, z, z], D[L, z, t]}, {D[L, t, x], D[L, t, y], D[L, t, z],  D[L, t, t]}} /. {x -> 3/2, y -> 1, z -> 1/2, t -> 1/4} 

Next, the result of

Resolve[ForAll[{a,b,c,d}, ({a,b,c,d}.M).{a,b,c,d} <= 0], PositiveReals] 

True

states this is the maximum point (local and global) since the quadratic form ({a,b,c,d}.M).{a,b,c,d} takes nonpositive values on the positive reals. Is there anoter way to symbolically solve this optimization problem?