test_iterativeschmidt.m
944f93f58245master

test_iterativeschmidt.m
View Options


			fprintf('Test 1: Random 2 site tensor ...');
			psi = rand(2,4,4,2);

			threshold = 1e-15;

			[Alist, slist, Vlist]= iterativeschmidt(psi, threshold);

			newpsi = ncon({Alist{1},Vlist{end}},{[-1 -2 1],[1 -3 -4]});
			error = max(abs(newpsi(:)-psi(:)));
			assert(error < 1e-12);
			fprintf('\tpassed\n');

			fprintf('Test 1: Random 3 site tensor ...');
			psi = rand(2,4,4,4,2);

			threshold = 1e-15;

			[Alist, slist, Vlist]= iterativeschmidt(psi, threshold);

			newpsi = ncon({Alist{1},Alist{2}, Vlist{end}},{[-1 -2 1],[1 -3 2],[2 -4 -5]});
			error = max(abs(newpsi(:)-psi(:)));
			assert(error < 1e-12);
			fprintf('\tpassed\n');