4 3 GfsPoisson GfsBox GfsGEdge {} {
  Time { iend = 10 }
  Refine LEVEL
  ApproxProjectionParams { nrelax = 4 tolerance = 1e-30 erelax = 2 }
  Init {} {
    Div = {
      int k = 3, l = 3;
      x = (x - 0.5)/2.;
      y = (y + 0.5)/2.;
      return -M_PI*M_PI*(k*k + l*l)*sin (M_PI*k*x)*sin (M_PI*l*y);
    }
  }
  OutputTime { istep = 1 } {
    awk '{print n++, $8}' > time
  }
  OutputProjectionStats { istep = 1 } {
    awk '{
      if ($1 == "niter:") printf ("%d ", $2);
      if ($1 == "residual.infty:") print $3 " " $4;
    }' > proj
  }
  OutputSimulation { start = end } sim-LEVEL { variables = P }
}
GfsBox {}
GfsBox {}
GfsBox {}
GfsBox {}
1 2 right
2 3 bottom
3 4 left