test_aklt_open.m
9ee3751175efmaster

test_aklt_open.m
View Options

			% 2 sites
			disp("2 sites")
			sites_lib = load('files3/aklt2sites_open.mat');
			psi_2 = sites_lib.psiaklt2;

			[Alist_2, slist_2, Vlist_2] = iterativeschmidt(psi_2, 1e-15)

			origpsi_2 = recpsi(Alist_2, Vlist_2{end});

			error = max(abs(origpsi_2(:)-psi_2(:)));
			assert(error < 1e-12, "aklt2sites_open failed the recompose");
			disp("Recomposition of aklt2sites_open ended up successfully")

			fprintf("Occupied size original psi: %d\n", occupation({psi_2}))
			fprintf("Occupied size compressed psi: %d\n", occupation({Alist_2{:}, Vlist_2{end}}))

			disp(" ")




			% more sites
			disp("More sites")
			aklt_lib = load('files3/aklt_open.mat');
			psi = aklt_lib.psiaklt;

			[Alist, slist, Vlist] = iterativeschmidt(psi, 1e-15)
			origpsi = recpsi(Alist, Vlist{end});

			error = max(abs(origpsi(:)-psi(:)));
			assert(error < 1e-12, "aklt_open failed the recompose");
			disp("Recomposition of aklt_open ended up successfully")

			fprintf("Occupied size original psi: %d\n", occupation({psi}))
			fprintf("Occupied size compressed psi: %d\n", occupation({Alist{:}, Vlist{end}}))

			function total = occupation(arr)
			total = 0;
			for k = 1:length(arr)
			total = total + prod(size(arr{k}));
			end
			end