edge fails because diff reacts go into edge. TODO: Make a shell of width
nneighborslimit around the site box which contain either pointer to
other end or a NULL ptr to avoid having to lookup boundaries during
solving. Will be groundwork for getting rid of neighborhoods as well.