double f(double x, double y, double z, double t) {
double r25809916 = x;
double r25809917 = 1.0;
double r25809918 = r25809916 + r25809917;
double r25809919 = sqrt(r25809918);
double r25809920 = sqrt(r25809916);
double r25809921 = r25809919 - r25809920;
double r25809922 = y;
double r25809923 = r25809922 + r25809917;
double r25809924 = sqrt(r25809923);
double r25809925 = sqrt(r25809922);
double r25809926 = r25809924 - r25809925;
double r25809927 = r25809921 + r25809926;
double r25809928 = z;
double r25809929 = r25809928 + r25809917;
double r25809930 = sqrt(r25809929);
double r25809931 = sqrt(r25809928);
double r25809932 = r25809930 - r25809931;
double r25809933 = r25809927 + r25809932;
double r25809934 = t;
double r25809935 = r25809934 + r25809917;
double r25809936 = sqrt(r25809935);
double r25809937 = sqrt(r25809934);
double r25809938 = r25809936 - r25809937;
double r25809939 = r25809933 + r25809938;
return r25809939;
}