test_cj.m
View Options

			% This script tests eig_power.m for correctness

			test_case_names = {
			'diag(1:2)';
			'diag(1:3)';
			'diag(1:10)';
			'rmg(1:2)';
			'rmg(1:3)';
			'rmg(1:10)';
			'rmg(-2.9:10.1)';
			'rmg(-9.9:3.1)';
			};

			N = length(test_case_names);

			for i=1:N

			fprintf('Test %d: %s ', i, test_case_names{i});

			rng(i);

			try
			[~, m_rl, m_cp] = evalc(test_case_names{i});
			catch e
			[~, m_rl] = evalc(test_case_names{i});
			m_cp = m_rl;
			end
			reference = sort(eig(m_rl));

			[P, D, ~] = eig_cj(m_rl); % cyclic

			val = sort(diag(D));

			% verify uniform convergence
			delta = norm(m_rl * P - P * D);

			assert(all(abs(val-reference)<1e-5), "Eigenvalues not corresponding");
			assert(delta<1e-5, "A * P - D * P != 0");

			fprintf('OK\n');
			end

compphys/exo3/test_cj.m33bba04159d6master