0001
0002
0003
0004
0005
0006
0007
0008 B = 3;
0009 zoomfactor = 1.4;
0010 J_min = 2;
0011 L = 128;
0012 J = s2let_jmax(L, B);
0013 type = 'colour';
0014 lighting = true;
0015
0016 ns = ceil(sqrt(2+J-J_min+1)) ;
0017 nx = ns - 1 + rem(2+J-J_min + 1, ns) ;
0018 ny = ns;
0019
0020 [kappa kappa0] = s2let_transform_axisym_tiling(B, L, J_min);
0021
0022 s2let_plot_axisym_tiling(B, L, J_min);
0023
0024 figure('Position',[100 100 1100 700])
0025
0026 h=subplot(nx, ny, 1);
0027 flm = zeros(L^2,1);
0028 for l = 0:L-1
0029 flm(l^2+l+1,1) = kappa0(l+1);
0030 end
0031 f = ssht_inverse(flm, L, 'Reality', true);
0032 ssht_plot_sphere(f, L, 'Type', type, 'Lighting', lighting);
0033 zoom(zoomfactor)
0034 v = caxis;
0035 temp = max(abs(v));
0036 caxis([-temp temp])
0037
0038 locate = get(h,'title');
0039 pos = get(locate,'position');
0040 pos(1,2) = pos(1,2)+0.7;
0041 pos(1,1) = pos(1,1)-0.7;
0042 set(locate,'pos',pos);
0043
0044 for j = J_min:5
0045 h=subplot(nx, ny, j-J_min+2);
0046 flm = zeros(L^2,1);
0047 for l = 0:L-1
0048 flm(l^2+l+1,1) = kappa(j+1,l+1);
0049 end
0050 f = ssht_inverse(flm, L, 'Reality', true);
0051 ssht_plot_sphere(f, L, 'Type', type, 'Lighting', lighting);
0052 v = caxis;
0053 temp = max(abs(v));
0054 caxis([-temp temp])
0055 colormap(jet)
0056
0057 locate = get(h,'title');
0058 pos = get(locate,'position');
0059 pos(1,2) = pos(1,2)+0.7;
0060 pos(1,1) = pos(1,1)-0.7;
0061 set(locate,'pos',pos);
0062 zoom(zoomfactor)
0063 end
0064
0065
0066