[184] | 1 | itload('../merger_debug.it'); |
---|
| 2 | |
---|
| 3 | ndat = 1000; % check if true! |
---|
| 4 | |
---|
| 5 | figure(1); |
---|
| 6 | for it=0:6 |
---|
| 7 | figure(it+1) |
---|
| 8 | subplot(2,2,1); |
---|
| 9 | hold off |
---|
| 10 | si = num2str(it); |
---|
| 11 | eval(['plot3(Smp' si '(1,:),Smp' si '(2,:),exp(Mpdf' si '),''r.'')']); |
---|
| 12 | if i==0, title('Proposal & mix'); end |
---|
| 13 | set(gca,'XLim',[-10,10]); |
---|
| 14 | |
---|
| 15 | subplot(2,2,2); |
---|
| 16 | hold off |
---|
| 17 | eval(['contour_2(Smp' si '(1,:),Smp' si '(2,:),exp(lW' si '(1,:)))']); |
---|
| 18 | % hold on |
---|
| 19 | % eval(['plot3(Smp' si '(1,:),Smp' si '(2,:),w' si ',''r.'')']); |
---|
| 20 | if i==0, title('Weighed Smp'); end |
---|
| 21 | |
---|
| 22 | subplot(2,2,3); |
---|
| 23 | hold off |
---|
| 24 | eval(['contour_2(Smp' si '(1,:),Smp' si '(2,:),exp(lW' si '(2,:)))']); |
---|
| 25 | % eval(['plot3(Smp' si '(1,:),Smp' si '(2,:),w' si ',''.'')']); |
---|
| 26 | % hold on |
---|
| 27 | % eval(['plot3(Smp' si '(1,:),Smp' si '(2,:),w_is_' si ',''r.'')']); |
---|
| 28 | if i==0, title('Weighed Smp'); end |
---|
| 29 | |
---|
| 30 | |
---|
| 31 | subplot(2,2,4); |
---|
| 32 | eval(['contour_2(Smp' si '(1,:),Smp' si '(2,:),w' si ')']); |
---|
| 33 | if i==0, title('Resampled hist'); end |
---|
| 34 | |
---|
| 35 | % keyboard |
---|
| 36 | end |
---|
| 37 | |
---|
| 38 | |
---|
| 39 | itload('../merger_iter_test.it'); |
---|
| 40 | XG = reshape(Grid(1,:),Npoints,Npoints); |
---|
| 41 | YG = reshape(Grid(2,:),Npoints,Npoints); |
---|
| 42 | |
---|
| 43 | figure(it+2); |
---|
| 44 | M1 = reshape(exp(Res1),Npoints,Npoints); |
---|
| 45 | contour(XG,YG,M1); |
---|
| 46 | |
---|
| 47 | figure(it+3); |
---|
| 48 | M2 = reshape(exp(Res2),Npoints,Npoints); |
---|
| 49 | contour(XG,YG,M2); |
---|
| 50 | |
---|
| 51 | figure(it+4); |
---|
| 52 | L1 = reshape(exp(lf1),Npoints,Npoints); |
---|
| 53 | contour(XG,YG,L1); |
---|