(1) rotation loop = generator of pi_1(SO(3)); its square is null-homotopic: OK
(2) tau well defined, restricts to x->Ax+b on the boundary, tau^k = id: 4000 random tests OK over 10 (lattice, A) cases
