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
|
function hist3d(f,T,A,leg,flags,ebox)
//!
// Copyright INRIA
nep=8
def=list(35,45,'X@Y@Z',[2 1 4])
[lhs,rhs]=argn(0)
if rhs<=0 then //demo
s_mat=['hist3d(10*rand(10,10));';]
write(%io(2),s_mat);execstr(s_mat);
return;
end
if typeof(f)=='list' then
[f,x,y]=f(1:3);
sx=prod(size(x));
sy=prod(size(y));
if [sx-1,sy-1]<>size(f) then
write(%io(2),'f: Incompatible length ');
return;
end
dx=(x(2)-x(1))/nep;
dy=(y(2)-y(1))/nep;
bnds=[x(1) x($), y(1) y($),mini(0,mini(f)) maxi(f)]
x=(x(1:sx-1)+x(2:sx))/2;
y=(y(1:sy-1)+y(2:sy))/2;
[nl,nc]=size(f);
else
[nl,nc]=size(f);
x=(1:nl)-(0.5)*ones(1,nl);
y=(1:nc)-(0.5)*ones(1,nc);
dx=1/nep; dy=1/nep;
bnds=[0 nl,0 nc,mini(0,mini(f)) maxi(f)]
end
x=x.*.[1,1] + dx*ones(x).*.[0,1] - dx*ones(x).*.[1,0];
y=y.*.[1,1] + dy*ones(y).*.[0,1] - dy*ones(y).*.[1,0];
a=[0;0;1;1]
b=[0;1;1;0]
c=[0;0;0;0]
d=[1;1;1;1]
ix=[b,b,a,a,c,d];
iy=[a,a,c,d,b,b];
indx=ones(1,nc) .*. (ones(1,nl).*.ix +(1:2:2*nl-1).*.ones(ix));
iy=matrix(iy,24,1);
//indy=(ones(1,nl).*.iy+(1:2:2*nl-1).*.ones(iy)) .*. ones(1,nc);
indy=(ones(1,nc).*.iy+(1:2:2*nc-1).*.ones(iy)) .*. ones(1,nl);
indy=matrix(indy,4,6*nc*nl);
[nnl,nnc]=size(indx);
xx=matrix(x(matrix(indx,1,nnl*nnc)),nnl,nnc);
yy=matrix(y(matrix(indy,1,nnl*nnc)),nnl,nnc);
zz=matrix(f,1,nl*nc).*.[c,d,b,b,a,a];
select rhs
case 1 then plot3d(xx,yy,zz,def(1),def(2),def(3),def(4),bnds)
case 2 then plot3d(xx,yy,zz,T,def(2),def(3),def(4),bnds)
case 3 then plot3d(xx,yy,zz,T,A,def(3),def(4),bnds)
case 4 then plot3d(xx,yy,zz,T,A,leg,def(4),bnds)
case 5 then plot3d(xx,yy,zz,T,A,leg,flags,bnds)
case 6 then plot3d(xx,yy,zz,T,A,leg,flags,ebox)
end
|
