%test sim_trans_t4()



global_defs;



example = 1;

switch example

	case (1)

		lam = 2*pi*rand();

		c_lam = cos(lam);

		s_lam = sin(lam);

		phi = 2*pi*rand();

		c_phi = cos(phi);

		s_phi = sin(phi);

		

		a0 = [c_lam;s_lam;0];

		a1 = a0;

		b0 = [1;0;0];

		b1 = b0;

		c0 = [ s_phi; 0; c_phi];

		c1 = [-s_phi; 0; c_phi];

	otherwise

		error("example number is out of range");

end

vecs_in=[c0,b0,a0;c1,b1,a1];

[vecs_out, U1, U0] = sim_trans_t4(vecs_in);

a0f=vecs_out(1:3,1);

a1f=vecs_out(4:6,1);



% c0   b0   a0         a0f   [U0]

% |    |    |    -->   |           

% c1   b1   a1         a1f   [U1]



%Let {sx,sy,sz} be the standard basis. Then 

%a0=a1=cos(lam)sx -sin(lam)sy

%b0=b1=sx

%c0=cos(phi)sz +sin(phi)sx

%c1=cos(phi)sz -sin(phi)sx





error = norm( 

dr11(c1,c0)*dr11(b1,b0)*dr11(a1,a0) - dr11(a1f,a0f)*kron(U1,U0))