From 8c3c49fd56de185396400a0e8a636ac03d8239ab Mon Sep 17 00:00:00 2001 From: Norbert Preining Date: Mon, 12 Feb 2018 14:52:54 +0900 Subject: update ketcindy --- .../ketlib/ketpicscifiles6/Phcutoffdata.sci | 184 --------------------- 1 file changed, 184 deletions(-) delete mode 100644 texmf-dist/scripts/ketcindy/ketlib/ketpicscifiles6/Phcutoffdata.sci (limited to 'texmf-dist/scripts/ketcindy/ketlib/ketpicscifiles6/Phcutoffdata.sci') diff --git a/texmf-dist/scripts/ketcindy/ketlib/ketpicscifiles6/Phcutoffdata.sci b/texmf-dist/scripts/ketcindy/ketlib/ketpicscifiles6/Phcutoffdata.sci deleted file mode 100644 index 879067ce..00000000 --- a/texmf-dist/scripts/ketcindy/ketlib/ketpicscifiles6/Phcutoffdata.sci +++ /dev/null @@ -1,184 +0,0 @@ -// 08.08.17 -// 10.08.15 near line 25 or so -// 13.10.21 ( __ added to varibles ) - -function Out2__=Phcutoffdata(varargin) - global PHCUTPOINTL PHVERTEXL PHFACEL; - Eps__=10^(-4); - Out2__=[]; - VL__=varargin(1); - FaceL__=varargin(2); - PlaneD__=varargin(3); - Sgnstr__=varargin(4); - Fugou__=1; - if Sgnstr__=='-' | Sgnstr__=='n' Fugou__=-1; end; - Out0__=Phcutdata(VL__,FaceL__,PlaneD__); - PtL__=PHCUTPOINTL; - PHVERTEXL=VL__; - PHFACEL=FaceL__; - if PtL__==[] - Out2__=Out0__; - return; - end; - N1__=Mixlength(VL__); - Face__=[]; - for I__=1:Mixlength(PtL__); - Tmp__=Mixop(I__,PtL__); - Tmp3__=Mixop(1,Tmp__); // 2010.08.15 - Tmp1__=Tmp3__(1); Tmp2__=Tmp3__(2); // 2010.08.15 - if Tmp1__~=Tmp2__ - VL__=Mixadd(VL__,Mixop(2,Tmp__)); - N1__=N1__+1; - Tmp1__=N1__; - end; - Face__=[Face__,Tmp1__]; - PtL__(I__)=list(Mixop(1,Tmp__),Tmp1__); // 2010.08.15 -// Tmp2__=Mixsub(1:I__-1,PtL__); -// Tmp2__=Mixadd(Tmp2__,Mix(Mixop(1,Tmp__),Tmp1__)); -// Tmp3__=Mixsub(I__+1:Mixlength(PtL__),PtL__); -// PtL__=Mixjoin(Tmp2__,Tmp3__); - end; - OutfL__=MixS(Face__); - if Mixtype(PlaneD__)~=1 - V1__=Mixop(1,PlaneD__); - Tmp__=Mixop(2,PlaneD__); - if length(Tmp__)>1 - d__=V1__(1)*Tmp__(1)+V1__(2)*Tmp__(2)+V1__(3)*Tmp__(3); - else - d__=Tmp__; - end; - elseif type(PlaneD__)==1 - V1__=PlaneD__(1:3); - d__=PlaneD__(4); - else - K__=mtlb_findstr(PlaneD__,'='); - if K__>0 - Tmp1__=part(PlaneD__,1:K__-1); - Tmp2__=part(PlaneD__,K__+1:length(PlaneD__)); - PlaneD__=Tmp1__+'-('+Tmp2__+')'; - end; - x=0; y=0; z=0; - d__=-evstr(PlaneD__); - x=1; y=0; z=0; - Tmp1__=evstr(PlaneD__)+d__; - x=0; y=1; z=0; - Tmp2__=evstr(PlaneD__)+d__; - x=0; y=0; z=1; - Tmp3__=evstr(PlaneD__)+d__; - V1__=[Tmp1__,Tmp2__,Tmp3__]; - end; - for I__=1:Mixlength(FaceL__) - Face__=Mixop(I__,FaceL__); - TmpL__=[]; - for J__=1:length(Face__) - N1__=Face__(J__); - if J__==length(Face__) - N2__=Face__(1); - else - N2__=Face__(J__+1); - end; - for K__=1:Mixlength(PtL__) - Pd__=Mixop(K__,PtL__); - Tmp__=Mixop(1,Pd__); - if Tmp__(1)==Tmp__(2) - if Tmp__(1)==N1__ - TmpL__=Mixadd(TmpL__,Mix(J__,[N1__,N2__],Mixop(2,Pd__))); - end - else - if Tmp__==[N1__,N2__] | Tmp__==[N2__,N1__] - TmpL__=Mixadd(TmpL__,Mix(J__,[N1__,N2__],Mixop(2,Pd__))); - end; - end; - end; - end; - if Mixlength(TmpL__)<2 - Flg__=0; - for J__=1:length(Face__) - Tmp__=Mixop(Face__(J__),VL__); - Tmp1__=Fugou__*(Naiseki(V1__,Tmp__)-d__); - if Tmp1__<-Eps__ - Flg__=1; - break; - end; - end; - if Flg__==0 - OutfL__=Mixadd(OutfL__,Face__); - end; - continue; - end; - Pd__=Mixop(1,TmpL__); - Qd__=Mixop(2,TmpL__); - Outf1__=[Mixop(3,Pd__)]; - Nf__=Mixop(1,Pd__)+1; - Tmp__=Mixop(2,Pd__); - JJ__=0; - while Tmp__~=Mixop(2,Qd__) - JJ__=JJ__+1; - if JJ__>20 - disp('bug'); - return; - end; - Tmp1__=Tmp__(2); - if Outf1__(length(Outf1__))~=Tmp1__ - Outf1__=[Outf1__,Tmp1__]; - end; - Tmp__=[Face__(Nf__)]; - Nf__=Nf__+1; - if Nf__>length(Face__) Nf__=1; end; - Tmp__=[Tmp__,Face__(Nf__)]; - end; - Tmp1__=Mixop(3,Qd__); - if Outf1__(length(Outf1__))~=Tmp1__ - Outf1__=[Outf1__,Tmp1__]; - end; - Outf2__=[Mixop(3,Pd__)]; - Nf__=Mixop(1,Pd__); - Tmp__=Mixop(2,Pd__); - JJ__=0; - while Tmp__~=Mixop(2,Qd__) - JJ__=JJ__+1; - if JJ__>20 - disp('bug'); - return; - end; - Tmp1__=Tmp__(1); - if Outf2__(length(Outf2__))~=Tmp1__ - Outf2__=[Outf2__,Tmp1__] - end; - Tmp__=[Face__(Nf__)]; - Nf__=Nf__-1; - if Nf__<1 Nf__=length(Face__); end; - Tmp__=[Face__(Nf__),Tmp__]; - end; - Tmp1__=Mixop(3,Qd__); - if Outf2__(length(Outf2__))~=Tmp1__ - Outf2__=[Outf2__,Tmp1__] - end; - if length(Outf1__)<3 | length(Outf2__)<3 - Face__=Outf1__; - if length(Outf1__)