double f(double x, double y, double z, double t) {
double r515215 = x;
double r515216 = 1.0;
double r515217 = r515215 + r515216;
double r515218 = sqrt(r515217);
double r515219 = sqrt(r515215);
double r515220 = r515218 - r515219;
double r515221 = y;
double r515222 = r515221 + r515216;
double r515223 = sqrt(r515222);
double r515224 = sqrt(r515221);
double r515225 = r515223 - r515224;
double r515226 = r515220 + r515225;
double r515227 = z;
double r515228 = r515227 + r515216;
double r515229 = sqrt(r515228);
double r515230 = sqrt(r515227);
double r515231 = r515229 - r515230;
double r515232 = r515226 + r515231;
double r515233 = t;
double r515234 = r515233 + r515216;
double r515235 = sqrt(r515234);
double r515236 = sqrt(r515233);
double r515237 = r515235 - r515236;
double r515238 = r515232 + r515237;
return r515238;
}