add test cases.  add refinement of the neighbor test, which is in-efficent and not very numpy like