double f(double x, double y, double z, double t) {
double r27816225 = x;
double r27816226 = 1.0;
double r27816227 = r27816225 + r27816226;
double r27816228 = sqrt(r27816227);
double r27816229 = sqrt(r27816225);
double r27816230 = r27816228 - r27816229;
double r27816231 = y;
double r27816232 = r27816231 + r27816226;
double r27816233 = sqrt(r27816232);
double r27816234 = sqrt(r27816231);
double r27816235 = r27816233 - r27816234;
double r27816236 = r27816230 + r27816235;
double r27816237 = z;
double r27816238 = r27816237 + r27816226;
double r27816239 = sqrt(r27816238);
double r27816240 = sqrt(r27816237);
double r27816241 = r27816239 - r27816240;
double r27816242 = r27816236 + r27816241;
double r27816243 = t;
double r27816244 = r27816243 + r27816226;
double r27816245 = sqrt(r27816244);
double r27816246 = sqrt(r27816243);
double r27816247 = r27816245 - r27816246;
double r27816248 = r27816242 + r27816247;
return r27816248;
}