double f(double x, double y, double z, double t) {
double r22925634 = x;
double r22925635 = 1.0;
double r22925636 = r22925634 + r22925635;
double r22925637 = sqrt(r22925636);
double r22925638 = sqrt(r22925634);
double r22925639 = r22925637 - r22925638;
double r22925640 = y;
double r22925641 = r22925640 + r22925635;
double r22925642 = sqrt(r22925641);
double r22925643 = sqrt(r22925640);
double r22925644 = r22925642 - r22925643;
double r22925645 = r22925639 + r22925644;
double r22925646 = z;
double r22925647 = r22925646 + r22925635;
double r22925648 = sqrt(r22925647);
double r22925649 = sqrt(r22925646);
double r22925650 = r22925648 - r22925649;
double r22925651 = r22925645 + r22925650;
double r22925652 = t;
double r22925653 = r22925652 + r22925635;
double r22925654 = sqrt(r22925653);
double r22925655 = sqrt(r22925652);
double r22925656 = r22925654 - r22925655;
double r22925657 = r22925651 + r22925656;
return r22925657;
}