private void loop() {
   float max = 0;
   int counter = 0;
   while (nodes.size() != 0) {
     Point2DInteger n = nodes.peek();
     ArrayList<Point2DInteger> children = FlowDirectionUtils.getUpstreamNodes(n, flowDirection);
     if (children.size() == 0) {
       flowAccumulation.set(n.getX(), n.getY(), 0);
       nodes.pop();
     } else {
       float sum = 0;
       boolean hasUndefined = false;
       for (Point2DInteger p : children) {
         float acc = flowAccumulation.get(p.getX(), p.getY());
         if (Float.isNaN(acc)) {
           hasUndefined = true;
           nodes.add(p);
         } else {
           sum += acc;
           sum += getArea(p);
         }
       }
       if (!hasUndefined) {
         max = Math.max(max, sum);
         flowAccumulation.set(n.getX(), n.getY(), sum);
         nodes.pop();
       }
     }
     counter++;
     if ((counter & 0xffff) == 0xffff) {
       System.out.println("FlowAcc " + counter + " " + max);
     }
   }
   System.out.println("Done");
 }
 private void clearFlowAccumulation() {
   for (Rectangle2DInteger rect : this.range) {
     int x0 = rect.getX();
     int y0 = rect.getY();
     int w = rect.getWidth();
     int h = rect.getHeight();
     int x1 = x0 + w;
     int y1 = y0 + h;
     for (int y = y0; y < y1; y++) {
       for (int x = x0; x < x1; x++) {
         flowAccumulation.set(x, y, Float.NaN);
       }
     }
   }
 }