Index: kernel/test/states/test_nonbonded_list.py
===================================================================
--- kernel/test/states/test_nonbonded_list.py	(revision 393)
+++ kernel/test/states/test_nonbonded_list.py	(working copy)
@@ -46,7 +46,6 @@
 
     def test_bl(self):
         """Test the bonded list"""
-        return
         m= IMP.Model()
         bds=[]
         for i in range(0,10):
Index: kernel/include/IMP/internal/ParticleGrid.h
===================================================================
--- kernel/include/IMP/internal/ParticleGrid.h	(revision 393)
+++ kernel/include/IMP/internal/ParticleGrid.h	(working copy)
@@ -43,7 +43,7 @@
   typedef Grid::VirtualIndex VirtualIndex;
   typedef Grid::Index Index;
   Grid::VirtualIndex get_virtual_index(Vector3D pt) const {
-    return grid_.get_index(pt);
+    return grid_.get_virtual_index(pt);
   }
 
   //! Apply the function F to all particles near to the center
Index: kernel/src/score_states/BipartiteNonbondedListScoreState.cpp
===================================================================
--- kernel/src/score_states/BipartiteNonbondedListScoreState.cpp	(revision 393)
+++ kernel/src/score_states/BipartiteNonbondedListScoreState.cpp	(working copy)
@@ -49,7 +49,8 @@
     XYZDecorator d= XYZDecorator::cast(p);
     internal::ParticleGrid::VirtualIndex index
       = grid_.get_virtual_index(Vector3D(d.get_x(), d.get_y(), d.get_z()));
-    //IMP_LOG(VERBOSE, "Index is " << index << std::endl);
+    IMP_LOG(VERBOSE, "Searching for " << p->get_index()
+            << " from cell " << index << std::endl);
     grid_.apply_to_nearby(f, index, cut, false);
   }
   IMP_LOG(VERBOSE, "Found " << P::size_nbl()