From dc71274f16f7ff65145e034900f582b62453a033 Mon Sep 17 00:00:00 2001 From: Norbert Preining Date: Sat, 14 Oct 2017 11:43:11 +0900 Subject: rename and update ketpic -> ketcindy --- .../ketlib/ketpicscifiles6/Wirerawparadata.sci | 86 ++++++++++++++++++++++ 1 file changed, 86 insertions(+) create mode 100644 texmf-dist/scripts/ketcindy/ketlib/ketpicscifiles6/Wirerawparadata.sci (limited to 'texmf-dist/scripts/ketcindy/ketlib/ketpicscifiles6/Wirerawparadata.sci') diff --git a/texmf-dist/scripts/ketcindy/ketlib/ketpicscifiles6/Wirerawparadata.sci b/texmf-dist/scripts/ketcindy/ketlib/ketpicscifiles6/Wirerawparadata.sci new file mode 100644 index 00000000..59c5ba54 --- /dev/null +++ b/texmf-dist/scripts/ketcindy/ketlib/ketpicscifiles6/Wirerawparadata.sci @@ -0,0 +1,86 @@ +// 08.09.15 +// 08.10.10 +// 08.10.21 +// 09.09.10 + +function FsL=Wirerawparadata(varargin) + timer(); + Nargs=length(varargin); + Fbdy=Projpara(varargin(1)); + Fd=varargin(2); N=3; + FdL=Fullformfunc(Fd); + Fxy=Mixop(1,FdL); + 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))); + DuL=varargin(N); + DvL=varargin(N+1); + N=N+2; + Np=[50,50]; + if Nargs>=N + Np=varargin(N); + if type(Np)==1 & length(Np)==1 + Np=[Np,Np]; + end; + end + Eps=0.05; + if Nargs>=N+1 + Eps=varargin(N+1); + end; + Eps2=0.2; + if Nargs>=N+2 + Eps2=varargin(N+2); + end; + Umin=Urg(1); Umax=Urg(2); + Vmin=Vrg(1); Vmax=Vrg(2); + if type(DuL)==10 + DuL=evstr(DuL); + elseif length(DuL)==1 + Tmp=(Umax-Umin)/(DuL+1); + Tmp1=Umin+Tmp; + Tmp2=Umin+DuL*Tmp + DuL=[Tmp1:Tmp:Tmp2]; + end; + if type(DvL)==10 + DvL=evstr(DvL); + elseif length(DvL)==1 + Tmp=(Vmax-Vmin)/(DvL+1); + Tmp1=Vmin+Tmp; + Tmp2=Vmin+DvL*Tmp + DvL=[Tmp1:Tmp:Tmp2]; + end; + Fbdxy=Makexybdy(Fd,Np); + FsL=[]; + for I=1:length(DuL) + U=DuL(I); + Tmp=Paramplot('[U,V]','V=[Vmin,Vmax]','N=Np(1)'); + FkL=Cuspsplitpara(Tmp,Fd,Eps); + for J=1:Mixlength(FkL) + Fk=Mixop(J,FkL); + FsL=Mixadd(FsL,Fk); + end; + end; + Nx=Mixlength(FsL); + if Nx>0 + Mixdisp(Mix(Uname,'-data=',1,':',Nx,' Time=',timer())); + end; + for I=1:length(DvL) + V=DvL(I); + Tmp=Paramplot('[U,V]','U=[Umin,Umax]','N=Np(2)'); + FkL=Cuspsplitpara(Tmp,Fd,Eps); + for J=1:Mixlength(FkL) + Fk=Mixop(J,FkL); + FsL=Mixadd(FsL,Fk); + end; + end; + Ny=Mixlength(FsL)-Nx; + if Ny>0 + Mixdisp(Mix(Vname,'-data=',Nx+1,':',Nx+Ny,' Time=',timer())); + end; +endfunction; + -- cgit v1.2.3