double f(double x, double y, double z, double t) {
double r20623857 = x;
double r20623858 = 1.0;
double r20623859 = r20623857 + r20623858;
double r20623860 = sqrt(r20623859);
double r20623861 = sqrt(r20623857);
double r20623862 = r20623860 - r20623861;
double r20623863 = y;
double r20623864 = r20623863 + r20623858;
double r20623865 = sqrt(r20623864);
double r20623866 = sqrt(r20623863);
double r20623867 = r20623865 - r20623866;
double r20623868 = r20623862 + r20623867;
double r20623869 = z;
double r20623870 = r20623869 + r20623858;
double r20623871 = sqrt(r20623870);
double r20623872 = sqrt(r20623869);
double r20623873 = r20623871 - r20623872;
double r20623874 = r20623868 + r20623873;
double r20623875 = t;
double r20623876 = r20623875 + r20623858;
double r20623877 = sqrt(r20623876);
double r20623878 = sqrt(r20623875);
double r20623879 = r20623877 - r20623878;
double r20623880 = r20623874 + r20623879;
return r20623880;
}