invert eigenvector V-

Resample/Flux/MatrixFlux.cpp

@@ -34,7 +34,7 @@ SpinMatrix PhaseRotationMatrix(const R3& B_direction, const Spinor& diagonal)
 
     if (std::abs(e.mag() - 1.) < eps) {
         const double factor1 = 2. * (1. + e.z());
-        SpinMatrix Q(1. + e.z(), -e.x() + I * e.y(), e.x() + I * e.y(), e.z() + 1.);
+        SpinMatrix Q(1. + e.z(), e.x() - I * e.y(), e.x() + I * e.y(), -e.z() - 1.);
         return Q * M * Q.adjoint() / factor1;
     }
     ASSERT(e.mag() < eps); // remaining case: no field