summaryrefslogtreecommitdiff
path: root/texmf-dist/scripts/ketpic/ketlib/ketpicscifiles6/Anglemark.sci
blob: 6deb302c07d4aaa59b4e0b6c05c2401bf04902f5 (plain)
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