summaryrefslogtreecommitdiff
path: root/texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Intersectcrvsf.sci
diff options
context:
space:
mode:
authorNorbert Preining <norbert@preining.info>2017-10-14 11:43:11 +0900
committerNorbert Preining <norbert@preining.info>2017-10-14 11:43:11 +0900
commitdc71274f16f7ff65145e034900f582b62453a033 (patch)
treea71452364cbe1a03feb6a84cdfd09eedae861590 /texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Intersectcrvsf.sci
parentaecda53b09be64698e9e5dcbf40f0fe289b37e2e (diff)
rename and update ketpic -> ketcindy
Diffstat (limited to 'texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Intersectcrvsf.sci')
-rw-r--r--texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Intersectcrvsf.sci47
1 files changed, 47 insertions, 0 deletions
diff --git a/texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Intersectcrvsf.sci b/texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Intersectcrvsf.sci
new file mode 100644
index 00000000..2a2e5b6c
--- /dev/null
+++ b/texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Intersectcrvsf.sci
@@ -0,0 +1,47 @@
+// 08.10.26
+// 10.02.16
+
+function Out=Intersectcrvsf(varargin)
+ Nargs=length(varargin);
+ Eps0=10^(-4);
+ Crv=varargin(1);
+ Fd=varargin(2);
+ N=3;Uveq=1;
+ if Nargs>=3
+ Tmp=varargin(3);
+ if type(Tmp)==10
+ Uveq=Tmp;
+ N=4;
+ end;
+ end;
+ Np=[50,50];
+ if Nargs>=N
+ Np=varargin(N);
+ if length(Np)==1
+ Np=[Np,Np];
+ end;
+ end;
+ Eps=[0.05,1]; //
+ if Nargs>=N+1
+ Eps=varargin(N+1);
+ end;
+ if length(Eps)==1 //
+ Eps=[Eps,1]; //
+ end; //
+ Out=[];
+ for I=1:Numptcrv(Crv)-1
+ PtA=Ptcrv(I,Crv);
+ PtB=Ptcrv(I+1,Crv);
+ if norm(PtB-PtA)<Eps0
+ continue;
+ end;
+ PL=Meetpoints(PtA,PtB,Fd,Uveq,Np,Eps(1));
+ for J=1:Mixlength(PL)
+ Pt=Mixop(J,PL);
+ Tmp=norm(Pt-PtA)/norm(PtB-PtA);
+ Tmp=min(max(Tmp,0),1);
+ Out=Mixadd(Out,[Pt,I+Tmp]);
+ end;
+ end;
+endfunction;
+