Daniel@0: function v = wrap(u,N)
Daniel@0: % WRAP Wrap a vector of indices around a torus.
Daniel@0: % v = wrap(u,N)
Daniel@0: %
Daniel@0: % e.g., wrap([-1 0 1 2 3 4], 3)   =   2 3 1 2 3 1
Daniel@0: 
Daniel@0: v = mod(u-1,N)+1;