1 files changed, 75 insertions, 0 deletions
diff --git a/texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/Rotate3pt.sci b/texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/Rotate3pt.sci
new file mode 100755
index 00000000..910aa734
--- /dev/null
+++ b/texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/Rotate3pt.sci
@@ -0,0 +1,75 @@
+// 08.08.15
+
+function Ans=Rotate3pt(varargin)
+  Eps=10^(-4);
+  Nargs=length(varargin);
+  P=varargin(1);
+  W1=varargin(2);
+  W2=varargin(3);
+  C=[0,0,0];
+  if Nargs>=4
+    C=varargin(4);
+  end;
+  if type(W2)==1 & length(W2)==1
+    Ct=cos(W2);
+    St=sin(W2);
+    V3=1/norm(W1)*W1;
+    if V3(1)==0
+      Tmp=[1,0,0];
+    else
+      Tmp=[0,1,0];
+    end;
+    W1=[Tmp(2)*V3(3)-Tmp(3)*V3(2),Tmp(3)*V3(1)-Tmp(1)*V3(3),Tmp(1)*V3(2)-Tmp(2)*V3(1)];
+    V1=1/norm(W1)*W1;
+    V2=[V3(2)*V1(3)-V3(3)*V1(2),V3(3)*V1(1)-V3(1)*V1(3),V3(1)*V1(2)-V3(2)*V1(1)];
+  else
+    Tmp=[W1(2)*W2(3)-W1(3)*W2(2),W1(3)*W2(1)-W1(1)*W2(3),W1(1)*W2(2)-W1(2)*W2(1)];
+    if norm(Tmp)<Eps
+      Ans=P;
+      return;
+    end;
+    V1=1/norm(W1)*W1;
+    Ns=V1(1)*W2(1)+V1(2)*W2(2)+V1(3)*W2(3);
+    Tmp=W2-Ns*V1;
+    V2=1/norm(Tmp)*Tmp;
+    Tmp=[V1(2)*V2(3)-V1(3)*V2(2),V1(3)*V2(1)-V1(1)*V2(3),V1(1)*V2(2)-V1(2)*V2(1)];
+    V3=1/norm(Tmp)*Tmp;
+    Ct=Ns/norm(W2);
+    St=sqrt(1-Ct^2);
+  end;
+  if norm(Tmp)<Eps
+    Ans=P;
+    return;
+  end;
+  V1x=V1(1); V1y=V1(2); V1z=V1(3);
+  V2x=V2(1); V2y=V2(2); V2z=V2(3);
+  V3x=V3(1); V3y=V3(2); V3z=V3(3);
+  if Mixtype(P)~=1
+    PtL=P;
+  else
+    PtL=MixS(P);
+  end;
+  Ans=[];
+  for N=1:Mixlength(PtL)
+    P=Mixop(N,PtL);
+    if P(1)==%inf
+      Ans=Mixadd(Ans,[%inf,%inf,%inf]);
+      continue;
+    end;
+    x=P(1)-C(1); y=P(2)-C(2); z=P(3)-C(3);
+    X=((V1x*Ct+V2x*St)*V1x+(-V1x*St+V2x*Ct)*V2x+V3x^2)*x;
+    X=X+((V1x*Ct+V2x*St)*V1y+(-V1x*St+V2x*Ct)*V2y+V3x*V3y)*y;
+    X=X+((V1x*Ct+V2x*St)*V1z+(-V1x*St+V2x*Ct)*V2z+V3x*V3z)*z; 
+    Y=((V1y*Ct+V2y*St)*V1x+(-V1y*St+V2y*Ct)*V2x+V3x*V3y)*x;
+    Y=Y+((V1y*Ct+V2y*St)*V1y+(-V1y*St+V2y*Ct)*V2y+V3y^2)*y;
+    Y=Y+((V1y*Ct+V2y*St)*V1z+(-V1y*St+V2y*Ct)*V2z+V3y*V3z)*z;
+    Z=((V1z*Ct+V2z*St)*V1x+(-V1z*St+V2z*Ct)*V2x+V3x*V3z)*x;
+    Z=Z+((V1z*Ct+V2z*St)*V1y+(-V1z*St+V2z*Ct)*V2y+V3y*V3z)*y;
+    Z=Z+((V1z*Ct+V2z*St)*V1z+(-V1z*St+V2z*Ct)*V2z+V3z^2)*z;
+    Ans=Mixadd(Ans,C+[X,Y,Z]);
+  end;
+  if Mixlength(Ans)==1
+    Ans=Mixop(1,Ans);
+  end;
+endfunction
+