Function main

Tests/Parser/main.cpp:237–621

}

int main (int argc, char* argv[])
{
    amrex::Initialize(argc, argv);

    {
        amrex::Print() << "\n";
        int nerror = 0;
        nerror += test3("if( ((z-zc)*(z-zc)+(y-yc)*(y-yc)+(x-xc)*(x-xc))^(0.5) < (r_star-dR), 0.0, if(((z-zc)*(z-zc)+(y-yc)*(y-yc)+(x-xc)*(x-xc))^(0.5) <= r_star, dens, 0.0))",
                        {{"xc", 0.1}, {"yc", -1.0}, {"zc", 0.2}, {"r_star", 0.73}, {"dR", 0.57}, {"dens", 12.}},
                        {"x","y","z"},
                        [=] (double x, double y, double z) -> double {
                            double xc=0.1, yc=-1.0, zc=0.2, r_star=0.73, dR=0.57, dens=12.;
                            double r = std::sqrt((z-zc)*(z-zc) + (y-yc)*(y-yc) + (x-xc)*(x-xc));
                            if (r >= r_star-dR && r <= r_star) {
                                return dens;
                            } else {
                                return 0.0;
                            }
                        },
                        {-1., -1., -1.0}, {1.0, 1.0, 1.0}, 100,
                        1.e-12, 1.e-15);

        nerror += test3("r=sqrt((z-zc)*(z-zc)+(y-yc)*(y-yc)+(x-xc)*(x-xc)); if(r < (r_star-dR), 0.0, if(r <= r_star, dens, 0.0))",
                        {{"xc", 0.1}, {"yc", -1.0}, {"zc", 0.2}, {"r_star", 0.73}, {"dR", 0.57}, {"dens", 12.}},
                        {"x","y","z"},
                        [=] (double x, double y, double z) -> double {
                            double xc=0.1, yc=-1.0, zc=0.2, r_star=0.73, dR=0.57, dens=12.;
                            double r = std::sqrt((z-zc)*(z-zc) + (y-yc)*(y-yc) + (x-xc)*(x-xc));
                            if (r >= r_star-dR && r <= r_star) {
                                return dens;
                            } else {
                                return 0.0;
                            }
                        },
                        {-1., -1., -1.0}, {1.0, 1.0, 1.0}, 100,
                        1.e-12, 1.e-15);

        nerror += test3("r=(z-zc)*(z-zc)+(y-yc)*(y-yc)+(x-xc)*(x-xc); r=sqrt(r); if(r < (r_star-dR), 0.0, if(r <= r_star, dens, 0.0))",
                        {{"xc", 0.1}, {"yc", -1.0}, {"zc", 0.2}, {"r_star", 0.73}, {"dR", 0.57}, {"dens", 12.}},
                        {"x","y","z"},
                        [=] (double x, double y, double z) -> double {
                            double xc=0.1, yc=-1.0, zc=0.2, r_star=0.73, dR=0.57, dens=12.;
                            double r = std::sqrt((z-zc)*(z-zc) + (y-yc)*(y-yc) + (x-xc)*(x-xc));
                            if (r >= r_star-dR && r <= r_star) {
                                return dens;
                            } else {
                                return 0.0;
                            }
                        },
                        {-1., -1., -1.0}, {1.0, 1.0, 1.0}, 100,
                        1.e-12, 1.e-15);

        nerror += test3("( ((( (z-zc)*(z-zc) + (y-yc)*(y-yc) + (x-xc)*(x-xc) )^(0.5))<=r_star) * ((( (z-zc)*(z-zc) + (y-yc)*(y-yc) + (x-xc)*(x-xc) )^(0.5))>=(r_star-dR)) )*dens",
                        {{"xc", 0.1}, {"yc", -1.0}, {"zc", 0.2}, {"r_star", 0.73}, {"dR", 0.57}, {"dens", 12.}},
                        {"x","y","z"},
                        [=] (double x, double y, double z) -> double {
                            double xc=0.1, yc=-1.0, zc=0.2, r_star=0.73, dR=0.57, dens=12.;
                            double r = std::sqrt((z-zc)*(z-zc) + (y-yc)*(y-yc) + (x-xc)*(x-xc));