double f(double x, double y, double z, double t) {
double r25828520 = x;
double r25828521 = 1.0;
double r25828522 = r25828520 + r25828521;
double r25828523 = sqrt(r25828522);
double r25828524 = sqrt(r25828520);
double r25828525 = r25828523 - r25828524;
double r25828526 = y;
double r25828527 = r25828526 + r25828521;
double r25828528 = sqrt(r25828527);
double r25828529 = sqrt(r25828526);
double r25828530 = r25828528 - r25828529;
double r25828531 = r25828525 + r25828530;
double r25828532 = z;
double r25828533 = r25828532 + r25828521;
double r25828534 = sqrt(r25828533);
double r25828535 = sqrt(r25828532);
double r25828536 = r25828534 - r25828535;
double r25828537 = r25828531 + r25828536;
double r25828538 = t;
double r25828539 = r25828538 + r25828521;
double r25828540 = sqrt(r25828539);
double r25828541 = sqrt(r25828538);
double r25828542 = r25828540 - r25828541;
double r25828543 = r25828537 + r25828542;
return r25828543;
}