I have a function that when I operate FullSimplify on, produces widly different results depending on if the parameters are left symbolic or values are substituted in before the simplification. I've narrowed down the part of my function that FullSimplify messes up.:

func = E^(-((-(x^2 + y^2 + z^2)^
     p + \[Rho])^2/\[Sigma])) + E^(-(((x^2 + y^2 + z^2)^
    p + \[Rho])^2/\[Sigma]))

This is just two gaussians translated symmetrically about the origin. Now, If I do FullSimplify on this function, I get the expected result which is that basically nothing happens. The function is already pretty simple. But if I input some parameter values:

E^(-33.3333 (0.8 - (x^2 + y^2 + 
      z^2)^0.26)^2) + E^(-33.3333 (0.8 + (x^2 + y^2 + z^2)^0.26)^2)

Then pass FullSimplify to this, Mathematica reduces it to zero.

FullSimplify[E^(-33.3333 (0.8 - (x^2 + y^2 + 
      z^2)^0.26)^2) + E^(-33.3333 (0.8 + (x^2 + y^2 + z^2)^0.26)^2)]

---> (*0.*)

This is really strange and almost introduced mistakes into my research since I didnt catch it at first. Luckily, I did some consistency checks and analysis and found the mistake. Unfortunately, the functions I work with are really large, so I would like to use simplification early on in the code so that it reduces computation time later.

Is there a way I can still use simplification without Mathematica killing the result? I see there have been some other issues with FullSimplify, but I think this one is unique enough to warrant a separate question, given the use of Exponentials.

I appreciate any insight and help!