double f(double x, double y, double z, double t) {
double r25771525 = x;
double r25771526 = 1.0;
double r25771527 = r25771525 + r25771526;
double r25771528 = sqrt(r25771527);
double r25771529 = sqrt(r25771525);
double r25771530 = r25771528 - r25771529;
double r25771531 = y;
double r25771532 = r25771531 + r25771526;
double r25771533 = sqrt(r25771532);
double r25771534 = sqrt(r25771531);
double r25771535 = r25771533 - r25771534;
double r25771536 = r25771530 + r25771535;
double r25771537 = z;
double r25771538 = r25771537 + r25771526;
double r25771539 = sqrt(r25771538);
double r25771540 = sqrt(r25771537);
double r25771541 = r25771539 - r25771540;
double r25771542 = r25771536 + r25771541;
double r25771543 = t;
double r25771544 = r25771543 + r25771526;
double r25771545 = sqrt(r25771544);
double r25771546 = sqrt(r25771543);
double r25771547 = r25771545 - r25771546;
double r25771548 = r25771542 + r25771547;
return r25771548;
}