Home > src > main > matlab > flaglet_perf_tests.m

flaglet_perf_tests

PURPOSE ^

flaglet_perf_tests

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

 flaglet_perf_tests

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 % flaglet_perf_tests
0002 
0003 ell = [ 4 8 16 32 64 128 256 512 ];
0004 accuracy = [ ...
0005     3.49e-15 ... % 4
0006     2.64e-14 ... % 8
0007     9.02e-14 ... % 16
0008     9.08e-13 ... % 32
0009     2.93e-12 ... % 64
0010     2.56e-11 ... % 128
0011     1.1e-10 ... % 256
0012     4.56e-10 ... % 512
0013     ];
0014 speed = [ ...
0015     0.0036 ...  % 4
0016     0.022 ... % 8
0017     0.3 ... % 16
0018     3.3 ... % 32
0019     50.0 ... % 64
0020     880 ... % 128
0021     10000 ... % 256
0022     110000 ... % 512
0023     ];
0024 accuracy_multires = [ ...
0025     3.31e-15 ... % 4
0026     2.27e-14 ... % 8
0027     8.42e-14 ... % 16
0028     8.15e-13 ... % 32
0029     2.86e-12 ... % 64
0030     2.35e-11 ... % 128
0031     1.1e-10 ... % 256
0032     4.56e-10 ... % 512
0033     ];
0034 speed_multires = [ ...
0035     0.00030 ...  % 4
0036     0.005 ... % 8
0037     0.08 ... % 16
0038     0.65 ... % 32
0039     8.0 ... % 64
0040     120 ... % 128
0041     2300 ... % 256
0042     23000 ... % 512
0043     ];
0044 
0045 nb_samples = 2*ell.^3;
0046 nb_samples_pow = str2mat('t04', 't07', 't10', 't13', 't16', 't19', 't22', 't25', 't28');
0047 
0048 figure('Position',[100 100 700 700])
0049 
0050 subplot(2,1,1)
0051 loglog( nb_samples, (1e-14)*ell.^2, 'red','LineWidth', 2 )
0052 hold on
0053 loglog( nb_samples, accuracy, '--oblack', 'LineWidth', 2, 'MarkerSize',10)
0054 loglog( nb_samples, accuracy_multires, '-oblack', 'LineWidth', 2, 'MarkerSize',10)
0055 grid on
0056 xlabel('x1','FontSize',20)
0057 ylabel('y1','FontSize',20)
0058 set(gca,'XTick',nb_samples,'XTickLabel',nb_samples_pow,'FontSize',20)
0059 set(gca,'YTick',[1e-16 1e-14 1e-12  1e-10 1e-8] ,'FontSize',20)
0060 axis([80 1.5*nb_samples(8) 1e-16 1e-8])
0061 
0062 %title('Accuracy of overall transform','FontSize',20)
0063 
0064 subplot(2,1,2)
0065 loglog( nb_samples, (5e-5)*ell.^4, 'red','LineWidth', 2 )
0066 hold on
0067 loglog( nb_samples, speed, '--oblack', 'LineWidth', 2, 'MarkerSize',10)
0068 loglog( nb_samples, speed_multires, '-oblack', 'LineWidth', 2, 'MarkerSize',10)
0069 grid on
0070 xlabel('x2','FontSize',20)
0071 ylabel('y2','FontSize',20)
0072 axis([80 1.5*nb_samples(8) 1e-5 2e7])
0073 set(gca,'YTick',[1e-6 1e-4 1e-2  1e0 1e2 1e4 1e6] ,'FontSize',20)
0074 set(gca,'XTick',nb_samples,'XTickLabel',nb_samples_pow,'FontSize',20)
0075 %title('Speed of overall transform','FontSize',20)
0076

Generated on Mon 24-Sep-2012 12:26:33 by m2html © 2005