jacobi_transition.m
33bba04159d6master

jacobi_transition.m
View Options

			% apply a rotation to A: J(p,q)' A J(p,q)
			function [A, P] = transition(A, P, p, q, c, s, N)
			% left multiplication by J'
			for j = 1:N
			t = A(p,j);
			A(p,j) = c * t - s * A(q,j);
			A(q,j) = s * t + c * A(q,j);
			end

			% right multiplication by J
			for i = 1:N
			t = A(i,p);
			A(i,p) = c * t - s * A(i,q);
			A(i,q) = s * t + c * A(i,q);

			tp = P(i,p);
			P(i,p) = c * tp - s * P(i,q);
			P(i,q) = s * tp + c * P(i,q);
			end
			end