SetDirectory["/www/user/fdahl/papers/Conjugation/"];
<<kappaLib.m
<<Petrov.m

Class V: (11 2 bar(2))

B=(({{0, 0, 0, 1, 0, 0}, {0, 0, 0, 0, 1, 0}, {0, 0, 0, 0, 0, 1}, {1, 0, 0, 0, 0, 0}, {0, 1, 0, 0, 0, 0}, {0, 0, 1, 0, 0, 0}}));

V=(({{ lambda1, 0, 0, 0, 0, 0}, {0, lambda2, 0, 0, 0, 0}, {0, 0, sigma1, tau1, 1, 0}, {0, 0, -tau1, sigma1, 0, 1}, {0, 0, 0, 0,  sigma1, tau1}, {0, 0, 0, 0, -tau1,  sigma1}}));

Eigenvalues[V]

W=(({{eps1, 0, 0, 0, 0, 0}, {0, eps2, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 1}, {0, 0, 0, 0, 1, 0}, {0, 0, 0, 1, 0, 0}, {0, 0, 1, 0, 0, 0}}));

Eigenvalues[W]

(* We know that NN should have signature (3,3) and eps1 ≤ eps2. *)

W=W/.{eps1->-1,eps2->1};
W//MatrixForm
Eigenvalues[W]

(* Permutation 6,6 with leading B *)
S=(({{0, 0, 0, 0, 0, 1}, {0, 0, 0, 0, 1, 0}, {(1/(Sqrt[2])), (1/(Sqrt[2])), 0, 0, 0, 0}, {0, 0, 1, 0, 0, 0}, {0, 0, 0, 1, 0, 0}, {-(1/(Sqrt[2])), (1/(Sqrt[2])), 0, 0, 0, 0}}));

Check that S is in set mathcal(S)

Transpose[S].B.S==W

Compute result

res=S.V.Inverse[S];
res//MatrixForm

Petrov[res]

Export notebook as .pdf

NotebookPrint[SelectedNotebook[],"/www/user/fdahl/papers/Conjugation/notebooks/ClassV.pdf"]