Efficient simultaneous permutation of two large matrices in R

Question

I have two objects, both matrices:

pop1[,1:5]
>            8:31:+    8:91:+     8:106:+      8:133:+        8:137:+
> [1,]        42        39