% Refinements
Ns = [10, 14, 20, 24, 30, 34, 40, 44];

methods_to_test = { '6s3k4p', '7s4k5p', '8s3k6p', '8s4k7p', '9s4k8p', '9s4k9p' };

c = 1.0;

diary('pde_convergence_log.txt')

for method=methods_to_test

	method = method{:};
	dudt = SSP_Tools.Integrators.MSRK('coefficients', method, 'initial_integrator', 'RK3', 'mini_dt_type', 'compatible', 'mini_dt_c', c);
	dudx = SSP_Tools.Discretizers.Spectral();
	
	problem = SSP_Tools.TestProblems.Advection('a', -1.0, ...
	                                           'integrator', dudt, ...
	                                           'discretizer', dudx, ...
	                                           'domain', [0, 1], ...
	                                           'initial_condition', @(x) sin(4*pi*x) );
	                                           
	Test = SSP_Tools.Tests.ConvergencePDE('problem', problem, ...
	                                      'refinements', Ns, ...
	                                      'refinement_type', 'time_space', ...
	                                      't', 1.0, ...
	                                      'cfl', 0.4 );

	Test.run_test();
	
	% Print the test buffer to the console 
	cellfun(@fprintf, Test.output_buffer)
	results = Test.results;
	
	format long;
	l2_errors = [results.l2error]'
	
end

diary('off');