-
Notifications
You must be signed in to change notification settings - Fork 4
/
sdfried.m
79 lines (72 loc) · 2.08 KB
/
sdfried.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
function sdfried
% SDFRIED
%
% Makes fried-egg plots for +- angular orders of the DOUBLE spherical
% polar cap.
%
% Second time round the color axis is OK.; I forgot why.
%
% Last modified by fjsimons-at-alum.mit.edu, 25.05.2005
TH=40;
L=18;
[lrnk,mrnk,lval]=sdelm(TH,L);
% Nearly double the amount of requested tapers
lrnk=gamini(lrnk,(mrnk~=0)+1);
lval=gamini(lval,(mrnk~=0)+1);
mrnk=gamini(mrnk,(mrnk~=0)+1);
% This is better than what we had before, in sdfried
% Because, e.g. the first and second eigenfunction of the same order can
% come right next to one another in the rank ordering.
mrnk((diff(mrnk)==0) & (diff(lrnk)==0))= ...
-mrnk((diff(mrnk)==0)& (diff(lrnk)==0));
[ah,ha]=krijetem(subnum(8,4));
for index=1:length(ah)
fnpl=sprintf('%s/SDFRIED-%i-%i-%i.mat',...
fullfile(getenv('IFILES'),'WIECZOREK'),TH,L,index);
if exist(fnpl,'file')==2
eval(sprintf('load %s',fnpl))
disp(sprintf('%s loaded by SDFRIED',fnpl))
axes(ah(index))
cbb=gray(10);
cbb=kelicol;
ispl=data/max(abs(data(:)));
ispl(abs(ispl)<1/100)=NaN;
imagefnan([-1 1],[1 -1],ispl,cbb,[-1 1])
colormap(cbb); hold on
ax2(1)=circ(1); ax2(2)=circ(sin(TH*pi/180));
set(ax2(1:2),'LineW',1)
set(ax2(2),'LineS','--')
axis([-1.0100 1.0100 -1.0100 1.0100])
axis off
if index==0
[axlim,handl,XYZ]=plotcont([0 90],[360 -90],4);
delete(handl)
hold on
p1=plot(XYZ(:,1),XYZ(:,2),'k-','LineWidth',1);
end
tl(index)=title(sprintf('%s_%i = %6.3f ; m = %i',...
'\lambda',...
lrnk(index),lval(index),mrnk(index)));
drawnow
else
[E,V,N,th,C]=grunbaum2(TH,L,mrnk(index),128,2);
E=E{lrnk(index)};
Emax=max(E(:)); thresh=100;
Ed=E;
Ed(abs(Ed)<Emax/thresh)=0;
axes(ah(index))
[data,ch,ph]=plotplm(Ed,[],[],5,[],TH);
tl(index)=...
title(sprintf('%s = %6.3f ; m = %i','\lambda',...
lval(index),mrnk(index)));
delete(ch)
drawnow
eval(sprintf('save %s data',fnpl))
end
end
seemax(ah,3)
fig2print(gcf,'tall')
set(ah,'CameraV',6.25)
figdisp
% Last-minute cosmetics
movev(tl,-0.25)