Question

I’ve got (binary) matrices represented by uint64_t (from C++11). And I’d like to be able to efficiently map from any column into first rank. For example

0 1 0 0 0 0 0 0
0 1 0 0 0 0 0 0
0 1 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0
0 1 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0

uint64_t matrice = 0x4040400040400040uLL;
uint64_t matrice_2 = map(matrice, ColumnEnum::Column2);

1 1 1 0 1 1 0 1
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0

matrice_2 contains 0xED00000000000000uLL;

Answer 1

Great question. I really enjoyed the hacking. Here is my solution:

uint64_t map(uint64_t x, int column)
{
  x = (x >> (7 - column)) & 0x0101010101010101uLL;
  x = (x | (x >> 7)) & 0x00FF00FF00FF00FFuLL;
  x = (x | (x >> 14))& 0x000000FF000000FFuLL;
  x = (x | (x >> 28))& 0x00000000000000FFuLL;
  return x << 56;
}

A working example can be found at ideone, where the call is really map(matrice, ColumnEnum::Column2).