# gerris2D star.gfs
# gfsview sim-6 solution.gfv

1 0 GfsPoisson GfsBox GfsGEdge {} {

    Time { iend = 1 }

    Refine LEVEL

    Solid ({
            double dr = 0.1;
            double theta = atan2 (y, x);
            double radius = 0.79*(0.45 - dr + dr*cos (6.*theta));
            return x*x + y*y - radius*radius;
          })

#    GModule hypre

    ApproxProjectionParams { erelax = 2 tolerance = 1e-30 nitermin = CYCLE nitermax = CYCLE }

    Init {} {
        Div = {
            int k = 3, l = 3;
            return -M_PI*M_PI*(k*k + l*l)*sin (M_PI*k*x)*sin (M_PI*l*y);
        }
    }

    OutputSimulation { start = end } sim-LEVEL { variables = P }
}

GfsBox {}