/** Surround all existing samples on the edge by a layer of closing samples. */
 public final void close() {
   if (closed) return;
   closed = true;
   List<ZSamplePoint<TZ>> processList = new ArrayList<ZSamplePoint<TZ>>(sampleGrid.size());
   for (ZSamplePoint<TZ> A : sampleGrid) {
     processList.add(A);
   }
   int n = 0;
   for (int i = 0; i < 2; i++) {
     List<ZSamplePoint<TZ>> newProcessList = new ArrayList<ZSamplePoint<TZ>>(processList.size());
     for (ZSamplePoint<TZ> A : processList) {
       if (A.right() == null) {
         newProcessList.add(closeSample(A.getX() + 1, A.getY()));
         n++;
       }
       if (A.left() == null) {
         newProcessList.add(closeSample(A.getX() - 1, A.getY()));
         n++;
       }
       if (A.up() == null) {
         newProcessList.add(closeSample(A.getX(), A.getY() + 1));
         n++;
       }
       if (A.down() == null) {
         newProcessList.add(closeSample(A.getX(), A.getY() - 1));
         n++;
       }
     }
     processList = newProcessList;
   }
   LOG.info("Added {} closing samples to get a total of {}.", n, sampleGrid.size());
 }
 public final void addSamplingPoint(Coordinate C0, double z) {
   if (closed) throw new IllegalStateException("Can't add a sample after closing.");
   int[] xy = sampleGrid.getLowerLeftIndex(C0);
   int x = xy[0];
   int y = xy[1];
   @SuppressWarnings("unchecked")
   ZSamplePoint<TZ>[] ABCD = new ZSamplePoint[4];
   ABCD[0] = sampleGrid.getOrCreate(x, y);
   ABCD[1] = sampleGrid.getOrCreate(x + 1, y);
   ABCD[2] = sampleGrid.getOrCreate(x, y + 1);
   ABCD[3] = sampleGrid.getOrCreate(x + 1, y + 1);
   for (ZSamplePoint<TZ> P : ABCD) {
     Coordinate C = sampleGrid.getCoordinates(P);
     P.setZ(metric.cumulateSample(C0, C, z, P.getZ()));
   }
 }
 private final ZSamplePoint<TZ> closeSample(int x, int y) {
   ZSamplePoint<TZ> A = sampleGrid.getOrCreate(x, y);
   A.setZ(
       metric.closeSample(
           A.up() != null ? A.up().getZ() : null,
           A.down() != null ? A.down().getZ() : null,
           A.right() != null ? A.right().getZ() : null,
           A.left() != null ? A.left().getZ() : null));
   return A;
 }