double f(double x, double y, double z, double t) {
double r325628 = x;
double r325629 = 1.0;
double r325630 = r325628 + r325629;
double r325631 = sqrt(r325630);
double r325632 = sqrt(r325628);
double r325633 = r325631 - r325632;
double r325634 = y;
double r325635 = r325634 + r325629;
double r325636 = sqrt(r325635);
double r325637 = sqrt(r325634);
double r325638 = r325636 - r325637;
double r325639 = r325633 + r325638;
double r325640 = z;
double r325641 = r325640 + r325629;
double r325642 = sqrt(r325641);
double r325643 = sqrt(r325640);
double r325644 = r325642 - r325643;
double r325645 = r325639 + r325644;
double r325646 = t;
double r325647 = r325646 + r325629;
double r325648 = sqrt(r325647);
double r325649 = sqrt(r325646);
double r325650 = r325648 - r325649;
double r325651 = r325645 + r325650;
return r325651;
}