test_solve_SD.m
e691a16f50fcmaster

test_solve_SD.m
View Options

			%This script tests solve_SD.m
			%Uses rmg.m

			fprintf('Test 1: Identity matrix ...');
			for n=1:1:10
			A=eye(n);
			X=rand(n,1);
			B=A*X;

			assert(all(abs(solve_SD(A,B) - X)<1e-8));
			end
			fprintf('\tpassed\n');

			fprintf('Test 2: Random matrices')
			for n=1:1:5
			A=rand(n);
			A=A'*A;
			while(cond(A)>90000)
			A=rand(n);
			A=A'*A;
			end
			X=rand(n,1);
			B=A*X;
			assert(norm(solve_SD(A,B)-X) <1e-8);
			end

			fprintf('\tpassed\n');

			fprintf('Test 3: Random matrices')
			for n=1:1:5
			A=rmg(rand(n));
			X=n*rand(n,1);
			B=A*X;


			assert(norm(solve_SD(A,B)-X)<1e-8);
			end

			fprintf('\tpassed\n');


			fprintf('All tests passed!\n');