aklt2sites.m
9ee3751175efmaster

aklt2sites.m
View Options

			matlib = load('files3/aklt2sites.mat');

			psi = matlib.psi

			[U,s,V] = svd(psi)
			V = V';
			chi = rankof(s);

			fprintf('Bond dimension: %d\n', chi);

			% contract
			U = U(:,1:chi);
			s = s(1:chi,1:chi);
			V = V(1:chi, :);

			origpsi = ncon({U, s, V},{[-1 1],[1 2],[2 -2]})