(a,b)=(1,1): det[V_R;L]/det[ReS;ImS] at 4 points: [0.33227572 0.33227572 0.33227572 0.33227572]  | (2):  [1.-0.j 1.-0.j 1.-0.j 1.+0.j]
(a,b)=(1,2): det[V_R;L]/det[ReS;ImS] at 4 points: [0.2578316 0.2578316 0.2578316 0.2578316]  | (2):  [1.+0.j 1.-0.j 1.+0.j 1.+0.j]
(a,b)=(2,2): det[V_R;L]/det[ReS;ImS] at 4 points: [-0.03788096 -0.03788096 -0.03788096 -0.03788096]  | (2):  [1.+0.j 1.+0.j 1.+0.j 1.+0.j]
(a,b)=(2,3): det[V_R;L]/det[ReS;ImS] at 4 points: [-0.03208408 -0.03208408 -0.03208408 -0.03208408]  | (2):  [1.-0.j 1.+0.j 1.-0.j 1.+0.j]
(a,b)=(3,3): det[V_R;L]/det[ReS;ImS] at 4 points: [-0.00155118 -0.00155118 -0.00155118 -0.00155118]  | (2):  [1.-0.j 1.-0.j 1.-0.j 1.+0.j]
j=1: det[[s,0,q],[0,sbar,qbar]] = (-0-1.76028294j)  tau2 conj(tau1) - tau1 conj(tau2) = -1.76028294j  ratio det[S;Sbar]/that = -2j
j=2: det[ReS;ImS] = -20.817314  |det S'|^2 = 20.817314
j=1: det[[s,0,q],[0,sbar,qbar]] = (-0+4.44588311j)  tau2 conj(tau1) - tau1 conj(tau2) = 4.44588311j  ratio det[S;Sbar]/that = (-0-2j)
j=2: det[ReS;ImS] = -8.14962017  |det S'|^2 = 8.14962017
