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
|
// 08.06.05
// 16.12.29 ( for too short segment )
function Out=Anglemark(varargin)
Nargs=length(varargin);
Eps=10^(-3);
PA=varargin(1); PB=varargin(2); PC=varargin(3);
r = 0.5;
if Nargs>=4
r = varargin(4)*r;
end
Out=[]; // 16.12.29
if r>min(norm(PA-PB),norm(PC-PB)) then
return;
end;
Cir=Circledata(PB,r);
Tmp=IntersectcrvsPp(Cir,Listplot(PA,PB));
P1=Op(2,Op(1,Tmp));
Tmp=IntersectcrvsPp(Cir,Listplot(PC,PB));
P2=Op(2,Mixop(1,Tmp));
if abs(P1-P2)<Eps
Out=[];
return;
end
if P1<P2-Eps
Out=Partcrv(P1,P2,Cir);
else
Tmp=Numptcrv(Cir);
Tmp1=Partcrv(P1,Tmp,Cir);
Tmp2=Partcrv(1,P2,Cir);
Out=Joincrvs(Tmp1,Tmp2);
end
endfunction
|