1 0 GfsPoisson GfsBox GfsGEdge {} { Time { iend = 1 } Refine LEVEL RefineSolid (LEVEL + 2) 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 { 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); } } OutputTime { istep = 1 } { awk '{if ($2 == 1) {print "CYCLE " $8}}' >> time } OutputProjectionStats { start = end } { awk '{ if ($1 == "residual.infty:") print "CYCLE "$3 " " $4; }' >> proj } OutputSimulation { start = end } sim-LEVEL-SOLVER { variables = P } } GfsBox {}