ground_state.m
09de8b211167master

ground_state.m
View Options

			A = zeros(2,3,2);
			A(:,1,:) = [0, 1; 0, 0] / sqrt(2); % index 1
			A(:,2,:) = -[1, 0; 0, -1] / 2.0; % index 0
			A(:,3,:) = [0, 0; -1, 0] / sqrt(2); % index -1

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

			recon_2 = ground_recover(A, 2);

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

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

			recon = ground_recover(A, 8);

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

			% function used to recover the original state starting
			% from the tensor A and the number of sites N
			function wavefct = ground_recover(A, N)
			Alist = cell(1,N-1);
			for k = 1:N-1
			Alist{k} = A;
			end
			wavefct = recpsi(Alist, A);
			end