double f(double x, double y, double z, double t) {
double r24625387 = x;
double r24625388 = 1.0;
double r24625389 = r24625387 + r24625388;
double r24625390 = sqrt(r24625389);
double r24625391 = sqrt(r24625387);
double r24625392 = r24625390 - r24625391;
double r24625393 = y;
double r24625394 = r24625393 + r24625388;
double r24625395 = sqrt(r24625394);
double r24625396 = sqrt(r24625393);
double r24625397 = r24625395 - r24625396;
double r24625398 = r24625392 + r24625397;
double r24625399 = z;
double r24625400 = r24625399 + r24625388;
double r24625401 = sqrt(r24625400);
double r24625402 = sqrt(r24625399);
double r24625403 = r24625401 - r24625402;
double r24625404 = r24625398 + r24625403;
double r24625405 = t;
double r24625406 = r24625405 + r24625388;
double r24625407 = sqrt(r24625406);
double r24625408 = sqrt(r24625405);
double r24625409 = r24625407 - r24625408;
double r24625410 = r24625404 + r24625409;
return r24625410;
}