s2let_perf_tests

 s2let_perf_tests

0001 % s2let_perf_tests
0002 
0003 ell = [ 4 8 16 32 64 128 256 512 1024 2048 4096 ];
0004 accuracy = [ ...
0005     2.5e-15 ... % 4
0006     3.97e-15 ... % 8
0007     8.0e-15 ... % 16
0008     2.1e-14 ... % 32
0009     4.10e-14 ... % 64
0010     7.32e-14 ... % 128
0011     1.38e-13 ... % 256
0012     2.31e-13 ... % 512
0013     5.2e-13 ... % 1024
0014     9.3e-13 ... % 2048
0015     2.0e-12 ... % 4096
0016     ];
0017 speed = [ ...
0018     0.00019 ...  % 4
0019     0.00038 ... % 8
0020     0.00160 ... % 16
0021     0.0060 ... % 32
0022     0.047 ... % 64
0023     0.32 ... % 128
0024     2.5 ... % 256
0025     20.5 ... % 512
0026     195 ... % 1024
0027     1750 ... % 2048
0028     15000 ... % 4096
0029     ];
0030 accuracy_multires = [ ...
0031     2.5e-15 ... % 4
0032     3.91e-15 ... % 8
0033     8.0e-15 ... % 16
0034     2.1e-14 ... % 32
0035     3.97e-14 ... % 64
0036     7.32e-14 ... % 128
0037     1.37e-13 ... % 256
0038     2.37e-13 ... % 512
0039     5.3e-13 ... % 1024
0040     9.3e-13 ... % 2048
0041     2.0e-12 ... % 4096
0042     ];
0043 speed_multires = [ ...
0044     0.00021 ...  % 4
0045     0.00040 ... % 8
0046     0.00105 ... % 16
0047     0.00300 ... % 32
0048     0.018 ... % 64
0049     0.1 ... % 128
0050     0.75 ... % 256
0051     5.2 ... % 512
0052     45.0 ... % 1024
0053     400 ... % 2048
0054     3000 ... % 4096
0055     ];
0056 nb_samples = ell;
0057 nb_samples_pow = str2mat('t04', 't08', 't16', 't32', 't64', 't28', 't56', 't12', 't24', 't48', 't96');
0058 
0059 figure('Position',[100 100 700 600])
0060 
0061 subplot(2,1,1)
0062 loglog( nb_samples, (1e-15)*ell, 'red','LineWidth', 2 )
0063 hold on
0064 loglog( nb_samples, accuracy, '--oblack', 'LineWidth', 2, 'MarkerSize',10)
0065 loglog( nb_samples, accuracy_multires, '-oblack', 'LineWidth', 2, 'MarkerSize',10)
0066 grid on
0067 xlabel('x1','FontSize',20)
0068 ylabel('y1','FontSize',20)
0069 set(gca,'XTick',nb_samples,'XTickLabel',nb_samples_pow,'FontSize',20)
0070 set(gca,'YTick',[1e-16 1e-15 1e-14 1e-13 1e-12 1e-11] ,'FontSize',20)
0071 axis([3e0 5300 1e-15 1e-11])
0072 grid minor
0073 %title('Accuracy of overall transform','FontSize',20)
0074 
0075 subplot(2,1,2)
0076 loglog( nb_samples, (1e-6)*ell.^3, 'red','LineWidth', 2 )
0077 hold on
0078 loglog( nb_samples, speed, '--oblack', 'LineWidth', 2, 'MarkerSize',10)
0079 loglog( nb_samples, speed_multires, '-oblack', 'LineWidth', 2, 'MarkerSize',10)
0080 grid on
0081 xlabel('x2','FontSize',20)
0082 ylabel('y2','FontSize',20)
0083 axis([3e0 5300 1e-5 1e5])
0084 set(gca,'YTick',[1e-6 1e-4 1e-2  1e0 1e2 1e4 1e6] ,'FontSize',20)
0085 set(gca,'XTick',nb_samples,'XTickLabel',nb_samples_pow,'FontSize',20)
0086 grid minor
0087 %title('Speed of overall transform','FontSize',20)
0088