staggered_NPR.py: Implement staggered matrices.  SLOW.