domain dirichlet_boundary (const geo& omega) {
  switch (omega.dimension()) {
    case 1:  return omega["left"];
    case 2:  return omega["left"]  + omega["bottom"];
    default: return omega["left"]  + omega["bottom"] 
                  + omega["right"] + omega["back"];
  }
}
