summaryrefslogtreecommitdiff
path: root/texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Cuspsplitpers.sci
diff options
context:
space:
mode:
Diffstat (limited to 'texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Cuspsplitpers.sci')
-rw-r--r--texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Cuspsplitpers.sci176
1 files changed, 0 insertions, 176 deletions
diff --git a/texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Cuspsplitpers.sci b/texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Cuspsplitpers.sci
deleted file mode 100644
index a8221f5b..00000000
--- a/texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Cuspsplitpers.sci
+++ /dev/null
@@ -1,176 +0,0 @@
-// 08.09.21
-// 08.10.10
-// 10.02.16 Eps__
-// 13.11.01 ( __ added to varibles )
-// 14.06.28 debug in the case of no cusp
-
-function OutkL__=Cuspsplitpers(varargin);
- global CUSPSPLITPT;
- Nargs__=length(varargin);
- Gdxy__=varargin(1);
- if Mixtype(Gdxy__)==1
- Gdxy__=Mix(Gdxy__);
- end;
- Fd__=varargin(2);
- FdL__=Fullformfunc(Fd__);
- Eps__0=10^(-5);
- Eps__=0.05;
- if Nargs__>2 Eps__=varargin(3); end;
- if length(Eps__)>1 //
- Eps__=Eps__(1);
- end; //
- Fxy__=Mixop(1,FdL__); N=2;
- Tmp__1=Mixop(2,FdL__);
- Tmp__2=Mixop(3,FdL__);
- Tmp__3=Mixop(4,FdL__);
- Xyzstr__=[Tmp__1,Tmp__2,Tmp__3];
- Tmp__=Mixop(5,FdL__);
- K__=mtlb_findstr(Tmp__,'=');
- Uname__=part(Tmp__,1:K__-1);
- Urg__=evstr(part(Tmp__,K__+1:length(Tmp__)));
- Tmp__=Mixop(6,FdL__);
- K__=mtlb_findstr(Tmp__,'=');
- Vname__=part(Tmp__,1:K__-1);
- Vrg__=evstr(part(Tmp__,K__+1:length(Tmp__)));
- CUSPSPLITPT=[];
- OutkL__=[];
- for Ng__=1:Mixlength(Gdxy__)
- PtxyL__=Mixop(Ng__,Gdxy__);
- PtkL__=[];
- PthL__=[];
- for I__=1:size(PtxyL__,1)
- Tmp__=PtxyL__(I__,:);
- Tmp__1=Uname__+'='+string(Tmp__(1));
- execstr(Tmp__1);
- Tmp__1=Vname__+'='+string(Tmp__(2))
- execstr(Tmp__1);
- Tmp__2=evstr(Xyzstr__);
- Tmp__3=Perspt(Tmp__2);
- if I__==1
- PtkL__=[Tmp__2];
- PthL__=[Tmp__3];
- else
- Tmp__4=PthL__(size(PthL__,1),:);
- if norm(Tmp__3-Tmp__4)>Eps__0
- PtkL__=[PtkL__;Tmp__2];
- PthL__=[PthL__;Tmp__3]
- end;
- end;
- end;
- if size(PthL__,1)==0 then // from 20140628
- OutkL__=[];
- return;
- end; // upto 20140628
- Ps__=PthL__(1,:); Pe=PthL__(size(PthL__,1),:);
- Cflg__=0;
- if norm(Ps__-Pe)<Eps__ Cflg__=1; end;
- CuspL__=[];
- Cva__=cos(10*%pi/180);
- for I__=1:size(PthL__,1)-1
- if size(PthL__,1)==2 break; end;
- Tmp__=PthL__(I__,:);
- if I__==1
- if Cflg__==0 continue; end;
- Tmp__1=PthL__(size(PthL__,1)-1,:);
- else
- Tmp__1=PthL__(I__-1,:);
- end;
- Tmp__2=PthL__(I__+1,:);
- V1__=Tmp__-Tmp__1;
- V2__=Tmp__2-Tmp__;
- Tmp__3=V1__(1)*V2__(1)+V1__(2)*V2__(2);
- Tmp__3=Tmp__3/(norm(V1__)*norm(V2__));
- Cuspflg__=0;
- if Tmp__3<Cva__
- P__=PthL__(I__,:);
- Kaku__=acos(Tmp__3)*180/%pi;
- if V1__(1)*V2__(2)-V1__(2)*V2__(1)<0
- Kaku__=-Kaku__
- end;
- Cuspflg__=0;
- for J__=I__+1:size(PthL__,1)-1
- if abs(Kaku__)>90
- Cuspflg__=1;
- break
- end;
- Q__=PthL__(J__,:);
- if norm(P__-Q__)>Eps__
- break
- end;
- V1__=Q__-PthL__(J__-1,:);
- V2__=PthL__(J__+1,:)-Q__;
- Tmp__3=V1__(1)*V2__(1)+V1__(2)*V2__(2);
- Tmp__3=Tmp__3/(norm(V1__)*norm(V2__));
- Tmp__=acos(Tmp__3)*180/%pi;
- if V1__(1)*V2__(2)-V1__(2)*V2__(1)<0
- Tmp__=-Tmp__
- end;
- Kaku__=Kaku__+Tmp__
- end;
- if Cuspflg__==1
- Tmp__=Trunc((I__+J__)*0.5);
- I__=J__;
- if length(CuspL__)==0
- CuspL__=[Tmp__]
- else
- CuspL__=[CuspL__,Tmp__];
- end;
- end;
- end;
- end;
- if Cflg__==0
- CuspL__=[1,CuspL__,size(PthL__,1)];
- elseif length(CuspL__)==0
- CuspL__=[1,size(PthL__,1)];
- elseif CuspL__(1)==1
- CuspL__=[CuspL__,size(PthL__,1)];
- else
- Tmp__=CuspL__(1);
- Tmp__1=size(PthL__,1);
- PthL__=[PthL__(Tmp__:Tmp__1,:);PthL__(2:Tmp__,:)];
- PtkL__=[PtkL__(Tmp__:Tmp__1,:);PtkL__(2:Tmp__,:)];
- CuspL__=CuspL__-Tmp__+1;
- CuspL__=[CuspL__,size(PthL__,1)];
- end;
- if length(CuspL__)==2
- //Outh__=[plot(PthL__)];
- Tmp__4=PthL__(size(PthL__,1),:);
- if size(PtkL__,1)>=2
- CUSPSPLITPT=Mixadd(CUSPSPLITPT,Tmp__4);
- OutkL__=Mixadd(OutkL__,PtkL__);
- end;
- continue
- end;
- //Outh__=[];
- Outk__=[];
- Is__=1;
- for I__=1:length(CuspL__)-1
- Tmp__1=CuspL__(Is__); Tmp__2=CuspL__(I__+1);
- Tmp__3=PthL__(Tmp__1,:); Tmp__4=PthL__(Tmp__2,:);
- if norm(Tmp__3-Tmp__4)>Eps__
- //Tmp__h=[op(Tmp__1..Tmp__2,PthL__)];
- Tmp__k=PtkL__(Tmp__1:Tmp__2,:);
- //Outh__=[op(Outh__),plot(Tmp__h)];
- Outk__=Mixadd(Outk__,Tmp__k);
- CUSPSPLITPT=Mixadd(CUSPSPLITPT,Tmp__4);
- Is__=I__+1;
- end;
- end;
- OutkL__=Mixjoin(OutkL__,Outk__);
- end;
- Tmp__1=Dropnumlistcrv(Projpers(OutkL__),Eps__*0.5);
- Tmp__=[];
- for I__=1:Mixlength(OutkL__)
- Tmp__2=Mixop(I__,OutkL__);
- Tmp__3=Mixop(I__,Tmp__1);
- Tmp__4=[];
- for J__=1:length(Tmp__3)
- Tmp__5=Tmp__2(Tmp__3(J__),:);
- Tmp__4=[Tmp__4;Tmp__5];
- end;
- if Tmp__4~=[]
- Tmp__=Mixadd(Tmp__,Tmp__4);
- end;
- end;
- OutkL__=Tmp__;
-endfunction;